<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:st1="urn:schemas-microsoft-com:office:smarttags" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv=Content-Type content="text/html; charset=iso-8859-1">
<meta name=Generator content="Microsoft Word 11 (filtered medium)">
<!--[if !mso]>
<style>
v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style>
<![endif]--><o:SmartTagType
namespaceuri="urn:schemas-microsoft-com:office:smarttags" name="City"/>
<o:SmartTagType namespaceuri="urn:schemas-microsoft-com:office:smarttags"
name="place"/>
<!--[if !mso]>
<style>
st1\:*{behavior:url(#default#ieooui) }
</style>
<![endif]-->
<style>
<!--
/* Font Definitions */
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman";}
a:link, span.MsoHyperlink
{color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{color:blue;
text-decoration:underline;}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:Arial;
color:navy;}
@page Section1
{size:595.3pt 841.9pt;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.Section1
{page:Section1;}
-->
</style>
</head>
<body lang=SV link=blue vlink=blue>
<div class=Section1>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'>Amro,<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'>Jose' is right, especially if you want agreement with
other data that have been transformed<br>
with those datum shift parameters. <o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'>But you wrote:<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=3 face="Times New Roman"><span lang=EN-GB
style='font-size:12.0pt'>> My need here is to get as accurate as possible
results (cm-level if possible).<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face="Times New Roman"><span lang=EN-GB
style='font-size:12.0pt'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'>Then, a better alternative could be to use the grid shift
file BETA2007.gsb. It is not distributed with Proj.4,<br>
but you can find a link at <a href="http://trac.osgeo.org/proj/">http://trac.osgeo.org/proj/</a>
under Download, Datum shift grids.<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'>If you download it to the directory where the other grid
shift files are in Proj.4 (the </span></font><font face="Courier New"><span
lang=EN-GB style='font-family:"Courier New"'>nad</span></font><font face=Arial><span
lang=EN-GB style='font-family:Arial'> directory),<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'>then you can do <o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=3 face="Courier New"><span lang=EN-GB
style='font-size:12.0pt;font-family:"Courier New"'>>cs2cs +init=epsg:4326
+to +ellps=bessel +proj=tmerc +lon_0=12 +x_0=4500000 +nadgrids=BETA2007.gsb<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face="Courier New"><span lang=EN-GB
style='font-size:12.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=3 face="Courier New"><span lang=EN-GB
style='font-size:12.0pt;font-family:"Courier New"'>13.755944 51.054555<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face="Courier New"><span lang=EN-GB
style='font-size:12.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=3 face="Courier New"><span lang=EN-GB
style='font-size:12.0pt;font-family:"Courier New"'>4623220.44 5659181.52
0.00<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face="Courier New"><span lang=EN-GB
style='font-size:12.0pt;font-family:"Courier New"'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'>This result differs by about 1.7 meters from Jose’s
example. <o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'>The BETA2007.gsb file is claimed to give submeter
accuracy, <o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'>“eine Koordinatengenauigkeit im Submeterbereich”.
(But maybe it <o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'>is still possible that the 7-parameter transformation
you cite is more <o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'>accurate locally around <st1:City w:st="on"><st1:place
w:st="on">Dresden</st1:place></st1:City>.)<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'>Best regards,<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span lang=EN-GB style='font-size:
12.0pt;font-family:Arial'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span style='font-size:12.0pt;
font-family:Arial'>Mikael Rittri<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span style='font-size:12.0pt;
font-family:Arial'>Carmenta<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span style='font-size:12.0pt;
font-family:Arial'>Sweden<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span style='font-size:12.0pt;
font-family:Arial'>http://www.carmenta.com<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face=Arial><span style='font-size:12.0pt;
font-family:Arial'><o:p> </o:p></span></font></p>
<div>
<div class=MsoNormal align=center style='text-align:center'><font size=3
face="Times New Roman"><span lang=EN-US style='font-size:12.0pt'>
<hr size=2 width="100%" align=center tabindex=-1>
</span></font></div>
<p class=MsoNormal><b><font size=2 face=Tahoma><span lang=EN-US
style='font-size:10.0pt;font-family:Tahoma;font-weight:bold'>From:</span></font></b><font
size=2 face=Tahoma><span lang=EN-US style='font-size:10.0pt;font-family:Tahoma'>
proj-bounces@lists.maptools.org [mailto:proj-bounces@lists.maptools.org] <b><span
style='font-weight:bold'>On Behalf Of </span></b>Jose Gonçalves<br>
<b><span style='font-weight:bold'>Sent:</span></b> den 7 maj 2012 15:59<br>
<b><span style='font-weight:bold'>To:</span></b> PROJ.4 and general Projections
Discussions<br>
<b><span style='font-weight:bold'>Subject:</span></b> Re: [Proj] Using EPSG
codes with cs2cs</span></font><span lang=EN-US><o:p></o:p></span></p>
</div>
<p class=MsoNormal><font size=3 face="Times New Roman"><span style='font-size:
12.0pt'><o:p> </o:p></span></font></p>
<p class=MsoNormal style='margin-bottom:12.0pt'><font size=3
face="Times New Roman"><span style='font-size:12.0pt'>In addition to my
previous e-mail, since you have your own datum shift parameters it is better to
write the full proj string, instead of using an epsg predefinition. Your
rotation angles are in radians but cs2cs expects arc-seconds. They become 1.04,
0.35, -3.08.<br>
According to the map projection definition you give, the command line should be
something like:<o:p></o:p></span></font></p>
<div>
<p class=MsoNormal><font size=3 face="Times New Roman"><span style='font-size:
12.0pt'>cs2cs +init=epsg:4326 +to +ellps=bessel +proj=tmerc +lon_0=12
+x_0=4500000 +towgs84=582,105,414,1.040,0.350,-3.080,8.3<br>
13.755944 51.054555 <br>
4623218.81 5659181.06 -46.16<br>
<br>
Regards<br>
<br>
Jose' Gonçalves<br>
<br>
<br>
<br>
<br>
2012/5/7 Jose Gonçalves <<a href="mailto:jagoncal@gmail.com" target="_blank">jagoncal@gmail.com</a>><o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face="Times New Roman"><span style='font-size:
12.0pt'>Hi<br>
<br>
Are you sure of EPSG code 31495? In my epsg file, installed with proj, I can
find only:<br>
<br>
epsg: 31461 # DHDN / 3-degree Gauss zone 1 (deprecated)<br>
epsg: 31462 # DHDN / 3-degree Gauss zone 2 (deprecated)<br>
epsg: 31463 # DHDN / 3-degree Gauss zone 3 (deprecated)<br>
epsg: 31464 # DHDN / 3-degree Gauss zone 4 (deprecated)<br>
epsg: 31465 # DHDN / 3-degree Gauss zone 5 (deprecated)<br>
epsg: 5520 # DHDN / 3-degree Gauss-Kruger zone 1<br>
epsg: 31466 # DHDN / 3-degree Gauss-Kruger zone 2<br>
epsg: 31467 # DHDN / 3-degree Gauss-Kruger zone 3<br>
epsg: 31468 # DHDN / 3-degree Gauss-Kruger zone 4<br>
epsg: 31469 # DHDN / 3-degree Gauss-Kruger zone 5<br>
<br>
If a projection ID is not in the file the expected answer is the one you got.
It is not a configuration problem.<br>
Check <a href="http://spatialreference.org" target="_blank">http://spatialreference.org</a>
to check the epsg code<br>
<br>
Regards<font color="#888888"><span style='color:#888888'><br>
<br>
<span class=hoenzb>Jose Gonçalves</span></span></font><o:p></o:p></span></font></p>
<div>
<div>
<p class=MsoNormal style='margin-bottom:12.0pt'><font size=3
face="Times New Roman"><span style='font-size:12.0pt'><br>
<br>
<br>
<br>
<br>
<o:p></o:p></span></font></p>
<div>
<p class=MsoNormal><font size=3 face="Times New Roman"><span style='font-size:
12.0pt'>2012/5/7 rps_rt <<a href="mailto:ec_eng01@yahoo.com" target="_blank">ec_eng01@yahoo.com</a>><o:p></o:p></span></font></p>
<p class=MsoNormal><font size=3 face="Times New Roman"><span style='font-size:
12.0pt'>Dear Warmerdam and all members,<br>
<br>
I have been browsing and reading for a while to convert my coordinates with<br>
no hope until finding this Proj. 4.<br>
<br>
Kindly, I seek your support here ... I am using Linux Ubuntu to convert from<br>
long/lat (WGS84) to Cartesian (x/y) and vice versa. My need here is to get<br>
as accurate as possible results (cm-level if possible).<br>
<br>
I was reading this thread because I had a problem of getting a cs2cs results<br>
for espg=31495, while I could get them for espg=31468. But the conversion<br>
results with the later one were not accurate ~ almost 100 meters away.<br>
<br>
e.g.,: echo "13.755944 51.054555" | cs2cs +init=epsg:4326 +to<br>
+init=epsg:31468<br>
<br>
but when I change it to: echo "13.755944 51.054555" | cs2cs
+init=epsg:4326<br>
+to +init=epsg:31495,<br>
I get the message "projection initialization failure cause: no options
found<br>
in 'init' file".<br>
<br>
I tried to change something based on this thread, to make sure that my<br>
PROJ_LIB path is in place, but it seems I ruined everything. Afterwards, all<br>
espg zones gave the below error:<br>
"<br>
pj_open_lib(epsg): call fopen(/usr/local/bin/cs2cs/epsg) - failed<br>
Using from definition: init=epsg:31467<br>
Rel. 4.8.0, 6 March 2012<br>
<cs2cs>:<br>
projection initialization cause: no system list, errno: 20<br>
"<br>
<br>
I tried to "make uninstall" and "make install" but nothing
was fixed!!<br>
<br>
Here are more details about my project:<br>
<br>
My locations are in Dresden city, Germany, (DHDN4), and I have (I assume)<br>
all what is needed to get my coordinates converted accurately, but no idea<br>
how to process them in Proj. 4 (I am totally new here):<br>
<br>
PROJCS["DHDN / Germany zone<br>
4",GEOGCS["DHDN",DATUM["Deutsche_Hauptdreiecksnetz",<br>
SPHEROID["Bessel 1841",6377397.155,299.1528128]],PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],<br>
PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",12],<br>
PARAMETER["scale_factor",1],PARAMETER["false_easting",4500000],<br>
PARAMETER["false_northing",0],UNIT["metre",1]]<br>
TOWGS84=582.0,105.0,414.0,5.042e-6,1.696e-6,-1.4932e-5,8.3<br>
EPSG=31495<br>
<br>
%Seven parameters:<br>
% 582.00 - Translation dx in meter to WGS84<br>
% 105.0 - Translation dy in meter to WGS84<br>
% 414.50 - Translation dz in meter to WGS84<br>
% 5.042e-6 - Rotation ro (omega) in Radiant to WGS84<br>
% 1.696e-6 - Rotation rf (phi) in radiant to WGS84<br>
% -1.4932e-5 - Rotation rk (Kappa) in Radiant to WGS84<br>
% 8.3e-06 - Mapscale factor in ppm (parts per million) to WGS84(m)<br>
<br>
Sorry for the long thread but did not want to miss any details, even in my<br>
first ever thread in this forum.<br>
<br>
I will be looking forward to reading form you.<br>
Best regards,<br>
Amro<br>
<br>
--<br>
View this message in context: <a
href="http://osgeo-org.1560.n6.nabble.com/Using-EPSG-codes-with-cs2cs-tp3841592p4957639.html"
target="_blank">http://osgeo-org.1560.n6.nabble.com/Using-EPSG-codes-with-cs2cs-tp3841592p4957639.html</a><o:p></o:p></span></font></p>
<div>
<p class=MsoNormal><font size=3 face="Times New Roman"><span style='font-size:
12.0pt'>Sent from the PROJ.4 mailing list archive at Nabble.com.<o:p></o:p></span></font></p>
</div>
<div>
<div>
<p class=MsoNormal><font size=3 face="Times New Roman"><span style='font-size:
12.0pt'>_______________________________________________<br>
Proj mailing list<br>
<a href="mailto:Proj@lists.maptools.org" target="_blank">Proj@lists.maptools.org</a><br>
<a href="http://lists.maptools.org/mailman/listinfo/proj" target="_blank">http://lists.maptools.org/mailman/listinfo/proj</a><o:p></o:p></span></font></p>
</div>
</div>
</div>
<p class=MsoNormal><font size=3 face="Times New Roman"><span style='font-size:
12.0pt'><o:p> </o:p></span></font></p>
</div>
</div>
</div>
<p class=MsoNormal><font size=3 face="Times New Roman"><span style='font-size:
12.0pt'><o:p> </o:p></span></font></p>
</div>
</body>
</html>