經度 phi，緯度 theta 處的坐標為：x =R* cos(phi) * cos(theta)y = Rsin(phi) * cos(theta)z = Rsin(theta)

經緯度采樣的采樣點是相同經緯度間隔的交點。但是采樣1000個點，如何劃分多少條經線，多少條緯線相交，才能使1000個采樣點最均勻的分布在球面上（雖然經緯度采樣本來就不均勻，但對于不同的采樣點個數應該有一種相對最均勻的經緯線劃分）？求大佬指教！我目前是將緯度每10度進行劃分。

import randomfrom mpl_toolkits import mplot3dimport numpy as npimport mathimport matplotlib.pyplot as plt%matplotlib inlineax=plt.axes(projection='3d')N=1000x=[]y=[]z=[]r=1#經度def longitude(lng): phi=(180+lng)*(math.pi/180) return phi#緯度def latitude(lat): theta=lat*(math.pi/180) return thetafor i in range(-80,90,10): for j in np.arange(-180,180,360/((N-2)/17)): #x.append(-r*math.sin(latitude(i))*math.cos(longitude(j))) #y.append(r*math.cos(latitude(i))) #z.append(r*math.sin(latitude(i))*math.sin(longitude(j)))x.append(r*math.cos(latitude(i))*math.cos(longitude(j))) z.append(r*math.sin(latitude(i))) y.append(r*math.cos(latitude(i))*math.sin(longitude(j)))x.append(r*math.cos(latitude(-90))*math.cos(longitude(0)))z.append(r*math.sin(latitude(-90)))y.append(r*math.cos(latitude(-90))*math.sin(longitude(0)))x.append(r*math.cos(latitude(90))*math.cos(longitude(0)))z.append(r*math.sin(latitude(90)))y.append(r*math.cos(latitude(90))*math.sin(longitude(0)))xline=np.array(x)yline=np.array(y)zline=np.array(z)print(xline.shape)ax.scatter3D(xline,yline,zline,s=2)plt.savefig('D:samples經緯度采樣.png')

效果

到此這篇關于python實現經緯度采樣的示例代碼的文章就介紹到這了,更多相關python 經緯度采樣內容請搜索好吧啦網以前的文章或繼續瀏覽下面的相關文章希望大家以后多多支持好吧啦網！