Dear list,

I am trying to change the projection of a SpatialPolygonsDataFrame but
there is something wrong going on that I can't figure. The projected map
show some lines that shouldn't be there.

The shape file was downloaded from here:
http://macroecology.ku.dk/resources/wallace/#Gis

Here is the very simple code I'm using:

setwd(PATH)
newrealms<-readShapePoly("newRealms.shp")
proj4string(newrealms) <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84
+no_defs")
moll_crs<-CRS("+proj=moll +ellps=WGS84")
newrealms_moll <- spTransform(newrealms, moll_crs)
plot(newrealms)
plot(newrealms_moll) #unexpected lines

Can anyone help me with this?

Thanks in advance!

Daniel

        [[alternative HTML version deleted]]

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Reply via email to