python将GRACE数据NC文件转为tiff

本文档详细介绍了如何在Python环境中,利用netCDF4和rasterio库将GRACE的NetCDF文件转换成地理TIFF格式。通过读取文件,设置投影信息,遍历时间序列数据并保存每个时间步的厚度数据,最终生成一系列TIFF文件。
摘要由CSDN通过智能技术生成

#我是在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,

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值