[GRASS-user] Unable to open the Corine land cover map using GrassGIS 7
Veronica Andreo
veroandreo at gmail.com
Tue Oct 24 08:14:41 PDT 2017
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 at 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 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/grass
>> 72/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_c
>> lc12_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/9270604e/attachment-0001.html>
More information about the grass-user
mailing list