淼的Python科研包|1.卫星降水数据下载

如果不明白怎么配置，可以参考往期付费推文https://mp.weixin.qq.com/s/FNDrmp9tZiYWWcVpAahsjw  ，一顿饭钱感谢各位老板！

from precipkit.download import configure_earthdataconfigure_earthdata()

from sysumkit.download.imerg import ImergDownloader# 1. 下载dl = ImergDownloader(    links_txt=下载链接地址,    output_folder=输出目标文件夹,    organize_by_date=True,)dl.download()# 2. 删除不完整的文件 这里设置为小于5MB的文件【这一步要小心，先确认正常文件大小是多少】dl.check_incomplete_files(    min_size_mb=5)# 3. 重试失败（可选）dl.retry_failed(max_retries=3)

GSMaP 数据下载是用固定的FTP账户密码端口，我已经在代码中写死，无需配置。

base_dir 指的是下载的产品，当前支持的是 V8最新版本 日尺度数据下载，如果需要其他版本的，可以参考源码进行更改。

start_mon, end_mon：起止年月

target: 输出文件夹

variable: 输出数据的变量名

from sysumkit.download.gsmap import GsmapDownloader# 1. 初始化下载器dl = GsmapDownloader(    base_dir='/standard/v8/daily_G/00Z-23Z/',    start_mon='202607',    end_mon='202607',    target=r'D:\GSMaP',    variable='precip')# 2. 执行下载dl.download()# 3. 如果有失败，重试dl.retry_failed(max_retries=3)# 4. 转换为 NCdl.gz_to_nc()

合并 IMERG 数据集

from sysumkit.download import IMERGAdapter, BasePrecipMerger# 数据适配器dataset = IMERGAdapter(    root=r"D:/IMERG")# 数据合并器merger = BasePrecipMerger(    adapter=dataset,    start_date="20200101",    end_date="20201231",    lat_range=(20, 30),    lon_range=(100, 120))# 合并并保存merger.merge(    r"D:/IMERG/IMERG.nc")

合并 GsMaP 数据集

from sysumkit.download import GSMaPAdapter, BasePrecipMergerimport xarray as xr# 数据适配器dataset = GSMaPAdapter(    root=r"D:/GSMaP")# 数据合并器merger = BasePrecipMerger(    adapter=dataset,    start_date="20260701",    end_date="20260727",    lat_range=(0, 60),    lon_range=(70, 140))# 合并并保存merger.merge(    r"D:/GSMaP/GSMaP.nc")ds = xr.open_dataset(r"D:/GSMaP/GSMaP.nc")ds.mean(dim='time').precip.plot()

欢迎加入 Python科研AI交流群：428657886

视频介绍