Ctrl+K

数据源

Contents

数据源#

cedarkit-maps 不限定数据加载的工具库，可以支持多种数据源。 只要加载后的数据对象满足下面给出的基本条件，就可以被正确识别。

基本条件

格点数据为 xarray.DataArray 格式，以纬度 (latitude)、经度 (longitude) 为坐标的二维数组，按行存储（纬度）。

示例#

使用 reki 加载 CMA-MESO 的 2 米温度场

import pandas as pd
from reki.data_finder import find_local_file
from reki.format.grib.eccodes import load_field_from_file

file_path = find_local_file(
    "cma_meso_3km/grib2/orig",
    start_time=pd.to_datetime("2024-04-01 00:00"),
    forecast_time=pd.to_timedelta("24h"),
)
  
t_2m_field = load_field_from_file(
    file_path,
    parameter="2t",  
) - 273.15
t_2m_field

<xarray.DataArray '2t' (latitude: 1671, longitude: 2501)> Size: 33MB
array([[-0.71931641, -0.75931641, -0.80931641, ..., -4.16931641,
        -2.87931641, -2.95931641],
       [-3.65931641, -3.66931641, -2.27931641, ..., -6.91931641,
        -6.31931641, -3.72931641],
       [-3.54931641, -3.58931641, -0.85931641, ..., -7.60931641,
        -7.47931641, -3.43931641],
       ...,
       [26.60068359, 26.60068359, 26.59068359, ..., 25.80068359,
        25.79068359, 25.97068359],
       [26.59068359, 26.60068359, 26.59068359, ..., 25.78068359,
        25.78068359, 26.00068359],
       [26.59068359, 26.59068359, 26.60068359, ..., 25.78068359,
        25.79068359, 25.97068359]])
Coordinates:
    time               datetime64[ns] 8B 2024-04-01
    step               timedelta64[ns] 8B 1 days
    valid_time         datetime64[ns] 8B 2024-04-02
    heightAboveGround  int64 8B 2
  * latitude           (latitude) float64 13kB 60.1 60.07 60.04 ... 10.03 10.0
  * longitude          (longitude) float64 20kB 70.0 70.03 70.06 ... 145.0 145.0

上述要素场包含两个维度属性

  • latitude：纬度

  • longitude：经度

同时包含其他非维度的坐标属性

  • time：起报时次

  • step：预报时效

  • valid_time：预报时间

  • heightAboveGround：层次类型，离地面高度，单位米

Contents