从域 ncfile 中提取数据

Question

我正在尝试从由纬度/经度边界定义的用户指定（非矩形）域中的每个网格单元中提取数据。我的输入文件位于曲线网格上。 python,cdo,ncks各种方法都试过了，还是想不通。我只想要输入 ncfile 的多边形域子集中每个网格单元的时间序列信息。我输入的 ncfile 信息在这里给出：

$ ncdump -h 1979_sfc_out.nc 

netcdf 79_sfc_out {
dimensions:
    x = 83 ;
    y = 94 ;
    time = UNLIMITED ; // (8736 currently)
    nv4 = 4 ;
variables:
    float time(time) ;
        time:axis = "T" ;
        time:long_name = "time" ;
        time:standard_name = "time" ;
        time:units = "hours since 1979-1-2 00:00:00" ;
        time:calendar = "standard" ;
    float x(x) ;
        x:axis = "x" ;
        x:long_name = "X-coordinate in Cartesian system" ;
        x:standard_name = "projection_x_coordinate" ;
        x:units = "meters" ;
    float y(y) ;
        y:axis = "y" ;
        y:long_name = "Y-coordinate in Cartesian system" ;
        y:standard_name = "projection_y_coordinate" ;
        y:units = "meters" ;
    float lon(y, x) ;
        lon:units = "degrees_east" ;
        lon:valid_range = -180., 180. ;
        lon:standard_name = "longitude" ;
        lon:bounds = "lon_bnds" ;
    float lat(y, x) ;
        lat:units = "degrees_north" ;
        lat:valid_range = -90., 90. ;
        lat:standard_name = "latitude" ;
        lat:bounds = "lat_bnds" ;
    float lon_bnds(y, x, nv4) ;
        lon_bnds:units = "degreesE" ;
    float lat_bnds(y, x, nv4) ;
        lat_bnds:units = "degreesN" ;
    char mapping ;
        mapping:false_easting = 0. ;
        mapping:false_northing = 0. ;
        mapping:grid_mapping_name = "polar_stereographic" ;
        mapping:latitude_of_projection_origin = 90. ;
        mapping:standard_parallel = 64. ;
        mapping:straight_vertical_longitude_from_pole = -152. ;
        mapping:semi_major_axis = 6370000. ;
        mapping:semi_minor_axis = 6370000. ;
    float SEAICE(time, y, x) ;
        SEAICE:_FillValue = -9999.f ;
        SEAICE:units = "fraction" ;
        SEAICE:long_name = "Ice concentration (ice=1;no ice=0)" ;
        SEAICE:grid_mapping = "mapping" ;
        SEAICE:coordinates = "lon lat" ;

我尝试过的一些东西是

lat1=71.2
lat2=72.9
lon1=-176.5
lon2=-160

cdo sellonlatbox,lon1,lon2,lat1,lat2 $ifile $box1_ofile
cdo sellonlatbox (Abort): Float parameter >lon1< contains invalid character at position 1!

我认为问题是我的输入文件有以米为单位的 x,y 维度（没有负号，可能是 'invalid character at position 1'），我要求 cdo 提取 lat/lon 以度为单位的维度。我的 ncfile 中有变量 'mapping'，它可能有助于将米转换为 lat/lon，但我不知道该怎么做。

而且 ncks 在这里对我不起作用，可能是由于我在 cdo 上看到的相同的仪表 <-> lon/lat 问题。

ncks -v SEAICE,U10,V10 -d latitude,71.2,72.9 -d longitude,-176.5,-160. $ifile -O $box1_ofile

ncks: ERROR dimension latitude is not in input file

虽然我想提取一个多边形，但我试过的这些例子只是矩形子集，我想得到多边形我可以做多个矩形子集来实现我最终的多边形形状，但如果有更好的方法如果有任何建议，我们将不胜感激。

谢谢

Answer 1

NCO 的辅助坐标 feature 设计用于在具有一维纬度和经度的非结构化网格上对曲线坐标进行超切片。试试这个：

ncks -X lon_min,lon_max,lat_min,lat_max in.nc out.nc
ncks -X -176.5,-160.,71.2,72.9 in.nc out.nc

您也可以菊花链式连接多个 -X 选项以获得奇怪的多边形。不幸的是，这可能对您不起作用，因为我刚刚注意到，您有一个带有 2D lat 和 lon 的曲线网格。为此，请尝试 feature 处的 ncap2。手册中的示例显示了如何沿矩形边界进行遮罩，您可以菊花链式链接 where() 语句中的条件以获得多边形。这不会改变输出文件的维度，但它允许您将多边形之外的所有内容设置为 _FillValue。

从域 ncfile 中提取数据

extracting data from a domain ncfile

python

coordinates

netcdf

nco

cdo-climate