Me gustaría convertir un raster a un vector de puntos utilizando grass:r.to.vect desde un script de pyqgis. Mi script no da ningún error ni salida. Ahora estoy atascado, lo que tengo que cambiar.
Aquí mi código:
input = raster.tif
output = points.shp
processing.runalg('grass:r.to.vect',{"input": input, "feature": 1,"GRASS_OUTPUT_TYPE_PARAMETER": 1, "GRASS_REGION_PARAMETER":"%f,%f,%f,%f"%(xmin, xmax, ymin, ymax), "output":output})
Hasta ahora lo he intentado:
Adelgazar la capa de entrada antes de convertirla en una trama. Ejecutar manualmente r.to.rast en QGIS sólo funciona bien en la versión adelgazada.
processing.runalg('grass:r.thin',{"input": input, "GRASS_REGION_PARAMETER":"%f,%f,%f,%f"%(xmin, xmax, ymin, ymax), "output":output})
Cargar la capa de entrada como una capa de QGIS:
fileInfo = QFileInfo(input)
baseName = fileInfo.baseName()
filePath = str(os.path.abspath(input))
rlayer = QgsRasterLayer(filePath, baseName)