二章空间数据.pptVIP

(X0’,Y0’)为该图幅左下角点的数字化坐标， (X0,Y0)为该点相应的投影坐标 (X’,Y’)为图幅中任意一点的数字化坐标， (X,Y)为该点相应的投影坐标 因图幅底边与数字化坐标系间会存在一交角?，则对任意点的投影坐标有下式成立： X=X0+dX Y=Y0+dY 据坐标系旋转变换理论，并考虑到投影坐标轴的X、Y与数字化坐标轴X’与Y’相反，有： dX＝[(X’-X0’)sin? +(Y’-Y0’)cos?]×100/M dY＝[(X’-X0’)cos? -(Y’-Y0’)sin?]×100/M 因而有数字化坐标与投影坐标的换算公式为： 将以厘米为单位的数字化坐标值换算成以米为单位的实际投影坐标值 X＝X0+[(X’-X0’)sin? +(Y’-Y0’)cos?]×100/M Y＝Y0+[(X’-X0’)cos? -(Y’-Y0’)sin?]×100/M 3、空间坐标计算 2、已知投影坐标（X,Y）求数字化坐标(X’,Y’) 当任意点的投影坐标（X,Y）已知时，有 dX =（X-X0）×M/100 dY=（X-X0）×M/100 则数字化坐标的计算公式为： X’ ＝ X0’ + dXsin? + dYcos? Y’ ＝ Y0’ + dXYcos? - dYsin? 3、空间坐标计算 二、地理坐标（?，?）与投影坐标（X,Y）的换算 （?，?）是全球共用的空间物体位置描述的统一的地理坐标系，任何物体在（?，?）地理坐标系中位置是唯一确定的，因为这是一个绝对坐标系统。 （X,Y）是投影平面直角坐标系，由于地图投影函数及投影区域不同，使得对投影坐标系选择会产生一定的差异，因此它是一个相对坐标系。 1、解析方法 包含两种情况，即正解和反解。所谓反解是指若设投影计算公式为： X=f1(?，?) ?=F1(X,Y) Y=f2(?，?) ? =F2(X,Y) 即（X,Y） (?，?) 反解方法一般用于投影方程已知且反解容易的情况，例如对于前面所介绍的各类正轴圆锥投影的地理坐标就可以用反解方法进行解算。 3、空间坐标计算 2、数值方法 对于以数字形式存在于计算机内的地图,其上任意一点的直角坐标（投影坐标）都应是可求的。 如果经纬线网格为相互垂直的两组平行线，我们就可以采用构造空间曲面的方法解求?＝F1（X,Y）, ? ＝F2（X,Y） 空间曲面的构造可以采用曲面插值方法或者曲面逼近方法，插值方法要求构造曲面在已知点处取已知点的值，即： ?i ＝ F1（Xi,Yi） ?i ＝ F2（Xi,Yi） i=1,2,… 曲面逼近方法采用其他一些约束条件，如离差平方和最小 ?（ ?i － F1（Xi,Yi））2＝min ?（ ?i － F2（Xi,Yi））2 ＝min 3、空间坐标计算 三、基于曲面插值的坐标变换 由于资料类型、来源以及获取时间的差异使得不同资料间的坐标系各不相同，因此通过坐标变换以求得坐标系的统一是空间分析的前期准备工作之一，也是其重要的内容之一。 利用曲面插值方法，就可以通过计算机来解决任意两幅图之间的拼贴和坐标的转换与统一。 1、两个三角形之间的线性变换 ?ABC和?A’B’C’为两个任意三角形，其形状和大小可以毫不相关。设两个三角形顶点坐标分别为： A(X1,Y1) A’(X’1,Y’1) B(X2,Y2) B’(X’2,Y’2) C(X3,Y3) C’(X’3,Y’3) 则我们可以建立?ABC和?A’B’C’之间的线性变换F，保证?ABC中任意一点P（X,Y）都唯一地对应于?A’B’C’中的一点P’(X’,Y’),该线性变换表示如下： F: X’=f1（X,Y）=a0+a1X+a2Y Y’=f2（X,Y）=b0+b1X+b2Y A C B P A’ B’ C’ P’ 3、空间坐标计算 这就意味着三角形边上的点经F变换后其象与相应顶点之间的位置比例关系不变。这说明了?ABC变换到?A’B’C’，同时?ACD变换到?A’C’D’尽管是两个不同的变换，但其公共边上保持了变换的连续性，这正好是空间数据坐标变换中所必须的。由于?ABC与?A’B’C’之间任意变形，可以保证任意变形的地图资料都可以通过三角形剖分来进行坐标归算。 A D B C P A’ D’ C’ B’ P’ 点P（X,Y）位于三角形的边AC上，并有： P＝?A＋（1－