<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
<title></title>
</head>
<body text="#000000" bgcolor="#ffffff">
Hi Jelmer,<br>
<br>
The correct parameters for the grid shift from epsg:28992 to WGS84
are described here:<br>
<br>
<a class="moz-txt-link-freetext" href="http://trac.osgeo.org/geotiff/ticket/22#comment:3">http://trac.osgeo.org/geotiff/ticket/22#comment:3</a><br>
<br>
There are several different formulas in the EPSG database, but the
most recent one is your number 1 (= 4). It has been put in the
PROJ4-database by Frank Warmerdam, as you can read in the posting
above. Numbers 2 and 3 are older definitions, although they deliver
about the same results. They should be discarded, but will probably
hang around yet for a long time in older software.<br>
<br>
Groeten,<br>
<br>
Jan<br>
<br>
On 30-7-2013 10:49, Jelmer Baas wrote:
<blockquote
cite="mid:BC9747AD8551184D9B674A8310FE337B0C99E10C@mailserver"
type="cite">
<meta http-equiv="Content-Type" content="text/html;
charset=ISO-8859-1">
<meta name="Generator" content="Microsoft Word 14 (filtered
medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.MsoPlainText, li.MsoPlainText, div.MsoPlainText
{mso-style-priority:99;
mso-style-link:"Plain Text Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
span.PlainTextChar
{mso-style-name:"Plain Text Char";
mso-style-priority:99;
mso-style-link:"Plain Text";
font-family:"Calibri","sans-serif";}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
<div class="WordSection1">
<p class="MsoPlainText"><span lang="NL">Hello,<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="NL"><o:p> </o:p></span></p>
<p class="MsoPlainText">For a new project we need to convert
coordinates from EPSG 28992 to EPSG 3044 with MapInfo, Oracle
AND a few Proj.4 programs (PostGIS, a custom DLL, Proj4js).
The problem we are having is that we're seeing some shifted
coordinates.<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">First test point (in the Netherlands)<o:p></o:p></p>
<p class="MsoPlainText"><span style="font-family: "Courier
New";">Origin 135745.924
517131.216<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family: "Courier
New";">Proj4js 236345.39474983042
5839505.550959114<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family: "Courier
New";">Postgis(1) 236345.390264227
5839505.55017783<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family: "Courier
New";">Postgis(2) 236345.394817119
5839505.55055913<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family: "Courier
New";" lang="NL">MapInfo
236345.47 5839504.058<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family: "Courier
New";" lang="NL">CS2CSN
236345.39 5839505.55<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family: "Courier
New";">DLL old(4) 236345.396465174
5839505.54496813<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family: "Courier
New";">DLL new(3) 236345.390259674
5839505.55009311<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family: "Courier
New";">Oracle 236345.49107424
5839504.12246656<o:p></o:p></span></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Test point (in Germany)<o:p></o:p></p>
<p class="MsoPlainText"><span style="font-family: "Courier
New";">Origin 254870.45
415123.56<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family: "Courier
New";">MapInfo 350265.960
5731670.259<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family: "Courier
New";">Postgis(1) 350263.4072329
5731667.880
<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family: "Courier
New";">Postgis(2) 350263.4132521
5731667.880
<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family: "Courier
New";">DLL new(3) 350263.407243136
5731667.88019826
<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family: "Courier
New";">DLL old(4) 350263.403514925
5731667.88565864
<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family: "Courier
New";">Oracle(5) 350265.967728329
5731670.33312477
<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family: "Courier
New";">CS2CS 350263.41
5731667.88 43.92 cs2cs +init=EPSG:28992 +to
+init=EPSG:3044<o:p></o:p></span></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">It seems MapInfo and Oracle more-or-less
agree, as do all Proj.4-based products. I don’t know which are
“right”, meaning I also don’t know which ones are wrong – let
alone find the fault.
<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Hopefully someone can help me with this
problem, or at least point me in the right direction.<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">(1): <span style="font-family:
"Courier New";">"+proj=sterea
+lat_0=52.15616055555555 +lon_0=5.38763888888889
+k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel
+towgs84=565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725
+units=m +no_defs "<o:p></o:p></span></p>
<p class="MsoPlainText">(2) Another server where I updated the
spatial_ref_system:
<span style="font-family: "Courier New";">"+proj=sterea
+lat_0=52.15616055555555 +lon_0=5.38763888888889
+k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m
+towgs84=565.2369,50.0087,465.658,-0.406857330322398,0.350732676542563,-1.8703473836068,4.0812
+no_defs no_defs"<o:p></o:p></span></p>
<p class="MsoPlainText">(3) A custom built DLL for usage in VB6.
"+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889
+k=0.999908 +x_0=155000 +y_0=463000 +ellps=bessel +units=m
+towgs84=565.2369,50.0087,465.658,-0.406857330322398,0.350732676542563,-1.8703473836068,4.0812
+no_defs <>"<o:p></o:p></p>
<p class="MsoPlainText">(4) Same DLL. “+proj=sterea
+lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079
+x_0=155000 +y_0=463000 +ellps=bessel
+towgs84=565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725
+units=m +no_defs <>”<o:p></o:p></p>
<p class="MsoPlainText">(5)<span style="font-family:
"Courier New";"> select
sdo_cs.transform(SDO_GEOMETRY('POINT(254870.45 415123.56)',
90112),25832) from dual;<o:p></o:p></span></p>
<p class="MsoPlainText">(5b) The Oracle definition contains some
tags I don’t see anywhere in Proj.4 (like the Bessel 1841
part): <o:p></o:p></p>
<p class="MsoPlainText">PROJCS["Netherlands National System",
GEOGCS [ "Netherlands Bessel", DATUM ["Netherlands Bessel",
SPHEROID ["Bessel 1841", 6377397.155, 299.1528128]], PRIMEM [
"Greenwich", 0.000000 ], UNIT ["Decimal Degree",
0.01745329251994330]], PROJECTION ["Stereographic"], PARAMETER
["Scale_Factor", 0.999908], PARAMETER ["Central_Meridian",
5.387639], PARAMETER ["Latitude_Of_Origin", 52.156161],
PARAMETER ["False_Easting", 155000.000000], PARAMETER
["False_Northing", 463000.000000], UNIT ["Meter",
1.000000000000]]
<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">A side note is that transformation via
my DLL doesn’t work properly (i.e., coordinates are shifted
tens of meters), unless I *<b>first</b>* transform to WGS84
and then to 3044. I assume this is caused by an easy-to-fix
error, but I don’t know what that is.<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">P.S.: PostGIS, Proj4JS and my epsg file
are all up-to-date.<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Thanks for any help,<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText"><span lang="NL">With kind regards,<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="NL">Jelmer Baas<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="NL"><o:p> </o:p></span></p>
</div>
<pre wrap="">
<fieldset class="mimeAttachmentHeader"></fieldset>
_______________________________________________
Proj mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Proj@lists.maptools.org">Proj@lists.maptools.org</a>
<a class="moz-txt-link-freetext" href="http://lists.maptools.org/mailman/listinfo/proj">http://lists.maptools.org/mailman/listinfo/proj</a></pre>
</blockquote>
</body>
</html>