在Python中使用GDAL为卫片制作缩略图

2014-11-24 13:12:27 · 作者: · 浏览: 2

卫片和航片文件一般都很大,


有时候需要制作一个缩略图,放在导航窗体中做底图


这个文件不用太清晰,只要反映整体面貌即可


在ArcMap,PCI和Erdas中,可以直接截屏,或者使用重投影等功能很方便地完成,


最近在用python和GDAL,于是乎,查了查资料,并不复杂,


首先要确定python安装了numpy库,代码



以上的这个方法,用到了numpy的矩阵运算,不过使用band的ReadAsArray


方法取得的数据,采用的是最邻近采样,虽然缩略图的要求不高,但是这个


采样方式总让人不爽,查了资料,使用gdal.ReprojectImage


可以设置采样方式,代码如下:



但这个生成的不是缩略图,而是从大图中裁剪出来的一部分