【python+GIS】python实现多个栅格值提取到点，写入EXCEL表

原创

已于 2023-01-27 16:50:29 修改 · 3.7k 阅读

66 ·

CC 4.0 BY-SA版权

文章标签：

#python

于 2022-10-02 19:36:14 首次发布

文章目录

一、应用背景
二、参考文献
三、代码
四、结果

一、应用背景

下载了FLUXNET的站点，也下载了地表温度数据，想要将所有的温度数据（tif）的值提取到点上。

二、参考文献

《python快速实现shp点提取栅格值 (两种方法对比)》
《Python+GDAL | 基于矢量（点）读取栅格数据》

三、代码

from osgeo import ogr
import os, sys
from osgeo import gdal
from osgeo.gdalconst import *
import csv
import xlwt 
inputSHP=r'D:\DATA\pointToraster\points.shp'  #点数据文件
InputRasterFolder=r'D:\DATA\pointToraster\tif'  #放栅格数据的文件夹
# 设置Excel编码
file = xlwt.Workbook('encoding = utf-8') 
# 创建sheet工作表
sheet1 = file.add_sheet('sheet1',cell_overwrite_ok=True) 
#改变工作空间
#############获取矢量点位的经纬度
#设置driver
driver = ogr.GetDriverByName('ESRI Shapefile')
#打开矢量
ds = driver.Open(inputSHP, 0)
#获取图层
layer = ds.GetLayer(