經度 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 經緯度采樣內容請搜索好吧啦網以前的文章或繼續(xù)瀏覽下面的相關文章希望大家以后多多支持好吧啦網！