实验八空间数据投影变换

实验八空间数据投影变换一、实验目的与要求

1．掌握图框生成和坐标换算方法。 2．掌握用户数据投影变换方法。 3．掌握图形投影变换方法。二、实验准备1．实验数据：本实验数据保存于文件夹Exercise-08中。2．预备知识：地图投影的基本原理和方法，特别应复习高斯—克吕格投影的相关知识。投影基础知识地图投影的基本问题：

地图投影是如何将地球椭球面表示在地图平面上，由于地球椭球面或圆球面是不可展开的曲面，即不可能展开成平面，而地图又必须是一个平面，所以需要将地球椭球面上的点投影到平面上，不同的投影方法，形成不同的投影坐标。投影基础知识——投影投影：

就是建立地球椭球面上点（λ，φ

）和平面上的点（ｘ，ｙ）之间的函数关系式的过程；投影的实质投影基础知识——投影变形将不可展的地球椭球面展开成平面，并且不能有断裂，则图形必将在某些地方被拉伸，某些地方被压缩，故投影变形是不可避免的。长度变形面积变形角度变形地图投影变形的图解示例

（摩尔维特投影－等积伪圆柱投影）地图投影变形的图解示例

（UTM－横轴等角割圆柱投影）投影基础知识地图投影的分类：

⑴、按地图投影的构成方法分：

a、几何投影：①、方位投影（见下一页示意图）②、圆柱投影（见下一页示意图）③、圆锥投影（见下一页示意图）b、非几何投影：用数学解析方法，求出投影公式，确定平面和球面之间点与点间的函数关系；（2）、按地图投影的变形性质分：

a、等角投影b、等积投影c、任意投影正、横、斜轴方位投影正、横、斜轴圆柱投影正、横、斜轴圆锥投影正轴投影经纬线形状投影基础知识——投影变换投影变换：就是将不同的地图投影函数关系式变换的过程；MAPGIS中的投影变换的定义：将当前地图投影坐标转换为另一种投影坐标，它包括坐标系的转换、不同投影系之间的变换以及同一投影系下不同坐标的变换等多种变换；地理坐标系用经（λ）纬度（φ）表示，记作M（λ，φ）投影基础知识——坐标系平面坐标系通过投影方式将椭球面投影到平面上。由曲变平平面直角坐标系采用直角坐标(笛卡尔坐标)来确定地面点的平面位置。可以通过投影将地理坐标转换成平面坐标。投影基础知识——坐标系投影基础知识——坐标系北京54坐标系：

解放后，为建立我国天文大地网，鉴于当时历史条件，在东北黑龙江边境上同苏联大地网联测，推算出其坐标作为我国天文大地网的起算数据；随后，通过锁网的大地坐标计算，推算出北京点的坐标，并定名为1954年北京坐标系。因此，1954年北京坐标系是苏联1942年坐标系的延伸，其原点不在北京，而在苏联普尔科沃。该坐标系采用克拉索夫斯基椭球作为参考椭球，高程系统采用正常高，以1956年黄海平均海水面为基准；

缺点：误差累计较大、参考椭球和国际不一致；西安80坐标系：

1978年4月召开的“全国天文大地网平差会议”上决定建立我国新的坐标系，称为1980年国家大地坐标系。其大地原点设在西安西北的永乐镇，简称西安原点。椭球参数选用1975年国际大地测量与地球物理联合会第16界大会的推荐值。简称IUUG-75地球椭球参数或IAG-75地球椭球；投影基础知识——高斯投影•由德国数学家高斯提出，后经克吕格扩充并推导出计算公式，故称为高斯-克吕格投影，简称高斯投影，为了控制变形，本投影采用分带的方法。6度分带从格林威治零度经线起，每6度分为一个投影带，全球共分为60个投影带；3度分带法从东经1度30分算起，每3度为一带。这样分带的方法在于使6度带的中央经线均为3度带的中央经线；我国1:2.5-1:50万地形图均采用6度分带；1:1万及更大比例尺地形图采用3度分带；

