<div dir="ltr">Hi, <div><br></div><div>I'm facing this same problem... </div><div><br></div><div>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: </div><div><div><br></div><div>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"</div><div>ERROR: Unable to convert input map projection to GRASS format; cannot create new location.</div><div>ERROR: Unable to read GDAL dataset <g250_clc06_V18_5.tif></div></div><div><br></div><div>while -o flag is supposed to override the check and import. This used to work before. Am I missunderstanding something here?</div><div><br></div><div>thanks a lot in advance, </div><div>Vero</div></div><div class="gmail_extra"><br><div class="gmail_quote">2017-10-09 22:07 GMT+02:00 Markus Metz <span dir="ltr"><<a href="mailto:markus.metz.giswork@gmail.com" target="_blank">markus.metz.giswork@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>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></div>while it should be<br><br>PROJCS["ETRS89 / LAEA Europe",<br>    GEOGCS["ETRS89",<br>        DATUM["European_Terrestrial_<wbr>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.<wbr>0174532925199433,<br>            AUTHORITY["EPSG","9122"]],<br>        AUTHORITY["EPSG","4258"]],<br>    PROJECTION["Lambert_Azimuthal_<wbr>Equal_Area"],<br>    PARAMETER["latitude_of_center"<wbr>,52],<br>    PARAMETER["longitude_of_<wbr>center",10],<br>    PARAMETER["false_easting",<wbr>4321000],<br>    PARAMETER["false_northing",<wbr>3210000],<br>    UNIT["metre",1,<br>        AUTHORITY["EPSG","9001"]],<br>    AUTHORITY["EPSG","3035"]]<br><br></div>Apparently the creators of the data tried to hide the Coordinate Reference System in the GIS data available for download.<br><br></div>Markus M<br><div><div><div><div><br>On Mon, Oct 9, 2017 at 7:13 PM, Hakim Benoudjit <<a href="mailto:h.benoudjit@gmail.com" target="_blank">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" target="_blank">markus.metz.giswork@gmail.com</a><wbr>> 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" target="_blank">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\<wbr>g100_clc12_V18_5a\subset.tif C:\Data\Corine-Land-Cover\<wbr>g100_clc12_V18_5a\subset-<wbr>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.<div><div class="h5"><br>>> ><br>>> > On 9 October 2017 at 17:23, Markus Metz <<a href="mailto:markus.metz.giswork@gmail.com" target="_blank">markus.metz.giswork@gmail.com</a><wbr>> wrote:<br>>> >><br>>> >><br>>> >><br>>> >> On Mon, Oct 9, 2017 at 4:37 PM, Stefan Blumentrath <<a href="mailto:Stefan.Blumentrath@nina.no" target="_blank">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" target="_blank">https://grass.osgeo.org/<wbr>grass72/manuals/r.import.html</a>)<wbr>.<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" target="_blank">grass-user-bounces@<wbr>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" target="_blank">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" target="_blank">http://land.copernicus.eu/pan-<wbr>european/corine-land-cover/<wbr>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_<wbr>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>>> >> > ______________________________<wbr>_________________<br>>> >> > grass-user mailing list<br>>> >> > <a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>>> >> > <a href="https://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/grass-user</a><br>>> >><br>>> ><br>><br>><br></div></div></div></div></div></div></div>
<br>______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/grass-user</a><br></blockquote></div><br></div>