problem with function project() in package rgdal

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

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?

Drs. Paul Hiemstra
Department of Physical Geography
Faculty of Geosciences
University of Utrecht
Heidelberglaan 2
P.O. Box 80.115
3508 TC Utrecht
Phone:  +3130 274 3113 Mon-Tue
Phone:  +3130 253 5773 Wed-Fri
http://intamap.geo.uu.nl/~paul