[R-pkg-devel] how to call PROJ.4 C code in a package?

Tue, May 26, 2015 2:53 PM

Roger,

I plan to make the next version of oce use rgdal::project() to do its work. The idea of using an external PROJ.4 is fine in my own work, but oce has a lot of users who are unaccustomed to building and installing external libraries. Also, many work at institutes that do not permit this kind of system modification. So, for oce I had to make a choice of including the C sources or using another package.

At first I felt the proj4 R package to be a natural candidate. However, it produces an error when any point in a lon-lat vector is ?offscale?, instead of inserting NA for the projected value. Indexing over a world coastline is too slow for practical work.

Over the weekend I set up a test version of oce to use rgdal::project(), and it?s early days but things seem to be working well.

I need to study the ?sp? package in more detail.  Maybe I?ll find the solution to a problem ?oce? has when lines in a coastal polygon ?cross over? the edge of the earth. This produces spurious domain-crossing lines when e.g. lon_0 is set to put the Pacific at the centre (a common choice of oceanographers).

Thanks VERY much for your detailed and helpful answer.  And thanks for your work on ?rgdal? and ?sp?, two very fine packages.

? Dan.

PS. glad to see more software going to github.

Dan Kelley
Oceanography Department
Dalhousie University
Halifax, NS, Canada

On May 26, 2015, at 4:14 PM, Roger Bivand <Roger.Bivand at nhh.no> wrote:

Dan,

I'm sorry to be late to the party (one is offline occasionally), but all the advice given so far is at best only partially correct.

Use of PROJ.4 is crucially conditional on access to PROJ_LIB a shell variable setting the address of metadata files not included in the C source. PROJ.4 without this access is very limited, not least because the EPSG listing of projections is not available, not to mention a small file of default values that recently floored rgdal and other packages using its PROJ.4 facilities (PROJ.4 4.9.1 omitted the file by mistake).

At present, oce does not ship with PROJ_LIB at all. I agree that it does not need to ship with the PROJ.4 C source, and should, like rgdal and proj4, use the system PROJ.4 (or the CRAN PROJ.4). Only rgdal checks in ./configure for settings related to PROJ_LIB. Using PROJ.4 externally makes sense because of version shifting and missed bug fixes (there were some bad bugs before the aborted 4.9.0 release), and crucially updated metadata files. Using a single PROJ.4 on a platform avoids different applications seeing possibly different metadata.

My guess would be that anything you might need to do can be done using function stubs exported by rgdal for use in raster (the package). You are free to choose not to use sp objects, but you'll find that the high-level support for these in rgdal is robust and well-motivated.

If you go for the proj4 rather than rgdal solution, do make sure that the CRAN Windows and OSX binaries are built with the metadata PROJ_LIB copied to the binary packages.

Again, apologies for not jumping in at the right time.

Roger

PS. Please also note that PROJ.4 is moving to

https://github.com/OSGeo/proj.4

--
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 91 00
e-mail: Roger.Bivand at nhh.no

[R-pkg-devel] how to call PROJ.4 C code in a package?

Thread (3 messages)