<div dir="ltr"><div><div><div><br><br>On Tue, Oct 24, 2017 at 4:53 PM, Veronica Andreo <<a href="mailto:veroandreo@gmail.com">veroandreo@gmail.com</a>> wrote:<br>><br>> Hi, <br>><br>> I'm facing this same problem... <br>><br>> 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: <br>><br>> 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"<br>> ERROR: Unable to convert input map projection to GRASS format; cannot create new location.<br>> ERROR: Unable to read GDAL dataset <g250_clc06_V18_5.tif><br>><br>> while -o flag is supposed to override the check and import. This used to work before. Am I missunderstanding something here?<br><br></div>r.import needs to be adjusted to import the raster map directly into the current location if the -o flag is given.<br><br></div>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<br></div><div><br></div>Markus M<br><div><div><div><br>><br>> thanks a lot in advance, <br>> Vero<br>><br>> 2017-10-09 22:07 GMT+02:00 Markus Metz <<a href="mailto:markus.metz.giswork@gmail.com">markus.metz.giswork@gmail.com</a>>:<br>>><br>>> for the record, the problem is that for Corine Land Cover 2012, the Coordinate Reference System embedded in the raster data is<br>>><br>>> LOCAL_CS["Geocoding information not available Projection Name = ETRS_1989_LAEA Units = meters GeoTIFF Units = ot",<br>>>     UNIT["unknown",1]]<br>>><br>>> while it should be<br>>><br>>> PROJCS["ETRS89 / LAEA Europe",<br>>>     GEOGCS["ETRS89",<br>>>         DATUM["European_Terrestrial_Reference_System_1989",<br>>>             SPHEROID["GRS 1980",6378137,298.257222101,<br>>>                 AUTHORITY["EPSG","7019"]],<br>>>             TOWGS84[0,0,0,0,0,0,0],<br>>>             AUTHORITY["EPSG","6258"]],<br>>>         PRIMEM["Greenwich",0,<br>>>             AUTHORITY["EPSG","8901"]],<br>>>         UNIT["degree",0.0174532925199433,<br>>>             AUTHORITY["EPSG","9122"]],<br>>>         AUTHORITY["EPSG","4258"]],<br>>>     PROJECTION["Lambert_Azimuthal_Equal_Area"],<br>>>     PARAMETER["latitude_of_center",52],<br>>>     PARAMETER["longitude_of_center",10],<br>>>     PARAMETER["false_easting",4321000],<br>>>     PARAMETER["false_northing",3210000],<br>>>     UNIT["metre",1,<br>>>         AUTHORITY["EPSG","9001"]],<br>>>     AUTHORITY["EPSG","3035"]]<br>>><br>>> Apparently the creators of the data tried to hide the Coordinate Reference System in the GIS data available for download.<br>>><br>>> Markus M<br>>><br>>> On Mon, Oct 9, 2017 at 7:13 PM, Hakim Benoudjit <<a href="mailto:h.benoudjit@gmail.com">h.benoudjit@gmail.com</a>> wrote:<br>>> ><br>>> > 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.<br>>> ><br>>> > Hakim.<br>>> ><br>>> ><br>>> > On 9 October 2017 at 17:54, Markus Metz <<a href="mailto:markus.metz.giswork@gmail.com">markus.metz.giswork@gmail.com</a>> wrote:<br>>> >><br>>> >> Hi Hakim,<br>>> >><br>>> >> On Mon, Oct 9, 2017 at 6:36 PM, Hakim Benoudjit <<a href="mailto:h.benoudjit@gmail.com">h.benoudjit@gmail.com</a>> wrote:<br>>> >> ><br>>> >> > Thanks Markus for your reponse.<br>>> >> ><br>>> >> > I tried the gdalwarp command to reproject the land cover raster to WGS84 using the following command:<br>>> >> > > 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<br>>> >> ><br>>> >> > But I'm not sure how to provide the csv file it's asking for to interpret the projection:<br>>> >> > 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.<br>>> >> ><br>>> >> > 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.<br>>> >> > 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)?<br>>> >><br>>> >> 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.<br>>> >><br>>> >> It seems that something is wrong with the GDAL installation. Are you using gdalwarp from OSGeo4W?<br>>> >><br>>> >> Best,<br>>> >><br>>> >> Markus<br>>> >><br>>> >><br>>> >> ><br>>> >> > Hakim.<br>>><br>>> >> ><br>>> >> > On 9 October 2017 at 17:23, Markus Metz <<a href="mailto:markus.metz.giswork@gmail.com">markus.metz.giswork@gmail.com</a>> wrote:<br>>> >> >><br>>> >> >><br>>> >> >><br>>> >> >> On Mon, Oct 9, 2017 at 4:37 PM, Stefan Blumentrath <<a href="mailto:Stefan.Blumentrath@nina.no">Stefan.Blumentrath@nina.no</a>> wrote:<br>>> >> >> ><br>>> >> >> > Hi Hakim,<br>>> >> >> ><br>>> >> >> >  <br>>> >> >> ><br>>> >> >> > Try r.import (<a href="https://grass.osgeo.org/grass72/manuals/r.import.html">https://grass.osgeo.org/grass72/manuals/r.import.html</a>).<br>>> >> >> ><br>>> >> >> > ETRS LAEA (EPSG:3035) is of course “supported”.<br>>> >> >><br>>> >> >> 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.<br>>> >> >><br>>> >> >> Markus M<br>>> >> >> ><br>>> >> >> >  <br>>> >> >> ><br>>> >> >> > Cheers<br>>> >> >> ><br>>> >> >> > Stefan<br>>> >> >> ><br>>> >> >> >  <br>>> >> >> ><br>>> >> >> > From: grass-user [mailto:<a href="mailto:grass-user-bounces@lists.osgeo.org">grass-user-bounces@lists.osgeo.org</a>] On Behalf Of Hakim Benoudjit<br>>> >> >> > Sent: mandag 9. oktober 2017 16.02<br>>> >> >> > To: <a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>>> >> >> > Subject: [GRASS-user] Unable to open the Corine land cover map using GrassGIS 7<br>>> >> >> ><br>>> >> >> >  <br>>> >> >> ><br>>> >> >> > Hi,<br>>> >> >> ><br>>> >> >> > I'm unable to open the Coring Land cover map (100m accuracy version) for Europe using GrassGIS, which can be downloaded from <a href="http://land.copernicus.eu/pan-european/corine-land-cover/clc-2012">http://land.copernicus.eu/pan-european/corine-land-cover/clc-2012</a>.<br>>> >> >> ><br>>> >> >> > 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><br>>> >> >> ><br>>> >> >> >  <br>>> >> >> ><br>>> >> >> > Is there a workaround in GrassGIS to open this image?<br>>> >> >> ><br>>> >> >> > Otherwise which tool should I use to reproject it?<br>>> >> >> ><br>>> >> >> > Thanks.<br>>> >> >> ><br>>> >> >> > Hakim.<br>>> >> >> ><br>>> >> >> ><br>>> >> >> > _______________________________________________<br>>> >> >> > grass-user mailing list<br>>> >> >> > <a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>>> >> >> > <a href="https://lists.osgeo.org/mailman/listinfo/grass-user">https://lists.osgeo.org/mailman/listinfo/grass-user</a><br>>> >> >><br>>> >> ><br>>> ><br>>> ><br>>><br>>> _______________________________________________<br>>> grass-user mailing list<br>>> <a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>>> <a href="https://lists.osgeo.org/mailman/listinfo/grass-user">https://lists.osgeo.org/mailman/listinfo/grass-user</a><br>><br>><br></div></div></div></div>