Nella documentazione per gdal_calc è indicato il calcolatore raster della riga di comando con sintassi numpy . Più avanti ci sono alcuni esempi in cui in uno di essi:
gdal_calc.py -A input.tif --outfile = result.tif --calc = "A * (A> 0)" --NoDataValue = 0 - significa impostare i valori di zero e sotto su null
Sfortunatamente non c'è esempio su operatori logici come:
--calc = "A * (A> 0 e A> B)" - significa mantenere A se A zero più grande e B più grande e impostare il resto su null
Sulla base delle funzioni logiche di Numpy / Scipy, mi aspetterei di scrivere operatori logici come:
--calc = "A * logical_and (A> 0, A> B)"
Ho provato questo e sembra funzionare, ma vorrei essere sicuro che sia corretto.
Allo stesso modo se si desidera un minimo di A e B:
--calc = "A * (A <= B) + B * (A> B)"
Puoi semplicemente scrivere:
--calc = "minima (A, B)"
Il mio problema è che non riesco a trovare alcun libro di cucina per assicurarmi di farlo bene. Esiste un buon libro di cucina con esempi avanzati di ciò che è e non è possibile con gdal_calc?