2016-10-12 3 views
0

pcolormesh 지도 위에. Python 2.7 용 코드 :홀수 파이썬 기본지도와 동작 및 I는 좌표가 X, Y는, 기본지도()를 이용하여 좌표를 매핑하도록 변환하는 numpy.meshgrid 제조하고 적용되는 pcolormesh함으로써 비닝 데이터 길이 스트립 플롯려고

import numpy as np 
from mpl_toolkits.basemap import Basemap, shiftgrid, addcyclic 
import matplotlib.pyplot as plt 

lon_tics = np.linspace(0, 360.0, 60) 
lat_tics = np.linspace(-90.0, 90.0, 30) 
map_bins = np.zeros((60,30), dtype = np.int) 
bin15 = np.random.randint(0,20,30) #we should see 2 strips 
bin45 = np.random.randint(0,20,30) #but we get lots of strange results 
map_bins[15] = bin15 
map_bins[45] = bin45 
m = Basemap(projection='moll',lon_0= -120,resolution='c') #NOTE changing lon_0 has weird results! 
lon_bins_2d, lat_bins_2d = np.meshgrid(lon_tics, lat_tics) 
xs, ys = m(lon_bins_2d, lat_bins_2d) 
plt.pcolormesh(xs, ys, np.transpose(map_bins)) 
plt.colorbar() 
m.drawparallels(np.arange(-90.,120.,30.), labels = [True]) 
m.drawmeridians(np.arange(0.,360.,60.), labels = [False]) 
plt.show() 

이것은 매우 이상한 행동입니다. 변경하여 중 무엇에, lon_0는베이스 맵() 인스턴스, 또는 설정되는 경우 데이터를 비닝 어떻게 위도/경도 쓰레기통있는 등, 우리가 행동을 다른 얻을 것이다 정의 :

  • 없음 종 방향 스트립
  • 1 종 스트립
  • 2 종 방향 스트립 (예상되는 동작) (사진 : lon_0 = -120 만든)
  • 맵의 가장자리 빈 '바탕색'(사진 : lon_0 만든 = 98)

나는이 문제를 해결하려고 노력 해왔다. 얼마 동안은 아프다. 아무도 내가 뭘 잘못 볼 수 있습니까?

읽어 주셔서 감사합니다.

Correct output

Smear

답변

0

일부를 둘러보고 내가이 글을 우연히 발견 추가 후 : Map projection and forced interpolation

이 지금은 내 문제를 해결했다, 그래서 나는 경우 다른 사람에 연결하고있어이 문제로 실행 .