<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=iso-8859-1"><meta name=Generator content="Microsoft Word 12 (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:0in;
        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;}
span.EmailStyle17
        {mso-style-type:personal-compose;
        font-family:"Calibri","sans-serif";
        color:windowtext;}
.MsoChpDefault
        {mso-style-type:export-only;}
@page WordSection1
        {size:8.5in 11.0in;
        margin:1.0in 1.0in 1.0in 1.0in;}
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=MsoNormal>Greetings to all,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I am looking for some assistance with a problem we have using the PROJ4 libraries, I am not a programmer neither am I geodesist – however I know enough to be dangerous! Please bear with me…<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>We have a customer in Shanghai China that is uses something he refers to as the Shanghai Local projection.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>It is built on Beijing 1954 and uses that datum, the 7 parameter shift he gives are as follows:<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>TX:          170.076<o:p></o:p></p><p class=MsoNormal>TY:          154.924<o:p></o:p></p><p class=MsoNormal>TZ:          97.308<o:p></o:p></p><p class=MsoNormal>RX:         1.666408<o:p></o:p></p><p class=MsoNormal>RY:          0.872624<o:p></o:p></p><p class=MsoNormal>RZ:          -8.648183<o:p></o:p></p><p class=MsoNormal>K:            0.00000186<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>The Ellipsoid is Krassovsky 1940<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>The local coordinate system parameters are TM:<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>CM: 121.46716061<o:p></o:p></p><p class=MsoNormal>Scale: 1.0<o:p></o:p></p><p class=MsoNormal>False Eastings: 0<o:p></o:p></p><p class=MsoNormal>False Northings: -34587140.589<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I constructed the following PROJ4 string:<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>+proj=tmerc +lat_0=0 +lon_0=121.46716061 +k=1 +x_0=0 +y_0=-3457140.589 +ellps=krass +towgs84=170.0760,154.9240,97.3080,1.6664,0.8726,-8.6482,0.00000186 +units=m +no_defs<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>However I am not getting the correct results, for example.<o:p></o:p></p><p class=MsoNormal>WGS84 coordinate<o:p></o:p></p><p class=MsoNormal>31° 17' 27.50965" N           121° 35' 54.41138" E         13.7729<o:p></o:p></p><p class=MsoNormal>Should give Shanghai Local coordinate<o:p></o:p></p><p class=MsoNormal>12501.0727E       6164.1228N<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>However the position from PROJ4 is:<o:p></o:p></p><p class=MsoNormal>12953.30E            6042.97N<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>We very rarely need to invoke this type of transformation so it’s quite possible something is awry with our code, but I sure would appreciate any thoughts on the problem.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Geoff Shipton,<br>Triton Imaging Inc.<br>Customer Support.<br>Capitola CA.<br>geoff@tritonimaginginc.com<br><br>The information contained in this message may be privileged or confidential information intended solely for the use of the addressee(s).  Any views or opinions presented in this communication do not necessarily represent those of Triton Imaging, Inc.  If the reader of this message is not the intended recipient, you are hereby notified that any dissemination, distribution or copying of this communication is strictly prohibited.  If you have received this communication in error, please immediately notify me and you are hereby instructed to delete all electronic copies and destroy all printed copies.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p></div></body></html>