Dear Paul,
Thank you for you reply. Everything works fine now with the current
versions of R, rgdal and sp:
> library(rgdal)
Le chargement a nécessité le package : sp
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 1.6.0, released 2008/12/04
Path to GDAL shared files: C:/R/RLIBS/rgdal/gdal
Loaded PROJ.4 runtime: Rel. 4.6.1, 21 August 2008
Path to PROJ.4 shared files: C:/R/RLIBS/rgdal/proj
> data(state)
> res <- project(cbind(state.center$x, state.center$y), "+proj=lcc
+lat_1=48 +lat_2=33 +lon_0=-100")
> sessionInfo()
R version 2.9.0 (2009-04-17)
i386-pc-mingw32
locale:
LC_COLLATE=French_France.1252;LC_CTYPE=French_France.1252;LC_MONETARY=French_France.1252;LC_NUMERIC=C;LC_TIME=French_France.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] rgdal_0.6-8 sp_0.9-36 fortunes_1.3-6
loaded via a namespace (and not attached):
[1] grid_2.9.0 lattice_0.17-22 tools_2.9.0
I don't know what happened...
Best regards,
Renaud
Paul Hiemstra a écrit :
Hi Renaud,
I don't get the same error, the piece of code runs fine. To which
version of proj did you build rgdal? My versions are:
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 1.4.4.0, released 2007/11/23
Path to GDAL shared files: /usr/share/gdal
Loaded PROJ.4 runtime: Rel. 4.6.0, 21 Dec 2007
Path to PROJ.4 shared files: (autodetected)
Googleing for "major axis or radius = 0 or not given" will give you some
useful things to start with.
As a remark, I always transform my spatial data to sp-classes (using
coordinates() or gridded(), see sp-docs) and then use spTransform to
perform the reprojection.
cheers,
Paul
This is my sessionInfo():
R version 2.9.0 (2009-04-17)
i486-pc-linux-gnu
locale:
LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=C;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_US.UTF-8;LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] rgdal_0.6-5 sp_0.9-29
loaded via a namespace (and not attached):
[1] grid_2.9.0 lattice_0.17-20 tools_2.9.0
Renaud Lancelot wrote:
Dear all,
I meet a problem with the function project() in package rgdal:
library(rgdal)
data(state)
res <- project(cbind(state.center$x, state.center$y),
+ "+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100")
Error in project(cbind(state.center$x, state.center$y), "+proj=lcc
+lat_1=48 +lat_2=33 +lon_0=-100") :
major axis or radius = 0 or not given
sessionInfo()
R version 2.8.1 Patched (2009-02-17 r47956)
i386-pc-mingw32
locale:
LC_COLLATE=French_France.1252;LC_CTYPE=French_France.1252;LC_MONETARY=French_France.1252;LC_NUMERIC=C;LC_TIME=French_France.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods
[7] base
other attached packages:
[1] rgdal_0.6-7 sp_0.9-34 fortunes_1.3-6
loaded via a namespace (and not attached):
[1] grid_2.8.1 lattice_0.17-20 tools_2.8.1
Is this a bug?
--
Renaud Lancelot
EDEN Project, coordinator
http://www.eden-fp6project.net/
UMR CIRAD-INRA "Contrôle des maladies animales exotiques et émergentes"
Joint research unit "Control of emerging and exotic animal diseases"
CIRAD
Campus International de Baillarguet TA A-DIR / B
F34398 Montpellier
http://www.cirad.fr http://bluetongue.cirad.fr/
Tel. +33 4 67 59 37 17 - Fax +33 4 67 59 37 95
Secr. +33 4 67 59 37 37 - Cell. +33 6 77 52 08 69
_______________________________________________
R-sig-Geo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
