Google Earth Engine——GFS全球天气预报模型数据集:384小时的预测,预测间隔为3小时,以6小时的时间分辨率进行(每天更新4次)

The Global Forecast System (GFS) is a weather forecast model produced by the National Centers for Environmental Prediction (NCEP). The GFS dataset consists of selected model outputs (described below) as gridded forecast variables. The 384-hour forecasts, with 3-hour forecast interval, are made at 6-hour temporal resolution (i.e. updated four times daily). Use the 'creation_time' and 'forecast_time' properties to select data of interest.

The GFS is a coupled model, composed of an atmosphere model, an ocean model, a land/soil model, and a sea ice model which work together to provide an accurate picture of weather conditions. See history of recent modifications to the global forecast/analysis system, the model performance statistical web page, and the documentation homepage for more information.

全球预测系统(GFS)是由美国国家环境预测中心(NCEP)制作的一个天气预报模型。GFS数据集由选定的模型输出(如下所述)组成,作为网格化的预测变量。384小时的预测,预测间隔为3小时,以6小时的时间分辨率进行(即每天更新4次)。使用 "创建时间 "和 "预报时间 "属性来选择感兴趣的数据。

GFS是一个耦合模型,由一个大气模型、一个海洋模型、一个土地/土壤模型和一个海冰模型组成,它们一起工作以提供一个准确的天气状况图。更多信息见全球预报/分析系统最近的修改历史、模型性能统计网页和文件主页。

Dataset Availability

2015-07-01T00:00:00 - 2021-09-30T00:00:00

Dataset Provider

NOAA/NCEP/EMC

Collection Snippet

ee.ImageCollection("NOAA/GFS0P25")

Resolution

27830 meters

Bands Table

NameDescriptionMin*Max*Units
temperature_2m_above_groundTemperature 2m above ground-69.1852.25°C
specific_humidity_2m_above_groundSpecific humidity 2m above ground00.03kg/kg
relative_humidity_2m_above_groundRelative humidity 2m above ground1100.05%
u_component_of_wind_10m_above_groundU component of wind 10m above ground-60.7359.28m/s
v_component_of_wind_10m_above_groundV component of wind 10m above ground-63.7859.39m/s
total_precipitation_surfaceUntil 2019-11-07 06:00:00, this field represents the precipitation at surface at the forecasted time. After that date, this field holds the cumulative precipitation at surface added together from all forecasts starting from hour 0 (only for assets with forecast_hours > 0)0626.75kg/m^2
precipitable_water_entire_atmospherePrecipitable water for entire atmosphere0100kg/m^2
total_cloud_cover_entire_atmosphereTotal cloud cover for entire atmosphere (only for assets with forecast_hours > 0)0100%
downward_shortwave_radiation_fluxDownward shortwave radiation flux (only for assets with forecast_hours > 0)01230W/m^2

* = Values are estimated

影像属性:

NameTypeDescription
creation_timeDoubleTime of creation
forecast_hoursDoubleForecast hours
forecast_timeDoubleForecast time

数据使用:

NOAA data, information, and products, regardless of the method of delivery, are not subject to copyright and carry no restrictions on their subsequent use by the public. Once obtained, they may be put to any lawful use. The forgoing data is in the public domain and is being provided without restriction on use and distribution.

数据引用:

Alpert, J., 2006 Sub-Grid Scale Mountain Blocking at NCEP, 20th Conf. WAF/16 Conf. NWP P2.4.

Alpert, J. C., S-Y. Hong and Y-J. Kim: 1996, Sensitivity of cyclogenesis to lower troposphere enhancement of gravity wave drag using the EMC MRF”, Proc. 11 Conf. On NWP, Norfolk, 322-323.

Alpert,J,, M. Kanamitsu, P. M. Caplan, J. G. Sela, G. H. White, and E. Kalnay, 1988: Mountain induced gravity wave drag parameterization in the NMC medium-range forecast model. Pre-prints, Eighth Conf. on Numerical Weather Prediction, Baltimore, MD, Amer. Meteor. Soc., 726-733.

Buehner, M., J. Morneau, and C. Charette, 2013: Four-dimensional ensemble-variational data assimilation for global deterministic weather prediction. Nonlinear Processes Geophys., 20, 669–682.

