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

Tue Oct 24 11:33:22 PDT 2017

On Tue, Oct 24, 2017 at 4:53 PM, Veronica Andreo <veroandreo at 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 at 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 at 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 at gmail.com>
wrote:
>> >>
>> >> Hi Hakim,
>> >>
>> >> On Mon, Oct 9, 2017 at 6:36 PM, Hakim Benoudjit <h.benoudjit at 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 at gmail.com> wrote:
>> >> >>
>> >> >>
>> >> >>
>> >> >> On Mon, Oct 9, 2017 at 4:37 PM, Stefan Blumentrath <
Stefan.Blumentrath at 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 at lists.osgeo.org] On
Behalf Of Hakim Benoudjit
>> >> >> > Sent: mandag 9. oktober 2017 16.02
>> >> >> > To: grass-user at 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 at lists.osgeo.org
>> >> >> > https://lists.osgeo.org/mailman/listinfo/grass-user
>> >> >>
>> >> >
>> >
>> >
>>
>> _______________________________________________
>> grass-user mailing list
>> grass-user at lists.osgeo.org
>> https://lists.osgeo.org/mailman/listinfo/grass-user
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20171024/f74bde4f/attachment-0001.html>