From b1ffd37980dc46a04912931401755c2f730cd9d7 Mon Sep 17 00:00:00 2001 From: march3 Date: Thu, 29 Jun 2023 14:42:26 +0800 Subject: [PATCH] =?UTF-8?q?Python=E8=B6=85=E4=BA=BA-=E5=AE=87=E5=AE=99?= =?UTF-8?q?=E6=A8=A1=E6=8B=9F=E5=99=A8?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- .../solar_system_habitable_zone.py | 67 +++++++++++++++--- textures/green_alpha_1.png | Bin 0 -> 2943 bytes 2 files changed, 57 insertions(+), 10 deletions(-) create mode 100644 textures/green_alpha_1.png diff --git a/sim_scenes/solar_system/solar_system_habitable_zone.py b/sim_scenes/solar_system/solar_system_habitable_zone.py index 3382186..0933701 100644 --- a/sim_scenes/solar_system/solar_system_habitable_zone.py +++ b/sim_scenes/solar_system/solar_system_habitable_zone.py @@ -6,29 +6,76 @@ # link :https://gitcode.net/pythoncr/ # python_version :3.8 # ============================================================================== -from bodies import Sun, Mercury, Venus, Earth, Moon, Mars, Jupiter, Saturn, Uranus, Neptune, Pluto, HabitableZone, Asteroids +from bodies import Sun, Mercury, Venus, Earth, Moon, Mars, Jupiter, Saturn, Uranus, Neptune, Pluto, HabitableZone, \ + Asteroids from common.consts import SECONDS_PER_WEEK, SECONDS_PER_DAY, SECONDS_PER_YEAR, AU from sim_scenes.func import mayavi_run, ursina_run +from simulators.ursina.entities.body_timer import TimeData +from simulators.ursina.ursina_config import UrsinaConfig +from simulators.ursina.ursina_event import UrsinaEvent +from ursina import Vec3 +import numpy as np if __name__ == '__main__': # 目前认为 太阳系 的宜居带范围是从距离太阳0.95个天文单位 (约1.42亿千米)到 2.4个天文单位（约3.59亿千米）的范围为宜居带， # 其宽度约为2.17亿千米， 按照这个标准，太阳系的宜居带中只有三个大型天体，分别是地球、 月球 以及火星（1.52天文单位）。 - sun = Sun(name="太阳", size_scale=0.5e2) # 太阳放大 80 倍，距离保持不变 + sun = Sun(name="太阳", size_scale=0.5e2) # 太阳放大 80 倍，距离保持不变 + earth = Earth(name="地球", size_scale=1.5e3) # 地球放大 1500 倍，距离保持不变 + moon_d = 10000000 + moon = Moon(name="月球", size_scale=1e2, + init_position=[moon_d, 0, 0], + init_velocity=[0, 0, 0], + distance_scale=0.2, + gravity_only_for_earth=True + ) bodies = [ sun, - Venus(name="金星", size_scale=1.5e3), # 金星放大 4000 倍，距离保持不变 - Earth(name="地球", size_scale=1.5e3), # 地球放大 4000 倍，距离保持不变 - Moon(name="月球", size_scale=2e3), # 地球放大 4000 倍，距离保持不变 - Mars(name="火星", size_scale=2e3), # 火星放大 4000 倍，距离保持不变 + Venus(name="金星", size_scale=1.5e3), # 金星放大 1500 倍，距离保持不变 + earth, moon, + # Moon(name="月球", size_scale=2e3, + # # init_position=[0, 0, 363104 + 1.12 * AU], + # # init_velocity=[-(29.79 + 1.03), 0, 0] + # init_position=[15000000, 0, AU], + # init_velocity=[-32.79, 0, 2.03], gravity_only_for_earth=True + # ), # 月球放大 2000 倍，距离保持不变 + Mars(name="火星", size_scale=2e3), # 火星放大 2000 倍，距离保持不变 Asteroids(name="小行星群", size_scale=3.2e2, - parent=sun), # 小行星群模拟(仅 ursina 模拟器支持) + parent=sun), # 小行星群模拟(仅 ursina 模拟器支持) HabitableZone(name="宜居带", size_scale=1e2, - parent=sun), # 小行星群模拟(仅 ursina 模拟器支持) - Jupiter(name="木星", size_scale=2e2), # 木星放大 800 倍，距离保持不变 + parent=sun), # 宜居带模拟(仅 ursina 模拟器支持) + Jupiter(name="木星", size_scale=2e2), # 木星放大 200 倍，距离保持不变 ] + + def on_ready(): + # 月球就会跟着地球自转而转 + moon.planet.parent = earth.planet + + + def on_timer_changed(time_data: TimeData): + # print(moon.planet.position) + # moon.planet.rotation = -Vec3(earth.planet.rotation_x, # - earth.planet.ring_rotation_x, + # earth.planet.rotation_y, + # earth.planet.rotation_z) + # 为了不让月球随着地球的周期旋转，则需要获取地球的旋转角度 + angle = earth.planet.rotation_y + # TODO:根据旋转的角度对月球的位置进行计算，保证月球公转和地球自转的关系 365天=12月 + # angle = np.array(angle * np.pi) + # + # px = moon_d * UrsinaConfig.SCALE_FACTOR * np.cos(angle) + # pz = moon_d * UrsinaConfig.SCALE_FACTOR * np.sin(angle) + # moon.planet.world_position = Vec3(0,0,0) + + + # 订阅事件后，上面2个函数功能才会起作用 + # 运行前会触发 on_ready + UrsinaEvent.on_ready_subscription(on_ready) + # 运行中，每时每刻都会触发 on_timer_changed + UrsinaEvent.on_timer_changed_subscription(on_timer_changed) + # 使用 ursina 查看的运行效果 # 常用快捷键： P：运行和暂停 O：重新开始 I：显示天体轨迹 # position = 左-右+、上+下-、前+后- - ursina_run(bodies, 1, position=(0, 2 * AU, -11 * AU), + ursina_run(bodies, SECONDS_PER_DAY, position=(0, 2 * AU, -11 * AU), + show_timer=True, bg_music="sounds/interstellar.mp3") diff --git a/textures/green_alpha_1.png b/textures/green_alpha_1.png new file mode 100644 index 0000000000000000000000000000000000000000..2a064f93fae150722c0057191ff9aa0acb04c14b GIT binary patch literal 2943 zcmchZ^;Z*o8^*sQbby4^2x&njq+R&~;Rry`Xr6JP zMJk0;b&}%9Y^FJ~UX$q&34}#(&@L9sB9<4x?g<5~Q@1k?LUYXakC@56%0tE&e;4jLD87Z3Z<1I!WUS;bNoirm0;WXw(F5rJd_ zcw2W<1QiE-ibjxSf!OLy+YZ1wH99C9RjqaidFu_RfX~p?i>#f~C(#XtV&*_F1xyLE zE;2wcNYHiFq>CMNae->%cMs`78UP;XP)%V_!3b8ibyaD>2sU*b21bf5V4)xp0IWeRKK6N86C0D+j^mw9(ko*ke5l_zj>eQHK%rq!cH(++lp#fffCI5#s# zgblt(jp#oI0M8!j^Rbc|W^rL|dckSYgS_2dbnL$FEvzcEF?rnhoQejp=ns%bZRY2P zseL*j)&4E+SSv7U23#ihtmVh(3hvQ=4c<#YD;=uy=NB4w-Q>GK&%m%MYG2qSGn=4# zB0udtl@Pvr+8&PuuY*U~>RFY9;^?DCc3uHh$fQLgji#&@Wl`j~Gv@J8(iPH-Ld|_i zj%OPdPw|GrdY$Nx$dp?`-RU!|OONF?1Tv23QUe(QF4lmc|KaR=q6Z;-4MLT`WHA^3 z)Wf%=p7SuI02XHjIWz@;)Rw(ubqNUQ%Pg4%pz$}?{qi&zq5}%Rt%6ALYK^n!KC|Np zOu3)=8VN8eyEC`c&s`+kWYIh`7jc2UkSjxj+k~JzLi4_wny27=_B2b~XIUl2#VOYN zk;^`?Gqzm$tz3<Yt+2ONB9oe+VS|jT3RQ~E8%LRCzbrkrKPzub}5gA49 zB%V5cm52{TWu&{m{ZSz7S&~BpNtVg!`A@t1&*QXhpIY#Ub#gk)TPT_(D?0F5{1d>z zOUSn)LoZS#M6j(uA3R}c^@ze?bSw}STIHfEB6X|{P;O_NC_u7K& zm1{5-w{c=?Li|Gcq7N=HvDu`;FEVw%nA4{h62-{rGU&R)$KX;W#wDtQEUvi-v1N>B z_pdoEGj4^m|(EVkuo?W>;eg;8TbP-vS#(GE5QEi^)^(#5etbz)oK5YJX zXhpp8t9!bh1p$WA#&K9S1^;Z(z8p?1Vu@7o(AR8NXIF}g^b)FpZX%jrJM(Mip!@H$ zON&d<+kxj8V;Q46;r5rg-MK@!%ed;frLz#FwBItbx=dt#$nj=5zhr*tZE7LDD;~;` z(!-jroNksLp1xRNV(MqAHiWKVFnMFj{N{mKSm|xkcZRcAy6>EChOm65pRm`e`jLUi zWWxa?o8o4hqo{oc#@g7J+Xb=kViSk@4(>!D&JH=p4GYa-Uo+}uN%xm}_ExofD~JHv z1Kt4X5QMRoUxgFSDa9#o;ZYxisiS6U{rb<`JN-iajQt@})KY3vHN$R|PL<)6|BMEY zHjmn``mS=WE{~G^ee%@vB=dHZJ(Vf`bPd-Uj2dDZ#{38VsD5);MlaFr^z1M#O^|eU zJhzK?V1LNNWRk!?G78XZ7;X<)C5QE1jKgpDm#%oy`jhGM^F)`1ls1=Ql5X)Oyx3Zt z^Q>YJOw|eJeZd~5@XK-B{idD0tK3*gTS@y|^Oe3d_iU#IQr$OHZ*Q06IBo>jy^?f& zDAzsTuK8!<*kNMKYwg_jI|9 zLtgWob2-X6#1G|_IhAbzIi8d3GwjJ%2CwWGdyUGj9Inu>WUriAy}oO=E4u5sJHXd_ zft7Fg2I@wrFUr3+z{96;xnx@53$_jWWy#!`*SXLyZ%YwYdKW`v%?!wd8GDvJEAx64 zOiZ8WS!)=4|81u@=TpvjE>|pnSjMw!tj%hwi5Oz5tg7O;-@iny+M!DID)?1gb7`~iq)M|$ z_;7eu_}Rn!gB$zAqvXFW)96jk4SiZ$$Y;nrq$!do@*TB1wE#^P&Fiy@w28Ez&PbnA z<!hM<|NqFmBHBz_}@DBN_Yrz^J_28Xj|__YtG$uQiekG+xB7IhZKO@*>8u~sqDQuT?p(G>2Fjsoe!S=-Y++o4Qqd96zUc{)HR(`RJ$A^|aiy++rl7Yet>K51IPQ z^#)E75~+-H=mpOn-qg4r%ACla-_{t1M+|~gT!*YOWa30fmrDd$aajf7;Z5*+-{mKr zj{oe5Fhp^NHesp_@zw7}@IjUvQ^GCUpEyMljA4OSWv!774})CS2g8Ui#M+D_!x;;t z>82Zcd@~6dJ^1>ymPC}fLOM3Bo9I=weea_Khh6$ZQ?gs;kd@hEaag(y&Ieas^|$&q zeqqL3-)1Ou2kkjJ`s+@8&+0OKxn$W3N3m>qd#A}wW9Vz_I%6g8M&dlLjpv((PvwjVm~HYf45SlN<}B!=V%)K*2;U? zkwRQGzgdOa)?2I{k87D-Ib z6A<8Kc@s#q