Trazado de contorno landsat con python

Question

Cómo hacerlo http://blog.remotesensing.io/2014/06/tracing-raster-outlines-with-gdal/ sin necesidad de guardar el shapefile. Me gustaría devolver sólo la geometría politonal.

gdal_trace_outline -ndv 0 -out-cs en -dp-toler 10 -ogr-out outline.shp

Answer 1

Desde el manual vinculado en la página a la que hace referencia, puede ver que hay varios métodos de salida diferentes:

Output:
  -report fn.ppm               Output graphical report of polygons found
  -mask-out fn.pbm             Output mask of bounding polygon in PBM format
  -out-cs [xy | en | ll]       Set coordinate system for following outputs
                               (pixel coords, easting/northing, or lon/lat)
                               Must be specified before -{wkt,wkb,ogr}-out options
  -llproj-toler val            Error tolerance for curved lines when
                               using '-out-cs ll' (in pixels, default is 1.0)
  -wkt-out fn.wkt              Output polygons in WKT format
  -wkb-out fn.wkb              Output polygons in WKB format
  -ogr-out fn.shp              Output polygons using an OGR format

La elección depende de cómo se planea interactuar con los resultados, pero WKT parece una solución sencilla. Sin embargo, el script parece destinado a crear archivos, no a informar a la línea de comandos.