矢量-栅格转换.ppt

空间数据结构转换 一、由矢量向栅格数据格式转换 重点：确定栅格单元的大小。 已知： Xmax、Xmin和Ymax、Ymin；I，J。 求：栅格单元的边长ΔX和ΔY。 它们之间的关系为： ΔY=（Xmax-Xmin）/J ΔX=（Ymax-Ymin）/I X Y （0，0） J I xmin xmax ymin ymax ΔX ΔY 设矢量坐标点（x,y)，转换后的栅格单元行列值（I,J),则有： I=1+INT[（Ymax-Y）/ΔY] J=1+INT[（X-Xmin）/ΔX] 式中INT表示取整函数。栅格点的值用点的属性表示。 点的栅格化 线的栅格化可以分解成对组成折线的每个线段的栅格化。 设线段的两端点坐标分别为（x1,y1）和（x2,y2）。其中有y2y1。 根据点栅格化方法，求得两端点的行列号为（I1,J1)和（I2,J2）。△I = | I2 - I1| ，△J= | J2 - J1| ，分两种情况处理：△I △J 和 △J △I 。 当△I △J ，平行于x轴的做每一行的中心扫描线，再求每一条扫描线与线段的交点，按点的栅格化方法将交点转化为栅格坐标。 当△J △I ，平行于y轴的做每一列的中心扫描线，再求每一条扫描线与线段的交点，按点的栅格化犯法将交点转化为栅格坐标。 线的栅格化（扫描线法） 如果△I △J ，扫描线与线段的交点的坐标为 Yi=Ymax-ΔY·(I-1/2) Xi=[(X2-X1)/(Y2-Y1)](Yi-Y1)+X1? 如果△J △I，扫描线与线段的交点的坐标为 Xi=Xmax-ΔX·(J-1/2) Xi=[(Y2-Y1)/(X2-X1)](Xi-X1)+Y1? 线的栅格化 例：设两个端点的行、列号已经求出，其行号为3和7，则中间网格的行号必为4、5、6。 线的栅格化 多边形（面域）栅格化 （1）射线算法 由待判点向图外某点引射线，判断该射线与某多边形所有边界相交的总次数，如果相交偶数次，则待判点在该多边形外部，如为奇数次，则待判点在该多边形内部，如下图所示。 二值化 将原来栅格单元的灰度值G（i，j）（256级）压缩为两个灰阶，即0和1两级。 方法：① 定义阈值T， RminTRmax; ②二值化，if R=T then R=1; else R=0; 二、由栅格向矢量数据格式转换 细化 细化的目的是消除线划横断面栅格数的差异，使得每一条线只保留代表其轴线或周围轮廓线（对面状符号而言）位置的单个栅格的宽度。 细化 1）N(p)为p的邻点的数值的和； 2）图像像素联接数T(p)，如果旋转着看像素周围的点，T(p)就是p周围8个点从0变成1的次数，它反映了像素邻点的联接的块数。 3）pW，pE，pS，pN分别指像素左侧、右侧、下边、上边邻点的数值。 像素周围的8个直接相邻像素 像素连接数 细化 算法步骤如下[Zhang-Suen，1984]： 1：对于栅格图像中的每个点p，进行如下操作：如果2?N(p)?6并且T(p)=1并且pNpSpE=0并且pWpEpS=0，则标志p点； 2：将所有被标志的栅格点赋值为0，如果没有被标志的点，则算法结束； 3：对于栅格图像中的每个点p，进行如下操作：如果2?N(p)?6并且T(p)=1并且pNpSpW=0并且pWpEpN=0，则标志p点； 4：将所有被标志的栅格点赋值为0，如果没有被标志的点，则算法结束； 5：转到第一步。 边界线跟踪 拓扑关系生成 去除多余点及曲线圆滑 (a) (b) (c) (d)