gdal对栅格数据进行重投影以及重采样

最新推荐文章于 2024-04-19 03:48:37 发布

陌云漫步

最新推荐文章于 2024-04-19 03:48:37 发布

阅读量3.3k

点赞数 3

CC 4.0 BY-SA版权

分类专栏：地理文章标签： python gdal

本文链接：https://blog.youkuaiyun.com/qq_39085138/article/details/111713295

地理专栏收录该内容

8 篇文章

订阅专栏

因服务器无arcpy，作者用GDAL进行Landsat数据重投影和重采样。目的是将数据从投影坐标转地理坐标并采样到0.00025度，借助GDAL的Warp函数，还介绍了目标坐标系EPSG及targetAlignedPixels选项的作用。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

最近要用服务器处理数据，但是服务器上没有arcpy，所以重投影以及重采样的任务只能交给gdal来做，目前网上的好多代码都是要有一个参考的栅格，但是像arcgis里往往不需要这么多参数

这里记录一下我的代码：

首先，我的目的是将landsat数据从投影坐标转成地理坐标，然后重采样到0.00025度。

GDAL库的Warp函数可以做到这一点，建议大家在使用的时候去仔细看看下面的官方文档（WarpOptions），看了之后才发现gdal真的可以完全取代arcpy（菜鸡发言）

https://gdal.org/python/osgeo.gdal-module.html#Warp

代码如下：

import gdal
pathin = ''
pathout = ''
gdal.Warp(pathout,pathin,dstSRS ='EPSG:4326',xRes=0.00025,yRes=0.00025,targetAlignedPixels=True)

代码其实就一行，解释一下dstSRS是destination Spatial Reference System（目标空间参考系统），就是要转成的目标坐标系。这里采用的是European Petroleum Survey Group [EPSG]（欧洲石油勘探组织），这个组织发布了全球不同坐标系的基准文件，只要我们输入EPSG：num（范围1024-32767），计算机就可以识别这个坐标系，下面是Wikipedia的解释：

EPSG Geodetic Parameter Dataset (also EPSG registry) is a public registry of geodetic datums, spatial reference systems, Earth ellipsoids, coordinate transformations and related units of measurement. Each entity is assigned an EPSG code between 1024-32767 along with a standard machine-readable well-known text (WKT) representation. （竟然是一个石油组织发布的）

EPSG：4326就是现在最常见的WGS84。

另外Warp还提供了一个非常好用的选项，targetAlignedPixels：whether to force output bounds to be multiple of output resolution，是否强制让输出的坐标的边界是分辨率的整数倍，这个功能真的太好用了，因为很多时候如果投影坐标转成地理坐标，不加上targetAlignedPixels=True的话出来的坐标都是无限不循环小数，如果后期想编程处理的话就没有一个很好的坐标基准。