依据矢量点批量计算遥感图像的植被指数（Python代码）

是仔仔啊

已于 2023-07-28 11:31:24 修改

阅读量345

点赞数

分类专栏：卫星影像特征文章标签： python

于 2023-07-28 11:24:42 首次发布

本文链接：https://blog.youkuaiyun.com/weixin_44576131/article/details/131976321

版权

依据矢量点批量计算遥感图像的植被指数（Python代码）

根据不同作物类型的矢量点计算植被指数：
注：目前代码所使用遥感图像为单景，若要计算时序的遥感数据，须修改下代码

from osgeo import ogr, osr
from osgeo import gdal
import os, sys
import glob
import numpy as np
import pandas as pd

#############获取矢量点位的经纬度
#设置driver
driver = ogr.GetDriverByName('ESRI Shapefile')
#打开矢量
ds = driver.Open('G:\\WS_XL.shp', 0)
if ds is None:
    print('Could not open ' +'sites.shp')
    sys.exit(1)
#获取图层
layer = ds.GetLayer(0)
#获取要素及要素地理位置
xValues = []
yValues = []
name=[]
feature = layer.GetNextFeature()
while feature:
    types = feature.GetField('Value')
    name.append(types)
    geometry = feature.GetGeometryRef()
    x = geometry