create polygon maskfile
Added by kunal bali over 4 years ago
Dear CDO users,
I need to create a one mask/or shape file of the given desired polygon box [see in the attached figure] from the attached .nc file.
The coordinates of the given polygon are
lonv = (/77.98,77.96,79.01,80.02,80.98,81.57,81.98,83.01,84.00,83.99,83.01,83.01,80.98,81.00,80.01,79.64,79.63,79.01,77.98/)
latv = (/31.00,29.70,29.15,28.82,28.43,28.01,27.98,27.47,27.38,28.64,28.99,29.40,29.56,29.99,30.05,30.01,30.61,30.64,31.00/)
Could you please provide some information on that.
Thank You
Replies (8)
RE: create polygon maskfile - Added by Brendan DeTracey over 4 years ago
RE: create polygon maskfile - Added by kunal bali over 4 years ago
Dear Brendan,
Thanks for your suggestion. I have tried this already.
This function creates only rectilinear box. But, in my case my area is different as you could see in the the attached figure. So, if you have any other suggestion then please let me as well.
Thank You
RE: create polygon maskfile - Added by Brendan DeTracey over 4 years ago
Welp, it would easier to help if you posted the command line you tried as well as upload your ascii region file.
Note that maskregion
documentation states:
The elements inside a region are untouched, the elements outside are set to missing value.
The output file has the same dimensions but the points outside the polygon are masked with missing values.
RE: create polygon maskfile - Added by kunal bali over 4 years ago
the function takes only 4 coordinates lat1,lat2,lon1,lon2. So it mask only square/rectangular box.
And in my case, I have more than 4 coordinates and not in rectangular box (as mentioned in attached figure). It's irregular shape. So this function can't be applicable as far as I understood.
RE: create polygon maskfile - Added by Brendan DeTracey over 4 years ago
maskbox takes box corners, maskregion takes a plain text file with coordinates....
RE: create polygon maskfile - Added by kunal bali over 4 years ago
I tried
cdo mskregion,myregion.txt file.nc out.nc
file created with NaN values only.
myregion.txt is
77.98,77.96,79.01,80.02,80.98,81.57,81.98,83.01,84.00,83.99,83.01,83.01,80.98,81.00,80.01,79.64,79.63,79.01,77.98 31.00,29.70,29.15,28.82,28.43,28.01,27.98,27.47,27.38,28.64,28.99,29.40,29.56,29.99,30.05,30.01,30.61,30.64,31.00
RE: create polygon maskfile - Added by Brendan DeTracey over 4 years ago
Carefully read the description for the polygon text file. You have not made it correctly. They have an example in the manual.