I’m trying to crop the region with area 109.865/109.885/-7.385/-7.365 then run the
grdcut function with a script like this
gmt grdcut -R109.865/109.885/-7.385/-7.365 Final.tif -GArea.tif -V
then it looks like this
but the Area.tif file doesn’t exist, please help
i am using gmt version 6.4.0
Hi, Nurul
Two things. Please don’t use shell screen captures. Just copy paste the text and wrap it under triple back tics.
If you didn’t get any errors then you have the result file. Sometimes you have to hit F5 for windows to update the file manager list.
Thanks for the advice Mr Joaquim, I think there is an error because the output file size is only 1 KB, maybe it’s like that
e:\Data di PC Windows 10\Coding\DEM\temp>gmt grdcut -R109.865/109.885/-7.385/-7.365 Final.tif -GArea.tif
grdcut [WARNING]: (w - x_min) must equal (NX + eps) * x_inc), where NX is an integer and |eps| <= 0.0001.
grdcut [WARNING]: w reset from 109.865 to 109.864166667
grdcut [WARNING]: (e - x_min) must equal (NX + eps) * x_inc), where NX is an integer and |eps| <= 0.0001.
grdcut [WARNING]: e reset from 109.885 to 109.885833333
grdcut [WARNING]: (s - y_min) must equal (NY + eps) * y_inc), where NY is an integer and |eps| <= 0.0001.
grdcut [WARNING]: s reset from -7.385 to -7.38583333333
grdcut [WARNING]: (n - y_min) must equal (NY + eps) * y_inc), where NY is an integer and |eps| <= 0.0001.
grdcut [WARNING]: n reset from -7.365 to -7.36416666667
Input file size is 3000, 3000
0
[Done] exited with code=0 in 0.764 seconds
Given that you cutting a region of about 2x2 km of an image with > 100 resolution, the result doesn’t seam unreasonable.