[Proj] What about datum shift via direct projection?
Mikael Rittri
Mikael.Rittri at carmenta.com
Thu Dec 11 07:43:08 PST 2008
Gerald wrote:
> While there seems to be lull in the hot debate about separation of
> church and state ... er ... datum and projection,
> [...]
> Thus, why is it so necessary to bind the two operations so tightly as done
> in the proj.4 distribution? I cannot find a precedence for this concept.
This post is not specifically about the PROJ.4 design (so I changed
the Subject line), but it is about how much datums and projections
can and should be separated.
There is method for datum shift that uses a direct projection.
As an example, the old Swedish Grid is traditionally defined
on the Swedish RT90 datum (ellipsoid: Bessel 1841) and using a
Transverse Mercator projection with
central meridian: 15° 48' 29.8" E
scale factor: 1
false easting: 1500000 m
false northing: 0 m
( http://www.lantmateriet.se/templates/LMV_Page.aspx?id=4766&lang=EN )
With this definition, one would need some datum shift method
to transform between RT90 lon/lat and WGS84 lon/lat.
However, a simpler method, now recommended by the Swedish Land Survey
instead of a 7-parameter shift, is to start from the WGS84 datum, and than
tweak the projection parameters a little: just use a Transverse Mercator
with
central meridian: 15° 48' 22.624306" E
scale factor: 1.00000561024
false easting: 1500064.274 m
false northing: -667.711 m
( http://www.lantmateriet.se/templates/LMV_Page.aspx?id=5197&lang=EN )
A paper describing this technique is
http://www.fig.net/pub/fig2006/papers/ps05_03/ps05_03_04_engberg_lilje_0670.pdf .
So, I have some rather vague questions to the readers of this list:
- What do you think of this technique?
- Is anyone else using it?
- Doesn't the technique imply that a projected coordinate system
may have an ambiguous geographic coordinate system? For the Swedish Grid,
I can think of the geographic coordinate system as RT90 lon/lat, if I use
the traditional projection parameters. Or I can think of it as WGS84 lon/lat,
if I use the direct projection instead.
- If the correct answer to the previous question is "No, you fool", then what?
If I wanted to express the Swedish Grid, datum-shifted by the direct projection,
in Well-Know Text, then I would be forced to say that the geographic coordinate
system is WGS84 lon/lat. But then the resulting CRS cannot be Swedish Grid,
because Swedish Grid has traditionally RT90 lon/lat as its geographic coordinate
system.
I think direct projections for datum shifts are efficient and easy to
use, and normally as accurate as a 7-parameter shift. But when I try
to fit this method into the traditional framework that separates datum
shifts and projections, and which insists that each projected CRS
has a unique geographic coordinate system, I run into problems.
Are these problems caused by inflexibility in the traditional framework?
Or is the method of direct projection just weird?
Or am I missing some good way to reconcile them?
Best regards,
--
Mikael Rittri
Carmenta AB
Box 11354
SE-404 28 Göteborg
Visitors: Sankt Eriksgatan 5
SWEDEN
mikael.rittri at carmenta.com
www.carmenta.com
More information about the Proj
mailing list