無法在 NetCDF 衛星文件 (Xarray) 中沿緯度和經度切片

Question

我的數據包含緯度和經度值，但我無法使用 Xarray 對它們進行切片。 當我嘗試運行這樣的代碼時：

SeaIce[0].sel(latitude = slice(72,80))

我收到以下錯誤

ValueError: dimensions or multi-index levels ['latitude'] do not exist

但維度確實存在。 這是我正在使用的數據集之一的示例：

Dimensions:                                 (time: 1, xgrid: 304, ygrid: 448)
Coordinates:
  * time                                    (time) datetime64[ns] 2002-01-01
  * ygrid                                   (ygrid) float32 5.838e+06 ... -5....
  * xgrid                                   (xgrid) float32 -3.838e+06 ... 3....
    latitude                                (ygrid, xgrid) float64 31.1 ... 3...
    longitude                               (ygrid, xgrid) float64 168.3 ... ...
Data variables:
    projection                              |S1 b''
    seaice_conc_monthly_cdr                 (time, ygrid, xgrid) float32 0.0 ...
    stdev_of_seaice_conc_monthly_cdr        (time, ygrid, xgrid) float32 ...
    melt_onset_day_seaice_conc_monthly_cdr  (time, ygrid, xgrid) float32 ...
    qa_of_seaice_conc_monthly_cdr           (time, ygrid, xgrid) float32 ...
    goddard_merged_seaice_conc_monthly      (time, ygrid, xgrid) float32 ...
    goddard_nt_seaice_conc_monthly          (time, ygrid, xgrid) float32 ...
    goddard_bt_seaice_conc_monthly          (time, ygrid, xgrid) float32 ...
Attributes: (12/41)
    Conventions:               CF-1.6, ACDD-1.3
    title:                     NOAA/NSIDC Climate Data Record of Passive Micr...
    references:                Comiso, J. C., and F. Nishio. 2008. Trends in ...
    comment:                   The variable melt_onset_day_seaice_conc_monthl...
    program:                   NOAA Climate Data Record Program
    cdr_variable:              seaice_conc_monthly_cdr
                       ...
    history:                   seaice_conc_monthly_nh_f13_200201_v03r01.nc cr...
    platform:                  DMSP 5D-2/F13 > Defense Meteorological Satelli...
    sensor:                    SSM/I > Special Sensor Microwave/Imager
    date_created:              2017-11-27T02:44:53Z
    time_coverage_start:       2002-01-01T00:00:00Z
    time_coverage_end:         2002-01-31T23:59:59Z

我嘗試過沿 xgrid 和 ygrid 維度進行切片，並且這些工作有效，但自然很難將注意力集中在沒有緯度和經度維度的正確區域上。

Answer 1

我設法解決了這個 using.where 和下面的一系列命令：

dr = te.where((te.latitude > 72) & (te.latitude < 80) & (te.longitude > -15) & (te.longitude < 10),drop=True)

這選擇了以 72 到 80 緯度和 -15 到 10 經度為界的區域。