由于驱动数据的要求,之前下载了北师大的GLASS LAI数据 1KM数据进行处理,具体是经过了
1.下载,地址http://glass-product.bnu.edu.cn/index.html
2.数据的格式转换及拼接,原始数据并非TIF格式,而为hdf 可以使用 arcpy.ExtractSubDataset_management工具进行批量转化,之后再进行拼接,这样就可以得到全球的TIF格式LAI数据
3数据的裁剪,根据自己研究区域 对数据进行裁剪(把这一步放在前面可以减少后续处理所需的处理量)
import arcpy
import os
Coordinate_System = "PROJCS['WGS_1984_Web_Mercator_Auxiliary_Sphere',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Mercator_Auxiliary_Sphere'],PARAMETER['False_Easting',0.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',0.0],PARAMETER['Standard_Parallel_1',0.0],PARAMETER['Auxiliary_Sphere_Type',0.0],UNIT['Meter',1.0]]"#这里使用的是Web84墨卡托投影坐标
inputdir=r"——输入文件夹——"
output_path = r"——输出文件夹——"
mask = "——掩膜文件——