<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> </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 lang="NL" style="font-family:"Courier New"">MapInfo 236345.47 5839504.058<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="NL" style="font-family:"Courier New"">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>
</body>
</html>