python语言GDAL读取所在点的像素值

由于项目需求，需要用到开源的软件开发，就开始研究GDAL，这个开源库，目前来说，调用GDAL库，最方便的语言还是python ，简单记录下，用python语言做GDAL开发的一些东西吧。一个在开发中经常遇见的问题，就是一个点数据所对应的像素值，这是很常见的，上代码：

#下面是GDAL打开shp数据，不细说了。

driver = ogr.GetDriverByName('ESRI Shapefile')

ds = driver.Open(sub_dir, 0)

sublayer = ds.GetLayer(0)

# 存放每个点的xy坐标的数组

xValues = []

yValues = []

#从这里开始

# 打开一个栅格数据Tif

dr = gdal.Open(Tif)

#存储着栅格数据集的地理坐标信息

transform = dr.GetGeoTransform()

# 影像左上角横坐标

x_origin = transform[0]

# 影像左上角纵坐标

y_origin = transform[3]

# 遥感图像的水平空间分辨率或者东西方向上的像素分辨率

pixel_width = transform[1]

# 遥感图像的垂直空间分辨率或者南北方向上的像素分辨率

pixel_height = transform[5]

x_offset = int((x - x_origin) / pixel_width)

y_offset = int((y - y_origin) / pixel_height)

pixelValue= dr.ReadAsArray(x_offset, y_offset,1, 1)

上面的代码就可以读取点对应的像素值pixelValue，重点要说，在GDAL中，影像都是左上角为起点，就是（0,0）点，看图：

备注：原创

开始时我也不理解，计算x_offset为啥要用x- x_origin,看见这个图就明白了。

再说下这个函数，ReadAsArray，这个函数就可以读取像素值，当然了，这个函数有很多的重载，感兴趣的可以再深入研究下。