weixin_45452300的博客

根据shp图对遥感影像进行裁剪，shp图有多个面的话将每个面裁剪出来并根据属性命名。

基于geopandas和rasterio实现区域统计分析，批量获取我们感兴趣的ROI的统计值。

利用GDAL库对多光谱影像进行分析，光谱5波段，顺序import numpy as npfrom osgeo import gdal, osrdef export_tif(out_tif_name, var_lat, var_lon, NDVI): # 分辨率计算 N_Lat = len(var_lat) N_Lon = len(var_lon) LonMin, LatMax, LonMax, LatMin = [min(var_lon), max(var_lat)

将tif转为DataFrame存为csv此数据为modis的LAI产品，仅有LAI一个波段from osgeo import gdalimport pandas as pdgdal.AllRegister()filePath = r'C:\Users\Administrator\遥感文件\MOD15A2H.006_Lai_500m_doy2020097_aid0001.tif'dataset = gdal.Open(filePath)adfGeoTransform = dataset.Ge