-
Notifications
You must be signed in to change notification settings - Fork 353
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Let mapproject also find encompassing rectangular regions for non-oblique projections #6669
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
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
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! |
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 -WE
and get his
-R-1.852421104674526/62.33269965938534/126/80.96267731894207+r