ArcGIS教程：细化 (Spatial Analyst)

1、摘要

    通过减少表示要素宽度的像元数来对栅格化的线状要素进行细化。

2、

用法

·      细化工具的典型应用是对已扫描的高程等值线地图进行处理。由于扫描仪分辨率和原始地图中线宽度的原因，等值线将在生成的栅格中表示为五至十个像元宽度的线状元素。运行细化工具后，各个等值线将表示为单个像元宽度的线状要素。

·      过滤器选项使用与边界清理相同的过滤算法，来移除从主要分支延伸出来的较短线状要素。它也可移除宽度小于三个像元的要素。

·      指定输入线状要素的最大线宽对于细化栅格必不可少，线状要素的宽度可超出或低于默认的最大线宽值。如果最大线宽适合要被细化的最宽线状要素，则可望获得最佳结果。

·      以下详细说明了细化工具中使用的常规算法：

Zhan, Cixiang，1993，《混合线细化方法》，自动制图会议记录 11，明尼阿波里斯，第 396-405 页

语法

    Thin (in_raster, {background_value}, {filter}, {corners}, {maximum_thickness})

3、代码实例

    “细化”示例 1（Python 窗口）

    本例对背景值为 NoData 的栅格进行细化，并在尝试保留拐角和交汇点的同时对边界进行平滑处理。

import arcpy

from arcpy import env

from arcpy.sa import *

env.workspace = "C:/sapyexamples/data"

thinOut = Thin("land","NODATA", "FILTER", "SHARP", 300)

thinOut.save("c:/sapyexamples/output/thinout")

    “细化”示例 2（独立脚本）

    本例对背景值为 NoData 的栅格进行细化，并在尝试保留拐角和交汇点的同时对边界进行平滑处理。

# Name: Thin_Ex_02.py

# Description: Thins rasterized linear features by

#              reducing the number of cells

#              representing the width of the features.

# Requirements: Spatial Analyst Extension

# Import system modules

import arcpy

from arcpy import env

from arcpy.sa import *

# Set environment settings

env.workspace = "C:/sapyexamples/data"

# Set local variables

inRaster = "land"

tolerance = 300

# Check out the ArcGIS Spatial Analyst extension license

arcpy.CheckOutExtension("Spatial")

# Execute Thin

thinOut = Thin(inRaster, "NODATA", "FILTER", "SHARP", tolerance)

# Save the output

thinOut.save("c:/sapyexamples/output/thinoutput")

4、环境

    像元大小, 当前工作空间, 掩膜, 输出坐标系, 范围, 临时工作空间, 捕捉栅格