python+gdal批量提取研究区光谱和植被指数信息

python+gdal批量提取研究区光谱和植被指数信息

最近在做一些遥感相关的图像处理项目，涉及到遥感影像反射率和植被指数的计算。待提取遥感数据为多时序影像，利用ENVI或arcgis提取显然不太现实，于是就自己实现了一个并分享在这里，供需要的朋友参考。
除提取影像各波段的反射率外，这里重点提取了常用的植被指数NDVI、RVI、OSVI、ENVI以及NDWI。

#!/usr/bin/env python3.7
# encoding: utf-8
"""
  @author: ISR
  @contact: [email protected]
  @file: get_ref_rvi.py
  @time: 2019/12/20 20:30
  @desc:
  利用shp图层提取实验小区反射率和植被指数
  """
import numpy
import gdal
import ogr
import osr
import sys
import shapefile
from pandas