Chun, H.-Y., and J.-J. Baik, 1998: Momentum Flux by Thermally Induced Internal Gravity Waves and Its Approximation for Large-Scale Models. J. Atmos. Sci., 55, 3299–3310.

Chun, H.-Y., Song, I.-S., Baik, J.-J. and Y.-J. Kim. 2004: Impact of a Convectively Forced Gravity Wave Drag Parameterization in NCAR CCM3. J. Climate, 17, 3530–3547.

Chun, H.-Y., Song, M.-D., Kim, J.-W., and J.-J. Baik, 2001: Effects of Gravity Wave Drag Induced by Cumulus Convection on the Atmospheric General Circulation. J. Atmos. Sci., 58, 302–319.

Clough, S.A., M.W. Shephard, E.J. Mlawer, J.S. Delamere, M.J. Iacono, K.Cady-Pereira, S. Boukabara, and P.D. Brown, 2005: Atmospheric radiative transfer modeling: A summary of the AER codes, J. Quant. Spectrosc. Radiat. Transfer, 91, 233-244. doi:10.1016/j.jqsrt.2004.05.058

Ebert, E.E., and J.A. Curry, 1992: A parameterization of ice cloud optical properties for climate models. J. Geophys. Res., 97, 3831-3836.

Fu, Q., 1996: An Accurate Parameterization of the Solar Radiative Properties of Cirrus Clouds for Climate Models. J. Climate, 9, 2058-2082.

Han, J., and H.-L. Pan, 2006: Sensitivity of hurricane intensity forecast to convective momentum transport parameterization. Mon. Wea. Rev., 134, 664-674.

Han, J., and H.-L. Pan, 2011: Revision of convection and vertical diffusion schemes in the NCEP global forecast system. Weather and Forecasting, 26, 520-533.

Han, J., M. Witek, J. Teixeira, R. Sun, H.-L. Pan, J. K. Fletcher, and C. S. Bretherton, 2016: Implementation in the NCEP GFS of a hybrid eddy-diffusivity mass-flux (EDMF) boundary layer parameterization with dissipative heating and modified stable boundary layer mixing. Weather and Forecasting, 31, 341-352.

Hou, Y., S. Moorthi and K. Campana, 2002: Parameterization of Solar Radiation Transfer in the NCEP Models, NCEP Office Note #441, pp46. Available here

Hu, Y.X., and K. Stamnes, 1993: An accurate parameterization of the radiative properties of water clouds suitable for use in climate models. J. Climate, 6, 728-74.

Iacono, M.J., E.J. Mlawer, S.A. Clough, and J.-J. Morcrette, 2000: Impact of an improved longwave radiation model, RRTM, on the energy budget and thermodynamic properties of the NCAR community climate model, CCM3, J. Geophys. Res., 105(D11), 14,873–14,890.2.

Johansson, Ake, 2008: Convectively Forced Gravity Wave Drag in the NCEP Global Weather and Climate Forecast Systems, SAIC/Environmental Modelling Center internal report.

Juang, H-M, et al. 2014:Regional Spectral Model workshop in memory of John Roads and Masao Kanamitsu, BAMS, A. Met. Soc, ES61-ES65.

Kim, Y.-J., and A. Arakawa (1995), Improvement of orographic gravity wave parameterization using a mesoscale gravity-wave model, J. Atmos. Sci.,52, 875–1902.

Kleist, D. T., 2012: An evaluation of hybrid variational-ensemble data assimilation for the NCEP GFS , Ph.D. Thesis, Dept. of Atmospheric and Oceanic Science, University of Maryland-College Park, 149 pp.

Lott, F and M. J. Miller: 1997, “A new subgrid-scale orographic drag parameterization: Its formulation and testing”, QJRMS, 123, pp101-127.

Mlawer, E.J., S.J. Taubman, P.D. Brown, M.J. Iacono, and S.A. Clough, 1997: Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave. J. Geophys. Res., 102, 16663-16682.

Sela, J., 2009: The implementation of the sigma-pressure hybrid coordinate into the GFS. NCEP Office Note #461, pp25.

Sela, J., 2010: The derivation of sigmapressure hybrid coordinate semi-Lagrangian model equations for the GFS. NCEP Office Note #462 pp31.

