Gmt spatial -Q[unit]

Q&A
1

Hi,

I’ve noticed that when I use gmt spatial -Q[unit], it seems to give me the area in km^2 (I think?) no matter what unit I choose. I’m wondering what the expected behaviour is supposed to be? If I use gmt spatial -Qd, I’m not sure I should expect an answer in degrees squared, should I?

cat << EOF > poly.tmp
0 0
0 1
1 1
1 0
0 0
EOF

gmt spatial poly.tmp -Q
gmt spatial poly.tmp -Qk
gmt spatial poly.tmp -Qe
gmt spatial poly.tmp -Qd

Gives the results

0.5	0.5	1
0.5	0.50000702007	12281.2810988
0.5	0.50000702007	12281.2810988
0.5	0.50000702007	12281.2810988

I’m using gmt v6.2.0.

Is this expected?

2

Hi tmerry. Yes, you are right. It is a bug and it was fixed 10 days ago. You should use the dev version (GMT 6.4)

1 Like
3

Ah okay, thanks! I didn’t think to check pull requests!