我知道之前曾有人問過這個問題,但是重點(diǎn)是Cartopy。我目前的工作是在日期線上觀察一個現(xiàn)象,并將在穿越日期線的太平洋中的一個盒子上工作。管理和繪制此數(shù)據(jù)一直是一個挑戰(zhàn),因?yàn)橐粋€索引的值從179變?yōu)?179。我可以輕松地將自己的緯度歸為一類lat_mask=[(lat<=20)&(lat>=-20)] #Tropics +-20* onlylat_masked=lat[lat_mask]哪個工作正常。lon=np.add(lon,180) #Reset longitudes to 0-360 from -180 to 180lon_mask=[(lon>=110)&(lon<=300)] 也可以正常工作...但是我正在尋找的完全相反的區(qū)域。但是我確實(shí)想要的區(qū)域:lon_mask=[(lon<=110)|(lon>=300)] 不起作用。 我可以分別做這些事,并為日期線的東部和西部的每個區(qū)域制作圖形,但我希望將它們組合為一個數(shù)據(jù)集/圖形。這樣我就可以為[300-360 0-110]設(shè)置一個蒙版,可以將其繪制為單個圖像。我正在尋找一種存儲信息的簡便方法,以提高可用性和可繪制性。我一直在嘗試尋找如何重新索引所有索引的方法,但這使我的頭部受傷。我正在尋找一種可靠的解決方案,因?yàn)樗鼘⒔?jīng)常使用。我本來以為這是matplotlib問題,但我認(rèn)為這是存儲數(shù)據(jù)的方式。編輯:更改為| -但是現(xiàn)在出現(xiàn)了這個問題。我感到lon= np.roll(lon,len(lon)/2)或可能有助于有效地將日期變更表移至其他地方,但我無法完全解決。
在日期線上管理和繪制數(shù)據(jù)
ibeautiful
2021-04-02 10:11:39