a minor bug in rgdal?
Dear Prof. Bivand,
I am running R-2.15.2 with rgdal version 0.7-22.
I noticed the following bug: when transforming a SpatialPoints object,
a white space is prepended to the proj4string of the target CRS. For
example:
crds <- matrix(data=c(9.05, 48.52), ncol=2)
spPoint <- SpatialPoints(coords=crds,
proj4string=CRS("+proj=longlat +ellps=sphere +no_defs"))
a <- spTransform(spPoint, CRS(paste("+proj=ob_tran +o_proj=longlat",
"+o_lon_p=-162 +o_lat_p=39.25 +lon_0=180 +ellps=sphere +no_defs")),
use_ob_tran=TRUE)
proj4string(a)
# " +proj=ob_tran +o_proj=longlat +o_lon_p=-162 +o_lat_p=39.25
+lon_0=180 +ellps=sphere +no_defs"
As You can see, a white space is prepended to the proj4string, from
"+proj=ob_tran ..."
it becomes " +proj=ob_tran ...". This leads to errors when one
interpolates using krige from gstat for example.
Best regards,
Dr. Martin Ivanov Eberhard-Karls-Universit?t T?bingen Mathematisch-Naturwissenschaftliche Fakult?t Fachbereich Geowissenschaften Water & Earth System Science (WESS) H?lderlinstra?e 12, 72074 T?bingen, Deutschland Tel. +4970712974213