Nearneighbor

Added by Francois 8 months ago

Hi,

I am using GMT 5.3.3. I have a file.xyz with boundaries of -800000/700000/-3400000/-600000 with a resolution of 5 km. I am trying to change its resolution by using:

blockmean $1 -R-800000/700000/-3400000/-600000 -I1000 -V > $2

nearneighbor $2 -R-800000/700000/-3400000/-600000 -I1000 -V -E0 -G$3 -S50000 -N4

However, I do have this error:

bash-3.2$ ./regrid.sh GrClimBamber.xyz GrClimBamberRg GrClimBamberRg.xyz
blockmean: Processing input table data
blockmean: W: -800000 E: 700000 S: -3400000 N: -600000 n_columns: 1501 n_rows: 2801
blockmean: Calculating block means
blockmean: N read: 168861 N used: 1681 N outside_area: 167180 N cells filled: 1681
nearneighbor: Grid dimensions are n_columns = 1501, n_rows = 2801
nearneighbor: Number of sectors = 4, minimum number of filled sectors = 2
nearneighbor: Processing input table data
nearneighbor: Processed record 1681
nearneighbor: Gridded row 2801
nearneighbor: 60001 nodes were assigned an average value
nearneighbor: 2091 nodes failed sector criteria and 4142209 nodes had no neighbor points (all set to 0)

I could change either -I or -S values, but I still have the same error...

What do I miss?!

Thank you!

François

GrClimBamber.xyz (3.76 MB)

regrid.sh Magnifier (666 Bytes)


Replies (2)

RE: Nearneighbor - Added by Paul 8 months ago

According to gmt info:
gmt info GrClimBamber.xyz

GrClimBamber.xyz: N = 168862    <-3400000/-600000>    <-800000/700000>    <-2693.357/4106.087>

so you are mixing up what is x and y here. Using -R-3400000/-600000/-800000/700000 is what you want.

RE: Nearneighbor - Added by Francois 8 months ago

Yep... it was a mixed up between x/y... Thanks.

(1-2/2)