<HTML dir=ltr><HEAD><TITLE>Re: [Proj] cs2cs usage</TITLE>
<META http-equiv=Content-Type content="text/html; charset=unicode">
<META content="MSHTML 6.00.6000.16608" name=GENERATOR><BASE href=http://mailstore.rossoalice.alice.it/exchange/Alice000000001264700/Bozze/R:%20[Proj]%20cs2cs%20usage.EML/1_text.htm></HEAD>
<BODY>
<DIV id=idOWAReplyText32082 dir=ltr>
<DIV dir=ltr><FONT face=Arial color=#000000 size=2>Hi,</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><FONT size=2>thanks you all for the precious help.</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><FONT size=2>I entered both *versions* of the command:</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><FONT face="Courier New" size=2>cs2cs +init=epsg:4326 +to +proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0 +ellps=intl +uni<BR>ts=m +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68 +no_defs<BR>10.73602778 45.75333333<BR>1635050.62 5068127.39 -45.30<BR></FONT></DIV>
<DIV dir=ltr><FONT face="Courier New" size=2>cs2cs +init=epsg:4326 +to +init=epsg:3003 +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68<BR>10.73602778 45.75333333<BR>1635050.62 5068127.39 -45.30</FONT></DIV>
<DIV dir=ltr><FONT face="Courier New" size=2></FONT> </DIV>
<DIV dir=ltr><FONT size=2>as you can see, the result it's the same, ...as it should be ;-)</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><FONT size=2>Markus, because I'll work on Lombardia DEM rasters only using WGS84 *input points* (for r.horizon and r.sun2 computations), do you think that I should transform *all rasters* from the current projection to WGS84 projection?</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><FONT size=2>About the *check-out procedure*, I fear I must trust the result as it is, I cannot see another solution :-( ....</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><FONT size=2>Thanks</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><FONT size=2>Marco</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><FONT size=2>PS: Frank, thank you very much for CRTL+Z; it was really frustrating to kill the process every time!... ;-)</FONT></DIV>
<DIV dir=ltr><FONT face="Courier New" size=2></FONT> </DIV></DIV>
<DIV dir=ltr><BR>
<HR tabIndex=-1>
<FONT face=Tahoma size=2><B>Da:</B> proj-bounces@lists.maptools.org per conto di Frank Warmerdam<BR><B>Inviato:</B> gio 20/03/2008 15.49<BR><B>A:</B> PROJ.4 and general Projections Discussions<BR><B>Oggetto:</B> Re: [Proj] cs2cs usage<BR></FONT><BR></DIV>
<DIV>
<P><FONT size=2>marco.pasetti@alice.it wrote:<BR>> Hi all,<BR>> <BR>> I read cs2cs manual page, but I still have big doubts about it! quickly:<BR>> <BR>> 0. Preamble: I built PROJ.4 4.6 on windows using MinGW<BR>> <BR>> 1. Question 1: when I launch cs2cs from command line it executes the requested operation and then waits for others inputs... how can I quit, without the need of killing the process in task manager?<BR><BR>Marco,<BR><BR>You normally need to indicate "end of input" to terminate cs2cs. On unix<BR>this is done with control-D, on windows with control-Z.<BR><BR>> 2. Question 2: I need to trasform single points from WGS84 (EPSG ID = 4326) to Gauss-Boaga Roma40 West (EPSG ID = 3003)<BR>> while ID=4326 determinates an univocal datum, to correctly define Gauss-Boaga Roma40 West I need to specify also +towgs84 parameters, retrieved by GRASS (creating a new location in GRASS using 3003 EPSG code, it asks for furthers options, suggesting the correct parameters to set):<BR>> the complete datum definition for Gauss-Boaga Roma40 West is as follows (from GRASS):<BR>> <BR>> EPSG ID = 3003<BR>> Option 2<BR>> Used in Italy (Peninsular Part)<BR>> towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68<BR>> Accuracy 3-4m<BR>> <BR>> how can I use cs2cs respecting the previous parameters?<BR>> I tried with cs2cs +init=epsg:4326 +to +init=epsg:3003<BR>> but I don't know how to set towgs84 parameters (and even if they are needed)<BR><BR>The epsg:3003 definition looks like:<BR><BR>+proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0 +ellps=intl<BR>+units=m +no_defs <><BR><BR>So you should use it directly and then add your towgs84 parameter.<BR><BR>eg.<BR><BR>cs2cs +init=epsg:4326 +to \<BR>+proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0 +ellps=intl \<BR>+units=m +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68 +no_defs<BR><BR>Best regards,<BR>--<BR>---------------------------------------+--------------------------------------<BR>I set the clouds in motion - turn up | Frank Warmerdam, warmerdam@pobox.com<BR>light and sound - activate the windows | <A href="http://pobox.com/~warmerdam">http://pobox.com/~warmerdam</A><BR>and watch the world go round - Rush | President OSGeo, <A href="http://osgeo.org/">http://osgeo.org</A><BR><BR>_______________________________________________<BR>Proj mailing list<BR>Proj@lists.maptools.org<BR><A href="http://lists.maptools.org/mailman/listinfo/proj">http://lists.maptools.org/mailman/listinfo/proj</A><BR></FONT></P></DIV></BODY></HTML>