<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<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]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoPlainText"><span lang="NL">Hello,<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="NL"><o:p>&nbsp;</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>&nbsp;</o:p></p>
<p class="MsoPlainText">First test point (in the Netherlands)<o:p></o:p></p>
<p class="MsoPlainText"><span style="font-family:&quot;Courier New&quot;">Origin&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 135745.924 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 517131.216<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family:&quot;Courier New&quot;">Proj4js&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 236345.39474983042&nbsp;&nbsp; 5839505.550959114<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family:&quot;Courier New&quot;">Postgis(1)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 236345.390264227&nbsp;&nbsp;&nbsp;&nbsp; 5839505.55017783<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family:&quot;Courier New&quot;">Postgis(2)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 236345.394817119&nbsp;&nbsp;&nbsp;&nbsp; 5839505.55055913<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="NL" style="font-family:&quot;Courier New&quot;">MapInfo&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 236345.47&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5839504.058<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="NL" style="font-family:&quot;Courier New&quot;">CS2CSN&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 236345.39&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5839505.55<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family:&quot;Courier New&quot;">DLL old(4) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 236345.396465174&nbsp;&nbsp;&nbsp;&nbsp; 5839505.54496813<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family:&quot;Courier New&quot;">DLL new(3)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 236345.390259674&nbsp;&nbsp;&nbsp;&nbsp; 5839505.55009311<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family:&quot;Courier New&quot;">Oracle&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 236345.49107424&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5839504.12246656<o:p></o:p></span></p>
<p class="MsoPlainText"><o:p>&nbsp;</o:p></p>
<p class="MsoPlainText">Test point (in Germany)<o:p></o:p></p>
<p class="MsoPlainText"><span style="font-family:&quot;Courier New&quot;">Origin&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 254870.45&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 415123.56<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family:&quot;Courier New&quot;">MapInfo&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 350265.960&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5731670.259<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family:&quot;Courier New&quot;">Postgis(1)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 350263.4072329&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5731667.880&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family:&quot;Courier New&quot;">Postgis(2)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 350263.4132521&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5731667.880&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family:&quot;Courier New&quot;">DLL new(3)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 350263.407243136&nbsp;&nbsp;&nbsp;&nbsp; 5731667.88019826
<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family:&quot;Courier New&quot;">DLL old(4)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 350263.403514925&nbsp;&nbsp;&nbsp;&nbsp; 5731667.88565864
<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family:&quot;Courier New&quot;">Oracle(5)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 350265.967728329&nbsp;&nbsp;&nbsp;&nbsp; 5731670.33312477&nbsp;&nbsp;&nbsp;&nbsp;
<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="font-family:&quot;Courier New&quot;">CS2CS&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 350263.41&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 5731667.88 43.92&nbsp;&nbsp;&nbsp;&nbsp; cs2cs &#43;init=EPSG:28992 &#43;to &#43;init=EPSG:3044<o:p></o:p></span></p>
<p class="MsoPlainText"><o:p>&nbsp;</o:p></p>
<p class="MsoPlainText">It seems MapInfo and Oracle more-or-less agree, as do all Proj.4-based products. I don&#8217;t know which are &#8220;right&#8221;, meaning I also don&#8217;t know which ones are wrong &#8211; let alone find the fault.
<o:p></o:p></p>
<p class="MsoPlainText"><o:p>&nbsp;</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>&nbsp;</o:p></p>
<p class="MsoPlainText">(1): <span style="font-family:&quot;Courier New&quot;">&quot;&#43;proj=sterea &#43;lat_0=52.15616055555555 &#43;lon_0=5.38763888888889 &#43;k=0.9999079 &#43;x_0=155000 &#43;y_0=463000 &#43;ellps=bessel &#43;towgs84=565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725 &#43;units=m
 &#43;no_defs &quot;<o:p></o:p></span></p>
<p class="MsoPlainText">(2) Another server where I updated the spatial_ref_system:
<span style="font-family:&quot;Courier New&quot;">&quot;&#43;proj=sterea &#43;lat_0=52.15616055555555 &#43;lon_0=5.38763888888889 &#43;k=0.9999079 &#43;x_0=155000 &#43;y_0=463000 &#43;ellps=bessel &#43;units=m &#43;towgs84=565.2369,50.0087,465.658,-0.406857330322398,0.350732676542563,-1.8703473836068,4.0812
 &#43;no_defs no_defs&quot;<o:p></o:p></span></p>
<p class="MsoPlainText">(3) A custom built DLL for usage in VB6. &quot;&#43;proj=sterea &#43;lat_0=52.15616055555555 &#43;lon_0=5.38763888888889 &#43;k=0.999908 &#43;x_0=155000 &#43;y_0=463000 &#43;ellps=bessel &#43;units=m &#43;towgs84=565.2369,50.0087,465.658,-0.406857330322398,0.350732676542563,-1.8703473836068,4.0812
 &#43;no_defs &lt;&gt;&quot;<o:p></o:p></p>
<p class="MsoPlainText">(4) Same DLL. &#8220;&#43;proj=sterea &#43;lat_0=52.15616055555555 &#43;lon_0=5.38763888888889 &#43;k=0.9999079 &#43;x_0=155000 &#43;y_0=463000 &#43;ellps=bessel &#43;towgs84=565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725 &#43;units=m &#43;no_defs&nbsp; &lt;&gt;&#8221;<o:p></o:p></p>
<p class="MsoPlainText">(5)<span style="font-family:&quot;Courier New&quot;"> select sdo_cs.transform(SDO_GEOMETRY('POINT(254870.45 415123.56)', 90112),25832)&nbsp; from dual;<o:p></o:p></span></p>
<p class="MsoPlainText">(5b) The Oracle definition contains some &nbsp;tags I don&#8217;t see anywhere in Proj.4 (like the Bessel 1841 part): &nbsp;<o:p></o:p></p>
<p class="MsoPlainText">PROJCS[&quot;Netherlands National System&quot;, GEOGCS [ &quot;Netherlands Bessel&quot;, DATUM [&quot;Netherlands Bessel&quot;, SPHEROID [&quot;Bessel 1841&quot;, 6377397.155, 299.1528128]], PRIMEM [ &quot;Greenwich&quot;, 0.000000 ], UNIT [&quot;Decimal Degree&quot;, 0.01745329251994330]], PROJECTION
 [&quot;Stereographic&quot;], PARAMETER [&quot;Scale_Factor&quot;, 0.999908], PARAMETER [&quot;Central_Meridian&quot;, 5.387639], PARAMETER [&quot;Latitude_Of_Origin&quot;, 52.156161], PARAMETER [&quot;False_Easting&quot;, 155000.000000], PARAMETER [&quot;False_Northing&quot;, 463000.000000], UNIT [&quot;Meter&quot;, 1.000000000000]]
<o:p></o:p></p>
<p class="MsoPlainText"><o:p>&nbsp;</o:p></p>
<p class="MsoPlainText">A side note is that transformation via my DLL doesn&#8217;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&#8217;t
 know what that is.<o:p></o:p></p>
<p class="MsoPlainText"><o:p>&nbsp;</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>&nbsp;</o:p></p>
<p class="MsoPlainText">Thanks for any help,<o:p></o:p></p>
<p class="MsoPlainText"><o:p>&nbsp;</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>&nbsp;</o:p></span></p>
</div>
</body>
</html>