使用python gdal模块打开相同投影但不同尺寸的图像
我不确定这是正确的站点,但是
我正在使用 python 2.7 上的 gdal.Open() 模块来打开 Landsat GeoTIFF 图像。它们具有相同的 UTM 地图投影,但图像大小不同。
如何将图像尺寸固定为相同?我想进行直接逐像素比较。
I'm not sure this is the correct site for this but,
I'm using the gdal.Open() module on python 2.7 to open Landsat GeoTIFF images. They have the same UTM map projection, but different image sizes.
How do I fix the image sizes to be the same? I'd like to make direct pixel-by-pixel comparisons.
如果你对这篇内容有疑问,欢迎到本站社区发帖提问 参与讨论,获取更多帮助,或者扫码二维码加入 Web 技术交流群。
绑定邮箱获取回复消息
由于您还没有绑定你的真实邮箱,如果其他用户或者作者回复了您的评论,将不能在第一时间通知您!
发布评论
评论(2)
gdal_translate 可以简单地使用 -outsize 参数来更改图像的大小,该参数分别采用两个整数值作为 xsize 和 ysize,或两个百分比值来缩放图像。 参见此处
gdal_translate can be used simply to change the size of an image using the -outsize parameter, which takes two integer values as the xsize and ysize respectively, or two percentage values to scale the image. see here
您需要将所有栅格重新投影为相同的栅格形状/大小。您可以使用 gdalwarp 对每个光栅文件执行此操作,例如:
您可以选择替换
- tr
(目标分辨率)与-ts width height
指定行/列数。-te
(目标范围)和-tr
/-ts
(目标分辨率/大小)选项强制每个源栅格具有相同的光栅形状叠加,因此您可以正确地进行逐像素比较。您还需要选择合适的-r
重采样方法,具体取决于您的分析试图实现什么目标。You need to reproject all rasters to the same raster shape/size. You can do this with gdalwarp to each raster file, e.g.:
You can optionally replace
-tr
(target resolution) with-ts width height
to specify the number of rows/columns.The
-te
(target extents) and-tr
/-ts
(target resolution/size) options forces each of your source raster to have the same raster shape overlay, so you can correctly do a pixel-by-pixel comparison. You also need to choose a suitable-r
resampling method, depending on what your analysis is trying to achieve.