[QGIS-it-user] QGIS 3.4.4 Poligonizzazione - da raster a vettore

Salve a tutti,
l'algoritmo in oggetto crea un vettore di output con un campo integer,
questo esclude tutti i casi in cui il raster ha valori con virgola.

Secondo voi è un bug da segnalare??

saluti

-----
https://pigrecoinfinito.wordpress.com/
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html

Forse l'utilizzo dei soli valori interi è una cosa voluta poiché in genere la poligonalizzazione si esegue su raster classificati, quindi contenenti un numero limitato di classi identificate da valori interi. I valori decimali sono spesso presenti in raster con molta variabilità nei valori che solitamente non ha senso vettorializzare (si rischia di avere un poligono per ogni pixel).

Alessandro

Il 11/02/2019 17:30, Totò ha scritto:

Salve a tutti,
l'algoritmo in oggetto crea un vettore di output con un campo integer,
questo esclude tutti i casi in cui il raster ha valori con virgola.

Secondo voi è un bug da segnalare??

saluti

-----
https://pigrecoinfinito.wordpress.com/
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html
_______________________________________________
QGIS-it-user mailing list
QGIS-it-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/qgis-it-user

non vedo la limitazione… a quale algoritmo ti riferisci? quello di gdal che usa gdal_plygonize.py ?

···

Luigi Pirelli



ginetto wrote

non vedo la limitazione... a quale algoritmo ti riferisci? quello di gdal
che usa gdal_plygonize.py ?

nella maschera che si avvia da processing compare:
console GDAL/OGR
gdal_polygonize.bat GPKG:C:/......./dati.gpkg:DTM_Palermo
C:/......./OUTPUT.shp -b 1 -f "ESRI Shapefile" DN

-----
https://pigrecoinfinito.wordpress.com/
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html

verifica usando direttameente gdal_polygonize da linea di comando (osgeo shell) per vedere se e’ una limitazione dello script originale o qualcosa introdotto dal wrapper di qgis/processing

···

Luigi Pirelli



e’ un limite dello script… nel senso che usando gdal_polygonize.py (o bat) si usa:
https://www.gdal.org/gdal__alg_8h.html#a3f522a9035d3512b5d414fb4752671b1
mentre semplicemente usando https://www.gdal.org/gdal__alg_8h.html#a7a789015334d677afcbef67e5a6b4a7c si risolverebbe

una nota di lutra di vari anni fa

https://trac.osgeo.org/gdal/ticket/4516

morale… fai prima a fare una PR a gdal con gdal_fpolygonize.py o aggiungere una opzione di gda_plygonize.py che gestisca la chiamata giusta in funzione del sourcetype

···

Luigi Pirelli



cambiando qui la chiamata in funzione del data source

https://github.com/OSGeo/gdal/blob/master/gdal/swig/python/scripts/gdal_polygonize.py#L263

···

Luigi Pirelli



grazie Luigi per le tue documentate risposte.

ho scaricato il file .py allegato alla risposta di Lutra
https://trac.osgeo.org/gdal/ticket/4516

il file l ho sostituito a quello che si trova qui:
C:\OSGeo4W64\apps\Python37\Scripts
oppure qui: C:\OSGeo4W64\bin

ma il risultato è sempre lo stesso, il valore p integer.

per il resto no ho capito se è un bug da segnalare o no.

saluti

-----
https://pigrecoinfinito.wordpress.com/
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html