ВНИМАНИЕ! На форуме начался конкурс - астрофотография месяца - ИЮНЬ!
0 Пользователей и 1 Гость просматривают эту тему.
вам желательно полнолуние в марте-апреле или сентябре-октябре, вблизи равноденствия, когда Луна находится в полнолунии близко к плоскости эклиптики - тогда можно поймать днем и ночью тени от Солнца и от Луны, так что они будут близки к совпадению, т.е. если выбрать время снимков, когда и азимут теней совпадет, то и длина теней будет почти совпадать
Подсказка на этом замечательном фотоhttp://www.astronet.ru/db/msg/1334803
from pylab import *import ephemfrom scipy.optimize import fsolveimport mathdef optifoo(args): t1,t2 = args d=ephem.Observer() d.lon='+83' d.lat='55' d.date='2015/04/22 00:00' d.date+=t1 m = ephem.Moon() S = ephem.Sun() m.compute(d) d.date+=t2 S.compute(d) lala = m.az-S.az baba = m.alt - S.alt while(lala<-math.pi/2): lala+=math.pi while(lala>math.pi/2): lala-=math.pi return lala,baba #return math.sin(lala),math.cos(lala),babadef print_optifoo(args): t1,t2 = args d=ephem.Observer() d.lon='+83' d.lat='55' d.date='2015/04/22 00:00' d.date+=t1 m = ephem.Moon() S = ephem.Sun() m.compute(d) S.compute(d) #print d print m.alt,m.az,m.phase d.date+=t2 S.compute(d) #print d print S.alt,S.az lala = m.az-S.az baba = m.alt - S.alt print lala**2 + baba**2def str_optifoo(args): r=[None, None, None, None,None] t1,t2 = args d=ephem.Observer() d.lon='+83' d.lat='55' d.date='2015/04/22 00:00' d.date+=t1 m = ephem.Moon() S = ephem.Sun() m.compute(d) S.compute(d) #if S.alt<0: return None r[0]=str(d) r[2]="{} {} {}".format(m.alt,m.az,m.phase) d.date+=t2 S.compute(d) r[1]=str(d) r[3]="{} {}".format(S.alt,S.az); lala = m.az-S.az baba = m.alt - S.alt #if(lala**2 + baba**2 >0.1): return None return rwith open("123.txt","w") as f: def woo(i,j): t1,t2 = fsolve (optifoo, (i, j)) k = str_optifoo((t1,t2)) if k != None: f.write("{} {}\n".format(t1,t2)) for i in xrange(4): f.write(k[i]) f.write("\n") f.write("\n") f.flush() for i in xrange(4,6): print i woo(i,-0.5) woo(i,0.5)#t1,t2 = fsolve (optifoo, (4.0, 0.0))#print_optifoo((t1,t2))
Пока Луна не закрыта и небо чистое, успевайте делать снимки.