- Messaggi: 488
- Ringraziamenti ricevuti 158
Accedi al sito per partecipare alle discussioni.
g=cos(thR)*cos(thI)+sin(thR)*sin(thI)*cos(phR-phI)
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.
m=cos(th_r)
g=m
m0=cos(th_i)=1
[m0/(m0+m)]
1/(1+cos(x))
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.
from matplotlib import pyplot as plt
import matplotlib.image as mpimg
import numpy as np
f=61.1
px_x,px_y=640,640
mm=53.
F=f/mm*px_x
xv,yv=np.meshgrid(np.linspace(int(-px_x/2),int(px_x/2),px_x+1),
np.linspace(int(-px_y/2),int(px_y/2),px_y+1))
zv=np.zeros((px_y+1,px_x+1))+F
raggio=340
r2=raggio*raggio
C_s=np.array([0,0,0]).reshape(3,1,1)
C_c=np.array([0,0,1000]).reshape(3,1,1)
a=0
i=np.array([0,np.sin(a),np.cos(a)]).reshape(3,1,1)
i=i/i.sum()
o=np.stack((xv,yv,zv))
r=(o/np.linalg.norm(o,axis=0))
o=o+C_c
L=C_s-o
tca=np.sum(L*r,axis=0)
d2=np.sum(L*L,axis=0)-tca*tca
thc=np.sqrt(r2-d2)
t1=tca+thc
P=o+t1*r
N=P-C_s
n=N/np.linalg.norm(N,axis=0)
th_r=np.sum(n*r,axis=0)
th_i=np.sum(n*i,axis=0)
m0=np.where(th_i<0,0,th_i)
m=th_r
R=np.nan_to_num((m0/(m0+m)))
i=raggio incidente/solare
r=raggio riflesso/entrante nella camera
P=punto della superficie
n=versore normale al punto P
th_r=prodotto scalare (n1*r1+n2*r2+n3*r3)
th_i=prodotto scalare (n1*i1+n2*i2+n3*i3)
m0=th_i
m=th_r
R=(m0/(m0+m))
i=[0,0,1] (Luna piena)
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.
plt.imshow(np.nan_to_num(R*m0))
plt.show()
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.
f=610
plt.imshow(np.nan_to_num(R))
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.
Ti ringrazio, e parlando per me stesso, si fa quel che si può. Rispetto a questi temi io sono un autodidatta, perciò la mia competenza è molto limitata, comunque in generale l'ignoranza è una malattia che si può combattere, anche se non la si sconfigge mai del tutto.Nichiren ha scritto: Sono passati tre anni credo dal mio ultimo post.
Colgo l'occasione per salutare l'utenza di LC,specialmente utenze da molto svanite.
Detto questo,LC è forse l'unico posto su internet, che mi faccia rendere conto di quanto si possa essere ignoranti su molte tematiche.
Grazie a doktorenko e kamiokade per divulgare e condividere il loro sapere.
Merce rara
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.
g=cos(th_r)*cos(th_i)+sin(th_r)*sin(th_i)*cos(ph_r-ph_i)
Accedi al sito per partecipare alle discussioni.
0.000,1.000,...x19
10.000,0.800,...
20.000,0.690,...
30.000,0.587,...
40.000,0.496,...
50.000,0.409,...
60.000,0.340,...
70.000,0.290,...
80.000,0.258,...
90.000,0.225,...
Accedi al sito per partecipare alle discussioni.
Dubito che lo possa fare (puoi misurare irradianza, radianza e fluence) ma credo sia possibile scrivere un plug-in per farlo (un integrator ad hoc, oppure più semplicemente un BSDF che invece della riflettanza ritorni g), e non dovrebbe essere nemmeno così complesso da fare ma occorrerebbe comunque del tempo per capire come funziona mitsuba e la sua interfaccia ai plug-in.Kamiokande, mi piacerebbe verificare la correttezza del mio modello riproducendolo con Mitsuba: con quest`ultimo programma, e` possibile mappare il valore di g sull`immagine?
Accedi al sito per partecipare alle discussioni.
0.000,0.116,0.116,0.116,0.116,0.116,0.116,0.116,0.116,0.116,0.116,0.116,0.116,0.116,0.116,0.116,0.116,0.116,0.116,0.116
10.000,0.142,0.139,0.136,0.133,0.130,0.127,0.124,0.121,0.118,0.116,0.114,0.111,0.109,0.108,0.107,0.106,0.104,0.102,0.100
20.000,0.181,0.173,0.166,0.159,0.153,0.146,0.139,0.131,0.125,0.116,0.112,0.106,0.100,0.096,0.092,0.090,0.088,0.086,0.085
30.000,0.232,0.223,0.212,0.204,0.195,0.187,0.166,0.150,0.134,0.116,0.109,0.103,0.097,0.092,0.087,0.082,0.078,0.074,0.070
40.000,0.290,0.280,0.270,0.260,0.237,0.214,0.182,0.160,0.146,0.115,0.100,0.085,0.071,0.058,0.049,0.043,0.040,0.038,0.036
50.000,0.380,0.350,0.325,0.300,0.273,0.240,0.200,0.170,0.147,0.113,0.093,0.075,0.060,0.054,0.048,0.042,0.039,0.037,0.035
60.000,0.530,0.460,0.403,0.350,0.296,0.252,0.220,0.185,0.150,0.112,0.086,0.063,0.056,0.051,0.046,0.042,0.038,0.035,0.032
70.000,1.000,0.745,0.537,0.460,0.360,0.290,0.230,0.174,0.138,0.104,0.087,0.075,0.057,0.048,0.043,0.035,0.032,0.030,0.028
80.000,0.818,0.713,0.593,0.487,0.393,0.286,0.215,0.158,0.122,0.091,0.076,0.063,0.052,0.045,0.041,0.036,0.034,0.032,0.031
90.000,0.665,0.594,0.510,0.440,0.362,0.280,0.200,0.141,0.107,0.087,0.065,0.056,0.045,0.036,0.030,0.026,0.024,0.022,0.021
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.
------,- 0 -,-10 -,-20 -,-30 -,-40 -,-50 -,-60 -,-70 -,-80 -,-90 -,-100-,-110-,-120-,-130-,-140-,-150-,-160-,-170-,-180-
80.000,0.818,0.713,0.593,0.487,0.393,0.286,0.215,0.158,0.122,0.091,0.076,0.063,0.052,0.045,0.041,0.036,0.034,0.032,0.031
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.
Accedi al sito per partecipare alle discussioni.