이 Stereographic 일도하기 matplotlib 극성 플롯 파이썬

Question

나는 다음과 비슷한 간단한 입체 태양 경로 다이어그램을 만들려고하고 있어요 y 축을 뒤집는 것에 대해? 현재 축이 0> 90으로 이동합니다. 방위각을 나타내려면 축을 90> 0으로 어떻게 바꿉니 까?

나는 시도했다 :

또한

ax.invert_yaxis() 
ax.yaxis_inverted() 

, 나는 등거리에 반대로 입체 프로젝션을 만드는 방법에 대해 갈 것 어떻게?

내 코드 :

import matplotlib.pylab as plt 
testFig = plt.figure(1, figsize=(8,8)) 
rect = [0.1,0.1,0.8,0.8] 
testAx = testFig.add_axes(rect,polar=True) 
testAx.invert_yaxis() 
testAx.set_theta_zero_location('N') 
testAx.set_theta_direction(-1) 

Azi = [90,180,270] 
Alt= [0,42,0] 
testAx.plot(Azi,Alt) 
plt.show() 

가 현재 내 코드도 정확하게 라인을 그릴하지 않는 것, 내가 뭔가에 각 또는 학위를 변환 할 필요가 필요합니까?

도움을 주시면 대단히 감사하겠습니다.

Answer 1

마침내 matplotlib로 놀 시간이있었습니다. 많은 검색을 한 후에 Joe Kington이 지적한 정확한 방법은 축을 서브 클래스 화하는 것입니다. 우수한 basemap 모듈을 사용하여 훨씬 빠른 방법을 발견했습니다.

아래는 내가 stackoverflow에 맞게 조정 한 코드입니다. 태양 고도와 방위각은 판다에서 생성 된 일련의 우표가있는 Pysolar로 계산되었습니다. 현재 난 단지 포인트를하려하고, 당신이 선 포인트 일출/일몰 고도 0 점을 가지고 있는지 확인해야합니다

import matplotlib.pylab as plt 
from mpl_toolkits.basemap import Basemap 
import numpy as np 

winterAzi = datafomPySolarAzi 
winterAlt = datafromPySolarAlt 

# create instance of basemap, note we want a south polar projection to 90 = E 
myMap = Basemap(projection='spstere',boundinglat=0,lon_0=180,resolution='l',round=True,suppress_ticks=True) 
# set the grid up 
gridX,gridY = 10.0,15.0 
parallelGrid = np.arange(-90.0,90.0,gridX) 
meridianGrid = np.arange(-180.0,180.0,gridY) 

# draw parallel and meridian grid, not labels are off. We have to manually create these. 
myMap.drawparallels(parallelGrid,labels=[False,False,False,False]) 
myMap.drawmeridians(meridianGrid,labels=[False,False,False,False],labelstyle='+/-',fmt='%i') 

# we have to send our values through basemap to convert coordinates, note -winterAlt 
winterX,winterY = myMap(winterAzi,-winterAlt) 

# plot azimuth labels, with a North label. 
ax = plt.gca() 
ax.text(0.5,1.025,'N',transform=ax.transAxes,horizontalalignment='center',verticalalignment='bottom',size=25) 
for para in np.arange(gridY,360,gridY): 
    x= (1.1*0.5*np.sin(np.deg2rad(para)))+0.5 
    y= (1.1*0.5*np.cos(np.deg2rad(para)))+0.5 
    ax.text(x,y,u'%i\N{DEGREE SIGN}'%para,transform=ax.transAxes,horizontalalignment='center',verticalalignment='center') 


# plot the winter values 
myMap.plot(winterX,winterY ,'bo') 

참고.