Azimuthal Equidistant projection problem

[Juan Diaz-Naveas]

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 ***@lists.hawaii.edu
Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead.

[Paul Wessel]

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

Post by Juan Diaz-Naveas
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.
pscoast -B10g10 -Je-70/-90/1:15000000 -R-95/-60/-75/-55 -Di -Ggreen -Sblue -P -Uc -V > az_eq_scale.ps
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 ***@lists.hawaii.edu
Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead.

[Juan Diaz-Naveas]

Dear Paul,

Thank you very much indeed. Very best wishes

Juan

Post by Paul Wessel
Hi Juan-
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

Post by Juan Diaz-Naveas
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.

Post by Juan Diaz-Naveas
The followind command (when -R is specififed as

-Rlonmin/lonmax/latmin/latmax) produces an output, but the -Sblue option

Post by Juan Diaz-Naveas
pscoast -B10g10 -Je-70/-90/1:15000000 -R-95/-60/-75/-55 -Di -Ggreen

-Sblue -P -Uc -V > az_eq_scale.ps

Post by Juan Diaz-Naveas
The next command (when -R is specified as

-Rlonmin/latmin/lonmax/latmaxr) behaves ALMOST correctly, because the 60èš–

Post by Juan Diaz-Naveas
pscoast -B10g10 -Je-70/-90/1:15000000 -R-95/-75/-60/-55r -Di -Ggreen

-Sblue -P -Uc -V > az_eq_scale_r.ps

Post by Juan Diaz-Naveas
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.

Post by Juan Diaz-Naveas
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
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
Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on
gmt.soest.hawaii.edu instead.

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 ***@lists.hawaii.edu
Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead.