<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Tahnk you Hamish,<div>running g.setproj in the PERMANENT mapset worked fine.</div><div><br></div><div>Thank you all for helping</div><div><br></div><div>Christian</div><div><br><div><div>Il giorno 27/feb/09, alle ore 03:29, Hamish ha scritto:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div><br>Christian Tiso wrote:<br><blockquote type="cite">Thank you markus, you're right, there's no datum<br></blockquote><blockquote type="cite">definition:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">GRASS 6.3.0 (fassaGB):~ > g.proj -w<br></blockquote><blockquote type="cite">PROJCS["Transverse Mercator",<br></blockquote><blockquote type="cite"> GEOGCS["international",<br></blockquote><blockquote type="cite"> DATUM["unknown",<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">SPHEROID["International_1924",6378388,297]],<br></blockquote><blockquote type="cite"> PRIMEM["Greenwich",0],<br></blockquote><blockquote type="cite"> UNIT["degree",0.0174532925199433]],<br></blockquote>....<br><blockquote type="cite">but this is beacause I created a new location from the GUI<br></blockquote><blockquote type="cite">start panel of GRASS by choosing the "EPSG code"<br></blockquote><blockquote type="cite">button and selectiong the 3003 code which should be the one<br></blockquote><blockquote type="cite">related to the GaussBoaga, Rome 1940 Monte Mario<br></blockquote><blockquote type="cite">projections.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Where am I wrong? It is not possible to create such a new<br></blockquote><blockquote type="cite">location only by choosing the EPSG code?<br></blockquote><br>The /usr/share/proj/epsg file's version of 3003 only defines the 1924<br>International ellipsoid, not the datum. so that's all grass knows about:<br><br># Monte Mario / Italy zone 1<br><3003> +proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0 +ellps=intl +units=m +no_defs <><br><br><br>try <a href="http://www.epsg-registry.org/">http://www.epsg-registry.org/</a> & go to "retrieve by code".<br><br>is "Monte Mario" aka "Rome 1940" datum? if so the epsg file may be faulty.<br>(don't know; ask on the proj4 mailing list)<br><br><br>grass's datum table has rome40:<br><br># Rome 1940 (Sardinia Parameters)<br>rome40 "Monte_Mario" international dx=-225 dy=-65 dz=9<br><br>so you could run g.setproj in the PERMANENT mapset, choosing only to<br>change the datum and set it that way, or hack the PROJ_INFO file manually.<br><br><br>Hamish<br><br><br><br><br><br></div></blockquote></div><br><div> <span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<br>Christian Tiso, collaborator<br><br>Department of Civil and Environmental Engineering<br>University of Trento<br>via Mesiano 77, I-38050 Trento, ITALY<br><a href="http://www.ing.unitn.it/dica">http://www.ing.unitn.it/dica</a><br>phone: +39 0461 882610<br>fax: +39 0461 882672<br>E-mail: <a href="mailto:christian.tiso@ing.unitn.it">christian.tiso@ing.unitn.it</a><br>~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~</div><div><br class="webkit-block-placeholder"></div></div></span><br class="Apple-interchange-newline"> </div><br></div></body></html>