Re: [otb-users] bug in Connected Components Segmentation/Vectorization?

[Guillaume Pasero]

Ok, then maybe it has something to do with vector simplification, or a bad rounding issue. I'll look into that. Guillaume On 09/29/2016 01:55 PM, André Stumpf wrote: 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 otb-...@googlegroups.com To unsubscribe from this group, send email to otb-users+...@googlegroups.com 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 otb-users+...@googlegroups.com. For more options, visit https://groups.google.com/d/optout. -- 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 - guillaum...@c-s.fr -- -- 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 otb-users@googlegroups.com To unsubscribe from this group, send email to otb-users+unsubscr...@googlegroups.com 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 otb-users+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout. -- 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 - guillaume.pas...@c-s.fr -- -- 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 otb-users@googlegroups.com To unsubscribe from this group, send email to otb-users+unsubscr...@googlegroups.com 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 otb-users+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.