Re: Azimuthal Equidistant projection problem
Paul Wessel <[email protected]> Mon, 26 May 2014 15:28:48 -1000
| Newsgroups | gmane.comp.gis.gmt.user |
|---|---|
| Message-ID | <[email protected]> |
Hi Juan- Yes, there were in fact two old bugs here, affecting both GMT 4 & 5: 1. The clip path for the entire region was hardwired to be a circle when -Je/E was used, regardless of -R. In effect, it assumed -Rg. This caused the blue circle. 2. The search that determines the min/max w/e/s/n for a region specified by -R...r did not check if either the S or N pole happened to be the antipode of the -Je/E pole [your case]. That caused a bad w/e/s/n setting and that gridline did not happen. Updated in both GMT 4 (r10233) and GMT 5 (r13201); please update and check if your other issues further down (psxy/grdimage) still need attention. If so please file a bug report on the tracker. Paul On May 24, 2014, at 5:56 PM, Juan Diaz-Naveas <[email protected]> wrote: > Dear GMT users, > > This might be a bug, but before reporting as that I want to share the strange behavior of the Azimuthal Equidistant projection. > > > The followind command (when -R is specififed as -Rlonmin/lonmax/latmin/latmax) produces an output, but the -Sblue option draws also a strange circle: > > pscoast -B10g10 -Je-70/-90/1:15000000 -R-95/-60/-75/-55 -Di -Ggreen -Sblue -P -Uc -V > az_eq_scale.ps > > The next command (when -R is specified as -Rlonmin/latmin/lonmax/latmaxr) behaves ALMOST correctly, because the 60°W meridian is NOT drawn: > > pscoast -B10g10 -Je-70/-90/1:15000000 -R-95/-75/-60/-55r -Di -Ggreen -Sblue -P -Uc -V > az_eq_scale_r.ps > > I attach both postscript files to let you verify this. > > Actually the problem becomes bigger when you add psxy or grdimage commands, but for now I want to solve this very basic issue. > > Does anyone have an idea on how to solve this problem? > > Any help will be very much appreciated. Best wishes > > Juan Diaz-Naveas > > Mailing list for GMT discussions of all kinds. If you are not sure you have found a bug, discuss it here first. To formally report bugs or request features, please register and add New Issue on gmt.soest.hawaii.edu To unsubscribe, send the message "signoff gmt-help" to [email protected] Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead. <az_eq_scale.ps><az_eq_scale_r.ps> Mailing list for GMT discussions of all kinds. If you are not sure you have found a bug, discuss it here first. To formally report bugs or request features, please register and add New Issue on gmt.soest.hawaii.edu To unsubscribe, send the message "signoff gmt-help" to [email protected] Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead.