[R-sig-Geo] Azimuthal equidistant projection and rgdal

[Johan Mazel]

Here is the code that I am using and the result.
The problem on the result is that, on the fourth and last page, it
looks like some parts of the map are not located where they are
supposed to be (Australia is thus right next to Canada).

My only modification is that I removed the xlim and ylim parameters in
the line "dta <- map("world", xlim=c(-120, -60), ylim=c(0,90),
plot=FALSE)" in the original script. It obviously has a impact on the
final output but I can`t figure out what is happening.

Thanks for your time.
Regards.
Johan

2012/10/25 Johan Mazel <johan.mazel at gmail.com>:
> Hi
> I am trying to use the solution rpesented here:
> https://stat.ethz.ch/pipermail/r-sig-geo/2007-December/002939.html.
>
> However it looks like the shifted (from the north pole) map is
> displayed on top on the original azimuthal equidistant projection
> centered on the north pole.
> Is there any way to correct this ?
>
> Thanks for your time.
> Regards.
> Johan
-------------- next part --------------
A non-text attachment was scrubbed...
Name: projection_azeq.r
Type: application/octet-stream
Size: 389 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20121112/88a1cf4d/attachment.obj>