From 74c8746464d1e19b1b94b92d0b439db5a8376e90 Mon Sep 17 00:00:00 2001 From: Zhang Xin Date: Thu, 29 Jun 2023 18:17:49 +0800 Subject: [PATCH] test --- libstraylight.dylib | Bin 112072 -> 112152 bytes pointing_50_1.dat | 227 +++++++++++++++++++++++++++++++++++++ pointing_50_2.dat | 229 +++++++++++++++++++++++++++++++++++++ pointing_50_3.dat | 262 +++++++++++++++++++++++++++++++++++++++++++ pointing_50_4.dat | 267 ++++++++++++++++++++++++++++++++++++++++++++ pointing_50_5.dat | 241 +++++++++++++++++++++++++++++++++++++++ straylight.py | 120 ++++++++++++++------ 7 files changed, 1310 insertions(+), 36 deletions(-) create mode 100644 pointing_50_1.dat create mode 100644 pointing_50_2.dat create mode 100644 pointing_50_3.dat create mode 100644 pointing_50_4.dat create mode 100644 pointing_50_5.dat diff --git a/libstraylight.dylib b/libstraylight.dylib index 74b2ca28177071ec4a050e811452cde080ed53c9..cf21e89dd40060dae994ca6bc920eab150144e0d 100755 GIT binary patch delta 21597 zcmbV!30#!b_xE#$859&4l|jIT9aj_)6mdgzbjsZ{#ati(McG1J`{|%KDvc<5GSt%0 zG|^HqT*(DYMOzeqWiFLcSq-@_6_&|+&b@bJ80-K3-}m`^=05j(wtMcm=Pu7P17~hJ zuDj)!zr5*acZFYHoD}!b?yEQSY2trHpsP&W{a07_ZTjS@0E;E@$yJDAeRSF13i7>% z;!aot7aFvHBIc{h_(@QYHxPHi4zNhm93HViUFdB=J=0X&2|IAzPzN;(Q7vb~;iBp@ zHWxNKv|wZ5s6&Tlb7MsNDxsnvMvYeBsrGF)U-CYbJj9$gP3^0T?l1JVN*-b1=OMMP zN*MrG)x+3wXr~#@3SpjRvdXg?T-A&k7|Ddp0a`^wBhYk16lNYSzrN1=2|WnD5pN z_g%Wg6oP=0x6x;ipxMIfoyy(Tl;zxGZU@RtWo~QB;7#2y)u@wjLpO;%fc_qt z>>(WRIN^5Wgrblwiy?1;*0LLxHJ_=ClC(-r>v#g3Jp=pv`586QxqP;h!L0cW!Od`L z?V%v+Ewh^(uf@xfjmk*|w}oXl-3?~F1*>UgmhB}eo{gE;)iOmfqRkM!6b?-N(vg4#N=90mjH|@{bT2wAJ$nLZgw=DawrWb%DdTU$!$DwVHt1(Skg(*F)BTeR!& zBNC61&~Xa7AXB34wGh}c#As~G6azyE$@-d1X;CXjZ(%o&J)>7x z`g7tfRy_(MgmiL;gc;J&$#MVe%c+frzn+UaS}8cdPXGU0F+5`YfJM&kHw^N z>$DkKNy%IAyhhU_aJuy?4fg~ww54qf3w)YGv`>WUl^~eo^NQ+eAawF=1^0b=s#Ayd0)F-L^VQ2aOlyQD zT7r8N2TxiseufVIFR`z{=C_6EV0QB&3F2is5izhmHeqQK)X^ zK(k6;DH`r-l``b&57m)Z9)Ny9?yVfTr?I9uoS5WUah=3OE}J9cOhxv}ynJPP(zdF%s^Ho>a6 zebA>(Kj*=?a*&lKN1;!0%a+{vg*Y!BxQ*k}fe1q}F}0d1C~xDS_j=~HY^UqqW;W~fkTC>4H1A&ffj)7S80`#{}! zY=L@s*oLzaT6XXov|e&vujA;Z{!fielM3VO6kIwTh9w;Wbe$!uoek@6IMSiD?#XV^ z{_ltCJo^;xcIe(exT$X&1@q-8f>Th=x)!P3fTeWih2ml-Ca^%zfORtuU*2jEu7 zP*wm@om!~O2TI0v+QL+GcEMkr+xmov6jiGz%gy@F4eJbf?zlG7m})Y1K}=|f%C7`o z3mv5LEP{9jT6NnM7p zl~B?p+U2_<>6anbG!D8%ghP1P9(Vor=K@g$e}wf{Dchk}*I4!!yw&w3_6z*lHK0Kb zZk(13!?WgH-nVc52&9JlI`>98HUCmE6od~^P22`Qh4*D2LQuqARq$5u=oYK;`v6|+ z){E7`?r!16L+_EDwZ(t)V?xzqIwtTf{d;1VwT7HuJq_lj23AR39F4>!Wo}F6;DTp! zaQOnsxpZ)mW}<7z3_mo@CNp_!FO%UslVEWBEjR<>yEpSz37Ip`$yCG2?p@eHDC-`n z8ovep>>i|Y+5&!&K2C3J5#20=7b06nrX#WryUU{wICz#FCM~YCvm-=;=xro5JK^2P zWL4y5aPHApwQO_A;2yo0XZJ!pjQqu_U!oIJo=HTc!ty)38#Q6fOPVl?j#3Th6IAqeF{-FR z+9zoH>5R!pEmtMEjm1yflhhhnuQ2s`rcE^;d~R}ERpf~No4GrrRc3KQW0cIl@T8{2 zFV-ZSiRGhe=4K32&X~D3hAEZMrf>71MUQLinJLD&keNH%=uT#unDR9hBt*`Lw#-~8 zmNRVaf?wZOLq|WXZ*icsxa5Db2>)&|Aedq<8XhKXc}rXNjk|sLbFCfC3--O+2&|0JQiXUJV~ zYJ~dg^X{no>D-Q@F7}9@%3NaJgU@5U@gSin)`n=1u%YIV+|RXT2@2DYtjDR(mtdQb z6T&L z<{Zbp#}<=eP9dYNNKvKjg*g+q7gBO&hfAjjz(Pl11`j( zm8Met5nns`fl*f_cL90Ae9p6LAbL}Yayo0^w{ink&%6e;NZumF-l;R~k|Cx3&9j@` zgj-2cH_k>k+)h^}wJ91l4)SEH;lQ8~s^~S)&@k0_d^L^&Zf&?^Q~$-v`TiO?VQH^1 zzRf|}W!!7<*`motxx!O;gm^WIW9Rq^Pwq;AI(skM+UqRsS zl|2K2gY9I?N?sk@hPfYlUbq-G4zUw1D*0vz3Ku>vybe5uh1B;VzhvaFHq5vjg()r# z`5bLx%0^tlAC@u0aW9F(IQ8UxUJzShN7^N2!LV8$6%PYJjiDOexPeR1+tc z;m9#_4&RWc@eTRq5nj$DVnZ4}vH=uNY0}7E^!k+M@Y4vLc3kd1*rC1tIXiYFvDYsd zsbhy>!^mK^4$4P*u=Q|xWUy-d@-?G6soF1x&{3V(pD=k;aQ82k(Hdu}#5IxDxI%>1 z)p#z3%!XyQVP(LPQNhO9q$^LEWqf8)){z?J^)!Ch@;gTOEG>X;OUJ84_K)YYKxBV~ zDMLom%xhv)-SRp~*_Wj(&5R0DmQ~ist753#@~%11>TI(nv|28^LZ0Cxb!?t>X8n>+ zKJYgXra0-Pn1xMD<0K}!Si)(q*GhV98R}Ud|A*4S^rTHxMo3*@W0Hw54uF3@46T|r855bH{aLO;;(bh+)|9;W|yLp$M{(LUbNkeIpkJObA-o_+qU0h(6mr24cW_llf>|IjHIH>Q(cG`U)S z7ax=!<2?PRhBlA!Z1zvur^j^a6Un2XHIWtD@e*`CUi4ysuO4noEBVS{H-!ej+^DTc z)>zX$Vb10AJidhT2ZRyU^MWOC2`qfsdqf5b*{>}1Li228;UO53`_xc=$AK~qhh#Y_ z+YJfq9o@``V@#u2!yUjpM4ZD%j#X>%0hC5udiG&r`3v8 zzm9d()o6K#@R>_=P*W&zTZso&t`M7s5TNZO5Lw63NUg;=hff$x{z)dGKzv{1n|Vna z2;WmYVfWY&S5hPr-&3?D)njq#Z}_~h0)eq1ZuNv~{w|4&4QB55@Uhy?v;{Vdw=-Q; za&bJGzFbe3H(CkW33bBUwKFAyC*ZD9hQij<#a!&0tn0HY^tfK&nyX*lllR;$?@bMT z2ZCWC>WZGcr*3((kV?zh1Gn}=ms4bW;k2f)r2~Zs9ptO%Ku~9b6K4pvJ5)OnG1)v)nM{zqm(TFt7(S?P;T`YxmUB|9iXq zKtkN<5BUFQSwFD9R@111^)^BCiE zJ(MG%uB_1G`b;s;w{8SIac!^A<8hX#Z(3LUmVCv}i=(*lP2R*8k1=oZ*yY>CPOGjQ zl3OjMc+VSJ{e%^hgIxGBW>th^lf7K(4FlYm?B`NR8T?G9X`Da%5+dUKT{7lB+cNd* zU~1eT_7Uuf^ENUYdZEFrUnh2f(s_7*;9rcsNw1)tWViL}P^3_Z6#R=wA*P#@lP_}0 zEY(NE;EMhCIK36J=utnZ8$E#IGu0ss_p`;ywA{86b`e7BF7H9VcpbiO92@UHq}NC4gQv+4s_vCpCBQV2b4j04T{(mo%tJx^mSuD~(@|3MiW{ogd9rkxEJH@;x=9fImP z`W3ADKlp8Gj($0-o`@BcG6|_keEWIe*0zZ$1A}=S#cM4+=g>mJzaQjO4DmB#&~eIj zb_5zH1hW~?HDOrW@LBw*6;1P)Mzsh>RhTJ84yR&sE)lRJ;XU>`3{4DX;~^`ti-ArD zbQ&u}1#*Z#>%}2bZ`|Mo$9Gl z&xDOreb^WR9eSZWExzP3g(+L;@F?*Yk2^AXi!KfRvu0XIvpV`!Nk0wI)7;r+7&9$! za2>ZLa`~W8!~E6HX?liqqqt}El z7oc~tKkEZ=$-0(t+?4n>t*m;DXk{pmn;uptx7T^^dw$GjLuGQ8=HKEn@%)g5gXeUA zHVb-A5A4y0OK?9iHMK*tpABJ-{zF!Ms)n+?h*cM11^F`D$LX_S?erkl5I&xMh8bbW zj2BotI5*=()xv$?okHKJ4o>ODcEIA4&P)#?s@tUj!zR22x$&RndG%N9j(wb>7HnRJ40!tR;rs?R?FpDZ8M2OmImRuEG|QdWD_l@DP})+VD0b!q0}$pfcx zx&VAm*GRx~^PvP~Bz8vvzQa^YfCGVF=}W)|tqGKo2Jfcoyw@Y)pE77sDx6F8Q%y>R`>7tPfvMn>7UTQtGMgltm137f2`Mlm zEz)(35E5C`3Kq_Z@r)59yDU060~)3GRoz?)L(+Yr$6S|Y@1RM!EQnl!6j+h&IeZas z3YSOf-hU)mPGF7%`2;d0SVtgLfhd2= zr+iU*Athd9&t+CBpg-WWApH*H((=G?(G(!-+_p%o_naUteInoy@HSN-KfH&Wre! zsOYjzz#wW0p`(CbO70&Mv3KFJcW)p-&0@*>M!+W$eom0Q)11JgnC#!+i3+^`pP?#; zN$}f$*0LHfFYslT;KKzijipnqsozAYmIDzv4K>Z0^u=VL!je7#?|BgsL`4`UV1|Sv z1#B&048mhWrK=ZT;hKM88>AU5MCVA+<_NHSj&%aPu|1y2GmJAm z*J2W%{bYUf%4K$D0GZW;cr`(p`;}e#f0suaU>;t z%qg3MKGO5bv6o30dv)E0t=Xq2O;NbM#4&JKyw`buf{7d1zO>iXFW;#@f*Z5Nztsrobdq3swU2O&_`Li03{F_kMw>rQR$L#w_j1 zA|ZEa2lfN(UmC;SEvZ?mcW~Ovrw0A{V%HlTA@hwer;(Df6e*Lm_bENmYo8Gm|5>3h z^bK!_F!yy@B^d|cd$S*GGDkUuNy;dsWZ^;lNYzNl^ zX3zi93LT5!!D$zFlD{Ksrw$@tZMHEb|! zcx$>-^HIWLFt`-Vc3Or*O7^$Iih_3o564i)Z=qFmOV5AznSmq>5;iPw5)5BAzKKRq zwN~n1P`b`+{BAfIlM8_6OM%AGi05g}3KS>Mae<}@v|peMfj$suu0WeP(r2TI4Whan zt01cRSY08i_}0utV-od3Brk*u>%+$k}S zjZI+Np$*Fx)D?O4Vl*t;;?9YW7h-dmF-pssNP zE^QR^wn)m=-=p0Yz1OV*_qToVXWkKSXEcgHRi2iXL@0jSLz~pMp6D0w_1l50BmDLD z%-a9HR%J{%VNZ;_cWAxt7-; z8>;EiQdHF-1tJE%L08OP?bQ3 z1+oaVN1%HGZ4-zd($LD=0?{KHqSXS?0~(@b0@33cq6J|B($g8jnF7&Tjwo3mzA@m> zTgW!f4HtN_y|>>y&Xd#kTp0ENe{Mt)at3_~emg>7SmQKcm< zPgHSH&@2;GTG5t>Y7thmMYR~K{}I(gSe+-TWmwG;)iZTe{Id{ECZ{^$LXJPsgrAGU zTnBN=c@eqo5WKg8YoIMD4JPmH$_V_M5xkW)Mj3v{J#T3M zVf#9SUql%^M0Ca-P(|wwF4X3=VmpbL*IgJ!9mO&b1;tX`l#DYSe7J9bN(qJs`=Z!M z2-)8OUwOsuAIic?w(sxB{C5Y)S57!47CK8~+(+QC0mCGIu?Z#WgV7FcgK?d}R?S1J z%wulgiglVY$hy4G@U=GD9a2h%sg(9`q;wp6zeM-(ZN{=o9FAtOxPT!3OiV8;#S1!r zB=8q>bDH6!C#ShQ4dk@)*)y@HoVlgqNSl&94TpX{1js+OD>gt%lr@a7S{42TA{UPzN|7`XVW|$@Knb`=eRcOpcY#1+kkbfPrPqq%3fE0>t6~2Kh z!@z4=h;~}yjikOsSCLE|a!Q5IP6px8`2NWjF~?GfPp`doA7d3CgJc|XEL!muZ9`-g z)??8(w-ab|(VET#!%lhN(J%2-Csp4RSbNI9a{x-qx3$!sJ3?B?$c;0tgd0@9HG{|Z zH~9Hf5WSK*-QvafkVJ3es7b8adecyNgY0p{Ui9y5>dSO^?XSx-nMW7F;(%>SFTd>TB?Apwvt;YyONQBP=W4thpNpuYWZWUcS)m#clW=kvp*l zw=C{NI?7dOa*<9KE;0C!fCoU$Sgx;fTXGV~W{p9Zb6c{X!ya)Eb#X5w{MY@(4fv43 zP3GYHYHl)NGMQY6bh`Mem8-8ZYla{l`{lM|Ca$(-O>YjF+mZwhJ54Sb_H`Ts|F3lY z>Fpr5c@n<~=QejtB%7Zgovw#CO1B>3l7)RM2ElN5WN!}J_bE5elf^1P5Sj4{$`T<^}{Pzh>V zHw&<*pCi?#r$;$_*3;c3L%+^q{o+%kv8WnptiyORrbjwmW=#vk+{|k{LX9~@Gp};^ zcZurqH1=MvtISv%)CphP)CPs}GK`nuyzIuyNM1(qvKKG=@G_d0F}xhW%YnQ!@bX1o z4&mi6UXI}9C|-`{<;zqi1&!liJTE8kauP4&c=;+X6L^`#%W1rv&dU^Drtva^mzliG z!g5GntKQQTrQb`nK_NjUFMKnoA)V&%x@#=S|M8|%>+hQ4pSEfnH)x1|Ppjantl(Cq zJP8&;!Ql7E4@z%)IQ<|gO*!RHr0kQFZbb3HeN#arSkgsyZs2F9l!@b7^Y3^#eOiwa zb*HHdCGW42nl9SA(^Rt+sorqmPEVH`SnJkC__LVvqwY3!2~;V{TuHqwtlw>_NkFO= zR^RQZ*~)7d?xH@{7Pk}n^Mw9A)W_$WBa%8&=-+GVn1IwT+QFK8Q7&7tdZ?pNYYyeV z`@x^6b_1!aItl77VeCE{Q)9>^6_2kNviV3&KMu0>i0gXBRMc1Zx5D0 z>F<5;#rWgj`?475@yA$I?NZoT6XrM?y}Y&*KK-K`er3XHx-)lpt){Ev6C~Fxg(#E= z$ASNrguAIe<~lqnHh|aT?oOYs5j|f6U7mPB=3@^Q2)RhVyiw5m!``P}@WW#d)sLIt z&&Sh?e zqmQ`J#ZNrgTd@8~cQzX?U>Oche|KlwpdXf@Fb~T-*!%Zzr_u{Tt0!!}?FIf%Jyb_8 zK-AN2PIn|@4`-}->cJ9W$J6et7kvM;yXwSO;9c8|^@EYM-SPP&7t3u>UOUG*a~_5U zuf+@Zk7RrLLLqBl^!P^53Koj{Xn9o1koOhx{+_(=mU4|H9)<@jx+>+}f;&51i*diE z>w&x%;UST(4R{iu%U{qNDJO9Erpq1YYJ-*S8%uwz|JvY_xq{u#U)-(b3U~b5dk^D$ zsqHJ(+RJ;my!VlJgS?ND_X+ZzB<~sWK40FK$ooopUnlQdJXh&@G#}X;j$_3tF3)<7n2k;yZYlGBp25)+hy zA6Kw%3d-_WcVm^RY4qgDlc&s$pB$f+I5#;<$yYKhf|wdVb4KQ@v_vA^6Qr#4l=L}? znMA3vQf6gj2uijC+KQi&K9ktXIhdK8Mj)Sqtn?%T>m*14Chue2iW$7PPvGZMKYRUBnxCZ}a_UyK$(m?c#2 zag>=DpFoCcB;xI7tI3_2iHS3awVb1@#FP|5`5ev3#NZNKM=&)$B_*BvQ^;X@0&lR0 zKt}o;k}Fm>HAYX)Oq`RP#@jnYJQ=YvUaxcX49Bl!^6{zUVCIzgG;Z@Uhgq3ZGUh30 zwo1Yj(PfK-X;V@&h<)!FMsb&Xo-u5)5~EgDzLn1c+`o#(f144nG)hm+z|hAkz0^u^ zLHB&t+wriogR*Q^!Tfv{HFkT9;-Kv7kHt^nSmgJ_{}1xgXyq(R$I!l%J}YxdqVfxy znWfxf4$A3aidxBDuEDbGq9cmm8-Q)-<1@3SDK%_F3TB?NO68!uKd_+DY8KFSdta2i zTFprQZyyvLq8wMHCnU#D!9(3em0CG+k+Ckyb*z=%B&omOEr?&uJRAI6=zwT#LH26a zwdYpaq#e41Mb1+WJ9G_8N>aXZ=oYS2IjEH%wv*iL-6)`%GGdx?)gdWCxs8nVJBV?h z;Ld8+o)r{$tzrIszj}&IsPAE^sFG(YTh&?V%5}9`y`q*Vd_E*ASJZqOBr6Y*^B!{4 z>fL`+xuRe`8fkdqT?fTMeXe7{)-`O1sxqYD!5S9yqNNKG7k9wo0To4^vDnZRi!J!H zt{Rb`9B@p)Z=KFKI;c;yL-L#9NbWi=cD%An(@mc-ZAyBoa?=sx@u(mkSbVRK{-9`l z_>dgyqMxmtYY>Lz4-M4Ha-^x1rH`x}3+FsRP9<_2)L*qNz`x1e>io%llu&n7mNhKs zy@L&R+v#kxZtQU`SiJ*R=6%i%>hb}owSt`bvn^UK+u@5zbk-T8T7Xe?#E}Y{o}hf< zjI$D_cY0>ROy#jNt_G$?8OoMMe2RV7NKBdH#u)NPDztLD8joYDg!EajrX(u0D4wOf z)kv*;{3|hQt`jppbEYwI@~b$vW+?A8N=d};UJs$jiF+inWgUqK=CMZ8a315zbP3s| ze5~SUPDxIltekB0YI0h9CVr`2gRH_mbvh5MI%aff-$pX|$)|Pf)s4ByL7B=<1ogwO2+TFZZtWaaxy zY$16Jv-Mvn{ULFj@RVsd_U7sVtTDl zS0sK&zl|S<%20`o`y}C2Nl5e-0@EbU$M=PFrQ(MS3~a-pGM_McF#w5p(LfR|+bQs3 z!dT59&+y8}58>bd4wX#;8PxbPmuV)MS(Lo z&Q`3?JqskZKL5NSvGqBqfUuEeE7k`hdcEL=*!p0!S7Pe}Qklfo2c=4m73GBN@pZ|! zK0w`-*!m#FRKmXXfvOo|BQg3q^5IB}FC84`m7bbJQqTK#W$frr%tyJJeR=(Wr z)=1n_>~>0_#5H?_{63>3#2gZY6B7Hc5csOZr8K|rdMI&WfxwM5B0~FC3LGqP|7`;I zmH5YX0>?`1DII7sQxdA*7Xm8_E^J~tqkVnT2Kb{g>}P&(59SX>;bxnR4$g6yJ8J zx#dup-zv^i;g@rjxW{>J*)iBX=?j5AVd9>&rG7}OLDvOXD~W@y!$3CJv+{-@o1Mg+ zuo^ye^n)@M;92&Apk8nkcfxkCM%@N}V?mxtKMQJpD{&`m$92Q)>>)~(&z^%%RqwNO z*lXvw5{K4Bkpmm1i9hN-AoR?%D;2AJN=#@$yVc}O>hUZaR;XvGYWl#J>WNddD6h?C$^*5cIIUw!1@@%U#FYN%z=O?9iE&UA z^2HlJwcD9eiTC*ZOi4n;7FdtQ4_yQ0M~8OFK7dm4A0S#MgUH3vYsSElyx zMM10Ov^^K0Z=0aOaqpuYP1wby&U&Ld7QwCF_cU7dR6*96m%3kQz{`9Om3e!&JoHAL znYU_gU0&R#CG&}@R202QSKF{e->_t8!;=0DOTzU9)nik_i_ zWPL`a)}B+89G#inI^RgIHlO6g+pK;fh63pUvfa#caM7!AT?A9JDyCaZ=(9*q$0$N&) z;E*C#ZSIHlTMP9C?_dad%)WmMRST{WZ#u6cNA|DU>`BsgTsp6}#rp3g*3_p@5K0-* z#HURbU0}U`kZ15E%A?pC%&&1ZB0EmK1(p8(o_}Lf+Hi8;f_whcRc{5s_;wRj^Ml~k zcGH;)G_;FV9qqWay+15#KcREV8IfG?=#9EMky<~#X-QV2?wTNtf;;W~JS!xl24y8< z|C`W0ptI^w1?U5MvcB+Qz)01aH{oJHJJtIY@J&E3)wYY^)uBVcK{*2+P>PpIa2J8j z`@Uwi>qMB^VGDj}xBMAe(8U_fAJ_ADd>b`2hQ3Imx_6~L#7 zzCiat#O6aJa}qVGhhu>qMkz=Vv(u=1AUZexG!LE5OmuQd?ToVuXQuA%Qipk(JC(Ou z5%KT^933aB?YjWp9Y?T*Fr#BQn+T?kk^UtoWf(=nyCtkdaaRxcp<`D=8P17%e>?Qk zsH+tXzhaR>k1LwUe<2xN$&06qH|i>da0aKwaw<7`v`OLJrz(^;Q-IAuo#G=;P`~b= zSGYV;PllsY1nLMpgMfNnr&-$%S6w7wMrcAX*(l+(>R>KB<06l!aB3j7!6A;J`lWN= z8650+6KCu_O*66kCt!SV7nf(Tc{Jo&at_u8M|iH0a4uI6ZH$4M;2}P5l}oeK3bjA3 zj2yUk76Lj2sMeOl@J^jnAC$v_PQzTz9T!HIoQ0E}I<)x~r_4P~%tcEV$6X`Mk??aT zU)BYjI(JfScniWin^-S6-+6@U&QY)n8E*gjQAL69u6}N`tKfn+^A;|6M<>G2UTtAh zh$kBdr6HXRkB&6YvNAE)fXdXZ@i>*xS8Q7AI+oYtE>b+pG)0xUm078JQ=clD&_;Fe zXL!S9VlHapl7I@LgZ{z-_SoH2_Gd$bu-{AAH+3S?JxyEv1Pqvu}9GFh<_z#Q7nGgWe? zS#f%RU)O$)OU*bQ?~do=v81afn+L|OotP&a?>biX&0%Qh+SB7DJeXLPfOkZ1!$H@r zr{*y-co?YTS*KuKx53>GS!fuVYqy*&bK;uN+iKyGhBMsl=A&6?MN+?g61=;Q3wE~9 z&`52Qky_MJ-w{(%mkn=(c$e?%p2Jk%9fr;!-OzD)s4fD8rEH9Iu5(TMmy$W z5*vG~Fa>>GF#RgO!|t{oi`|WII$B1uX% zcfx>(J07PGJQ0P7(6Q$TRmK5G>KV)C!`nS4IZncXvubX09X9vvXx^M*zg2TR;N@Nc zE=Q0=&7D31@AQgNT`dLo-a}YBn9%!3}TDn%RaqXC}{e6kD2r; z8Zv72IX`*pjji;omJ=)7?BG2i5eK4?5InD;??-!{EW zsHiu+RQu41&v*GO!KitHF&)nIjld~?yKjHT5BBm-Z5_}S!u#8+Ufc`g`nB)#2&vRB zdm2YgSdEybE;x^HrS`)WP-Me5DAt%iF|fB^Tes6zB()p7-!EBps08}XUGp4pdCJiZd7Ky$HZ|TI*K=d4EG1264!J=HZ zr-^3Ky@@n==4cJ82O3pB?uHhT!&R^DhVhZZF<)K=?g1@*v3bI)*gW#aNM;-x z^mA&Ad3yL}yw3#p+)L4A4{-~&veH&snaJk^S#i%RR9BM_`68P7b0;+w%Ng9Pl~s83 zjLT?>mJmJ8MR6+Qqm$S>o?@@z|6nhiyi|T}?&TS2b+xqmgcC;wJG)S<8hh&t(a}!L zN&aAnRyQZLfrudiF0Z`I)Az1pH6)Dm_j!(>Oqq+8mhMszw+tEXbHDh%bK7d{zZ7w| zb|hk`hb9bdVt`gx(uf_{2a|^e*mx4EH)6+Q>SP?-qXNl3#^nIxh-28Sdu^ay+ZH}Q zCNUDSS3O0-zcmgujew1&PiqSB(T{EYYy)8VY1n15XVLHq=qkiI?Th(SY*{R z75q7@V~}MoZNgQBJj0waQBK*~6hg?#H{1%h zbpgAkX+99Xtu{U_JpQdwpUz=67JT5`h!&AFRglT&on;#5@DXG6dR!{8TXIE-k#kh6 zFgci1b|bsF7j0`wDm(+_dBf0=gAMV5I+Q04QnlF11^hIaZy<-VqZn698e^(&<8^*Q z%kxGdvHCB^laxydS~sc};$e)I=hw+BKctyR4t3S>!VVv)bjNDZik*4swynCh278h> z^{K}doYbcfH103P0j>W&$g;Hck{;dJ-RP0$Ev$?4cHDau0cs8RrHdK8E?}*}+yJ5% zA2;c^Uf}hnUMzJqRA5Rfx72usqP;e2h){kf4NVS5OQqgTO0G__P3#!&;8I5$e+P*MJ*2K;rTYnAe8mRzG=FdP zMcpwNscTzBq^h`r81G^mWqIb~9mSJHHVz~1$#?ndh0FSa&&_u5?Wkzg@h#AM^!%WR zEjWs}$l*FqEsGi03Y+1{M*GvP9dKc^R||@p;7)8lHaYe zAn(yXLUS0mWX*Yi9Ah4pdgo%jDIHHKC+~uKOt>5A=#jJ>N$1QkV$4I2x+Y@|pf z)<(9reABa?nMci&!jpnXKKi6^84MrSy}1`9<*UYZX5))dnBtJaCl+p^ zh^zC`vG+6&jJy`f`@AZ)Ixnk{PPIth*H!(a3Lfe&*4ycgCyQb6WWVV$MeHx9KPi)#{-5x^qAw(!j_i! z1Wl)XFaHT;v0syYD|@tGHzT!y{m}ueD_oB5?h!}^m^zv zp@(YMdYCbx2TK9tgiig2tVab?UM)=wT*q1=_uvy=P52yyIPmc$4i$j21{eQjgJa6u$5AZ624lOE2MJ8 z&}mHjQL5S8d@-rI94e7=F>=+#^uhto(dDxGs7HLcl5g{x1Sk)4blI$4jRHA* zabWfJ52#Md46AogtCZX~bc@ev3;AVL)H2A-27UdaXGB%c=i-pILtMjFgjYe_i< zvpG)66F;2#!b@7>`KliRfm6G?ks^`!svl53b1E**-cJf^p*W_yrkSwQKjojqbYdRL zlhO~sZJLegyX9l1p=tADYd?6Sb&xx~NtnBKt^C6D&dgAY!j{)J`IbsywGJOUFa>nQ zk-Xa)(=;_M%{n{-5Or51@218SiBwvU?rB2D6;l*%;UvY#@TWr&Sv(41QV7eTx?Pxi z_zFAr9EBNI7U$LSCGu`qJk!&_qcUZee$GAPlw$pyk8tXKsAuee=EC+q*A9efoX zds4Qu`YtVFmFD7=iV{^&PZNuBnOwtlgohWKjrqT}`#8_Gf3Z9Rq6)sT*raZwO+E^i zuE2A%edbfBjTTB%=VZ|BYpZo7BD^J*wc8O_R};K+Ph(h8K1Wx|>YcH_Il2;7{}V20|HUBu zli=gpJmYpHc(;nZh|z7<03{-y^G``c?7cm`d>xlBX0W?1KS#HV)t|y|A9HlYtiA{< z$k!o1j&BsS)Z%; z+SAFa64WF=Vqno+U;h*oC2jdeE?GiutPV=&j`R+zcycfwXX2L{!{!7)z&tlL1;XZe zvbR7tFEHd$hDaOgI(CvetShv{tg!xoJaPR1o9B7=xQ#0gdZepBDfN#CaF^9pXx44V zSxlr4B%hp7SHq|Cy0tB&wUHn3Jor@SZ#u=~hyefjZ4K+tGfV7ogsHQF{Wte4R)*g_ zFCRH5aI{e%>c46B;*2KEY@r#0@f4;@M5kTEOkr04udJIFO`sc3RA1=2{8AB?Y$W z&KDc;Rn<|9-Z&n@@z6clm-PmHa;W$I6Hkt6W_r`8s@`7#;qyFL1YAuHWEpTTImov= zE;Z;qxkjhOaZLjY;R`$+Dfwe`g)n|WFjK*j1vShE1}_}LPDB2}XH_jv!Y2#E*+-yG z>C1M)pp-E7F{Gr-R%uT_ZOUkdUwkg=6IifRQI_)WM8j~(U}HMm65XlL^1?L+(iXL0 zLt*`*Kwp0p7BdOY7`jhL#*BrldMOkx9o~LIp(C=jE%BY4R^1zC37NYl;;L15grL*` zZVxwe84oUV0~V$ZZ*x$P>{U2eIOk)hpo&(V168R*@f%rnS`Wi?G2`gy&of*-ksB95 zsT9$1*I`Ag9?CiHd3P~y3HO)yfyChwTb`vuBzEGHRjUq=SSztNVpqyrBTr3H{W+8pnO#SBtZp%A0*%rtCIkO zt5si@0LMYA{!9Ws9v@0jN1P8hh@)@V@!RKW-6L94b5=p}l4BjuF0y3Q?9``blyMQH zFY;1lE`t0;9ekou{AoGWA2u%?(LQ&*RZ7iIu}P`t;n$@DG~I-ZNGVMUq@)gOb7LLX zv&pGbuw~g$Ro*&qNDDA%W{ZWG2VoPMrl2=wSRQZdX&IKuIqte_3Ph>|`2>BE}k#mkEgxFbHh4ZE=_&z*YK>##WN;l z%}vgH#&sut;WG_mj)Dn7@UINwza*;oa0HLdqlx7$o5o0evMLh48$$Z0zeeuKDKzVJK9H-j~pi zV8PpZqwX&K1PNcqX*gvhlGGRQ>0eSep_N#=p4L z@O3g*xQWLl>PS1`&JCg>)Ckx{x}ykqNpgRki5=-JJF)=*I|*j;L4i zmUv&r&Qu_N`SC*fO2Lhr_MwQBQzSh;$l%A5_g zt9!5zXjmPpG7<2Bk%m~d8TJ^w*lIXw2=v!vlcJ&@kXM20Z9?@fbQxjy(M zwp(_4x7c|oS5RGnolwj>ATm3wbCyA<7F;~jAF1aV7Wj3#CrX;zVd`HXg84FR$@`SR8;HyF*0lozo=@avh0A|WC4(>VzL|6%xKa|>w2-CFk)R-b_3GZ4QDTxA77`lbAFZ2g=+JUTQ|fxkB}gH>xOOs8+tl_M_Gp* z_JJ!K;+@w>bQE?PJ3x$apmUI^9)-tVzZxbvU#w`7TRS>Mu-^MvWqzA59u4_dV{<6i)n# z$l7um%85uuVhw}Bgcfu$d9YygrYEz9#U?*C8BC`6?o+MX_&ADAbDKIqr%ijk+T-z% zXTHvy)e zFA7jv(QLzoipDe$Nd>yQv_+zyo?YO|jsSKW?(9g9=!~lTdo}Z{A*h#Ag)z10#*uD3 zlC5;-L*dRK)(KAT%_`&MVw{133JG_*~;HR)e{ zMBnS@!pz;ij^lAIo6Gueqcc&?oR55pfP;oON$|}6F>a5r+SsL^VcY(jZlB{+Z*+MK zWFLrf3+5#0@@w$rfv$#kaa!vO@@gr>CiMk@t_gHfpgMtG7sxEoUV-ST4E2fxdLWQV zAbKK0UXDQYJcek6K=d?*D6Or4^el#Ou0R0-O^1@#Jeeb$dhMAG3pq#5>KfcHacIM- zH96h>4cE)M!?3dDf&a#dCuWG6#Ki29qI`BaplBKX%)ufS37Zajvw`s1!D6?uNNo&b zZy0yzIkyIyqn3$&4Gtda?(i;>%B)lt2nQY9a&&7Ql=L0^)6pjX}!x4Ra);1qKfN?dWERc za<^1e%dnausz1q zceJaUmo+H~e2(>S`(=PtNV>t4V~gC@Sd)%H?Xf-{&s&pHLui}13)7Evgu;Ce<-?CN z=GF^m2ZfelM9kqZp*(@T4P#IE;<4;Xxj+6C?sj?C-d8)(zLOP(o8wtpIdB~}z9jaP zJ`4Lo^Z2@Yiw6TGJB)J*T_g4dS$mOEa-MFH%m6Z3)w)@j$$>J4Pe8Em`Al zOSPC8M)h;B<$O1{709DGt)_^bDu4g{CB{ybZ>hY&0vF?oP1a0(G-$0cV;{A#_+A5> z-~Qh4oq@us$;DsV<5BW9=H$gX{H70|u~1N7poT^G)M-?2q#Md@qaR1p3A@)IVIEJb z?W~+_+XWr@%FEJ59@d&!tw91>)d{PE7wK)xeueA5jDfY;|6sb9*EgCC!&&ozeg=>@ju zloaT8*$1cp&sDzO?U8QOJ(5A8Z*{d!DB}d}FZED;u@J^z3Jg29&=H4~|4g6kzJZVy zA!_c7zi5u?>xgbtZ-$af!K&niaJb61t$b6{Q~%< z$_r**o{f)Dm6zXBl_tZ*%l)CM`dPLdG&RrRQ$<`&Pds05!7>NSHgLVhk400h-0j^U zyY3!7im$82gA#Ny zzBSRRyF$>_5o`=3T@CO(-(B?d#?_$!611-SVc*pdwh`XHI=tUFJgVbV^){E` z0L?9~ z1sYFg>R)GZCp5WLWMBO)hfI@ulEaT@!J-c@J=+_P{oJn{d}qzg{S!xic_W=JG>3;~ zU7b<=AEaZSG`X3$o*C7jbBN{=Ijo5*U-{7zHeyAZv{&`Fv^QcV*_(-Uy6)prar-{K zu*H`L$e}1bC}1eiY1F*Bb_dz`Xpj*YcY@00S?jDYz~Xc&wif73|r?ZN<;8AUT)`Q zF)w%Uau+Z6@Un!L`*>N(%LBYD{guWuTg}?GO0kmq>?FthEXO-#aLV)yW|+?an}TkV;W*!^?Lf)rq|J#T}HyT3jgj z-bG2Q7Kg$#tOnsZ|6nNWz-l5^w{?dvv094Ng<;_PE6QT^S{O{jYLGor!(j(j6R}!B z)xj!t1o++?sQw-)zd`K1NQWTYw-!Xej(aGE)o1Y!P45kKE5+(qJx+VPUJHZJ--G-% z+ikMPr$2`rNBgZ|ih`k_Fx*DLLy&MkQr*o#QO?3%tfgRW4}5b!()BRbobdG)U3HM~ za|d|#w^r)gNG*ee-%!Up? z1_>Q#mkX*t2CF=BVfY{TBkf$#*izVgrvn`N!^?SSo}jFUZ|-)0pZ@S->w*0_1i#cC z_UBa9xb?8(&%Vyx)(fFHSbX0H?)~Y7&ppl$`mxzC49ihi_HyWr0~Si+Jsnc8>bxF) zdeE1dp~J&|DsO;k5BuU5N?QFR)w-?vMMi`!EK} z)vyN3>+t5|@hYc}z~!I5&MiI?9ZBXLnf;HKYS2gU!asdkIUM_^AKL)+|MYX#ekz1c zbD`l4UidmOwxJ)(f-P9yhj$y6K?!qmsaS#W#j9pD{*CXEAxLvHJbzu#zQS#ut{w6K za7x}klJ|Of|4qtGlsFKNG<229`&rxz>3SZI6Lejc_m#LO(v|8f?tggldkbZ|s3>M% z{A=n>zERSj#=jkxK6`dVb~{6^l#Z8T3*$a_0^50Uqt@;*r3N6CA%yib$&M0rn@_vP}wR^B)A zyPu+LmxK~|KP2xbiWU!NpP!B&!Ab@}onyfEY}T^Hi1;AIedw%N)1Dh!H0)oj=jjRA z>{I)q{wB7j=wKe3k1Nq>?K<|YLy^6SeNyz9i49=ainN=UhoMf@3KeH9O2}N8u_Se_ zQmAB@1t}{nC2iT<45Hk%P?n^p3(5ncB+N<6BzCqPfz0Go0{I+drOhKyC_yrTVggd7 zgz9sbCmXVeQ$}dsoRqXo?!ZwF=L!=Q0%T;QaaSrioIf{f4k^}H5ci{&qqOw7+>vV> zrDr6kW^rHY1YROE%^YPI<|ZVPp}PXkaw>1`C_xHX9P%nSXmYfY~G&LM# z&Phn+HfuS|%9xYBLP6DQ5~heQ*GZT7koZ|QEG~D);!!vjzx2l94}4ppicVBcIV2_~Dm4yv_S?H6+1#fn0hpI+ z?e9e^fF(FS_>($%ujm_Ktpl&zC-A{<Dv-E(3e+{T4OTxByh=2!u zi8s^;yj$Wwt_XZeV&heTKOl^U!%DeLm?gd_@jnv3Epba%Azw$B_IPhe2>41E43aor z;{F!?*MdJx;-L~xl34jh@aIWf`K`cN7XEhv=UDU+Q~OGxB*aUB9TF!>Z0TY855j<@ zhl?e)^e{u>QqsqhsO5jXSTJs%+`tLqOg~3=vK|B%>uEUE(B(mrGng7(dw(mr9%~@db(VCH_+4jS~MMaiPSm zw8EZYnjx>&Pp#;x#1A+44?bSuRYe&2s)ZaXS0q-J3jBq{BWeY{ zCGla2A4@#_L%}Co8p^)~9|=Mr5%3f8iNJj%F8^HMF%o|*`7fLNde0vO<##Ek2X;fTOM%|B$kglN}9x$C!W<3o8=?XMvgHz zD0%SuScT(HA%F0qAmDz=@jba2F*|{;NbD|glr+#mVlRm!IL=lq&sKvZc9#0rB=+nm z+UsV(je{a!dFmQXn6lVItPA5gUr~B}AQa*xUU^2~#TGu_&9G+@uaUS=;=xCSe5u6w z6#^TMOG45{L3mGMC0pR{B>qp4z>g&UYm>lUY7wD={Q`$c+~s9~$4DG^MBq6RhrPzJ zL0OFi903fxxVGBheH>?_h`7GB!EX3!gegt@wl+A_2KQ}-as2ho2#PYv2FKapQX72D z2K(W1VeQaZ8*H-1v})B@6A;_s&|4c=Y=e(;j3V-S-v)nggYVercXe*EPwgu&SrG_d zvBAA?FSj0|VFEWgFx7@X#|Ed`;N?x2x3A 300: +# break +# # if i != 300: +# # continue +# # if i != 693438: +# # continue +# ju = surveylist[i,0] +# pointing = transRaDec2D((c_equtor[i].ra*u.degree).value, (c_equtor[i].dec*u.degree).value) +# # print(ju, pointing, surveylist[i,3:9]) +# sl = StrayLight(jtime = ju, sat = surveylist[i,3:6], radec = np.array([(c_equtor[i].ra*u.degree).value, (c_equtor[i].dec*u.degree).value])) +# # strayl_i,s_zoldical ,s_earth, s_earth1 = sl.caculateStrayLightFilter(filter = 'i') +# # print(i,strayl_i,s_zoldical,s_earth, s_earth1) +# p_cart= transRaDec2D((c_equtor[i].ra*u.degree).value, (c_equtor[i].dec*u.degree).value) +# sky_pix[i,6] = getAngle132(x1 = surveylist[i,6], y1 = surveylist[i,7], z1 = surveylist[i,8], x2 = p_cart[0], y2 = p_cart[1], z2 = p_cart[2], x3 = 0, y3 = 0, z3 = 0) + +# earthZenithAngle,isInSunSide = calculateAnglePwithEarth(sat = surveylist[i,3:6], pointing = pointing, sun = surveylist[i,6:9]) +# sky_pix[i,4] = earthZenithAngle +# sky_pix[i,5] = isInSunSide + +# e1_,py = sl.caculateEarthShineFilter(filter = 'i') +# # e2, e2_ = sl.calculateZodiacalFilter1(filter = 'i') +# e3, v_mag = sl.calculateZaodiacalFilter2(filter = 'i', sun_pos=surveylist[i,6:9]) +# # e4 = sl.caculateStarLightFilter(filter = 'i',pointYaxis = py) +# # e4 = 0 + +# e1,py = sl.caculateEarthshineGrating(grating = 'GI', pixel_size_phy = 10, normFilter = 'g') + +# # e2 = sl.caculateStarLightGrating(grating = 'GV', pointYaxis = py) +# e2 = sl.caculateStarLightGrating(grating = 'GI', pointYaxis = py) +# e4 = sl.caculateStarLightFilter(filter = 'i',pointYaxis = py) + +# e5=sl.caculateStrayLightFilter(filter = 'i', pixel_size_phy = 10, pixel_scale = 0.074, sun_pos = surveylist[i,6:9]) +# e6,_=sl.caculateStrayLightGrating(grating = 'GI', normFilter_es = 'g', sun_pos = surveylist[i,6:9]) + +# sky_pix[i,0] = e1 +# sky_pix[i,1] = e2 +# sky_pix[i,2] = e3 +# sky_pix[i,3] = e4 +# print(e1+e2,e1_+e3+e4,e5,e6) + + # print(e1,e2,e3,e4) + + +import sys +fn = sys.argv[1] surveylist = np.loadtxt(fn) -sky_pix = np.zeros([surveylist.shape[0],7]) +sky_pix = np.zeros([surveylist.shape[0],4]) -i = 693438 -c_eclip = coord.SkyCoord(surveylist[:,2]*u.degree, surveylist[:,1]*u.degree,frame='barycentrictrueecliptic') -c_equtor = c_eclip.transform_to('icrs') +# i = 693438 +# c_eclip = coord.SkyCoord(surveylist[:,2]*u.degree, surveylist[:,1]*u.degree,frame='barycentrictrueecliptic') +# c_equtor = c_eclip.transform_to('icrs') # pointing = transRaDec2D((c_equtor[i].ra*u.degree).value, (c_equtor[i].dec*u.degree).value) @@ -489,47 +552,32 @@ c_equtor = c_eclip.transform_to('icrs') # sl = StrayLight(jtime = ju, sat = surveylist[i,3:6], pointing = pointing) # sl.caculateStrayLightGrating(grating = 'GI', pixel_size_phy = 10, normFilter = 'g') +strRes = "" for i in np.arange(surveylist.shape[0]): - print(i) - if i > 300: - break - # if i != 300: - # continue - # if i != 693438: - # continue - ju = surveylist[i,0] - pointing = transRaDec2D((c_equtor[i].ra*u.degree).value, (c_equtor[i].dec*u.degree).value) + + + ju = surveylist[i,5] + pointing = transRaDec2D(surveylist[i,0],surveylist[i,1]) # print(ju, pointing, surveylist[i,3:9]) - sl = StrayLight(jtime = ju, sat = surveylist[i,3:6], radec = np.array([(c_equtor[i].ra*u.degree).value, (c_equtor[i].dec*u.degree).value])) - # strayl_i,s_zoldical ,s_earth, s_earth1 = sl.caculateStrayLightFilter(filter = 'i') - # print(i,strayl_i,s_zoldical,s_earth, s_earth1) - p_cart= transRaDec2D((c_equtor[i].ra*u.degree).value, (c_equtor[i].dec*u.degree).value) - sky_pix[i,6] = getAngle132(x1 = surveylist[i,6], y1 = surveylist[i,7], z1 = surveylist[i,8], x2 = p_cart[0], y2 = p_cart[1], z2 = p_cart[2], x3 = 0, y3 = 0, z3 = 0) + sl = StrayLight(jtime = ju, sat = surveylist[i,6:9], radec = np.array([surveylist[i,0],surveylist[i,1]])) - earthZenithAngle,isInSunSide = calculateAnglePwithEarth(sat = surveylist[i,3:6], pointing = pointing, sun = surveylist[i,6:9]) - sky_pix[i,4] = earthZenithAngle - sky_pix[i,5] = isInSunSide - e1_,py = sl.caculateEarthShineFilter(filter = 'i') - # e2, e2_ = sl.calculateZodiacalFilter1(filter = 'i') - e3, v_mag = sl.calculateZaodiacalFilter2(filter = 'i', sun_pos=surveylist[i,6:9]) - # e4 = sl.caculateStarLightFilter(filter = 'i',pointYaxis = py) - # e4 = 0 + e1,py = sl.caculateEarthShineFilter(filter = 'i') - e1,py = sl.caculateEarthshineGrating(grating = 'GI', pixel_size_phy = 10, normFilter = 'g') + e2, v_mag = sl.calculateZaodiacalFilter2(filter = 'i', sun_pos=surveylist[i,9:12]) - # e2 = sl.caculateStarLightGrating(grating = 'GV', pointYaxis = py) - e2 = sl.caculateStarLightGrating(grating = 'GI', pointYaxis = py) - e4 = sl.caculateStarLightFilter(filter = 'i',pointYaxis = py) + e3 = sl.caculateStarLightFilter(filter = 'i',pointYaxis = py) - e5=sl.caculateStrayLightFilter(filter = 'i', pixel_size_phy = 10, pixel_scale = 0.074, sun_pos = surveylist[i,6:9]) - e6,_=sl.caculateStrayLightGrating(grating = 'GI', normFilter_es = 'g', sun_pos = surveylist[i,6:9]) + print(i, e1, e2, e3) sky_pix[i,0] = e1 sky_pix[i,1] = e2 sky_pix[i,2] = e3 - sky_pix[i,3] = e4 - print(e1+e2,e1_+e3+e4,e5,e6) + sky_pix[i,3] = e1+e2+e3 - # print(e1,e2,e3,e4) + strRes = str(e1) + ' ' + str(e2) + ' ' + str(e3) + ' ' + str(e1+e2+e3) + '\n' +outF = open(fn[0:-4]+'stralight.dat',"w") +outF.write(strRes) +outF.flush() +outF.close() -- GitLab