#!/usr/bin/env sh
gmt begin aditya1 png
gmt basemap -R73.4/81/28/35 -JM7i -Ba2f2
gmt grdimage best_1s.tif -Itopo_adjust.grd -Cgray99 -E600 -V
gmt grdmask digitizedpolyjksh.dat -R73.4/81/28/36.7 -I1s -NNaN/1 -Gmask.grd=nb -fg -V
gmt grdmath resampled_siteupdate1.grd mask.grd MUL = resampled_siteupdate1_masked.grd
gmt grdmath coupling_updatedmap.grd mask.grd MUL = coupling_updatedmap_masked.grd
gmt grdimage coupling_updatedmap_masked.grd -Ccommon.cpt -E600 -V
gmt end show
digitizedpolyjksh.dat (1.1 KB)
i have uploaded the .dat file. All i want is to remove the polygon or mask it but i dont understand why the polygon is coming like this. please help
The image you posted first cannot have been produced by your script (that should have been wrapped in triple back ticks ``` to facilitate reading) that has no calls to text or velo). So, I don’t even understand what the problem is.
So far it appears the topic starter wasn’t the one who wrote the code, and the code wasn’t the one that produced the plot.
I think this is a really difficult situation.
#!/usr/bin/env sh
gmt begin test png
gmt basemap -R73.4/81/28/35 -JM7i -Ba2f2
gmt grdgradient best_1s.tif -A135 -Ne0.8+o0.04+a0.05 -Gtopo_adjust.grd -V
gmt grdimage best_1s.tif -Itopo_adjust.grd -Cgray99 -E600 -V
gmt grdmask digitizedpolyjksh.dat -R73.4/81/28/36.7 -I1s -N1/0/NaN -Gmask.grd=nb -fg -V
gmt grdmath resampled_siteupdate1.grd mask.grd MUL = resampled_siteupdate1_masked.grd
gmt grdmath coupling_updatedmap.grd mask.grd MUL = coupling_updatedmap_masked.grd
gmt grdimage resampled_siteupdate1.grd -R73.4/81/28/35 -Icoupling_updatedmap.grd -Ccommon.cpt -E600 -Q -V
gmt end show
sorry for the confusion but this is the script for recent figure uploaded.
the above line created mask. If you don’t want the mask, you don’t need that code line.
the above code lines applied the mask to the resampled_siteupdate1.grd and coupling_updatedmap.grd
if you don’t want to do masking, don’t run these gmt math code lines.
NB In those gmt math code lines that apply the mask both source and destination file names are identical. Every time these lines are executed, the resulting masked grids overwrite the source non-masked files. If you want to do masking and to keep the source grids, give the output files different names.
