Use of CORINE land use over a large domain

Dataset specific topics involving WPS.

Use of CORINE land use over a large domain

Postby blheveder35 » Fri May 25, 2018 6:51 am

Dear all,
I want to use SRTM-1Arc and 100m CORINE Landuse (mapped to 24 USGS categories) for a large domain (Europe).

A have begun with CORINE Landuse at 100m resolution and downloaded the CLC12 file g100_clc12_V18_5.tif.
I have followed the steps described in the different posts:
- reproject from "Lambert azimutal Equal Area" to WGS84
gdalwarp -s_srs etrs_1989_laea_l52_m10.prj -t_srs EPSG:4326 g100_clc12_V18_5.tif g100_clc12_WGS84.tif

- reclassify CORINE classes to USGS classes
python reclassify.py g100_clc12_WGS84.tif g100_clc12_WGS84_USGS.tif

- file format conversion : tif => binary for MET
gdal_translate -of ENVI -ot Int16 g100_clc12_WGS84_USGS.tif 00001-185328.00001-107317

I have two problems :
=> my python script end with an error

Traceback (most recent call last):
File "reclassify.py", line 73, in <module>
reclassify(sys.argv[1],sys.argv[2])
File "reclassify.py", line 13, in reclassify
data[i,j]=nclass(data[i,j])
File "reclassify.py", line 70, in nclass
return newclass[oldclass]
IndexError: list index out of range

I use the classical script available on different threads ...
I get an output file g100_clc12_WGS84_USGS.tif which looks OK ???

=> It ends with a very large binary file of 38Gb 00001-185328.00001-107317 and I am afraid that WPS is not able to work with a so large file.
Do you know how to cut the TIF file into tiles with gdal ?

Thanks a lot for the help

Blandine
blheveder35
 
Posts: 2
Joined: Tue Oct 18, 2016 12:39 pm

Re: Use of CORINE land use over a large domain

Postby Lucen » Thu Jun 21, 2018 12:22 am

Hi blheveder35,
I am playing with the similar problem recently, and I just found lots of helpful threads here, yet I didn't see one with the code you are using here. I am wondering would you please post here with some links of where these codes come from and we could see what these errors really are.

Thanks,
Lucen
Lucen
 
Posts: 1
Joined: Thu Jun 21, 2018 12:09 am


Return to Working with Various Datasets

Who is online

Users browsing this forum: No registered users and 11 guests