[Proj] problem in trunk/4.9.0 in R/rgdal

Frank Warmerdam warmerdam at pobox.com
Thu Feb 6 15:30:28 PST 2014


Roger,

I have reproduced a seemingly related problem with the cs2cs commandline
program and filed a ticket:

  http://trac.osgeo.org/proj/ticket/229

I'll work on the problem this afternoon.

Best regards,
Frank



On Thu, Feb 6, 2014 at 3:11 PM, Roger Bivand <Roger.Bivand at nhh.no> wrote:

> On Thu, 6 Feb 2014, Roger Bivand wrote:
>
> > We see a problem in 4.9.0 that is not in 4.8.0 in R. The code provoking
> the
> > problem in rgdal 0.8-15 built against 4.9.0 (trunk):
> >
> > http://cran.r-project.org/src/contrib/rgdal_0.8-15.tar.gz
>
> A simpler example is:
>
> minimal install is install.packages(c("sp", "rgdal"))
>
> library(rgdal)
> checkCRSArgs("+init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84
> +no_defs +towgs84=0,0,0")
>
> failing with:
>
> > library(rgdal)
> Loading required package: sp
> rgdal: version: 0.8-15, (SVN revision 496M)
> Geospatial Data Abstraction Library extensions to R successfully loaded
> Loaded GDAL runtime: GDAL 1.11dev, released 2013/04/13
> Path to GDAL shared files: /usr/local/share/gdal
> Loaded PROJ.4 runtime: Rel. 4.9.0, 27 October 2013, [PJ_VERSION: 490]
> Path to PROJ.4 shared files: (autodetected)
> > checkCRSArgs("+init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84
> +no_defs +towgs84=0,0,0")
> [[1]]
> [1] FALSE
>
> [[2]]
> [1] "no options found in 'init' file"
>
> Roger
>
> >
> > minimal install is install.packages(c("sp", "spacetime", "rgdal"))
> >
> > library(sp)
> > library(spacetime)
> > data(air)
> > library(rgdal)
> > checkCRSArgs(proj4string(rural at sp))
> >
> > which shows:
> >
> >> checkCRSArgs(proj4string(rural at sp))
> > [[1]]
> > [1] FALSE
> >
> > [[2]]
> > [1] "no options found in 'init' file"
> >
> > first time round. The object rural is read in from a saved R object, and
> has
> > a stored projection of:
> >
> > "+init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs
> > +towgs84=0,0,0"
> >
> > So when this string goes via checkCRSArgs() to C code (same function
> name)
> > calling:
> >
> >
> > if (!(pj = pj_init_plus(CHAR(STRING_ELT(args, 0))))) {
> >       SET_STRING_ELT(VECTOR_ELT(res, 1), 0,
> >        COPY_TO_USER_STRING(pj_strerrno(*pj_get_errno_ref())));
> >       UNPROTECT(1);
> >       return(res);
> > }
> >
> > where res is a list with first element FALSE by default, second contains
> an
> > error message on failure (or what pj_get_def() returns on success. This
> > around line 292 in:
> >
> >
> https://r-forge.r-project.org/scm/viewvc.php/pkg/src/projectit.cpp?view=markup&root=rgdal
> >
> > However, if we then do:
> >
> > checkCRSArgs("+init=epsg:4326")
> >
> > we tickle something in pj_init_plus() and it works, and subsequently:
> >
> > checkCRSArgs(proj4string(rural at sp))
> >
> > works. I cannot see any changes in pj_init.c in the section looking for
> +init
> > between 4.8.0 and 4.9.0 that might affect this. Under rgdal 0.8-15 and
> 4.8.0,
> > there is no trouble, but the issue arises predictably with the same rgdal
> > 0.8-15 and current 4.9.0 trunk.
> >
> > I'd be grateful for any advice.
> >
> > Best wishes,
> >
> > Roger
> >
> >
>
> --
> Roger Bivand
> Department of Economics, Norwegian School of Economics,
> Helleveien 30, N-5045 Bergen, Norway.
> voice: +47 55 95 93 55; fax +47 55 95 95 43
> e-mail: Roger.Bivand at nhh.no
>
> _______________________________________________
> Proj mailing list
> Proj at lists.maptools.org
> http://lists.maptools.org/mailman/listinfo/proj
>



-- 
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam,
warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Software Developer
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20140206/85f9a701/attachment.html>


More information about the Proj mailing list