Reclassificar raster usando GDAL?

Estou procurando uma ferramenta do GDAL capaz de reclassificar uma varredura que contém vários valores discretos. Até agora, encontrei rasladores Reclassify usando GDAL e Python ; no entanto, parece que apenas valores únicos podem ser manipulados aqui.

Existe algo como uma ferramenta nativa do conjunto?

gdal_calc pode ser usado para uma reclassificação de muitas classes.

Por exemplo, você pode alterar os valores abaixo (e iguais) de 12 a 10, os valores de 20, 30, 40, 50 permanecem os mesmos e os valores entre acima de 50 e 62 são alterados para 60:

  python gdal_calc.py -A input.tif --outfile=output.file --calc="10*(A<=12)+20*(A==20)+30*(A==30)+40*(A==40)+50*(A==50)+60*((A>50)*(A<=62))" --NoDataValue=0

gdal_reclassify é uma ferramenta não oficial do Python, baseada em ligações do Python GDAL, capaz de reclassificar de acordo com várias classes de valores.

Requisitos:

python
numpy
gdal binaries
python-gdal bindings

Exemplo:

python gdal_reclassify.py source_dataset.tif destination_dataset.tif -c "<30, <50, <80, ==130, <210" -r "1, 2, 3, 4, 5" -d 0 -n true -p "COMPRESS=LZW"

Se você estiver trabalhando em um script python, use o .ReadAsArraymétodo Você pode reclassificar usando numpy.

import numpy as np
sample = np.random.randint(low = 0, high = 9, size =(5,5))
print(sample)
sample[sample == 4] = 40
sample[sample <= 2] = -20
print(sample)