首页 > 解决方案 > 停止重新采样栅格数据的栅格

问题描述

我正在尝试使用 rasterio 绘制栅格,但我认为数据正在以某种方式重新采样。

创建的地图未显示数据中的详细信息。见这里: 蟒蛇情节

...与原始数据(用 GIS 绘制)相比: 地形数据

知道如何阻止 rasterio 对数据进行插值/重采样吗?

这是我的代码:

import rasterio 
from rasterio.plot import show
import matplotlib.pyplot as plt

topo = rasterio.open('../topo_raster/nz100dem2ihs21.tif')
nz_mask_file  = '../gis_data/NZ_mask.shp'
nzmask  = gpd.read_file(nz_mask_file)

fig= plt.figure(dpi=300. )
ax = plt.axes( projection=ccrs.epsg(27200))

show(topo.read(1,masked=True), ax=ax,transform=topo.transform,interpolation='none',zorder=1,cmap='gist_gray')
nzmask.plot(ax=ax,facecolor="white", edgecolor='black', lw=linewidths_rr,zorder=0.5)

plt.savefig('../newzealand.png')
plt.clf()

标签: pythonrasterio

解决方案


推荐阅读