[GRASS-user] Problems with a GeoTIFF referenced in Gauss Boaga
	system
    Giacomo Piva 
    piva at meeo.it
       
    Thu Aug  5 10:14:48 EDT 2010
    
    
  
Hi all,
I'm fighting with r.in.gdal function and a GeoTiff gereferenced in the 
Gauss Boaga system.
I have this GeoTiff that is originlly georeferenced in the WGS84 system 
and I need to convert it into Gauss Boaga system using the "gdalwarp" 
GDAL command tool.
Moreover I need to convert this raster in a shapefile using the 
r.to.vect GRASS command.
This is the (correct) gdalinfo output of mu Gauss Boaga GeoTiff:
Coordinate System is:
PROJCS["Monte Mario / Italy zone 1",
     GEOGCS["Monte Mario",
         DATUM["Monte_Mario",
             SPHEROID["International 1924",6378388,297.000000000005,
                 AUTHORITY["EPSG","7022"]],
             AUTHORITY["EPSG","6265"]],
  PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433],
         AUTHORITY["EPSG","4265"]],
     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["metre",1,
         AUTHORITY["EPSG","9001"]],
     AUTHORITY["EPSG","3003"]]
When I run the command
r.in.gdal -o -e input=<inputfile.tif> output=inputraster location=location
I get this error about the region
ERROR: region for current mapset is not set
        run "g.region"
The g.region -p command return:
projection: 0 (x,y)
zone:       0
north:      2030
south:      0
west:       0
east:       1354
nsres:      1
ewres:      1
rows:       2030
cols:       1354
cells:      2748620
How should I do?
The same image georeferenced in the WGS84 system is properly loaded, so 
I'm wondering about to convert the GeoTiff AFTER loading int GRASS, but 
I don't know how.
Does someone can help me?
Thank you
-- 
Giacomo
    
    
More information about the grass-user
mailing list