[GRASS-user] Unable to open the Corine land cover map using GrassGIS 7

GRASS GRASS user
1

Hi,

I’m unable to open the Coring Land cover map (100m accuracy version) for Europe using GrassGIS, which can be downloaded from http://land.copernicus.eu/pan-european/corine-land-cover/clc-2012.

Apparentely the coordinate system of this land cover map isn’t supported (The projection is ETRS89 (EPSG:3035)), while the projection of the location I'm trying to import the raster to is WGS84 (EPSG:4326). Here is the error shown in GrassGIS: *ERROR: Coordinate reference system not available for input <C:\Data\Corine Land Cover\g100_clc12_V18_5a\g100_clc12_V18_5.tif>*

Is there a workaround in GrassGIS to open this image?
Otherwise which tool should I use to reproject it?

Thanks.

Hakim.

2

Hi Hakim,

Try r.import (https://grass.osgeo.org/grass72/manuals/r.import.html).

ETRS LAEA (EPSG:3035) is of course “supported”.

Cheers

Stefan

···

Hi,

I’m unable to open the Coring Land cover map (100m accuracy version) for Europe using GrassGIS, which can be downloaded from http://land.copernicus.eu/pan-european/corine-land-cover/clc-2012.

Apparentely the coordinate system of this land cover map isn’t supported (The projection is ETRS89 (EPSG:3035)), while the projection of the location I’m trying to import the raster to is WGS84 (EPSG:4326). Here is the error shown in GrassGIS: ERROR: Coordinate reference system not available for input <C:\Data\Corine Land Cover\g100_clc12_V18_5a\g100_clc12_V18_5.tif>

Is there a workaround in GrassGIS to open this image?

Otherwise which tool should I use to reproject it?

Thanks.

Hakim.

3

Hakim Benoudjit:

Hi,

I'm unable to open the Coring Land cover map (100m accuracy version) for
Europe using GrassGIS, which can be downloaded from
http://land.copernicus.eu/pan-european/corine-land-cover/clc-2012.
Apparentely the coordinate system of this land cover map isn't supported
(The projection is ETRS89 (EPSG:3035)), while the projection of the
location I'm trying to import the raster to is WGS84 (EPSG:4326). Here is
the error shown in GrassGIS: *ERROR: Coordinate reference system not
available for input <C:\Data\Corine Land
Cover\g100_clc12_V18_5a\g100_clc12_V18_5.tif>*

Is there a workaround in GrassGIS to open this image?
Otherwise which tool should I use to reproject it?

Dear Hakim, you can

1) create an ETRS89-based Location and then import the map in this
Location. If it is required to land in WGS84, there is r.proj that will
transform the raster map from one Location (say the "new" one based on
EPSG 3035) to the other (say the WGS84).

2) transform the map to wgs84 (use, for example GDAL tools -- perhaps
from inside QGIS)

It all depends what the further processing steps are. Better to avoid
transformations if precision is of interest.

If only quick viewing is of concern, consider using QGIS.

Some links to our Wiki:

https://grasswiki.osgeo.org/wiki/GRASS_Location_Wizard
https://grasswiki.osgeo.org/wiki/Map_Reprojection

Nikos

4

On Mon, Oct 9, 2017 at 4:47 PM, Nikos Alexandris <nik@nikosalexandris.net> wrote:

Hakim Benoudjit:

Hi,

I’m unable to open the Coring Land cover map (100m accuracy version) for
Europe using GrassGIS, which can be downloaded from
http://land.copernicus.eu/pan-european/corine-land-cover/clc-2012.
Apparentely the coordinate system of this land cover map isn’t supported
(The projection is ETRS89 (EPSG:3035)), while the projection of the
location I’m trying to import the raster to is WGS84 (EPSG:4326). Here is
the error shown in GrassGIS: ERROR: Coordinate reference system not
available for input <C:\Data\Corine Land
Cover\g100_clc12_V18_5a\g100_clc12_V18_5.tif>

As the error says, there is no coordinate reference system defined with the input. The coordinate reference system is only defined in non-standard metadata which are not readable for GDAL. Therefore you need to specify the CRS (EPSG:3035) manually, no matter what tool you are going to use (GRASS GIS/GDAL/QGIS).

Is there a workaround in GrassGIS to open this image?
Otherwise which tool should I use to reproject it?

Dear Hakim, you can

  1. create an ETRS89-based Location and then import the map in this
    Location.

You will need the -o flag to override projection check.

If it is required to land in WGS84, there is r.proj that will
transform the raster map from one Location (say the “new” one based on
EPSG 3035) to the other (say the WGS84).

  1. transform the map to wgs84 (use, for example GDAL tools – perhaps
    from inside QGIS)

That will only work if you manually assign a CRS to the input, e.g. with gdalwarp -s_srs EPSG3035

HTH,

Markus M

It all depends what the further processing steps are. Better to avoid
transformations if precision is of interest.

If only quick viewing is of concern, consider using QGIS.

Some links to our Wiki:

https://grasswiki.osgeo.org/wiki/GRASS_Location_Wizard
https://grasswiki.osgeo.org/wiki/Map_Reprojection

Nikos

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

5

On Mon, Oct 9, 2017 at 4:37 PM, Stefan Blumentrath <Stefan.Blumentrath@nina.no> wrote:

Hi Hakim,

