How to use grdcut ? I have tried but the doesn’t work

Q&A
#1

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

gambar
gambar1102×390 20.8 KB

but the Area.tif file doesn’t exist, please help
i am using gmt version 6.4.0

#2

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.

#3

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

#4

Given that you cutting a region of about 2x2 km of an image with > 100 resolution, the result doesn’t seam unreasonable.