Let mapproject also find encompassing rectangular regions for non-oblique projections#6669
Let mapproject also find encompassing rectangular regions for non-oblique projections#6669PaulWessel merged 3 commits intomasterfrom
Conversation
See this forum post for background. What is needed is -We or -WE, hereby implemented. Now, one can run gmt mapproject -JS36/90/30c -R-15/60/68/90 -WE and get -R-1.852421104674526/62.33269965938534/126/80.96267731894207+r
anbj
left a comment
anbj
left a comment
There was a problem hiding this comment.
It's like Christmas, just in May.
|
This is great Paul, thanks alot. I find it strange that I'm the first one(?) to nag enough about this to get it implemented. Guess confusion can sometimes lead to great things. |
|
Note to self and a 'fun' observation; I see that the rounding off is different between the two methods |
|
Sometimes one finds oneself on a train with nothing to do but hack away. |
This is because the -Wn method involves a multiplication of the rectangular extent by the factor "1" which in a computer is not 1 but perhaps 1 + 2e-12 and hence that extra multiplication will mess with insignificant digits down the line. |
|
Thanks! |
2 participants
See this forum post for background. What is needed is -We or -WE, hereby implemented. Now, @anbj can run
gmt mapproject -JS36/90/30c -R-15/60/68/90 -WEand get his
-R-1.852421104674526/62.33269965938534/126/80.96267731894207+r