[Proj] GitHub "project" for the next release
richard.greenwood at gmail.com
Thu Aug 10 08:31:46 EST 2017
Computationally datum transformations are completely separable from map
projections, but from the end users perspective they are inextricably
combined in the modern era. A user should reasonably expect to get the same
results in QGIS and ArcMap given the same ESPG codes, but in many cases
they do not because of datum differences. I know Gerald disagreed with me
on this in the past and obviously others still disagree with me. I do not
believe that moving Proj into the 21st century threatens users who do not
care about datum transformations. But ignoring datum transformations will
make Proj irrelevant in the future.
On Thu, Aug 10, 2017 at 5:44 AM, Thomas Knudsen <knudsen.thomas at gmail.com>
> Duncan (and daan),
> Yes, map projections are mathematical functions, while reference frame
> shifts are geophysical/geodetical transformations.
> The latter are, however, also *in practice* represented by sequential
> application of parameterized mathematical functions, drawn from a rather
> limited gamut - basically consisting of Helmert, Horner, Molodensky, Grid
> interpolation and ellipsoidal/cartesian transformation.
> All of these fundamental operations are architecturally similar to the map
> projection functions of the original libproj, and are implemented in a
> similar way.
> Hence, what Kristian Evers and I have done over the last 18 months is
> technically to add 5 (or so) new projections to the original libproj.
> To string the stuff together, we designed the pipeline metaphor, which is
> also materialized as just another projection.
> So we have added what amounts to a few more projections to the library.
> These can be used or ignored, in the same way you may ignore the Transverse
> Mercator implementation if you are only interested in the Lambert Conformal
> Conic (You can probably even compile a version of the library leaving most
> of the “pseudoprojections” out).
> All parameters for implementing actual reference frame transformations
> *are stored in parameter files, not in the library per se*.
> If you are *not interested* in reference frame shifts, the code supporting
> reference frame shifts will *not affect you*. And the parameter files will
> only affect you through the (limited) disk space required. But you are not
> forced to install them anyway.
> So if you are *not interested* in reference frame shifts, please *ignore
> the matter* and stop arguing against stuff that will *not affect you*. This
> stuff is certainly useful for a range of other people.
> Should we want to follow daan’s recipe of splitting the reference frame
> handling entirely from the projection library, we would have to duplicate
> most of the well tested libproj infrastructural architecture into another
> library, which would, in most cases, be used in applications also depending
> on libproj.
> Simply leveraging the existing infrastructure and adding a handful of
> pseudoprojections was definitely the way to go. Especially since the
> libproj architecture makes it possible to do so in a neat and non-intrusive
> Regarding the new API:
> The new API emerged as a by-product of attempts to make the libproj
> internals more maintainable, by improving namespacing and introducing new
> internal data types in order to reduce type-punning. This makes it more
> clear exactly what kind of coordinate is expected at each spot in the code.
> The new API exposes these data types, simplifies a wide field of
> operations, and provides a unified interface to 2D, 3D and 4D
> The new API does not replace the old API, but obviously the designers hope
> that it will become the preferred one for new applications. The old API may
> some day gain a “deprecated” status, but will only rot entirely away if no
> one is interested in maintaining it. So if anyone needs it badly enough to
> sponsor its maintenance, it will live on forever.
> In a limited number of cases, software pre-dating the original API,
> described in proj-api.h, may experience breakage, since the original
> projects.h file does not describe an actual interface (API), but rather
> exposes every single function and data structure in the library.
> So there is only one way to ensure that each-and-every program, using
> projects.h to do calls to libproj internals, would keep working with new
> versions of libproj: That would be to freeze both development and
> maintenance of the PROJ system. And anyone wishing for that may get theirs
> by just sticking with proj 4.3.0.
> I personally believe, though, that we should make efforts not to hamper
> more well disciplined programs using projects.h just to do calls to the
> “obvious API” of projects.h (pj_init, pj_fwd, pj_inv, pj_free, and probably
> a few more).
> The proj and cs2cs filters will live on as long as anyone wants to
> maintain them. But they definitely need some tlc, so test cases involving
> some of your more esoteric use cases of proj and cs2cs would be very
> welcome, and would greatly reduce the risk of regressions during
> 2017-08-10 6:56 GMT+02:00 Duncan Agnew <dagnew at ucsd.edu>:
>> I pretty much agree with Dan. For clarity (lacking in my previous entry),
>> there are two things, easily conflated, to talk about:
>> A. proj the program that people use (often on the command line)
>> to go between lat-long and projection coordinates. Most of these users
>> don't care about which lat/long they are using: certainly not important
>> if making a map, and probably not if using it to set up a grid system
>> for your own use (I've used it for both).
>> B. PROJ the package, which includes not just projections
>> but all the additional machinery to handle different datums--and
>> in the future heights and temporal changes.
>> I completely agree that PROJ needs to keep evolving; but at the
>> same time I (and probably many other users) would like to have proj
>> stay the same. This ought to be relatively easy since it is almost all
>> mathematics, the only parameters being the size and shape of the
>> ellipsoid--and I don't imagine there will be new values of this being
>> produced. But it depends on allowing the functionality of proj (not
>> PROJ) to not require the user who just wants a projection to have to
>> master the additional features needed to do datum transformations.
>> Of course grids are just projections, so it was easy for Gerald
>> to add things like SPCS, blurring the line a bit; but, with Daan (and
>> Gerald) I'd argue that it would be helpful to many users to have things
>> somewhat distinct.
>> And, I'd like to raise a question (I'm genuinely curious) about how
>> time-dependent motion will be added. Having been professionally engaged in
>> measuring and modeling crustal motion for the last 30+ years, I can say
>> that it
>> is going to be a lot more complicated to keep track of. proj (as above)
>> only the -le parameters (a fixed set); once you add the rigid-body
>> motions for
>> datum transformations there are a lot more parameters, many conventional
>> so unchanging), which EPSG fortunately collects; adding grids for
>> (HARN) needs more information yet. But once you get into temporal changes,
>> things are even worse: models will keep changing either because of better
>> or earthquakes (cf HTDP). This seems like a progression from a package
>> dominated by algebra to a small algebraic component attached to a large
>> and ever-growing database of parameters. I'm not arguing against it, at
>> all, just saying that it is going to be quite a challenge beyond modifying
>> the programs. (I completely lack the expertise to contribute to the code,
>> but if anyone has any geophysics questions I'm happy to try to answer
>> Duncan Agnew
>>> Proj mailing list
>>> Proj at lists.maptools.org
>> Proj mailing list
>> Proj at lists.maptools.org
> Proj mailing list
> Proj at lists.maptools.org
Richard W. Greenwood, PLS
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Proj