GIS中栅格数据转换矢量数据算法

在上一篇文章GIS中将矢量数据转换删格数据算法中提到了矢量数据转换删格数据的方法；同样，在GIS中删格数据转换为矢量数据也是同样重要的操作，在本文中将介绍GIS系统中把删格数据转换为矢量数据的具体算法。下面是操作步骤：

1.二值化

由于扫描后的图像是以不同灰度级存储的，为了进行栅格数据矢量化的转换，需压缩为两级(0和1)，称为二值化。

2、二值图像的预处理

对于扫描输入的图幅，由于原稿不干净等原因，总是会出现一些飞白、污点、线划边缘凹凸不平等。

3.细化（剥皮法、骨架法）

所谓细化就是将二值图像象元阵列逐步剥除轮廓边缘的点，使之成为线划宽度只有一个象元的骨架图形。细化后的图形骨架既保留了原图形的绝大部分特征，又便于下一步的跟踪处理。

细化的基本过程是：

1. 确定需细化的象元集合；
2. 移去不是骨架的象元；
3. 重复，直到仅剩骨架象元。

如果是对扫描后的地图图像进行细化处理，应符合下列基本要求：

1. 保持原线划的连续性；
2. 线宽只为一个象元；
3. 细划后的骨架应是原线划的中心线；
4. 保持图形的原有特征。

4、追踪

细化后的二值图像形成了骨架图，追踪就是把骨架转换为矢量图形的坐标序列。其基本步骤为：

1. 从左向右，从上向下搜索线划起始点，并记下坐标。
2. 朝该点的8个方向追踪点，若没有，则本条线的追踪结束，转(1)进行下条线的追踪；否则记下坐标。
3. 把搜索点移到新取的点上，转（2）

注意的是，已追踪点应作标记，防止重复追踪。

5.拓扑化

为了进行拓扑化，需找出线的端点和结点，以及孤立点。

1. 孤立点：8邻城中没有为1的象元。如图(1)。
2. 端点：8邻城中只有一个为1的象元。如图(2)。
3. 结点：8邻城中有三个或三个以上为1的象元。如图（3）。