第4章地理信息系统数据库的组织.pptVIP

P1点的行列值： i=1+Integer(300-100)/1=201 j=1+Integer(51-0)/1=52 P2点的行列值： i=1+Integer(300-103)/1=198 j=1+Integer(53-0)/1=54 直线经过的行范围：198-201 P1 P2 199行中心处同直线相交的y值： y=300-199*1+1/2=101.5 该行y值对应的x值： x=(53-51)(101.5-100)/(103-100)+51=52 x坐标对应的栅格列值为： j=1+Integer(52-0)/1=53 P1 P2 200行中心处同直线相交的y值： y=300-200*1+1/2=100.5 该行y值对应的x值： x=(53-51)(100.5-100)/(103-100)+51=51.3 x坐标对应的栅格列值为： j=1+Integer(51.3-0)/1=52 P1、P2所连成的直线的栅格坐标为 （201，52），（200，52），（199，53），（198，54） （3）区域填充 区域矢量数据转成栅格数据是通过矢量边界轮廓的转换实现的。矢量边界线段转成栅格数据后，还要进行面域的填充。 1）射线法 ①判断疑问点p(x,y)是否在多边形内，从点向左引水平扫描线（射线）与区域边界相交次数为c，若c为奇数，在多边形内，c为偶数，则在多边形外。 ②作一系列水平扫描线，求出扫描线和区域边界交点，将交点按x值大小排序，其相邻坐标点之间的射线在区域内。 产生奇异性 I1 O X Y X3 X1 X2 X1 X2 X3 X4 X2 X1 I2 I3 O X Y P P1 ③奇异性：射线遇到极值点出现判断失误。 a.采用邻点分析法区分出极值点。对极值点看作2个同值交点，对非极值点看作一个斜交点，从而解决奇异性。 判断极值点：极值点为两直线交点，若两直线在扫描线的同一侧，则为极值点，否则为非极值点。 P1 P2 P3 P4 b.简化方法：它对组成多边形的每条直线的高端点y值坐标进行负修正，避开了奇异性。又称‘上闭下开’法。即在二直线的交点处扫描线上面的边与扫描线交点有效，扫描线下面的边与该扫描线交点无效，扫描线与多边形重合时不作求交运算。 Y O X I1 I1 相交零次 I2 I2 相交两次 I3 I3 相交两次 I4 I4 相交两次 I5 I5 相交四次 I6 I6 相交两次 2）边界点跟踪法 从边界上某一点栅格单元开始按顺时针方向跟踪边界上各栅格，对多边形中岛则按逆时针方向跟踪，将跟踪的每个栅格分别赋予R，L，N。 R--该栅格同相邻元素的行数不同，且行数增加的单元 L--该栅格同相邻元素的行数不同，且行数减少的单元 N--该栅格值单元同相邻单元行数相同 最后，逐行扫描，根据填充字符值，填充L-R之间的栅格。 N N N N L R L R L N N R L R L R L R L N N N N L R L R L R L R L N N R L N N R L R L R L R L R L R L R N N N N N N 3）边界代数法 前提：已知组成多边形边界（弧段）的拓扑关系，即边界前进方向的左右多边形号。 方法：假定沿边界前进方向y值下降时称下行，y值上升时称上行。上行时填充值为左多边形号减右多边形号，下行时为右多边形号减左多边形号，将每次填充值同该处的原始值作代数运算得到最终填充值。 N2 N4 N1 N3 1 2 3 N1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 N2 （1） N3 -1 -1 -1 2 2 2 2 1 1 2 2 2 2 1 1 2 2 2 2 1 1 2 2 2 1 1 1 N2 （2） N3 3 3 2 2 2 2 2 2 2 1 1 2 2 2