Use of CORINE land use over a large domain

[blheveder35]

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

[Lucen]

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