On Tue, Oct 24, 2017 at 4:53 PM, Veronica Andreo <[email protected]> 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 <[email protected]>: >> >> 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 <[email protected]> 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 <[email protected]> wrote: >> >> >> >> Hi Hakim, >> >> >> >> On Mon, Oct 9, 2017 at 6:36 PM, Hakim Benoudjit <[email protected]> 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 < [email protected]> wrote: >> >> >> >> >> >> >> >> >> >> >> >> On Mon, Oct 9, 2017 at 4:37 PM, Stefan Blumentrath < [email protected]> 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:[email protected]] On Behalf Of Hakim Benoudjit >> >> >> > Sent: mandag 9. oktober 2017 16.02 >> >> >> > To: [email protected] >> >> >> > 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 >> >> >> > [email protected] >> >> >> > https://lists.osgeo.org/mailman/listinfo/grass-user >> >> >> >> >> > >> > >> > >> >> _______________________________________________ >> grass-user mailing list >> [email protected] >> https://lists.osgeo.org/mailman/listinfo/grass-user > >
_______________________________________________ grass-user mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/grass-user
