小男孩‘自慰网亚洲一区二区,亚洲一级在线播放毛片,亚洲中文字幕av每天更新,黄aⅴ永久免费无码,91成人午夜在线精品,色网站免费在线观看,亚洲欧洲wwwww在线观看

分享

一招搞定各大機(jī)構(gòu)數(shù)值預(yù)報產(chǎn)品下載,包括EC最新公開預(yù)報數(shù)據(jù)

 江海博覽 2022-03-14

不知不覺停更了一個多月,除了寫論文之外,還參加了上海氣象學(xué)會主辦的 第二屆人工智能天氣預(yù)報創(chuàng)新大賽,主要是對數(shù)值模式預(yù)報的2m氣溫和降水進(jìn)行訂正。

這是第二次參加上海氣象學(xué)會主辦的比賽,第一屆的短臨預(yù)報比賽和小伙伴一起最后獲得了二等獎。這一次只有我一個人,最后因?yàn)楦鞣N原因只參加了前半程,提交了大概3次結(jié)果,從最初的雙賽道第一最后定格到第四。

這幾年國內(nèi)外和氣象相關(guān)的AI大賽比較多,改天再細(xì)聊一下AI+氣象大賽

剛好近兩天ECMWF開放了預(yù)報數(shù)據(jù),那就來看一下如何獲取到全球各大氣象機(jī)構(gòu)發(fā)布的數(shù)值預(yù)報產(chǎn)品吧。

全球有不少機(jī)構(gòu)對外發(fā)布了數(shù)值預(yù)報產(chǎn)品,但是大多都比較分散,通常需要自己去尋找再寫不少代碼來下載這些產(chǎn)品。前幾天發(fā)現(xiàn)了一個工具,可以非常方便的下載常用的數(shù)值預(yù)報產(chǎn)品。ECMWF剛對外發(fā)布預(yù)報數(shù)據(jù),開發(fā)者就更新了下載接口,非常給力了哈哈!

廢話不多說了,下面就來介紹一下這個工具。

圖片

Herbie 屬于 NOAA 的大數(shù)據(jù)項目,旨在更方便的從各云存儲中獲取數(shù)值天氣預(yù)報結(jié)果。最初,是在 Brian K.Blaylock 博士期間開發(fā)用來獲取 HRRR數(shù)據(jù)的。2020年此工具成為 NOAA 大數(shù)據(jù)項目的一部分,于是開始開發(fā)適用更多的數(shù)據(jù)源,目前支持獲取如下數(shù)據(jù):

  • High Resolution Rapid Refresh (HRRR)
  • Rapid Refresh (RAP)
  • Global Forecast System (GFS)
  • National Blend of Models (NBM)
  • Rapid Refresh Forecast System - Prototype (RRFS)
  • ECMWF Open Data Forecast Prodducts (? new in Herbie 0.0.8)

其數(shù)值天氣預(yù)報結(jié)果的數(shù)據(jù)源包括:

  • NOMADS
  • Big Data Program Partners (AWS, Google, Azure)
  • ECMWF Open Data Azure storage
  • University of Utah CHPC Pando archive

安裝

最好使用 Python3.8+ 版本,安裝方式如下:

pip install herbie-data

pip install git+https://github.com/blaylockbk/Herbie.git

食用方法

Herbie 目前具有如下功能:

  • 從不同數(shù)據(jù)源搜索模式輸出
  • 下載完整的GRIB2文件
  • 從GRIB2文件中提取子集并下載
  • 使用 xarray 讀取數(shù)據(jù)
  • 使用 Cartopy 繪圖

總的來說,Herbie 使用方便,而且兼容目前常用的工具。

from herbie.archive import Herbie

# Herbie object for the HRRR model 6-hr surface forecast product
H = Herbie('2021-01-01 12:00',
           model='hrrr',
           product='sfc',
           fxx=6)

# Download the full GRIB2 file
H.download()

# Download a subset, like all fields at 500 mb
H.download(':500 mb')

# Read subset with xarray, like 2-m temperature.
H.xarray('TMP:2 m')

有時可能需要下載大量的數(shù)據(jù), Herbie 也提供了非常方便的工具函數(shù):

from herbie.tools import bulk_download
import pandas as pd

# Download HRRR analysis for every hour of a day
DATES = pd.date_range('2021-01-01 00:00''2021-01-01 06:00', freq='1H')
bulk_download(DATES, model='hrrr', product='sfc', fxx=0, searchString='TMP:2 m')

可以通過如下命令下載ECMWF最新公開的數(shù)據(jù)集:

from datetime import datetime
from herbie.archive import Herbie

now = datetime.now()
today = datetime(now.year, now.month, now.day)

H = Herbie(today, model='ecmwf', product='oper', save_dir='./')

H.download() # 下載所有文件

H.download(':t:'# 下載所有層的溫度

就是下載速度有點(diǎn)慢,好了就介紹到這,感興趣的可以去 GitHub (https://github.com/blaylockbk/Herbie) 網(wǎng)頁查看。目前主要是獲取數(shù)值預(yù)報產(chǎn)品,數(shù)據(jù)處理和繪圖部分正處于開發(fā)中,如果有什么問題也可以去反饋。

圖片
END

    本站是提供個人知識管理的網(wǎng)絡(luò)存儲空間,所有內(nèi)容均由用戶發(fā)布,不代表本站觀點(diǎn)。請注意甄別內(nèi)容中的聯(lián)系方式、誘導(dǎo)購買等信息,謹(jǐn)防詐騙。如發(fā)現(xiàn)有害或侵權(quán)內(nèi)容,請點(diǎn)擊一鍵舉報。
    轉(zhuǎn)藏 分享 獻(xiàn)花(0

    0條評論

    發(fā)表

    請遵守用戶 評論公約

    類似文章 更多