Al rasterizar una capa vectorial me gustaría conseguir una salida que contenga 1 para las áreas cubiertas por la forma y 0 para las áreas no cubiertas. Hasta ahora GDAL siempre convierte las áreas que no están cubiertas en nódulos. He intentado cambiar eso con el siguiente comando de GDAL pero no funciona correctamente. Supongo que quema nodata como cero pero también la metainformación de que cero = nodata por lo que no puedo utilizar la capa para el álgebra raster después (todas las zells con nodata se excluyen, sin embargo me gustaría que fueran cero en la operación de álgebra).
Este es mi comando gdal:
gdal_rasterize -l GMW_clipped_2010 -a pxlval -tr 0.000269959 0.000269959 -init -999999999.0 -a_nodata 0.0 -te -94.9017777777779 10.918977291452409 -78.1278774106337 23.2768888888896 -ot Byte -of GTiff /home/GMW_clipped_2010.shp /tmp/processing_46ba07b600a54c5581ad63a22017470a/b2dc5cbedbf746a2b9dbbdb406bd35e8/OUTPUT.tif