1、影像边界提取步骤说明1 二值化将影像转换为二值化的影像。1.1 处理波段从理论上来说,比之前处理 3 个波段至少要快 3 倍。单波段处理。1.2 分割值二值化分隔值为 254,修改为 250。1.3 分割左右两侧数据的处理之前只处理一侧是因为分割值为 254.修改为 250 之后,0-249 和 250-255 都是需要处理的。1.4 原始数据情况说明原始数据大小 780M.图幅格式 tiff。图幅个数 9。单幅影像大小约为 87M.1.5 处理工具和时间处理工具:二值化.fmw处理时间:共 2 minutes 49.0 seconds,每幅影像约 18seconds。2 转化为多边形将二值化
2、后的影像转换为多边形,这一步是处理过程中最耗时间的一步,也是最关键的一步。2.1 成果数据格式设置为 FGDB.2.2 成果数据格式的选择原方案中成果数据的格式为 shp 或者 PGDB。Shp 不支持嵌套查询,没有面积字段不利于进行下一步的操作。PGDB 数据要求比较严格,有的数据输出到 PGDB 会出错。这些数据如果导出 shp 格式之后再转换到 PGDB 也会出错。2.3 处理工具和时间处理工具: ARCTOOL 里面的 raster to polygon.处理时间:共 27 minutes,每幅影像约 3 minutes。3 从多边形中提取影像边界从影像背景和影像范围构成的多边形数据集
3、中提取影像范围多边形。3.1 成果数据格式PGDB 或 FGDB 均可。FGDB 比较安全不会出问题。转换失败之后一般重新导入即可。上一步的数据格式为 FGDB,转换到 PGDB 中可能会出现问题。例如:Warning empty output generated;System Error (0): Unexpected Error.Failed to execute (FeatureClassToFeatureClass_8).转换完成之后要进行检查。3.2 处理工具ARCTOOL 里面的 select.3.3 SQL 语句File geodatabase 和 personal geodat
4、abase 查询的语句是有区别的。PGDB 中的语句AREA =(select max(AREA) from a66 where GRIDCODE=0 )FGDB 中的语句“Shape_Area“ = (select MAX( “Shape_Area“ ) from ta18 where “GRIDCODE“ =0 )3.4 处理过程中的意外情况的处理有时发生意外的错误 重新处理即可。数据文件名太长,修改为简短的文件名。检查 ARCGIS 中是否打开正在处理的数据。3.5 处理之后的检查处理完成之后要检查是否有图形要素。4 多边形整饰填充影像范围多边形中的洞;修改多边形的边界。可以通过修改步骤
5、提高处理速度。填充影像范围多边形中的洞,之后提取最大的多边形,在对多边形的边界进行修整。4.1 成果数据格式PGDB 或 FGDB 均可。4.2 处理工具和时间处理工具:多边形整饰.fmw处理时间:8 minutes 16.0 seconds其中的 perserve shared boundaries 设置为 yes 之后处理的速度变慢。4.3 处理成果检查5 提取影像范围多边形从影像范围及其中的洞构成的多边形集合中提取影像范围。5.1 成果数据格式PGDB 或 FGDB 均可。5.2 处理 SQLShape_Area = (select max(Shape_Area) from a01 )这
6、里没有对 GRIDCODE 的选择。5.3 处理工具和时间ARCTOOL 里面的 select.6 修复多边形存在以下两种情况需要修复部分边界缺失边界不规整部分边界缺失步骤 时间二值化 2 minutes 49.0 seconds转为多边形 26 minutes 59 seconds提取影像边界 2 minutes 44 seconds 有时发生意外的错误 重新处理即可。处理完成之后要检查是否有图形要素。多边形整饰 8 minutes 16.0 seconds 检查处理成果二次提取影像边界 6.00 seconds 检查处理成果Shape_Area = (select max(Shape_Area) from ta18 )合计 40 minutes 16 seconds遗留问题处理之后的影像存在以下问题:暂无方法解决