第七色在线视频,2021少妇久久久久久久久久,亚洲欧洲精品成人久久av18,亚洲国产精品特色大片观看完整版,孙宇晨将参加特朗普的晚宴

為了賬號(hào)安全,請(qǐng)及時(shí)綁定郵箱和手機(jī)立即綁定
已解決430363個(gè)問(wèn)題,去搜搜看,總會(huì)有你想問(wèn)的

繪制具有周期性邊界且無(wú)線段的尼斯曲線

繪制具有周期性邊界且無(wú)線段的尼斯曲線

尚方寶劍之說(shuō) 2021-10-10 13:46:18
使用庫(kù)pyephem,我想找到一種方法來(lái)繪制表示衛(wèi)星地面軌跡的漂亮繪圖(經(jīng)度、緯度)曲線。我已經(jīng)計(jì)算了 (longitude,latitude) 但是當(dāng)經(jīng)度超過(guò) +180° 時(shí),計(jì)算的下一個(gè)值是例如 -178°,這樣就繪制了一個(gè)線段:這就是糟糕的數(shù)字。在我目前擁有的代碼和獲得的數(shù)字下方。currentDate = date.datetime(2018,12,1,0,0,0);for i in range(nPoints):  iss.compute(currentDate)  # compute latitude  posLat[i] = iss.sublat*(180/math.pi)  # compute longitude  posLong[i] = iss.sublong*(180/math.pi)  currentDate = currentDate + date.timedelta(seconds = (time3orbits/nPoints*3600))                  print currentDate, posLong[i], posLat[i]plt.plot(posLong,posLat)plt.show()這是我得到的圖(繞地球 3 圈):如何防止繪制的線段?問(wèn)候iss.sublat和iss.sublong分別是 iss 位置的當(dāng)前緯度和經(jīng)度(從一個(gè)currentDate固定的)。
查看完整描述

1 回答

?
牛魔王的故事

TA貢獻(xiàn)1830條經(jīng)驗(yàn) 獲得超3個(gè)贊

您可以將符號(hào)從正變?yōu)樨?fù)的數(shù)據(jù)數(shù)組拆分,然后繪制單個(gè)子數(shù)組:


orbits = np.argwhere(

    (data[:-1, 0] * data[1:, 0] < 0)

    & (data[1:, 0] < 0)

).ravel() + 1


for orbit in np.split(data, orbits, axis=0):

    plt.plot(orbit[:, 0], orbit[:, 1])

http://img1.sycdn.imooc.com//61627e4e0001d4ce06140432.jpg

查看完整回答
反對(duì) 回復(fù) 2021-10-10
  • 1 回答
  • 0 關(guān)注
  • 260 瀏覽
慕課專欄
更多

添加回答

舉報(bào)

0/150
提交
取消
微信客服

購(gòu)課補(bǔ)貼
聯(lián)系客服咨詢優(yōu)惠詳情

幫助反饋 APP下載

慕課網(wǎng)APP
您的移動(dòng)學(xué)習(xí)伙伴

公眾號(hào)

掃描二維碼
關(guān)注慕課網(wǎng)微信公眾號(hào)