python GDAL的应用

最新推荐文章于 2024-07-09 09:42:52 发布

似曾相识_燕归来

最新推荐文章于 2024-07-09 09:42:52 发布

阅读量505

点赞数

分类专栏： python

本文链接：https://blog.csdn.net/u010624166/article/details/80140023

版权

这篇博客介绍了如何在Python中运用GDAL库来处理地理数据。通过遍历图像像素，计算每个像素的(X, Y)坐标，并存储在二维列表`arrSlope`中，展示了GDAL在地理信息系统的应用。" 80726504,7405643,二维几何变换与齐次坐标的原理,"['图形学', '线性代数', '几何变换', '计算机图形学']

摘要由CSDN通过智能技术生成

from osgeo import gdal
gdal.AllRegister()
filePath="D://data3//J47G085041.tif"

dataset=gdal.Open(filePath)
#if dataset is None:
# sys.exit(1)

#从一个栅格数据集中取出Geotransform作为一个list，然后读取其中的数据
adfGeoTransform=dataset.GetGeoTransform()
#读取栅格数据集的x方向像素数，y方向像素数，和波段数
cols=dataset.RasterXSize
rows=dataset.RasterYSize
bands=dataset.RasterCount
driver=dataset.GetDriver()#.LongName # 文件类型
h=dataset.RasterIO

#左上角地理坐标
originX=adfGeoTransform[0]
originY=adfGeoTransform[3]
pixelWidth=adfGeoTransform[1]
pixelHeight=adfGeoTransform[5]
#计算某一坐标对应像素的相对位置(pixel offset)，也就是该坐标与左上角的像素的相对位置，按像素数计算
X=readfile
xOffset=int((x-originX)/pixelWidth)
yOffset=int((Y-originY)/pixelHeight)
#首先读取一个波段