извлечение данных на основе диапазона широты

>>> dfForecastMSL
<xarray.Dataset>
Dimensions:     (latitude: 1536, longitude: 3072, valid_time: 29)
Coordinates:
    meanSea     int64 0
    time        datetime64[ns] 2020-05-30
  * longitude   (longitude) float64 0.0 0.1172 0.2344 ... 359.6 359.8 359.9
  * latitude    (latitude) float64 89.91 89.79 89.68 ... -89.68 -89.79 -89.91
    step        (valid_time) timedelta64[ns] 0 days 00:00:00 ... 7 days 00:00:00
  * valid_time  (valid_time) datetime64[ns] 2020-05-30 ... 2020-06-06
Data variables:
    msl         (valid_time, latitude, longitude) float32 102246.484 ... 103826.81
Attributes:
    GRIB_edition:            2
    GRIB_centre:             kwbc
    GRIB_centreDescription:  US National Weather Service - NCEP
    GRIB_subCentre:          0
    Conventions:             CF-1.7
    institution:             US National Weather Service - NCEP
    history:                 2020-09-19T07:31:46 GRIB to CDM+CF via cfgrib-0....
>>> dfForecastMSL.history
'2020-09-19T07:31:46 GRIB to CDM+CF via cfgrib-0.9.8.4/ecCodes-2.17.0 with {"source": "/home/NCMRWFTEMP/vsprasad/EXP_HY2B/data/gdasv14/gdas/prodCNTL/gdas.20200530/gdas.t00z.master.grb2f00", "filter_by_keys": {"cfVarName": "msl", "typeOfLevel": "meanSea"}, "encode_cf": ["parameter", "time", "geography", "vertical"]}'

Я пытаюсь извлечь интересующую меня область вот так.

india = dfForecastMSL.sel(longitude=slice(60,100),latitude=slice(0,40))

И результат выглядит так ....

>>> india = dfForecastMSL.sel(longitude=slice(60,100),latitude=slice(0,40))
>>> india
<xarray.Dataset>
Dimensions:     (latitude: 0, longitude: 342, valid_time: 29)
Coordinates:
    meanSea     int64 0
    time        datetime64[ns] 2020-05-30
  * longitude   (longitude) float64 60.0 60.12 60.23 60.35 ... 99.73 99.84 99.96
  * latitude    (latitude) float64
    step        (valid_time) timedelta64[ns] 0 days 00:00:00 ... 7 days 00:00:00
  * valid_time  (valid_time) datetime64[ns] 2020-05-30 ... 2020-06-06
Data variables:
    msl         (valid_time, latitude, longitude) float32
Attributes:
    GRIB_edition:            2
    GRIB_centre:             kwbc
    GRIB_centreDescription:  US National Weather Service - NCEP
    GRIB_subCentre:          0
    Conventions:             CF-1.7
    institution:             US National Weather Service - NCEP
    history:                 2020-09-19T07:31:46 GRIB to CDM+CF via cfgrib-0....

** Почему не хватает широты? ** Или я что-то изношен?


python-xarray
person M. Sateesh    schedule 19.09.2020    source источник
comment
stackoverflow.com/help/minimal-reproducible-example   -  person Maximilian    schedule 19.09.2020

Ответы (1)


arrow_upward
0
arrow_downward

Похоже, что координата широты сохраняется от +90 до -90, то есть в порядке убывания. Срез должен быть взят в том же порядке, чтобы он был .sel(latitude=slice(40,0)).

person observingClouds    schedule 11.10.2020