ESA CCI 土壤水分数据处理结果分享

ESA 土壤水分下载网址

https://www.esa-landcover-cci.org/

批量下载

可用Motrix-1.8.19-x64

Motrix 是一个功能全面的下载管理器，支持包括HTTP、FTP、BitTorrent及Magnet链接等多种下载协议。此项目是开源的，由agalwood维护，并在GitHub上拥有活跃的社区支持。它以其简洁的界面和强大的后台处理能力而受到用户的喜爱，为用户提供免费且高效的文件下载解决方案。

项目网址：https://github.com/agalwood/Motrix

motrix支持

1. 断点续传: Motrix支持断点续传，即使下载中断，也能从上次中断处继续。
2. 批量下载: 利用任务列表特性，你可以一次性添加多个下载任务，特别适合资源包下载。
3. 计划下载: 结合系统调度工具，可以安排在低峰时段自动下载，节省带宽资源。

或者尝试火狐插件DownThemAll

nc下载与处理

nc转tif代码

import xarray as xr
import rasterio
from rasterio.transform import from_origin

# 打开 NetCDF 文件
nc_file = r'C:\Users\iQier\Desktop\ESACCI.nc'
ds = xr.open_dataset(nc_file)

# 假设需要提取变量 'variable_name'（根据具体变量名进行调整）
variable = ds['sm']  # 例如：'temperature'，'precipitation' 等

# 如果是4D数据（例如：time, lat, lon），选择一个时间点（例如第0个时间点）
# 如果你的数据是3D的（例如：lat, lon），可以直接使用
data = variable.isel(time=0).values  # 选择第0个时间点的数据

# 获取坐标轴信息（假设为纬度、经度信息）
lons = ds['lon'].values
lats = ds['lat'].values

# 创建栅格数据集的转换信息
transform = from_origin(lons.min(), lats.max(), abs(lons[1] - lons[0]), abs(lats[1] - lats[0]))

# 设置输出文件路径
tif_file = r'C:\Users\iQier\Desktop\output.tif'

# 使用 rasterio 创建 GeoTIFF 文件
with rasterio.open(tif_file, 'w', driver='GTiff',
                   height=data.shape[0], width=data.shape[1],
                   count=1, dtype=data.dtype,
                   crs='+proj=latlong', transform=transform) as dst:
    dst.write(data, 1)

print(f"转换完成，文件保存为 {tif_file}")

运行结果

nc转tif结果显示

2000年逐月结果

2000-2020逐年结果

练习数据与python代码如下

超链接：完整代码 测试数据 2000年逐月结果 2000-2020逐年结果