二维数组中的纬度/经度索引答案

【问题标题】：lat/lon indexing in 2d array二维数组中的纬度/经度索引
【发布时间】：2018-08-16 08:14:33
【问题描述】：

我有一个.nc 文件，其中纬度和经度在二维数组中。我想提取某些已知点的值。我没有找到找到那些已知点的索引的方法。

数据结构如下：

OrderedDict([(u'lon', <type 'netCDF4._netCDF4.Variable'>
              float32 lon(y, x)
                  standard_name: longitude
                  long_name: longitude
                  units: degrees_east
                  _CoordinateAxisType: Lon
              unlimited dimensions: 
              current shape = (146, 209)
              filling off), (u'lat', <type 'netCDF4._netCDF4.Variable'>
              float32 lat(y, x)
                  standard_name: latitude
                  long_name: latitude
                  units: degrees_north
                  _CoordinateAxisType: Lat
              unlimited dimensions: 
              current shape = (146, 209)
              filling off), (u'time', <type 'netCDF4._netCDF4.Variable'>
              float64 time(time)
                  standard_name: time
                  units: seconds since 1949-12-1 00:00:00
                  calendar: proleptic_gregorian
                  axis: T
              unlimited dimensions: time
              current shape = (55152,)
              filling off), (u'TOT_PREC', <type 'netCDF4._netCDF4.Variable'>
              float32 TOT_PREC(time, y, x)
                  standard_name: precipitation_amount
                  long_name: total precipitation amount
                  units: kg m-2
                  coordinates: lat lon
                  cell_methods: time: sum
              unlimited dimensions: time
              current shape = (55152, 146, 209)
              filling off)])

假设如果我想提取坐标为：lat=27.25、lon=85.75 的点的值，那么如何获取这些坐标的纬度和经度索引。

任何帮助或建议将不胜感激！

【问题讨论】：

标签： python netcdf python-xarray

【解决方案1】：

首先，你读取 .nc 文件和里面的变量

ds1 = NetCDFFile('file.nc')
ds2 = ...

lon = ds.variables['lon'][:]
lat = ds.variables['lat'][:]

然后使用 'for' 循环 {}{}{} lat,lon,data 写入文件然后你很容易读取数据。希望这有点道理。

【讨论】：