[GRASS-user] Problems with a GeoTIFF referenced in Gauss Boaga system

Giacomo_Piva · August 5, 2010, 2:14pm

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

Paul_Kelly · August 5, 2010, 2:39pm

On Thu, 5 Aug 2010, Giacomo Piva wrote:

[...]

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"

Just a suggestion to try - does it work without the -e flag? The "location" option automatically sets the extents of the new location to the same as your data, so as far as I can see you should not need the -e flag. Although that may not be related to the error; I'm not sure...

However if that turns out to be the problem, a check should be added to warn that -e and location= are incompatible. Also I think the message to just "run g.region" is confusing; I guess it refers to the old interactive g.region?

Paul

Daniel_Victoria · August 5, 2010, 2:44pm

You need to create a Location in Grass were the coordinate system is
the Gauss Boaga. Your current location has no reference system, it's a
XY location, as g.region -p output shows.

When you start Grass ask for it to create a new location based on a
georeferenced file. Then point to your GeoTiff and the location should
be set correctly. After that, r.in.gdal should work.

Daniel

On Thu, Aug 5, 2010 at 11:14 AM, Giacomo Piva <piva@meeo.it> wrote:

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

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Giacomo_Piva · August 5, 2010, 3:07pm

No... it don't work.

On 08/05/2010 04:39 PM, Paul Kelly wrote:

On Thu, 5 Aug 2010, Giacomo Piva wrote:

[...]

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"

Just a suggestion to try - does it work without the -e flag? The "location" option automatically sets the extents of the new location to the same as your data, so as far as I can see you should not need the -e flag. Although that may not be related to the error; I'm not sure...

However if that turns out to be the problem, a check should be added to warn that -e and location= are incompatible. Also I think the message to just "run g.region" is confusing; I guess it refers to the old interactive g.region?

Paul

--
Giacomo Piva

MEEO S.r.l.
Via Saragat 9
I-44122 Ferrara, Italy
Tel: +39 0532 1861501
Fax: +39 0532 1861637
http://www.meeo.it

========================================
"ATTENZIONE:le informazioni contenute in questo messaggio sono
da considerarsi confidenziali ed il loro utilizzo è riservato unicamente
al destinatario sopra indicato. Chi dovesse ricevere questo messaggio
per errore è tenuto ad informare il mittente ed a rimuoverlo
definitivamente da ogni supporto elettronico o cartaceo."

"WARNING:This message contains confidential and/or proprietary
information which may be subject to privilege or immunity and which
is intended for use of its addressee only. Should you receive this
message in error, you are kindly requested to inform the sender and
to definitively remove it from any paper or electronic format."

Giacomo_Piva · August 5, 2010, 3:19pm

On 08/05/2010 04:44 PM, Daniel Victoria wrote:

You need to create a Location in Grass were the coordinate system is
the Gauss Boaga. Your current location has no reference system, it's a
XY location, as g.region -p output shows.

When you start Grass ask for it to create a new location based on a
georeferenced file. Then point to your GeoTiff and the location should
be set correctly. After that, r.in.gdal should work.

This is the point.
I'm using GRASS without UI, only running command by a bash script.
Actually I'm using a location (I don't remember how I made it) by setting up the grassrc file with the information about it:

[...]
echo "LOCATION_NAME: start" > .grassrc6
echo "MAPSET: PERMANENT" >> .grassrc6
echo "DIGITIZER: none" >> .grassrc6
echo "GISDBASE: $GISDBASE" >> .grassrc6

export GISRC=.grassrc6
[...]

then I read the input data ant then I change the location "jumping" into the region just created with the r.in.gdal command.
How can I set the region by scratch or convert the location to a georeferenced location?

I also tied to use the g.proj command in order to get the georefernce of the input image and try to set the region but the result is the same:

g.proj georef=<input file> location=<location>

Trying to open with OGR...
Trying to open with GDAL...
...succeeded.

g.proj -d

XY location (unprojected)

g.region -p

projection: 0 (x,y)
zone: 0
north: 2030
south: 0
west: 0
east: 1354
nsres: 1
ewres: 1
rows: 2030
cols: 1354
cells: 2748620

Daniel

On Thu, Aug 5, 2010 at 11:14 AM, Giacomo Piva<piva@meeo.it> wrote:


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

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

--
Giacomo Piva

MEEO S.r.l.
Via Saragat 9
I-44122 Ferrara, Italy
Tel: +39 0532 1861501
Fax: +39 0532 1861637
http://www.meeo.it

