#我是在mini conda环境运行的
import netCDF4 as nc
import rasterio
from rasterio.transform import from_origin
from rasterio.crs import CRS
# 定义输入和输出文件路径
input_file = r'G:\xxx\xxx\xxxxxxx.nc'#grace的文件位置和名称
output_folder = r'G:\xxxx'#输出文件位置
output_file_format = r'{}/{}.tiff'
#定义投影信息
projection = CRS.from_epsg(4326)
# 打开nc文件
dataset = nc.Dataset(input_file)
# 读取数据
time_data = dataset.variables['time'][:]
lon_data = dataset.variables['lon'][:]
lat_data = dataset.variables['lat'][:]
lwe_thickness_data = dataset.variables['lwe_thickness'][:]#我用的csr mascon的数据
# 获取数据的形状和维度
time_len, lon_len, lat_len = len(time_data), len(lon_data), len(lat_data)
# 日期列表
dates = [
200204, 200205, 200208, 200209, 200210, 200211, 200212, 200301, 200302, 200303, 200304, 200305,