利用python读取点矢量对应栅格值

 

每行代码都有详细注解

所需库
# GDAL是栅格和矢量地理空间数据格式的转换器库
# 旧版本加载库的方法
#import gdal,gdalconst
# from 语句让你从模块中导入一个指定的部分到当前命名空间中
# GDAL数据驱动，与OGR数据驱动类似 gdal是读取栅格数据 ogr是读取矢量数据
# 查看支持数据格式https://gdal.org/drivers/raster/index.html
# import ogr或者
# from osgeo import ogr或者
try:
    from osgeo import ogr
except:
    import ogr
#gdalconst中的常量都加了前缀，力图与其他的module冲突最小
from osgeo import gdal,gdalconst
from osgeo import gdal_array as ga  # 用于引入一个模块的同时为该模块取一个别名
from osgeo.gdalconst import GA_ReadOnly


# shp矢量的绝对路径 栅格的绝对路径
def getPixels(shp, img):
    # 要读取某种类型的数据，必须要先载入数据驱动，也就是初始化一个对象，让它“知道”某种数据结构。
    # GDAL数据驱动，与OGR数据驱动类似，需要先创建某一类型的数据驱动，再创建响应的栅格数据集。
    # 一次性注册所有的数据驱动，但是只能读不能写：gdal.AllRegister()
    # 单独注册某一类型的数据驱动，这样的话可以读也可以写，可以新建数据集：
    driver = ogr.GetDriverByName(&#