MAPGIS图幅拼接..doc

4.3图幅拼接 4.3.1问题和数据分析 1．问题提出 我国地理信息系统中都采用与基本比例尺地形图系列一致的地图投影系统，大于等于1：50万采用高斯--克吕格投影，1：100万采用正轴等角割圆锥投影，这两种投影均实现了将大区域空间进行分割，对于一些大比例尺专题图件还可选择矩形分幅，一幅图最多可以与邻近的8个图幅相接，如图4.3-1a、b所示。对于大幅面地图数字化，常常按矩形分块，相当于图4.3-1c。无论是那种情况，数字化时每幅图的坐标系均不一样，不能反映出图幅间位置邻接关系，若要拼接，必须建立统一的坐标系，将数字化数据通过平移、旋转操作转换到统一坐标系上拼接成一幅完整的大图。MapGIS 提供了平移、旋转、缩放等图形变换功能，可实现图幅坐标变换。 a高斯—克吕格投影 b正轴等角圆锥投影 c矩形分幅 图4.3-1 图幅邻接关系 2．数据准备 MapGIS 图形编辑系统提供了平移、旋转、缩放等图形变换功能，可实现图幅坐标变换。表4.3-1给出了相邻九幅1:5万标准图幅的左下角起始经度和起始纬度，通过利用MapGIS投影变换和图形编辑系统对下表九幅标准图框进行拼接。 表4.3-1 图幅起始经纬度值 图幅1 图幅2 图幅3 图幅4 图幅5 图幅6 图幅7 图幅8 图幅9 起始 经度 10630′ 106°30′ 106°30′ 106°45′ 106°45′ 106°45′ 107°00′ 107°00′ 107°00′ 起始 纬度 2930′ 29°40′ 29°50′ 29°30′ 29°40′ 29°50′ 29°30′ 29°40′ 29°50′ 4.3.2拼图基本原理 高斯投影和圆锥投影在小范围内，图框可看成梯形，现以九幅高斯投影图框为例，说明图幅拼接过程。高斯投影在小范围内，图框可看成梯形，如图4.3-2所示，九幅图可有两种简单拼接方法。 图4.3-2 梯形图幅邻接关系 1．（1）九幅图经平移、旋转操作，分别将左下角平移至坐标原点，上下图框平行X轴。 （2）Ⅰ固定，Ⅱ沿+X平移（l1-l2）/2 ，沿+Y平移h1，Ⅲ沿+X平移（l1-l3）/2，沿+Y平移h1+h2，三个文件加添到一起，则Ⅰ、Ⅱ、Ⅲ三幅图相拼得A ，同样Ⅳ、Ⅴ、Ⅵ相拼得B，Ⅶ、Ⅷ、Ⅸ 相拼得C。 （3）B固定，A 沿-X平移l1， 则A 右下角点同B的左下角点重合，如图4.3-3a所示，A 顺时针旋转θ角，将A文件添加到B中，A与B 相拼，如图4.3-3b所示。 a b c d e 图4.3-3 ABC拼接过程 （4）B固定，C逆时针旋转θ/角，如图4.3-3c所示，再沿+X平移l1，将C文件添加到A、B合并的文件中，A、B、C相拼，如图4.3-3d所示，。 （5）拼接后的图沿+X平移l1cosθ，沿-Y平移-l1sinθ，拼接完毕，如图4.3-3e所示。 2． （1）同 （2）同 （3）A固定，B逆时针旋转θ角，如图4.3-4a所示，再沿+X平移l1，将B文件添加到A中，B与A相拼，如图4.3-4b所示。 （4） C 逆时针旋转θ+θ/角，如图4.3-4c所示，再沿+X平移l1+ l1/cosθ，沿+Y平移l1/sinθ，将C文件添加到A、B合并的文件中，A、B、C相拼，如图4.3-4d所示。 （5）拼接后的图顺时针旋转α角，拼接完毕，如图4.3-4e所示。 a b c d e 图4.3-4 ABC拼接过程 4.3.2图幅拼接过程 1．系列标准图框生成 在投影变换子系统中，利用“系列标准图框”生成功能，如图4.3-5所示。按照表4.3-1给定的起始经、纬度坐标生成九幅1：5万的标准图框Frame1-Frame9。在图框参数输入，比例尺、接图表、外图廓线均不选择，如图4.3-6所示，因为图形的平移、旋转、缩放操作都是相对于原点的，所以生成的图框应将左下角平移为原点，同时要保证图框底边处于水平状态。 图4.3-5 1：5万图框 图4.3-6 1：5万图框参数 2．图幅拼接过程 （1）按照方法一，先完成Frame1、Frame2、Frame3的拼接，在图形编辑系统中打开Frame1，用鼠标读取Frame1图幅右上角坐标，并记录下来(484.58，369.60)，关闭Frame1。打开Frame2，利用“其它整图变换”功能，输入x、y平移量为