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