高斯投影采用分带投影。将椭球面按一定经差分带，分别进行投影。NSc中央子午线赤道高斯投影平面赤道中央子午线高斯—克吕格投影赤道中央子午线平行圈子午线Oxy高斯投影的特性（等角横切椭圆柱投影）中央子午线投影后为直线，且长度不变。除中央子午线外，其余子午线的投影均为凹向中央子午线的曲线，并以中央子午线为对称轴。投影后有长度变形。赤道线投影后为直线，但有长度变形。赤道中央子午线平行圈子午线Oxy高斯投影的特性除赤道外的其余纬线，投影后为凸向赤道的曲线，并以赤道为对称轴。经线与纬线投影后仍然保持正交，其他角度投影后也保持不变。所有长度变形的线段，其长度变形比均大于l。离中央子午线愈远，长度变形愈大。高斯投影带划分我国规定按经差6°和3°进行投影分带。6°带格林威治零度经线开始，按6°的经差自西向东分成60个带。3°带自1.5°开始，按3°的经差自西向东分成120个带。高斯投影——投影带的划分高斯投影——投影带的划分6°带与3°带中央经线之间的关系如图:这样分带的方法在于3°带中央经线在奇数带时与6°带中央经线重合；偶数为6°带分带经线高斯投影——投影带的划分6°带与3°带中央经线之间的关系如图:由于高斯-克吕格投影每一个投影带的坐标都是对本带坐标原点的相对值，所以各带的坐标完全相同，使用时只需变一个带号即可；高斯投影——投影带的划分6°带与3°带中央经线之间的关系如图:1:2.5-1:50万的地形图均采用6°分带1:1万及更大比例尺地形图采用3°分带MAPGIS投影变换为满足不同的投影方法实现不同图件的需要，在进行图形数据处理中很可能要从一个地图投影坐标系统转换到另一个投影坐标系统，MAPGIS投影变换是为实现这一功能服务的，该系统共提供了20种不同投影间的相互转换及经纬网生成功能。三、实验步骤与内容1．数据准备 将实验数据复制，粘贴至各自文件夹内。启动MAPGIS主程序。在主菜单界面中，点击参数按钮，在弹出的对话框中，设置工作目录最终指向Exercise-08（盘符依据各人具体情况设置）。三、实验步骤与内容2．标准图框生成 执行如下命令：实用服务投影变换系列标准图框生成1：10万图框（或者点击图框工具栏中的10按钮）。三、实验步骤与内容三、实验步骤与内容三、实验步骤与内容三、实验步骤与内容3．非标准图框生成 执行如下命令：系列标准图框生成1：1000图框（或者点击图框工具栏中的1000千按钮）。三、实验步骤与内容三、实验步骤与内容4．单点数据投影变换 执行如下命令：投影转换输入单点投影转换。三、实验步骤与内容三、实验步骤与内容三、实验步骤与内容三、实验步骤与内容三、实验步骤与内容依据上述方法，完成下表空白单元格坐标值的输入。原始数据投影结果说明103°25′00″E25°30′00″N转换为6度带平面直角坐标Y=17460kmX=2758km转换为地理坐标，直角平面坐标为6度带坐标Y=17460kmX=2758km转换为3度带平面直角坐标Y=35580kmX=3460km转换为地理坐标，直角平面坐标为3度带坐标高斯投影——投影带的计算6°带中央经线经度:L0=6°N-3°某点大地经度为L，该点6°带带号：N=[L/6](取整数)+13°带中央经线经度:L0=3°n式中n为3°带带号高斯投影我国幅员辽阔，南从北纬4°，北至北纬54°，西从东经75°，东至东经135°。中央经度从75°起共计11个六度带，21个3度带。高斯投影由于我国领土全部位于赤道以北，因此x值为正值，y值有正有负。Oxy500km为使y坐标为正，讲坐标轴西移500km，并在坐标前冠以代号。m点m的坐标：Xm=4346216.985mYm=19634527.165m其中19表示第19带xO三、实验步骤与内容5．空间数据投影变换 5.1．数据加载 执行如下命令：文件打开文件，加载投影变换所需数据。三、实验步骤与内容三、实验步骤与内容5.2．采集控制点 执行如下命令：投影变换MAPGIS文件投影选择点文件，选择控制点数据文件。三、实验步骤与内容三、实验步骤与内容在视图窗口点击鼠标右键，选择复位窗口命令。三、实验步骤与内容执行如下命令：投影变换当前文件TIC点输入TIC点，输入TIC点信息，在窗口中按顺时针或逆时针方向依次采集四个控制点。三、实验步骤与内容三、实验步骤与内容三、实验步骤与内容三、实验步骤与内容三、实验步骤与内容5.3．拷贝控制点信息 执行如下命令：投影变换当前文件TIC点浏览编辑TIC点，检查所采集的TIC点的信息是否有误。三、实验步骤与内容三、实验步骤与内容执行如下命令：投影变换文件间拷贝TIC点输入TIC点，将所采集的控制点信息拷贝给其他文件。三、实验步骤与内容三、实验步骤与内容 5.4．投影变换 执行如下命令：投影变换进行投影变换，完成图形的投影变换。三、实验步骤与内容三、实验步骤与内容三、实验步骤与内容三、实验步骤与内容三、实验步骤与内容三、实验步骤与内容执行如下命令：文件另存文件，将投影变换的结果文件换名保存。三、实验步骤与内容三、实验步骤与内容三、实验步骤与内容执行如下命令：投影变换文件间拷贝投影参数，将投影变换的结果文件换名保存。三、实验步骤与内容三、实验步骤与内容执行如下命令：投影变换进行投影变换，按上述Step7和Step11的步骤依次完成大竹蓬组.WT和大竹蓬组.WP的投影变换。投影结果文件分别保存为大竹蓬组-1.WT和大竹蓬组-1.WP。三、实验步骤与内容5.5．图框生成 按3所示步骤完成非标准矩形图框的生成，图框主要技术参数为：左下角坐标：X＝2537.7km，Y＝18346km右上角坐标：X＝2540.8km，Y＝18350.5km 方里网：实线公里坐标网，网间距为2km 图框名：演示用图框三、实验步骤与内容三、实验步骤与内容5.6．用户数据投影变换 执行如下命令：投影变换用户文件投影变换，将所提供的数据转换为图形