Yang, F., 2009: On the Negative Water Vapor in the NCEP GFS: Sources and Solution. 23rd Conference on Weather Analysis and Forecasting/19th Conference on Numerical Weather Prediction, 1-5 June 2009, Omaha, NE.

Yang, F., K. Mitchell, Y. Hou, Y. Dai, X. Zeng, Z. Wang, and X. Liang, 2008: Dependence of land surface albedo on solar zenith angle: observations and model parameterizations. Journal of Applied Meteorology and Climatology.No.11, Vol 47, 2963-2982.

代码:

var dataset = ee.ImageCollection('NOAA/GFS0P25')
                  .filter(ee.Filter.date('2018-03-01', '2018-03-02'));
var temperatureAboveGround = dataset.select('temperature_2m_above_ground');
var visParams = {
  min: -40.0,
  max: 35.0,
  palette: ['blue', 'purple', 'cyan', 'green', 'yellow', 'red'],
};
Map.setCenter(71.72, 52.48, 3.0);
Map.addLayer(temperatureAboveGround, visParams, 'Temperature Above Ground');

好的,您的问题是如何编写基于TensorFlow和LSTM模型的Python自定义类ModelLstm预测全球风速代码,读取720个gfs预报模型文件,表示720小时全球格点数据,每个文件是等经纬投影0.25度的分辨率,作为训练数据(batch_size, time_steps, features) 维度数据样例,用24小时数据预测未来24小时风速。是吗? 以下是一个示例代码,可以实现您的需求: ```python import tensorflow as tf import numpy as np import os class ModelLstm: def __init__(self, batch_size, time_steps, features): self.batch_size = batch_size self.time_steps = time_steps self.features = features self.model = self.build_model() def build_model(self): model = tf.keras.Sequential([ tf.keras.layers.LSTM(64, input_shape=(self.time_steps, self.features), return_sequences=True), tf.keras.layers.Dropout(0.2), tf.keras.layers.LSTM(32), tf.keras.layers.Dropout(0.2), tf.keras.layers.Dense(1) ]) model.compile(loss='mse', optimizer='adam') return model def train(self, data, labels, epochs): self.model.fit(data, labels, epochs=epochs, batch_size=self.batch_size) def predict(self, data): return self.model.predict(data) def read_data(): data_path = 'your_path_to_gfs_data' file_list = os.listdir(data_path) data_list = [] for file in file_list: file_data = np.load(os.path.join(data_path, file)) data_list.append(file_data) data = np.concatenate(data_list, axis=0) return data def prepare_data(data): data_x = [] data_y = [] for i in range(data.shape[0] - 48): data_x.append(data[i:i+24]) data_y.append(data[i+24:i+48, :, :, 0]) data_x = np.array(data_x) data_y = np.array(data_y) return data_x, data_y if __name__ == '__main__': data = read_data() data_x, data_y = prepare_data(data) model = ModelLstm(batch_size=32, time_steps=24, features=2) model.train(data_x, data_y, epochs=10) test_data = np.random.rand(1, 24, 2) prediction = model.predict(test_data) print(prediction) ``` 这个代码的大致流程如下: 1. 定义了一个名为 ModelLstm 的类,它包含了 LSTM 模型的构建、训练和预测方法。 2. 定义了一个名为 read_data 的函数,用来读取数据。 3. 定义了一个名为 prepare_data 的函数,用来将读取到的数据进行处理,得到训练数据和标签。 4. 在 main 函数中,读取数据并进行处理,然后使用 ModelLstm 类进行训练和预测。 其中,prepare_data 函数中的代码需要根据您的数据格式进行修改。具体来说,data_x 表示训练数据,data_y 表示标签,它们的维度应该为 (样本数, 时间步数, 特征数) 和 (样本数, 预测时间步数, 特征数)。在这个示例中,特征数为 2,第一维表示经度,第二维表示纬度。您需要根据您的数据格式进行修改。 另外,如果您的数据量很大,可能需要使用 TensorFlow 中的 Dataset API 进行数据读取和处理。这个示例中只是简单地使用了 Python 中的 for 循环进行数据读取和处理。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

此星光明

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值