学习笔记——【python】GetGeoTransform()使用，gdal截取图像，使用GDAL进行影像投影坐标、地理坐标、图上坐标的转换-CSDN博客

本文链接：https://blog.csdn.net/weixin_43570470/article/details/121555964

本文介绍了如何使用Python的GDAL库进行图像处理，包括通过GetGeoTransform()获取图像坐标信息，利用ReadAsArray()进行图像截取，以及进行投影坐标、地理坐标和图上坐标的转换。详细步骤和实例演示帮助理解GDAL在地图处理中的应用。

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

1. GetGeoTransform()使用、gdal截取图像

GetGeoTransform()
GeoTransform[0]，左上角横坐标（应该是投影坐标）
GeoTransform[2]，行旋转
GeoTransform[1]，像元宽度(影像在水平空间的分辨率)
GeoTransform[3]，左上角纵坐标（应该是投影坐标）
GeoTransform[4]，列旋转
GeoTransform[5]，像元高度(影像在垂直空间的分辨率)，
如果影像是指北的，GeoTransform[2]和GeoTransform[4]这两个参数的值为0，GeoTransform[5]为负；
如果图像不含地理坐标信息，默认返回值是：(0,1,0,0,0,1)

ReadAsArray(xoff,yoff,xsize,ysize)
#xoff,yoff是取值窗口的左上角在实际数据中所处象元的xy位置。
#xsize,ysize是取值窗口覆盖的区域大小

import gdal  #导入库

dataset = gdal.open("filename")  #文件名，如*.tif
geotransform = dataset.GetGeoTransform()
print('数据投影：')
print(dataset.GetProjection())
width = int(dataset.RasterXSize) #读取图像的宽度，x方向上的像素个数
height = int(dataset.RasterYSize)  #读取图像的高度，y方向上的像素个数
RasterCount = int(dataset.RasterCount)
print('数据的大小（行/height，列/weight, RasterCount）：')
print('(%s %s %s)' % (height, width, RasterCount))

if geotransform[0] == 0.0 and geotransform[1] == 1.0:
	print("EMPTY IMAGE!")
img_tu = dataset.ReadAsArray(xoff=xp_1