Thank you markus, you're right, there's no datum definition:

GRASS 6.3.0 (fassaGB):~ > g.proj -w
PROJCS["Transverse Mercator",
    GEOGCS["international",
        DATUM["unknown",
            SPHEROID["International_1924",6378388,297]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",9],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",1500000],
    PARAMETER["false_northing",0],
    UNIT["Meter",1]]

but this is beacause I created a new location from the GUI start panel of GRASS by choosing the "EPSG code" button and selectiong the 3003 code which should be the one related to the GaussBoaga, Rome 1940 Monte Mario projections.

Where am I wrong? It is not possible to create such a new location only by choosing the EPSG code?

Thank you for helping me.

Christian Tiso



Il giorno 26/feb/09, alle ore 00:19, Markus Neteler ha scritto:

On Wed, Feb 25, 2009 at 5:02 PM, Christian Tiso
<[email protected]> wrote:
Hi,
using command r.proj to reproject a raster topographic map from a Gauss Boaga location (EPSG:3003) to a UTM location (EPSG:32632) I obtain a shift error of the reprojected map of about 40-50m in south direction. It is
related to the towgs84 bursa wolf parameters?
Is this a known bug?
I'm using Grass 6.3.0-2 on a Mac Intel with Moretti's Build binaries.

please post

g.proj -w

I assume that the datum definition is lacking or wrong (towgs84).

Markus

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Christian Tiso, collaborator

Department of Civil and Environmental Engineering
University of  Trento
via Mesiano 77, I-38050 Trento, ITALY
http://www.ing.unitn.it/dica
phone: +39 0461 882610
fax:   +39 0461 882672
E-mail: [email protected]
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~



_______________________________________________
grass-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-dev

Reply via email to