2

我有一个分辨率为 0.25 度的全局数据，我想屏蔽它，以便它只包含陆地上的数据值。

数据涵盖 lon 维度的完整 360 度和 lat 维度的 -60 到 60 度。

下面列出了文件头以及汇总的经纬度坐标值：

netcdf cmorph_global_daily {
dimensions:
    lat = UNLIMITED ; // (480 currently)
    lon = 1440 ;
    time = 7305 ;
variables:
    float lat(lat) ;
        lat:units = "degrees_north" ;
        lat:long_name = "Latitude" ;
    float lon(lon) ;
        lon:units = "degrees_east" ;
        lon:long_name = "Longitude" ;
    float prcp(lat, lon, time) ;
        prcp:_FillValue = NaNf ;
        prcp:units = "mm" ;
        prcp:standard_name = "precipitation" ;
        prcp:long_name = "Precipitation" ;
        prcp:description = "CMORPH Version 1.0BETA Version, daily precip from 00Z-24Z" ;
    int time(time) ;
        time:units = "days since 1900-01-01" ;
        time:long_name = "Time" ;
        time:calendar = "gregorian" ;

// global attributes:
        :history = "Mon Mar 26 10:44:42 2018: ncpdq -a lat,lon,time cmorph_adjusted_daily.nc latlontime/cmorph_adjusted_daily.nc\nThu Mar 15 10:21:10 2018: ncks -4 cmorph_adjusted_daily.nc cmorph_adjusted_daily.nc" ;
        :nco_openmp_thread_number = 1 ;
        :title = "CMORPH Version 1.0BETA Version, daily precip from 00Z-24Z" ;
        :NCO = "4.7.2" ;
data:

 lat = -59.875, -59.625, -59.375, -59.125, ..., 59.125, 59.375, 59.625, 59.875 ;
 lon = 0.125, 0.375, 0.625, 0.875, 1.125, ..., 359.125, 359.375, 359.625, 359.875 ;

我更喜欢使用 Python/numpy 和/或 NCO 来执行此操作，因为这是我的典型工具集。

4

2 回答 2

5

一旦你在同一个网格上有一个变量来屏蔽，你可以使用 ncap2 where，例如，

ncap2 -s 'where(LANDMASK != 1) prcp=prcp@_FillValue' in.nc out.nc

如果您的掩码与您的数据位于不同的网格上，您可以使用（在 Linux/Mac 上）ncremap 的掩码功能，例如，将您的数据重新映射到您的掩码（反之亦然），例如

ncremap --msk_dst=LANDMASK -d mask.nc prcp_in.nc prcp_out.nc
于 2018-03-29T12:01:22.497 回答
4

CDO 的解决方案

# make landseamask on 0.25 grid with 1 for land and missing for sea
cdo -P 8 -f nc2 setctomiss,0 -gtc,0 -remapcon,r1440x720 -topo seamask.nc

# not sure if you need to do this but lets put your data on same grid
cdo remapcon,r1440x720 yourdata.nc yourdataremap.nc

# now mask the data
cdo mul yourdataremap.nc seamask.nc yourdatamasked.nc

希望有效...

编辑：您也可以通过这种方式直接在数据网格上制作landsea mask：

# make landseamask on 0.25 grid with 1 for land and missing for sea
cdo -P 8 -f nc setctomiss,0 -gtc,0 -remapcon,yourdata.nc -topo seamask.nc

# now mask the data
cdo mul yourdataremap.nc seamask.nc yourdatamasked.nc
于 2018-03-29T10:51:56.647 回答