GDAL(python) 之GeoTransform

该文章已生成可运行项目,
演示

使用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()
 
本文章已经生成可运行项目
评论 4
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值