2017-09-18 29 views

答えて

0

私はあなたのプロットを生成するために使用しているコードなしでこれをテストすることはできませんが、あなたが使用している場合、それはおそらく

plt.xlim(112.4, 114.3) 
plt.ylim(35.8, 37.7) 

のようなものですaxインターフェイスではなく、

ax.set_xlim(112.4, 114.3) 
ax.set_ylim(35.8, 37.7) 
+0

fr matplotlibのインポートCMからPLT としてNPインポートmatplotlib.pyplotとしてOM OSGEOインポートGDAL インポートnumpyの gdal.UseExceptions() DS = gdal.Open( 'D:\ PyGDAL_Case \ DEM_ZH \ Dem_Zh.tif') バンド= ds.GetRasterBand(1) elev1 = band.ReadAsArray() NROWS、ncolsの= elev1.shape YY = np.where(elev1> = 0)[0] XX = np.where(elev1> = 0 )[1] ZZ = elev1 [np.where(elev1> = 0)] elev_map = np.nan * np.empty((NROWS、ncolsの)) elev_map [YY、XX] ZZ X0、DXを= 、dxdy、y0、dydx、dy = ds.GetGeoTransform() x 1 = X0 + DX * ncolsの Y1 = Y0 + DY * NROWS plt.imshow(elev_map、CMAP = 'gist_earth_r'、程度は= [X0、X1、Y1、Y0]) plt.show() –

+0

以上でありますコード(私はapproriateスタイルを編集する方法を知らない)。優しい指示のためにジェレミーに感謝します。やってみます。 –

関連する問題