========================================
"ATTENZIONE:le informazioni contenute in questo messaggio sono
da considerarsi confidenziali ed il loro utilizzo è riservato unicamente
al destinatario sopra indicato. Chi dovesse ricevere questo messaggio
per errore è tenuto ad informare il mittente ed a rimuoverlo
definitivamente da ogni supporto elettronico o cartaceo."

"WARNING:This message contains confidential and/or proprietary
information which may be subject to privilege or immunity and which
is intended for use of its addressee only. Should you receive this
message in error, you are kindly requested to inform the sender and
to definitively remove it from any paper or electronic format."

Daniel_Victoria · August 5, 2010, 5:14pm

Well, you are running Grass in a different way then... not sure if I
can help. Maybe try the location parameter in r.in.gdal. It's suposed
to create an appropriate location before importing the data...

Daniel

On Thu, Aug 5, 2010 at 12:19 PM, Giacomo Piva <piva@meeo.it> wrote:

On 08/05/2010 04:44 PM, Daniel Victoria wrote:

You need to create a Location in Grass were the coordinate system is
the Gauss Boaga. Your current location has no reference system, it's a
XY location, as g.region -p output shows.

When you start Grass ask for it to create a new location based on a
georeferenced file. Then point to your GeoTiff and the location should
be set correctly. After that, r.in.gdal should work.

This is the point.
I'm using GRASS without UI, only running command by a bash script.
Actually I'm using a location (I don't remember how I made it) by setting up
the grassrc file with the information about it:

[...]
echo "LOCATION_NAME: start" > .grassrc6
echo "MAPSET: PERMANENT" >> .grassrc6
echo "DIGITIZER: none" >> .grassrc6
echo "GISDBASE: $GISDBASE" >> .grassrc6

export GISRC=.grassrc6
[...]

then I read the input data ant then I change the location "jumping" into the
region just created with the r.in.gdal command.
How can I set the region by scratch or convert the location to a
georeferenced location?

I also tied to use the g.proj command in order to get the georefernce of the
input image and try to set the region but the result is the same:

g.proj georef=<input file> location=<location>

Trying to open with OGR...
Trying to open with GDAL...
...succeeded.

g.proj -d

XY location (unprojected)

g.region -p

projection: 0 (x,y)
zone: 0
north: 2030
south: 0
west: 0
east: 1354
nsres: 1
ewres: 1
rows: 2030
cols: 1354
cells: 2748620

Daniel

On Thu, Aug 5, 2010 at 11:14 AM, Giacomo Piva<piva@meeo.it> wrote:

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

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

--
Giacomo Piva

MEEO S.r.l.
Via Saragat 9
I-44122 Ferrara, Italy
Tel: +39 0532 1861501
Fax: +39 0532 1861637
http://www.meeo.it

========================================
"ATTENZIONE:le informazioni contenute in questo messaggio sono
da considerarsi confidenziali ed il loro utilizzo è riservato unicamente
al destinatario sopra indicato. Chi dovesse ricevere questo messaggio
per errore è tenuto ad informare il mittente ed a rimuoverlo
definitivamente da ogni supporto elettronico o cartaceo."

"WARNING:This message contains confidential and/or proprietary
information which may be subject to privilege or immunity and which
is intended for use of its addressee only. Should you receive this
message in error, you are kindly requested to inform the sender and
to definitively remove it from any paper or electronic format."

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Giacomo_Piva · August 6, 2010, 8:55am

Thank you to all,
I solved the problem by avoiding the location option.

On 08/05/2010 04:14 PM, Giacomo Piva wrote:

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 Piva

MEEO S.r.l.
Via Saragat 9
I-44122 Ferrara, Italy
Tel: +39 0532 1861501
Fax: +39 0532 1861637
http://www.meeo.it

========================================
"ATTENZIONE:le informazioni contenute in questo messaggio sono
da considerarsi confidenziali ed il loro utilizzo è riservato unicamente
al destinatario sopra indicato. Chi dovesse ricevere questo messaggio
per errore è tenuto ad informare il mittente ed a rimuoverlo
definitivamente da ogni supporto elettronico o cartaceo."

"WARNING:This message contains confidential and/or proprietary
information which may be subject to privilege or immunity and which
is intended for use of its addressee only. Should you receive this
message in error, you are kindly requested to inform the sender and
to definitively remove it from any paper or electronic format."