Try r.import (https://grass.osgeo.org/grass72/manuals/r.import.html).

ETRS LAEA (EPSG:3035) is of course “supported”.

True, but a coordinate reference system not included in the input data, therefore r.import will not work and you need to define the CRS (EPSG:3035) manually.

Markus M

Cheers

Stefan

From: grass-user [mailto:grass-user-bounces@lists.osgeo.org] On Behalf Of Hakim Benoudjit
Sent: mandag 9. oktober 2017 16.02
To: grass-user@lists.osgeo.org
Subject: [GRASS-user] Unable to open the Corine land cover map using GrassGIS 7

Hi,

I’m unable to open the Coring Land cover map (100m accuracy version) for Europe using GrassGIS, which can be downloaded from http://land.copernicus.eu/pan-european/corine-land-cover/clc-2012.

Apparentely the coordinate system of this land cover map isn’t supported (The projection is ETRS89 (EPSG:3035)), while the projection of the location I’m trying to import the raster to is WGS84 (EPSG:4326). Here is the error shown in GrassGIS: ERROR: Coordinate reference system not available for input <C:\Data\Corine Land Cover\g100_clc12_V18_5a\g100_clc12_V18_5.tif>

Is there a workaround in GrassGIS to open this image?

Otherwise which tool should I use to reproject it?

Thanks.

Hakim.

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

6

You’re right I wasn’t using the gdal commands that come with OSGeo4W (with that version of gdal it’s working well), I was working with a version of gdal that I’ve compiled from source. So probably some flags had to be set during the compilation, but I forgot to do so.

Hakim.

···

On 9 October 2017 at 17:54, Markus Metz <markus.metz.giswork@gmail.com> wrote:

Hi Hakim,

On Mon, Oct 9, 2017 at 6:36 PM, Hakim Benoudjit <h.benoudjit@gmail.com> wrote:

Thanks Markus for your reponse.

I tried the gdalwarp command to reproject the land cover raster to WGS84 using the following command:

gdalwarp -s_srs EPSG:3035 -t_srs EPSG:4326 C:\Data\Corine-Land-Cover\g100_clc12_V18_5a\subset.tif C:\Data\Corine-Land-Cover\g100_clc12_V18_5a\subset-wgs84.tif

But I’m not sure how to provide the csv file it’s asking for to interpret the projection:
ERROR 4: Unable to open EPSG support file gcs.csv. Try setting the GDAL_DATA environment variable to point to the directory containing EPSG csv files.

The only csv file provided with this land cover archive is inside the Legend folder (clc_legend.csv), which gives the index of each type of land cover.
Do you have any idea where to locate this gcs.csv file it’s asking for, and how to provide it on Windows (should I add an environment variable set to the path of this file)?

this gcs.csv file is part of a standard GDAL installation. GDAL uses by default its own gcs.csv file, there should be no need to specify this file or GDAL_DATA explicitely.

It seems that something is wrong with the GDAL installation. Are you using gdalwarp from OSGeo4W?

Best,

Markus

Hakim.

On 9 October 2017 at 17:23, Markus Metz <markus.metz.giswork@gmail.com> wrote:

On Mon, Oct 9, 2017 at 4:37 PM, Stefan Blumentrath <Stefan.Blumentrath@nina.no> wrote:

Hi Hakim,

Try r.import (https://grass.osgeo.org/grass72/manuals/r.import.html).

ETRS LAEA (EPSG:3035) is of course “supported”.

True, but a coordinate reference system not included in the input data, therefore r.import will not work and you need to define the CRS (EPSG:3035) manually.

Markus M

Cheers

Stefan

From: grass-user [mailto:grass-user-bounces@lists.osgeo.org] On Behalf Of Hakim Benoudjit
Sent: mandag 9. oktober 2017 16.02
To: grass-user@lists.osgeo.org
Subject: [GRASS-user] Unable to open the Corine land cover map using GrassGIS 7

Hi,

I’m unable to open the Coring Land cover map (100m accuracy version) for Europe using GrassGIS, which can be downloaded from http://land.copernicus.eu/pan-european/corine-land-cover/clc-2012.

Apparentely the coordinate system of this land cover map isn’t supported (The projection is ETRS89 (EPSG:3035)), while the projection of the location I’m trying to import the raster to is WGS84 (EPSG:4326). Here is the error shown in GrassGIS: ERROR: Coordinate reference system not available for input <C:\Data\Corine Land Cover\g100_clc12_V18_5a\g100_clc12_V18_5.tif>

Is there a workaround in GrassGIS to open this image?

Otherwise which tool should I use to reproject it?

Thanks.

Hakim.

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

7

for the record, the problem is that for Corine Land Cover 2012, the Coordinate Reference System embedded in the raster data is

LOCAL_CS[“Geocoding information not available Projection Name = ETRS_1989_LAEA Units = meters GeoTIFF Units = ot”,
UNIT[“unknown”,1]]

while it should be

PROJCS[“ETRS89 / LAEA Europe”,
GEOGCS[“ETRS89”,
DATUM[“European_Terrestrial_Reference_System_1989”,
SPHEROID[“GRS 1980”,6378137,298.257222101,
AUTHORITY[“EPSG”,“7019”]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY[“EPSG”,“6258”]],
PRIMEM[“Greenwich”,0,
AUTHORITY[“EPSG”,“8901”]],
UNIT[“degree”,0.0174532925199433,
AUTHORITY[“EPSG”,“9122”]],
AUTHORITY[“EPSG”,“4258”]],
PROJECTION[“Lambert_Azimuthal_Equal_Area”],
PARAMETER[“latitude_of_center”,52],
PARAMETER[“longitude_of_center”,10],
PARAMETER[“false_easting”,4321000],
PARAMETER[“false_northing”,3210000],
UNIT[“metre”,1,
AUTHORITY[“EPSG”,“9001”]],
AUTHORITY[“EPSG”,“3035”]]

Apparently the creators of the data tried to hide the Coordinate Reference System in the GIS data available for download.

Markus M

On Mon, Oct 9, 2017 at 7:13 PM, Hakim Benoudjit <h.benoudjit@gmail.com> wrote:

You’re right I wasn’t using the gdal commands that come with OSGeo4W (with that version of gdal it’s working well), I was working with a version of gdal that I’ve compiled from source. So probably some flags had to be set during the compilation, but I forgot to do so.

Hakim.

On 9 October 2017 at 17:54, Markus Metz <markus.metz.giswork@gmail.com> wrote:

Hi Hakim,

On Mon, Oct 9, 2017 at 6:36 PM, Hakim Benoudjit <h.benoudjit@gmail.com> wrote:

Thanks Markus for your reponse.

I tried the gdalwarp command to reproject the land cover raster to WGS84 using the following command:

gdalwarp -s_srs EPSG:3035 -t_srs EPSG:4326 C:\Data\Corine-Land-Cover\g100_clc12_V18_5a\subset.tif C:\Data\Corine-Land-Cover\g100_clc12_V18_5a\subset-wgs84.tif

But I’m not sure how to provide the csv file it’s asking for to interpret the projection:
ERROR 4: Unable to open EPSG support file gcs.csv. Try setting the GDAL_DATA environment variable to point to the directory containing EPSG csv files.

The only csv file provided with this land cover archive is inside the Legend folder (clc_legend.csv), which gives the index of each type of land cover.
Do you have any idea where to locate this gcs.csv file it’s asking for, and how to provide it on Windows (should I add an environment variable set to the path of this file)?

this gcs.csv file is part of a standard GDAL installation. GDAL uses by default its own gcs.csv file, there should be no need to specify this file or GDAL_DATA explicitely.

It seems that something is wrong with the GDAL installation. Are you using gdalwarp from OSGeo4W?

Best,

Markus

Hakim.

On 9 October 2017 at 17:23, Markus Metz <markus.metz.giswork@gmail.com> wrote:

On Mon, Oct 9, 2017 at 4:37 PM, Stefan Blumentrath <Stefan.Blumentrath@nina.no> wrote:

Hi Hakim,

Try r.import (https://grass.osgeo.org/grass72/manuals/r.import.html).

ETRS LAEA (EPSG:3035) is of course “supported”.

True, but a coordinate reference system not included in the input data, therefore r.import will not work and you need to define the CRS (EPSG:3035) manually.

Markus M

Cheers

Stefan

From: grass-user [mailto:grass-user-bounces@lists.osgeo.org] On Behalf Of Hakim Benoudjit
Sent: mandag 9. oktober 2017 16.02
To: grass-user@lists.osgeo.org
Subject: [GRASS-user] Unable to open the Corine land cover map using GrassGIS 7

Hi,

I’m unable to open the Coring Land cover map (100m accuracy version) for Europe using GrassGIS, which can be downloaded from http://land.copernicus.eu/pan-european/corine-land-cover/clc-2012.

Apparentely the coordinate system of this land cover map isn’t supported (The projection is ETRS89 (EPSG:3035)), while the projection of the location I’m trying to import the raster to is WGS84 (EPSG:4326). Here is the error shown in GrassGIS: ERROR: Coordinate reference system not available for input <C:\Data\Corine Land Cover\g100_clc12_V18_5a\g100_clc12_V18_5.tif>

Is there a workaround in GrassGIS to open this image?

Otherwise which tool should I use to reproject it?

Thanks.

Hakim.

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

8

Hi,

I’m facing this same problem…

I have my EPSG:3035 location and when using r.import with -o flag to override projection check and import corine land cover, I get:

GRASS 7.3.svn (eu_laea):~> r.import -o input=g250_clc06_V18_5.tif output=clc06_V18_5_250m resample=nearest extent=region resolution=value resolution_value=250 title=“Corine Land Cover 2006 - 250m - V18.5”
ERROR: Unable to convert input map projection to GRASS format; cannot create new location.
ERROR: Unable to read GDAL dataset <g250_clc06_V18_5.tif>

while -o flag is supposed to override the check and import. This used to work before. Am I missunderstanding something here?

thanks a lot in advance,
Vero

···

2017-10-09 22:07 GMT+02:00 Markus Metz <markus.metz.giswork@gmail.com>:

for the record, the problem is that for Corine Land Cover 2012, the Coordinate Reference System embedded in the raster data is

LOCAL_CS[“Geocoding information not available Projection Name = ETRS_1989_LAEA Units = meters GeoTIFF Units = ot”,
UNIT[“unknown”,1]]

while it should be

PROJCS[“ETRS89 / LAEA Europe”,
GEOGCS[“ETRS89”,
DATUM[“European_Terrestrial_Reference_System_1989”,
SPHEROID[“GRS 1980”,6378137,298.257222101,
AUTHORITY[“EPSG”,“7019”]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY[“EPSG”,“6258”]],
PRIMEM[“Greenwich”,0,
AUTHORITY[“EPSG”,“8901”]],
UNIT[“degree”,0.0174532925199433,
AUTHORITY[“EPSG”,“9122”]],
AUTHORITY[“EPSG”,“4258”]],
PROJECTION[“Lambert_Azimuthal_Equal_Area”],
PARAMETER[“latitude_of_center”,52],
PARAMETER[“longitude_of_center”,10],
PARAMETER[“false_easting”,4321000],
PARAMETER[“false_northing”,3210000],
UNIT[“metre”,1,
AUTHORITY[“EPSG”,“9001”]],
AUTHORITY[“EPSG”,“3035”]]

Apparently the creators of the data tried to hide the Coordinate Reference System in the GIS data available for download.

Markus M

On Mon, Oct 9, 2017 at 7:13 PM, Hakim Benoudjit <h.benoudjit@gmail.com> wrote:

You’re right I wasn’t using the gdal commands that come with OSGeo4W (with that version of gdal it’s working well), I was working with a version of gdal that I’ve compiled from source. So probably some flags had to be set during the compilation, but I forgot to do so.

Hakim.

On 9 October 2017 at 17:54, Markus Metz <markus.metz.giswork@gmail.com> wrote:

Hi Hakim,

On Mon, Oct 9, 2017 at 6:36 PM, Hakim Benoudjit <h.benoudjit@gmail.com> wrote:

Thanks Markus for your reponse.

I tried the gdalwarp command to reproject the land cover raster to WGS84 using the following command:

gdalwarp -s_srs EPSG:3035 -t_srs EPSG:4326 C:\Data\Corine-Land-Cover\g100_clc12_V18_5a\subset.tif C:\Data\Corine-Land-Cover\g100_clc12_V18_5a\subset-wgs84.tif

But I’m not sure how to provide the csv file it’s asking for to interpret the projection:
ERROR 4: Unable to open EPSG support file gcs.csv. Try setting the GDAL_DATA environment variable to point to the directory containing EPSG csv files.

The only csv file provided with this land cover archive is inside the Legend folder (clc_legend.csv), which gives the index of each type of land cover.
Do you have any idea where to locate this gcs.csv file it’s asking for, and how to provide it on Windows (should I add an environment variable set to the path of this file)?

this gcs.csv file is part of a standard GDAL installation. GDAL uses by default its own gcs.csv file, there should be no need to specify this file or GDAL_DATA explicitely.

It seems that something is wrong with the GDAL installation. Are you using gdalwarp from OSGeo4W?

Best,

Markus

Hakim.

On 9 October 2017 at 17:23, Markus Metz <markus.metz.giswork@gmail.com> wrote:

On Mon, Oct 9, 2017 at 4:37 PM, Stefan Blumentrath <Stefan.Blumentrath@nina.no> wrote:

Hi Hakim,

Try r.import (https://grass.osgeo.org/grass72/manuals/r.import.html).

ETRS LAEA (EPSG:3035) is of course “supported”.

True, but a coordinate reference system not included in the input data, therefore r.import will not work and you need to define the CRS (EPSG:3035) manually.

Markus M

Cheers

Stefan

From: grass-user [mailto:grass-user-bounces@lists.osgeo.org] On Behalf Of Hakim Benoudjit
Sent: mandag 9. oktober 2017 16.02
To: grass-user@lists.osgeo.org
Subject: [GRASS-user] Unable to open the Corine land cover map using GrassGIS 7

Hi,

I’m unable to open the Coring Land cover map (100m accuracy version) for Europe using GrassGIS, which can be downloaded from http://land.copernicus.eu/pan-european/corine-land-cover/clc-2012.

Apparentely the coordinate system of this land cover map isn’t supported (The projection is ETRS89 (EPSG:3035)), while the projection of the location I’m trying to import the raster to is WGS84 (EPSG:4326). Here is the error shown in GrassGIS: ERROR: Coordinate reference system not available for input <C:\Data\Corine Land Cover\g100_clc12_V18_5a\g100_clc12_V18_5.tif>

Is there a workaround in GrassGIS to open this image?

Otherwise which tool should I use to reproject it?

Thanks.

Hakim.

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

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

9

I tested with:

r.in.gdal -o -r input=g250_clc06_V18_5.tif output=clc_test

-o works there and map gets imported in the epsg 3035 location…

vero

···

2017-10-24 16:53 GMT+02:00 Veronica Andreo <veroandreo@gmail.com>:

Hi,

I’m facing this same problem…

I have my EPSG:3035 location and when using r.import with -o flag to override projection check and import corine land cover, I get:

GRASS 7.3.svn (eu_laea):~> r.import -o input=g250_clc06_V18_5.tif output=clc06_V18_5_250m resample=nearest extent=region resolution=value resolution_value=250 title=“Corine Land Cover 2006 - 250m - V18.5”
ERROR: Unable to convert input map projection to GRASS format; cannot create new location.
ERROR: Unable to read GDAL dataset <g250_clc06_V18_5.tif>

while -o flag is supposed to override the check and import. This used to work before. Am I missunderstanding something here?

thanks a lot in advance,
Vero

2017-10-09 22:07 GMT+02:00 Markus Metz <markus.metz.giswork@gmail.com>:

for the record, the problem is that for Corine Land Cover 2012, the Coordinate Reference System embedded in the raster data is

LOCAL_CS[“Geocoding information not available Projection Name = ETRS_1989_LAEA Units = meters GeoTIFF Units = ot”,
UNIT[“unknown”,1]]

while it should be

PROJCS[“ETRS89 / LAEA Europe”,
GEOGCS[“ETRS89”,
DATUM[“European_Terrestrial_Reference_System_1989”,
SPHEROID[“GRS 1980”,6378137,298.257222101,
AUTHORITY[“EPSG”,“7019”]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY[“EPSG”,“6258”]],
PRIMEM[“Greenwich”,0,
AUTHORITY[“EPSG”,“8901”]],
UNIT[“degree”,0.0174532925199433,
AUTHORITY[“EPSG”,“9122”]],
AUTHORITY[“EPSG”,“4258”]],
PROJECTION[“Lambert_Azimuthal_Equal_Area”],
PARAMETER[“latitude_of_center”,52],
PARAMETER[“longitude_of_center”,10],
PARAMETER[“false_easting”,4321000],
PARAMETER[“false_northing”,3210000],
UNIT[“metre”,1,
AUTHORITY[“EPSG”,“9001”]],
AUTHORITY[“EPSG”,“3035”]]

Apparently the creators of the data tried to hide the Coordinate Reference System in the GIS data available for download.

Markus M

On Mon, Oct 9, 2017 at 7:13 PM, Hakim Benoudjit <h.benoudjit@gmail.com> wrote:

You’re right I wasn’t using the gdal commands that come with OSGeo4W (with that version of gdal it’s working well), I was working with a version of gdal that I’ve compiled from source. So probably some flags had to be set during the compilation, but I forgot to do so.

Hakim.

On 9 October 2017 at 17:54, Markus Metz <markus.metz.giswork@gmail.com> wrote:

Hi Hakim,

On Mon, Oct 9, 2017 at 6:36 PM, Hakim Benoudjit <h.benoudjit@gmail.com> wrote:

Thanks Markus for your reponse.

I tried the gdalwarp command to reproject the land cover raster to WGS84 using the following command:

gdalwarp -s_srs EPSG:3035 -t_srs EPSG:4326 C:\Data\Corine-Land-Cover\g100_clc12_V18_5a\subset.tif C:\Data\Corine-Land-Cover\g100_clc12_V18_5a\subset-wgs84.tif

But I’m not sure how to provide the csv file it’s asking for to interpret the projection:
ERROR 4: Unable to open EPSG support file gcs.csv. Try setting the GDAL_DATA environment variable to point to the directory containing EPSG csv files.

The only csv file provided with this land cover archive is inside the Legend folder (clc_legend.csv), which gives the index of each type of land cover.
Do you have any idea where to locate this gcs.csv file it’s asking for, and how to provide it on Windows (should I add an environment variable set to the path of this file)?

this gcs.csv file is part of a standard GDAL installation. GDAL uses by default its own gcs.csv file, there should be no need to specify this file or GDAL_DATA explicitely.

It seems that something is wrong with the GDAL installation. Are you using gdalwarp from OSGeo4W?

Best,

Markus

Hakim.

On 9 October 2017 at 17:23, Markus Metz <markus.metz.giswork@gmail.com> wrote:

On Mon, Oct 9, 2017 at 4:37 PM, Stefan Blumentrath <Stefan.Blumentrath@nina.no> wrote:

Hi Hakim,

Try r.import (https://grass.osgeo.org/grass72/manuals/r.import.html).

ETRS LAEA (EPSG:3035) is of course “supported”.

True, but a coordinate reference system not included in the input data, therefore r.import will not work and you need to define the CRS (EPSG:3035) manually.

Markus M

Cheers

Stefan

From: grass-user [mailto:grass-user-bounces@lists.osgeo.org] On Behalf Of Hakim Benoudjit
Sent: mandag 9. oktober 2017 16.02
To: grass-user@lists.osgeo.org
Subject: [GRASS-user] Unable to open the Corine land cover map using GrassGIS 7

Hi,

I’m unable to open the Coring Land cover map (100m accuracy version) for Europe using GrassGIS, which can be downloaded from http://land.copernicus.eu/pan-european/corine-land-cover/clc-2012.

Apparentely the coordinate system of this land cover map isn’t supported (The projection is ETRS89 (EPSG:3035)), while the projection of the location I’m trying to import the raster to is WGS84 (EPSG:4326). Here is the error shown in GrassGIS: ERROR: Coordinate reference system not available for input <C:\Data\Corine Land Cover\g100_clc12_V18_5a\g100_clc12_V18_5.tif>

Is there a workaround in GrassGIS to open this image?

Otherwise which tool should I use to reproject it?

Thanks.

Hakim.

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

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

10

Veronica Andreo wrote

Hi,

I'm facing this same problem...

I have my EPSG:3035 location and when using r.import with -o flag to
override projection check and import corine land cover, I get:

GRASS 7.3.svn (eu_laea):~> r.import -o input=g250_clc06_V18_5.tif
output=clc06_V18_5_250m resample=nearest extent=region resolution=value
resolution_value=250 title="Corine Land Cover 2006 - 250m - V18.5"
ERROR: Unable to convert input map projection to GRASS format; cannot
create new location.
ERROR: Unable to read GDAL dataset
<g250_clc06_V18_5.tif>
while -o flag is supposed to override the check and import. This used to
work before. Am I missunderstanding something here?

confirmed

r.import input=D:\dl\ne\g250_clc12_V18_5a\g250_clc12_V18_5.tif
output=g250_clc12_V18_5 -o
ERROR: Unable to convert input map projection to GRASS format; cannot create
new location.
ERROR: Unable to read GDAL dataset
<D:\dl\ne\g250_clc12_V18_5a\g250_clc12_V18_5.tif>

where does 'cannot create new location' comes from?

-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

11

Helmut Kudrnovsky wrote

Veronica Andreo wrote

Hi,

I'm facing this same problem...

I have my EPSG:3035 location and when using r.import with -o flag to
override projection check and import corine land cover, I get:

GRASS 7.3.svn (eu_laea):~> r.import -o input=g250_clc06_V18_5.tif
output=clc06_V18_5_250m resample=nearest extent=region resolution=value
resolution_value=250 title="Corine Land Cover 2006 - 250m - V18.5"
ERROR: Unable to convert input map projection to GRASS format; cannot
create new location.
ERROR: Unable to read GDAL dataset

<g250_clc06_V18_5.tif>

while -o flag is supposed to override the check and import. This used to
work before. Am I missunderstanding something here?

confirmed

r.import input=D:\dl\ne\g250_clc12_V18_5a\g250_clc12_V18_5.tif
output=g250_clc12_V18_5 -o
ERROR: Unable to convert input map projection to GRASS format; cannot
create
new location.
ERROR: Unable to read GDAL dataset
&lt;D:\dl\ne\g250_clc12_V18_5a\g250_clc12_V18_5.tif&gt;

where does 'cannot create new location' comes from?

grass7.2.2

r.in.gdal -o input=D:\dl\ne\g250_clc12_V18_5a\g250_clc12_V18_5.tif
output=fhfghgfh
WARNING: Over-riding projection check
Proceeding with import of 1 raster bands...
Importing raster map <fhfghgfh>...

grass7.3.svn

r.in.gdal -o input=D:\dl\ne\g250_clc12_V18_5a\g250_clc12_V18_5.tif
output=sgdfg
Input contains an invalid SRS. WKT definition:
LOCAL_CS["Geocoding information not available Projection Name =
ETRS_1989_LAEA Units = meters GeoTIFF Units = ot",UNIT["unknown",1]]
WARNING: Unable to convert input map projection information to GRASS format.
Over-riding projection check
Importing raster map <sgdfg>...

-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

12

Helmut Kudrnovsky wrote

Veronica Andreo wrote

Hi,

I'm facing this same problem...

I have my EPSG:3035 location and when using r.import with -o flag to
override projection check and import corine land cover, I get:

GRASS 7.3.svn (eu_laea):~> r.import -o input=g250_clc06_V18_5.tif
output=clc06_V18_5_250m resample=nearest extent=region resolution=value
resolution_value=250 title="Corine Land Cover 2006 - 250m - V18.5"
ERROR: Unable to convert input map projection to GRASS format; cannot
create new location.
ERROR: Unable to read GDAL dataset

<g250_clc06_V18_5.tif>

while -o flag is supposed to override the check and import. This used to
work before. Am I missunderstanding something here?

confirmed

r.import input=D:\dl\ne\g250_clc12_V18_5a\g250_clc12_V18_5.tif
output=g250_clc12_V18_5 -o
ERROR: Unable to convert input map projection to GRASS format; cannot
create
new location.
ERROR: Unable to read GDAL dataset
&lt;D:\dl\ne\g250_clc12_V18_5a\g250_clc12_V18_5.tif&gt;

where does 'cannot create new location' comes from?

DEBUG=3

r.import input=D:\dl\ne\g250_clc12_V18_5a\g250_clc12_V18_5.tif
output=debugimport -o
D2/3: filename =
C:\OSGEO4~1\apps\grass\grass-7.3.svn/scripts/r.import.py
D1/3: G_set_program_name(): r.import
D2/3: G_file_name(): path = D:\grassdata\test_loc3035\data
D2/3: G_file_name(): path =
D:\grassdata\test_loc3035\data\cell\debugimport
D1/3: G_set_program_name(): g.gisenv
D3/3: G_option_to_separator(): key = separator -> sep = '
'
D1/3: grass.script.core.start_command(): g.gisenv -n
D1/3: G_set_program_name(): g.gisenv
D3/3: G_option_to_separator(): key = separator -> sep = '
'
D1/3: grass.script.core.start_command(): g.tempfile
pid=13236
D1/3: G_set_program_name(): g.tempfile
D2/3: G_file_name(): path = D:\grassdata\test_loc3035\data
D2/3: G_file_name(): path = D:\grassdata\test_loc3035\data
D2/3: G__temp_element(): .tmp/unknown (tmp=0)
D2/3: G_file_name(): path =
D:\grassdata\test_loc3035\data\.tmp/unknown\13236.0
D2/3: G_tempfile_pid():
D:\grassdata\test_loc3035\data\.tmp/unknown\13236.0
D1/3: grass.script.core.start_command(): g.proj --q -j
D1/3: G_set_program_name(): g.proj
D2/3: G_file_name(): path =
D:\grassdata\test_loc3035\PERMANENT\DEFAULT_WIND
D2/3: G_file_name(): path =
D:\grassdata\test_loc3035\PERMANENT\DEFAULT_WIND
D2/3: file open: read (mode = r)
D2/3: G__read_Cell_head
D2/3: G__read_Cell_head_array
D3/3: region item: proj: 99
D3/3: region item: zone: 0
D3/3: region item: north: 1
D3/3: region item: south: 0
D3/3: region item: east: 1
D3/3: region item: west: 0
D3/3: region item: cols: 1
D3/3: region item: rows: 1
D3/3: region item: e-w resol: 1
D3/3: region item: n-s resol: 1
D3/3: region item: top: 1.000000000000000
D3/3: region item: bottom: 0.000000000000000
D3/3: region item: cols3: 1
D3/3: region item: rows3: 1
D3/3: region item: depths: 1
D3/3: region item: e-w resol3: 1
D3/3: region item: n-s resol3: 1
D3/3: region item: t-b resol: 1
D2/3: G_file_name(): path =
D:\grassdata\test_loc3035\PERMANENT\PROJ_INFO
D2/3: G_file_name(): path =
D:\grassdata\test_loc3035\PERMANENT\PROJ_UNITS
D2/3: G_file_name(): path =
D:\grassdata\test_loc3035\PERMANENT\PROJ_EPSG
D3/3: GPJ__get_datum_params: datumname: <etrs89>
D3/3: set_datumtrans(): GPJ__get_datum_params() status=2
D3/3: set_datumtrans(): datum transform terms found with 1
options
D3/3: GPJ__get_datum_params: datumname: <etrs89>
D1/3: grass.script.core.start_command(): g.message -v
message=Creating temporary location for
<D:\dl\ne\g250_clc12_V18_5a\g250_clc12_V18_5.tif>...
D1/3: G_set_program_name(): g.message
D2/3: G_file_name(): path = D:\grassdata\test_loc3035\data
D1/3: grass.script.core.start_command(): r.in.gdal --q -c
output=debugimport
input=D:\dl\ne\g250_clc12_V18_5a\g250_clc12_V18_5.tif title=
location=temp_import_location_13236 memory=300
D1/3: G_set_program_name(): r.in.gdal
D2/3: G_file_name(): path = D:\grassdata\test_loc3035\data
D3/3: GDAL Driver: GTiff
D3/3: GDAL size: row = 34360, col = 50992
ERROR: Unable to convert input map projection to GRASS format; cannot create
new location.
D1/3: grass.script.core.start_command(): g.message -e
message=Unable to read GDAL dataset
<D:\dl\ne\g250_clc12_V18_5a\g250_clc12_V18_5.tif>
D1/3: G_set_program_name(): g.message
D2/3: G_file_name(): path = D:\grassdata\test_loc3035\data
ERROR: Unable to read GDAL dataset
<D:\dl\ne\g250_clc12_V18_5a\g250_clc12_V18_5.tif>

-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

13

On Tue, Oct 24, 2017 at 4:53 PM, Veronica Andreo <veroandreo@gmail.com> wrote:

Hi,

I’m facing this same problem…

I have my EPSG:3035 location and when using r.import with -o flag to override projection check and import corine land cover, I get:

GRASS 7.3.svn (eu_laea):~> r.import -o input=g250_clc06_V18_5.tif output=clc06_V18_5_250m resample=nearest extent=region resolution=value resolution_value=250 title=“Corine Land Cover 2006 - 250m - V18.5”
ERROR: Unable to convert input map projection to GRASS format; cannot create new location.
ERROR: Unable to read GDAL dataset <g250_clc06_V18_5.tif>

while -o flag is supposed to override the check and import. This used to work before. Am I missunderstanding something here?

r.import needs to be adjusted to import the raster map directly into the current location if the -o flag is given.

The underlying reason is that g250_clc06_V18_5.tif has projection information that is invalid. This invalid projection information is printed out as WKT before the error messages. You can also check with gdalinfo g250_clc06_V18_5.tif

Markus M

thanks a lot in advance,
Vero

2017-10-09 22:07 GMT+02:00 Markus Metz <markus.metz.giswork@gmail.com>:

for the record, the problem is that for Corine Land Cover 2012, the Coordinate Reference System embedded in the raster data is

LOCAL_CS[“Geocoding information not available Projection Name = ETRS_1989_LAEA Units = meters GeoTIFF Units = ot”,
UNIT[“unknown”,1]]

while it should be

PROJCS[“ETRS89 / LAEA Europe”,
GEOGCS[“ETRS89”,
DATUM[“European_Terrestrial_Reference_System_1989”,
SPHEROID[“GRS 1980”,6378137,298.257222101,
AUTHORITY[“EPSG”,“7019”]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY[“EPSG”,“6258”]],
PRIMEM[“Greenwich”,0,
AUTHORITY[“EPSG”,“8901”]],
UNIT[“degree”,0.0174532925199433,
AUTHORITY[“EPSG”,“9122”]],
AUTHORITY[“EPSG”,“4258”]],
PROJECTION[“Lambert_Azimuthal_Equal_Area”],
PARAMETER[“latitude_of_center”,52],
PARAMETER[“longitude_of_center”,10],
PARAMETER[“false_easting”,4321000],
PARAMETER[“false_northing”,3210000],
UNIT[“metre”,1,
AUTHORITY[“EPSG”,“9001”]],
AUTHORITY[“EPSG”,“3035”]]

Apparently the creators of the data tried to hide the Coordinate Reference System in the GIS data available for download.

Markus M

On Mon, Oct 9, 2017 at 7:13 PM, Hakim Benoudjit <h.benoudjit@gmail.com> wrote:

You’re right I wasn’t using the gdal commands that come with OSGeo4W (with that version of gdal it’s working well), I was working with a version of gdal that I’ve compiled from source. So probably some flags had to be set during the compilation, but I forgot to do so.

Hakim.

On 9 October 2017 at 17:54, Markus Metz <markus.metz.giswork@gmail.com> wrote:

Hi Hakim,

On Mon, Oct 9, 2017 at 6:36 PM, Hakim Benoudjit <h.benoudjit@gmail.com> wrote:

Thanks Markus for your reponse.

I tried the gdalwarp command to reproject the land cover raster to WGS84 using the following command:

gdalwarp -s_srs EPSG:3035 -t_srs EPSG:4326 C:\Data\Corine-Land-Cover\g100_clc12_V18_5a\subset.tif C:\Data\Corine-Land-Cover\g100_clc12_V18_5a\subset-wgs84.tif

But I’m not sure how to provide the csv file it’s asking for to interpret the projection:
ERROR 4: Unable to open EPSG support file gcs.csv. Try setting the GDAL_DATA environment variable to point to the directory containing EPSG csv files.

The only csv file provided with this land cover archive is inside the Legend folder (clc_legend.csv), which gives the index of each type of land cover.
Do you have any idea where to locate this gcs.csv file it’s asking for, and how to provide it on Windows (should I add an environment variable set to the path of this file)?

this gcs.csv file is part of a standard GDAL installation. GDAL uses by default its own gcs.csv file, there should be no need to specify this file or GDAL_DATA explicitely.

It seems that something is wrong with the GDAL installation. Are you using gdalwarp from OSGeo4W?

Best,

Markus

Hakim.

On 9 October 2017 at 17:23, Markus Metz <markus.metz.giswork@gmail.com> wrote:

On Mon, Oct 9, 2017 at 4:37 PM, Stefan Blumentrath <Stefan.Blumentrath@nina.no> wrote:

Hi Hakim,

Try r.import (https://grass.osgeo.org/grass72/manuals/r.import.html).

ETRS LAEA (EPSG:3035) is of course “supported”.

True, but a coordinate reference system not included in the input data, therefore r.import will not work and you need to define the CRS (EPSG:3035) manually.

Markus M

Cheers

Stefan

From: grass-user [mailto:grass-user-bounces@lists.osgeo.org] On Behalf Of Hakim Benoudjit
Sent: mandag 9. oktober 2017 16.02
To: grass-user@lists.osgeo.org
Subject: [GRASS-user] Unable to open the Corine land cover map using GrassGIS 7

Hi,

I’m unable to open the Coring Land cover map (100m accuracy version) for Europe using GrassGIS, which can be downloaded from http://land.copernicus.eu/pan-european/corine-land-cover/clc-2012.

Apparentely the coordinate system of this land cover map isn’t supported (The projection is ETRS89 (EPSG:3035)), while the projection of the location I’m trying to import the raster to is WGS84 (EPSG:4326). Here is the error shown in GrassGIS: ERROR: Coordinate reference system not available for input <C:\Data\Corine Land Cover\g100_clc12_V18_5a\g100_clc12_V18_5.tif>

Is there a workaround in GrassGIS to open this image?

Otherwise which tool should I use to reproject it?

Thanks.

Hakim.

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

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

14

Markus Metz-3 wrote

On Tue, Oct 24, 2017 at 4:53 PM, Veronica Andreo &lt;

veroandreo@

&gt;
wrote:

Hi,

I'm facing this same problem...

I have my EPSG:3035 location and when using r.import with -o flag to

override projection check and import corine land cover, I get:

GRASS 7.3.svn (eu_laea):~> r.import -o input=g250_clc06_V18_5.tif

output=clc06_V18_5_250m resample=nearest extent=region resolution=value
resolution_value=250 title="Corine Land Cover 2006 - 250m - V18.5"

ERROR: Unable to convert input map projection to GRASS format; cannot

create new location.

ERROR: Unable to read GDAL dataset

<g250_clc06_V18_5.tif>

while -o flag is supposed to override the check and import. This used to

work before. Am I missunderstanding something here?

r.import needs to be adjusted to import the raster map directly into the
current location if the -o flag is given.

The underlying reason is that g250_clc06_V18_5.tif has projection
information that is invalid. This invalid projection information is
printed
out as WKT before the error messages. You can also check with gdalinfo
g250_clc06_V18_5.tif

just reported the invalid SRS to the data provider.

-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

15

On Tue, Oct 24, 2017 at 5:15 PM, Helmut Kudrnovsky <hellik@web.de> wrote:

Veronica Andreo wrote

Hi,

I’m facing this same problem…

I have my EPSG:3035 location and when using r.import with -o flag to
override projection check and import corine land cover, I get:

GRASS 7.3.svn (eu_laea):~> r.import -o input=g250_clc06_V18_5.tif
output=clc06_V18_5_250m resample=nearest extent=region resolution=value
resolution_value=250 title=“Corine Land Cover 2006 - 250m - V18.5”
ERROR: Unable to convert input map projection to GRASS format; cannot
create new location.
ERROR: Unable to read GDAL dataset
<g250_clc06_V18_5.tif>
while -o flag is supposed to override the check and import. This used to
work before. Am I missunderstanding something here?

confirmed

r.import input=D:\dl\ne\g250_clc12_V18_5a\g250_clc12_V18_5.tif
output=g250_clc12_V18_5 -o
ERROR: Unable to convert input map projection to GRASS format; cannot create
new location.
ERROR: Unable to read GDAL dataset
<D:\dl\ne\g250_clc12_V18_5a\g250_clc12_V18_5.tif>

where does ‘cannot create new location’ comes from?

r.import / v.import tried to create a temporary location even if it is not needed and direct import works. Fixed in trunk r71590.

Markus M

best regards
Helmut

Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

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

16

Helmut Kudrnovsky wrote

Markus Metz-3 wrote

On Tue, Oct 24, 2017 at 4:53 PM, Veronica Andreo &lt;

veroandreo@

&gt;
wrote:

Hi,

I'm facing this same problem...

I have my EPSG:3035 location and when using r.import with -o flag to

override projection check and import corine land cover, I get:

GRASS 7.3.svn (eu_laea):~> r.import -o input=g250_clc06_V18_5.tif

output=clc06_V18_5_250m resample=nearest extent=region resolution=value
resolution_value=250 title="Corine Land Cover 2006 - 250m - V18.5"

ERROR: Unable to convert input map projection to GRASS format; cannot

create new location.

ERROR: Unable to read GDAL dataset

<g250_clc06_V18_5.tif>

while -o flag is supposed to override the check and import. This used to

work before. Am I missunderstanding something here?

r.import needs to be adjusted to import the raster map directly into the
current location if the -o flag is given.

The underlying reason is that g250_clc06_V18_5.tif has projection
information that is invalid. This invalid projection information is
printed
out as WKT before the error messages. You can also check with gdalinfo
g250_clc06_V18_5.tif

just reported the invalid SRS to the data provider.

for the record, just get feedback by the data provider:

the invalid/incorrect embedded SRS will be corrected in the CLC2018 update.

-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html