在计算过程中,需要根据gis里对于土地利用的不同类型对应不同的计算公式,所以我先导出为tif文件再转为nc文件
from osgeo import gdal
input_file = "F:/LUCC/NJ_LU2000.tif"
output_file = "F:/LUCC/lucc_2000.nc"
gdal.Translate(output_file,input_file,format='netCDF')
osgeo文件可以用conda安装
pip install osgeo