opposite of "enlarge"
Added by Florian Rothe about 9 years ago
Hello everyone,
is there an operator with the opposite function of "enlarge"?
I have global satellite Data but want to reduce my file to the extend of Europe.
At first I thought it can be done with "sellonlatbox" but all other variables than lon and lat remain untouched
I'm not sure if the "remap" operators are the right ones for this job...
For a short hint I would be very thankful.
Best regards
Florian
Replies (4)
RE: opposite of "enlarge" - Added by Ralf Mueller about 9 years ago
pls upload the output of ncdump -h if you file is netcdf. seems strange to me, that only coordinates get changes, this is pretty unusual for CDO
RE: opposite of "enlarge" - Added by Florian Rothe about 9 years ago
Hi Ralf,
looks like I was wrong: the other variables are also touched by the calculation.
When i plot the variable "sm" with panoply I see values outside of my selected lonlatbox. Inside the box there are no values at all.
Maybe I made a mistake with the syntax?
Anyway here is the ncdump:
netcdf file:/C:/Users/Florian/Documents/UNI/Geographie/8SemesterSS15/BA/ESA_CCI/CCI_SM_monmeans/SM_1991_eur.nc {
dimensions:
lon = 1044;
lat = 188;
time = UNLIMITED; // (12 currently)
nb2 = 2;
variables:
float lon(lon=1044);
:standard_name = "longitude";
:long_name = "longitude";
:units = "degrees_east";
:axis = "X";
float lat(lat=188);
:standard_name = "latitude";
:long_name = "latitude";
:units = "degrees_north";
:axis = "Y";double time(time=12);
:standard_name = "time";
:bounds = "time_bnds";
:units = "days since 1991-01-01 00:00:00";
:calendar = "standard";double time_bnds(time=12, nb2=2);short sm(time=12, lat=188, lon=1044);
:long_name = "Volumetric Soil Moisture";
:units = "m3 m-3";
:add_offset = 0.0f; // float
:scale_factor = 1.0E-4f; // float
:_FillValue = -9999S; // short
:missing_value = -9999S; // short
:_CoordinateAxes = "lat lon time";short sm_uncertainty(time=12, lat=188, lon=1044);
:long_name = "Volumetric Soil Moisture Uncertainty";
:units = "m3 m-3";
:add_offset = 0.0f; // float
:scale_factor = 1.0E-4f; // float
:_FillValue = -9999S; // short
:missing_value = -9999S; // short
:_CoordinateAxes = "lat lon time";byte dnflag(time=12, lat=188, lon=1044);
:long_name = "Day(100) or Night(110) Flag";
:_FillValue = 127B; // byte
:missing_value = 127B; // byte
:flag_meanings = "day night";
:_CoordinateAxes = "lat lon time";byte flag(time=12, lat=188, lon=1044);
:long_name = "Flag";
:_FillValue = 127B; // byte
:missing_value = 127B; // byte
:flag_meanings = "no_flag snow_coverage_or_temperature_below_zero dense_vegetation others_no_convergence_in_the_model_thus_no_valid_sm_estimates";
:_CoordinateAxes = "lat lon time";byte freqband(time=12, lat=188, lon=1044);
:long_name = "Frequency Band";
:_FillValue = 127B; // byte
:missing_value = 127B; // byte
:flag_meanings = "c_band c_band c_band x_band ku_band";
:_CoordinateAxes = "lat lon time";byte mode(time=12, lat=188, lon=1044);
:long_name = "Satellite Mode";
:_FillValue = 127B; // byte
:missing_value = 127B; // byte
:flag_meanings = "ascending descending";
:_CoordinateAxes = "lat lon time";byte sensor(time=12, lat=188, lon=1044);
:long_name = "Sensor";
:_FillValue = 127B; // byte
:missing_value = 127B; // byte
:flag_meanings = "SMMR SSMI TMI AMSR_E AMI_WS ASCAT SSMI+AMI_WS TMI+AMI_WS AMSR_E+AMI_WS AMSR_E+ASCAT WindSat AMSR2 WindSat+ASCAT AMSR2+ASCAT";
:_CoordinateAxes = "lat lon time";double t0(time=12, lat=188, lon=1044);
:long_name = "Observation Timestamp";
:units = "days since 1970-01-01 00:00:00 UTC";
:_FillValue = -9999.0; // double
:missing_value = -9999.0; // double
:_CoordinateAxes = "lat lon time";// global attributes:
:CDI = "Climate Data Interface version 1.6.4 (http://code.zmaw.de/projects/cdi)";
:Conventions = "CF-1.5";
:history = "Tue May 05 10:46:47 2015: cdo sellonlatbox,67,-32,81,34 ESA_CCI_1991_mm.nc SM_1991_eur.nc\nWed Apr 29 11:43:13 2015: cdo -f nc monmean -mergetime ESACCI*.nc ESA_CCI_1991_mm.nc\n2014-11-13 11:28:12 - product produced";
:source = "WARP 5.5R1.1/AMI-WS/ERS12 Level 2 Soil Moisture; WARP 5.4R1.0/AMI-WS/ERS2 Level 2 Soil Moisture; WARP 5.5R2.2/ASCAT/MetOp-A Level 2 Soil Moisture; LPRMv05/SMMR/Nimbus 7 L3 Surface Soil Moisture, Ancillary Params, and quality flags; LPRMv05/SSMI/F08, F11, F13 DMSP L3 Surface Soil Moisture, Ancillary Params, and quality flags; LPRMv05/TMI/TRMM L2 Surface Soil Moisture, Ancillary Params, and QC; LPRMv05/AMSR-E/Aqua L2B Surface Soil Moisture, Ancillary Params, and QC; LPRMv05/WINDSAT/CORIOLIS L2 Surface Soil Moisture, Ancillary Params, and QC; LPRMv05/AMSR2/GCOM-W1 L3 Surface Soil Moisture, Ancillary Params";
:institution = "Vienna University of Technology, VU University Amsterdam";
:title = "ESA CCI Surface Soil Moisture Combined Product";
:contact = "ecv_sm_contact@ipf.tuwien.ac.at";
:references = "http://www.esa-soilmoisture-cci.org; Liu, Y.Y., Dorigo, W.A., Parinussa, R.M., de Jeu, R.A.M. , Wagner, W., McCabe, M.F., Evans, J.P., van Dijk, A.I.J.M. (2012). Trend-preserving blending of passive and active microwave soil moisture retrievals, Remote Sensing of Environment, 123, 280-297, doi: 10.1016/j.rse.2012.03.014; Liu, Y.Y., Parinussa, R.M., Dorigo, W.A., De Jeu, R.A.M., Wagner, W., van Dijk, A.I.J.M., McCabe, M.F., & Evans, J.P. (2011): Developing an improved soil moisture dataset by blending passive and active microwave satellite based retrievals. Hydrology and Earth System Sciences, 15, 425-436; Wagner, W., W. Dorigo, R. de Jeu, D. Fernandez, J. Benveniste, E. Haas, M. Ertl (2012): Fusion of active and passive microwave observations to create an Essential Climate Variable data record on soil moisture. ISPRS Annals of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Volume I-7, 2012. XXII ISPRS Congress, 25 August – 01 September 2012, Melbourne, Australia";
:tracking_id = "d2708d55-5503-49ba-bc74-7d7b1961f8d3";
:product_version = "02.1";
:summary = "The data set was produced during the CCI project.";
:keywords = "Soil Moisture/Water Content";
:id = "ESACCI-SOILMOISTURE-L3S-SSMV-COMBINED-19910101000000-fv02.1.nc";
:naming_authority = "TU Wien";
:keywords_vocabulary = "NASA Global Change Master Directory (GCMD) Science Keywords";
:cdm_data_type = "Grid";
:comment = "These data were produced as part of the ESA soil moisture CCI projects. Contract No. ESRIN 22086/08/I-EC, ESRIN 4000104814/11/I-NB";
:date_created = "20141126T003957Z";
:creator_name = "Department of Geodesy and Geoinformation, Vienna University of Technology";
:creator_url = "http://rs.geo.tuwien.ac.at";
:creator_email = "ecv_sm_developer@ipf.tuwien.ac.at";
:project = "Climate Change Initiative - European Space Agency";
:geospatial_lat_min = "-90.0";
:geospatial_lat_max = "90.0";
:geospatial_lon_min = "-180.0";
:geospatial_lon_max = "180.0";
:geospatial_vertical_min = "0.0";
:geospatial_vertical_max = "0.0";
:time_coverage_start = "19781101T000000Z";
:time_coverage_end = "20131231T000000Z";
:time_coverage_duration = "P35Y";
:time_coverage_resolution = "P1D";
:standard_name_vocabulary = "NetCDF Climate and Forecast (CF) Metadata Convention";
:license = "data use is free and open for all registered users";
:platform = "Nimbus 7, DMSP, TRMM, ERS-1, ERS-2, METOP-A, AQUA, Coriolis, GCOM-W1";
:sensor = "SMMR, SSM/I, TMI, AMI-WS, ASCAT, AMSR-E, WindSat, AMSR2";
:spatial_resolution = "25km";
:geospatial_lat_units = "degrees_north";
:geospatial_lon_units = "degrees_east";
:geospatial_lat_resolution = "0.25 degree";
:geospatial_lon_resolution = "0.25 degree";
:CDO = "Climate Data Operators version 1.6.4 (http://code.zmaw.de/projects/cdo)";
}Thanks for your help
RE: opposite of "enlarge" - Added by Ralf Mueller about 9 years ago
I guess you mixed up the order of parameters
sellonlatbox,lon1,lon2,lat1,lat2 ifile ofile
is needed. Try
cdo -f nc sellonlatbox,-20,120,-60,-10 -topo selTopo.ncand have a look at it with ncview
RE: opposite of "enlarge" - Added by Florian Rothe about 9 years ago
You are right.
I really mixed up the order.
Thanks for your help!