GDAL(python) 之GeoTransform

最新推荐文章于 2025-09-07 19:36:40 发布

原创

最新推荐文章于 2025-09-07 19:36:40 发布 · 1.4w 阅读

65 ·

CC 4.0 BY-SA版权

该文章已生成可运行项目，

演示

使用GDAL库读出的dataset带有两个重要的地理参数，分别是Projection和GeoTransform。有了这两个参数，就确定了影像的地理位置。

再GDAL for Python中，GeoTransform是一个六个元素的元组。
例如，我找了一个影像，读取并显示它的GeoTransform，则为如下形式：
在这里插入图片描述

形式

(486892.5, 15.0, 0.0, 4105507.5, 0.0, -15.0)

六个参数分别为：

左上角x坐标，水平分辨率，旋转参数，左上角y坐标，旋转参数，竖直分辨率。如另一篇博客中所说，满足如下关系式：
在这里插入图片描述
一般来说，旋转参数都为0。

实践

有了上述知识，不妨利用一张图片进行实际的验证。

总体的思路是，首先读取我下载的某Landsat8数据集的全色波段，然后截取一个2000 * 2000的部分。

然后我们手动的把图片分成上下两部分，然后给上下两部分按照他们的情况，定义新的GeoTransform，再将图片写出，看看我们定义的GeoTransform是否正确，以此验证我们的知识。

下面是代码验证：

import numpy as np
from osgeo import gdal, gdal_array
import os

filename = "F:\SRGAN_program\dataset\LC81290352019095LGN00\LC08_L1TP_129035_20190405_20190422_01_T1_B8.tif"
dataset = gdal.Open(filename, gdal.GA_ReadOnly)
if dataset == None:
    raise Exception("Image name error.")
else:
    datatype = np.float16
    height = dataset.RasterYSize
    width = dataset.RasterXSize
    projection = dataset.GetProjection()