Hi Guillaume,
Thanks for your reply. The binary raster I try to vectorize is derived from
a series of orthorectified images. So it is ground geometry (UTM).
I guess using a DEM in the process would not solve the problem in this
case, or?
Here the gdalinfo on the binary GeoTiff I try to vectorize:
gdalinfo VC_5_bin_close.tif
Driver: GTiff/GeoTIFF
Files: VC_5_bin_close.tif
VC_5_bin_close.tif.aux.xml
Size is 10980, 10980
Coordinate System is:
PROJCS["WGS 84 / UTM zone 31N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",3],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AXIS["Easting",EAST],
AXIS["Northing",NORTH],
AUTHORITY["EPSG","32631"]]
Origin = (699960.000000000000000,5000040.000000000000000)
Pixel Size = (10.000000000000000,-10.000000000000000)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 699960.000, 5000040.000) ( 5d32'33.47"E, 45d 7'31.89"N)
Lower Left ( 699960.000, 4890240.000) ( 5d29'59.46"E, 44d 8'16.83"N)
Upper Right ( 809760.000, 5000040.000) ( 6d56' 9.95"E, 45d 5' 9.45"N)
Lower Right ( 809760.000, 4890240.000) ( 6d52'11.95"E, 44d 5'59.18"N)
Center ( 754860.000, 4945140.000) ( 6d12'43.71"E, 44d36'51.95"N)
Band 1 Block=10980x1 Type=Byte, ColorInterp=Gray
Min=0.000 Max=1.000
Minimum=0.000, Maximum=1.000, Mean=0.341, StdDev=0.474
Metadata:
STATISTICS_MAXIMUM=1
STATISTICS_MEAN=0.34128088806407
STATISTICS_MINIMUM=0
STATISTICS_STDDEV=0.47413947684861
Cheers,
André
On Wednesday, September 28, 2016 at 12:34:55 PM UTC+2, Guillaume Pasero
wrote:
>
> Hi,
>
> It looks like errors when reprojecting the output vector data.
>
> The input image is in sensor geometry ?
>
> Have you tried to set an DEM (elev parameter) ?
> Regards,
> Guillaume
>
> On 09/19/2016 11:09 AM, André Stumpf wrote:
>
> Hello,
>
> I'm currently working on some simple Python wrappers around some OTB
> modules and got some strange results when trying to use
> the connected components application to vectorize a binary image.
>
> Using...
>
> subprocess.call([otbpath + "/otbcli_ConnectedComponentSegmentation",
> "-in", inputimage,
> "-expr", "distance<1",
> "-minsize", str(minsize),
> "-out", output_shapefile])
>
> ... I get the attached results with strong artefacts at what seems to be the
> borders between individually process strips. Furthermore, the application also
> seems to have the undocumented feature to abstract the polygons.
>
> As an alternative I tried to use otbcli_Segmentation which provides the
> expected result (see attached as well) but unfortunately takes > 6h on a
> 10000x10000 image.
>
> subprocess.call([otbpath + "/otbcli_Segmentation",
> "-in", inputimage,
> "-filter", "cc",
> "-filter.cc.expr", "distance<1",
> "-mode", "vector",
> "-mode.vector.minsize", str(minsize),
> "-mode.vector.out", output_shapefile,
> "-mode.vector.tilesize", str(4096)])
> Since such type of vectorization seems like a common task I wonder if anybody
> came across such issues before and has solved it?
>
> Many thanks in advance for any advice,
> André
>
> --
> --
> Check the OTB FAQ at
> http://www.orfeo-toolbox.org/FAQ.html
>
> You received this message because you are subscribed to the Google
> Groups "otb-users" group.
> To post to this group, send email to [email protected]
> <javascript:>
> To unsubscribe from this group, send email to
> [email protected] <javascript:>
> For more options, visit this group at
> http://groups.google.com/group/otb-users?hl=en
> ---
> You received this message because you are subscribed to the Google Groups
> "otb-users" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to [email protected] <javascript:>.
> For more options, visit https://groups.google.com/d/optout.
>
>
> --
> <http://www.c-s.fr> *Guillaume PASERO*
> Responsable technique
> *Business Unit E-SPACE & Geo Information - Département Image &
> Applications*
>
> *CS Systèmes d'Information*
> Parc de la Grande Plaine - 5, Rue Brindejonc des Moulinais - BP 15872
> 31506 Toulouse Cedex 05 - FRANCE
> +33 561 17 64 21 - [email protected] <javascript:>
>
--
--
Check the OTB FAQ at
http://www.orfeo-toolbox.org/FAQ.html
You received this message because you are subscribed to the Google
Groups "otb-users" group.
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/otb-users?hl=en
---
You received this message because you are subscribed to the Google Groups
"otb-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
