ERDAS西安80坐标系创建
1980西安坐标系
1980西安坐标--1980年国家大地坐标系我国最初覆盖全国的坐标系是1954年北京坐标系，采用了克拉索夫椭球元素（a=6378245ｍ，α=1/298.3）。
1954年北京坐标系的建立方法是，依照1953年我国东北边境内若干三角点与前苏联境内的大地控制网联接，将其坐标延伸到我国，并在北京市建立了名义上的坐标原点，并定名为1954年北京坐标系。
以后经分区域局部平差，扩展、加密而遍及全国。
因此，1954年北京坐标系，实际上是前苏联1942年坐标系，原点不在北京，而在前苏联的普尔科沃。
几十年来，我国按1954年北京坐标系建立了全国大地控制网，完成了覆盖全国的各种比例尺地形图，满足了经济、国防建设的需要。
由于各种原因，1954年北京坐标系存在如下主要缺点和问题：（1）克拉索夫斯基椭球体长半轴（a=6378245ｍ）比1975年国际大地测量与地球物理联合会推荐的更精确地球椭球长半轴（a=63781 40ｍ）大105m；（2）1954年北京坐标系所对应的参考椭球面与我国大地水准面存在着自西向东递增的系统性倾斜，高程异常（大地高与海拔高之差）最大为+65m（全国范围平均为29m），且出现在我国东部沿海经济发达地区。
（3）提供的大地点坐标，未经整体平差，是分级、分区域的局部平差结果。
使点位之间（特别是分别位于不同平差区域的点位）的兼容性较差，影响了坐标系本身的精度。
针对1954年北京坐标系的缺点和问题，1978年我国决定建立新的国家大地坐标系，该坐标系统取名为1980年国家大地坐标系。
大地坐标系原点设在处于我国中心位置的陕西省泾阳县永乐镇，它位于西安市西北方向约30km处，简称西安原点。
该坐标系的主要优点是：（1）地球椭球体元素，采用1975年国际大地测量与地球物理联合会推荐的更精确的参数，其中主要参数为：长半轴ａ=6378140m；短半轴ｂ=6356755.29；扁率α=1：298.257。
（2）椭球定位以我国范围高程异常值平方和最小为原则求解参数，椭球面与我国大地水准面获得了较好的吻合。
ERDAS西安80坐标系创建
ERDAS中设置西安80坐标系1.打开ERDAS目录一般情况下在C:\Leica Geosystems\Geospatial Imaging 9.2\etc中找到spheroid.tab文件选择记事本打开在记事本中输入以下红色内容"Xian_1980"{75 6378140 6356755.2882"D_xian_1980" 0 0 0 0 0 0 0 }保存并关闭记事本2.在ERDAS中配置投影参数打开ERDAS选择DataPrep选择Reproject Images（重投影影像）在Categories处选择自定义在Custom选项卡中，首先选择投影类型Gauss Kruger （高斯克鲁格投影）Spheriod（椭球）中即可选择刚才添加的“Xian_1980”椭球设置中央经线为东经99°，向东偏移33500000米设置完成后保存在Save as中输入Xian_1980_3_Degree_GK_Zone_33 In Category GCS_Xian_1980双击回车点击Cancel 退出自定义面板便可以在Categories中选择刚刚创建的“GCS_Xian_1980”系列，以及33度带的投影了。
3.影像的重投影选择需要重投影的源文件，并设置导出路径和文件名，点击左下角的OK即可对影像进行重投影图中为将名为“mz-pinjie.img”的影像重投影到“GCS_Xian_1980”系列下的“Xian_1980_3_Degree_GK_Zone_33”投影下。
附：（摘自网络）RDAS中定义西安80坐标系的方法ERDAS中添加椭球体和基准面：ERDAS包含了一个能够自定义椭球体、基准面、投影方式的扩展库，通过这个扩展库，可以在ERDAS中添加任何可能存在的投影系统。
基于以上2点，就可以将IAG75椭球的参数添加到ERDAS中，并且应用这个椭球对栅格数据进行投影变换。
Erdas实习报告
ERDAS遥感影像处理综合实习报告一：自定义坐标系（北京54、西安80、2000坐标系）操作步骤：1.1 添加椭球体修改文件为ellipse.txt，语法为<椭球体名称>,<长半轴>,<短半轴>。
这里的逗号为英文半角输入状态下的逗号，建议直接复制文件中已有的椭球体进行修改。
这里将下面三行加在ellipse.txt文件的末尾，保存关闭即可。
最终效果如下图所示。
• Krasovsky,6378245.0,6356863.0• IAG-75,6378140.0,6356755.3• CGCS2000,6378137.0,6356752.3注：ellipse.txt文件中已经有了克拉索夫斯基椭球，由于翻译原因，这里的英文名称是Krassovsky，为了让其他软件平台识别，这里新建一个Krasovsky椭球体。
1.2 添加基准面修改文件为datum.txt，语法为<基准面名称>,<椭球体名称>,<平移三参数>。
这里将下面三行添加在datum.txt文件末尾，保存关闭即可。
最终效果如下图所示。
• D_Beijing_1954, Krasovsky, -12, -113, -41• D_Xian_1980,IAG-75,0,0,0• D_China_2000,CGCS2000,0,0,01.3 定义坐标系操作步骤如下：（1）打开ENVI Classic，选择Map > Customize Map Projection工具；（2）在弹出的CustomizedMap ProjectionDefinition对话框内填写如图所示参数，其中Projection Name保持与ArcGIS中的名称一致；（3）选择Projection > Add New Project ion…，保存投影坐标系；（4）选择File > Save Projections…，在弹出对话框中点击OK，将新建坐标系保存在map_proj.txt文件内，以便下次启动ENVI后依然可以使用。
一种建立与1980西安坐标系极为接近的独立坐标系的方法
一种建立与1980西安坐标系极为接近的独立坐标系的方法摘要：利用AtuoCAD软件建立与1980西安坐标系极为接近的独立坐标系的方法，解决在小范围工程测量中，满足投影变形≤2.5cm/km的要求，保证在1980西安坐标系下初步概略设计及图文在独立坐标系下的有效性。
关键词：1980西安坐标系；独立坐标系；AtuoCADAbstract: A kind Method that establishing extremely independent coordinate system which clos e to the 1980 Xi’an coordinate system by AtuoCAD . on a small scale engineering survey, projection distortion ≤ 2.5cm/km meet the requirements to ensure that the 1980 Xi’an coordinate system schematic design and preliminary Text in the effectiveness of the independent coordinate system.Keywords: 1980 Xi’an coordinate system; independent coordinate system; AtuoCAD中图分类号: TD17引言在工程测量中，当作业区域的可能位于投影带的边缘，且测区高差起伏太大，平均高程面与国家大地基准参考椭球面有较大的差距时，必须要建立地方独立系以控制投影变形≤2.5cm/km[1] 。
有些工程施工方在1980西安坐标系下进行了工程的初步概略设计,为了确保原有初步概略设计成果在地方独立坐标系中的有效性,这就需要建立一套能够与1980西安坐标系极为接近的地方独立坐标系。
LPS操作流程
目录使用LPS处理无人机数据操作流程 (3)一、处理流程 (3)二、数据准备 (4)2.1相机参数 (4)2.2POS与相片数据 (5)2.3其他数据 (6)三、建立工程 (6)四、添加数据 (9)五、内定向与外方位元素导入 (10)5.1内定向 (10)5.2外方位元素导入 (12)六、自动生成同名点与添加控制点 (15)6.1自动生成同名点 (15)6.2添加控制点 (18)七、空中三角测量 (19)八、提取DEM (24)九、正射校正 (25)十、影像镶嵌 (27)附录1、无人机处理常见问题 (31)附录2、ERDAS IMAGINE中添加西安80坐标系 (33)使用LPS处理SPOT推扫式数据 (35)使用LPS处理无人机数据操作流程一、处理流程使用LPS处理无人机数据流程如下：二、数据准备2.1相机参数无人机搭载的相机一般为数码相机，在LPS中至少需要的相机参数有：焦距长、CCD尺寸。
同时还支持Australis校验参数：像主点偏移x0、像主点偏移y0、焦距长c、径向畸变系数k1、径向畸变系数k2、偏心畸变系数p1、偏心畸变系数p2、CCD非正方形比例系数b1、CCD非正交性畸变系数b2。
通常拿到的相机参数与下例类似：相机参数实例参数说明2141.8223 1421.2097 4677.4097 x y c0.000000004973527526 -0.000000000000000202 k1k20.000000042747151103 -0.000000007783852979 p1p20.000258582239 -0.000041822938 b1b2需将焦距长的单位换算成毫米，并计算像主点偏移，其他参数可直接使用(本例CCD尺寸为5.2μm)：焦距长=c×CCD=4677.4097×5.2=24.323mm像主点偏移x0=(像主点x-相片宽/2) ×CCD=(2141.8223-4272/2)×5.2=0.03027596mm像主点偏移y0=(像主点y -相片高/2) ×CCD=(1421.2097-2848)/2×5.2=-0.01450956mm2.2POS与相片数据目前无人机大多都搭载有GPS/IMU，可获取飞行POS数据，LPS可利用该数据对影像进行相对定向。
如何设置北京54或西安80坐标系
如何设置北京54或西安80坐标系？方法很简单，但是需要5个参数，这5个参数可以从当地测绘部门获取，也可以从设置好的GPS中读出。
第一步，在任意画面按MENU键，进入“系统设置”，然后进入”选择坐标系“，选”主位置画面“，在出现的列表中选择最下边的”自定义“或”User Grid“，选第一行的“横向墨卡托”，用箭头键和ENTER键输入当地的中央纬线（00.00000N）和中央经线(XXX.00000E)，进入下一页，然后用同样的方法设置投影比例（1.00000000）、尺度比（1.00000000）、东西偏差(00500000.0)和南北偏差(00000000.0)，然后按“完成”按钮。
中央纬线的计算方法用GOOGLE 查找吧，网上有很多。
第二步，在任意画面按MENU键，进入“系统设置”，然后进入”选择大地基准“，选”主位置画面“，在出现的列表中选择倒数第四个的”自定“或”User“，设置Delta A, Detla F, Delta X, Delta Y, Delta Z等5个参数，这5个参数需要找当地测绘部门查询。
如何计算当地的中央经线？子午线GPS要设置北京54或者西安80坐标系的时候，除了五个Delta 参数需要输入外，还有一个中央经线参数需要输入，与Delta参数不同，中央经线是可以自己计算的，计算方法如下：1、六度带中央经线经度的计算：当地中央经线经度＝６°×当地带号－３°，例如：地形图上的横坐标为２０３４５，其所处的六度带的中央经线经度为：６°×２０－３°＝１１７°（适用于１∶２．５万和１∶５万地形图）。
2、三度带中央经线经度的计算：中央经线经度＝３°×当地带号（适用于１∶１万地形图）。
附ARCGIS中定义的坐标系：北京54坐标系：Beijing 1954 3 Degree GK CM 75E.prj（三度分带法，中央经线东经75度，横坐标前不加带号）Beijing 1954 3 Degree GK Zone 25.prj（三度分带法，带号25，横坐标前加带号）Beijing 1954 GK Zone 13.prj（六度分带法，中央经线东经75度，横坐标前加带号）Beijing 1954 GK Zone 13N.prj（六度分带法，中央经线东经75度，横坐标前不加带号）西安80坐标系：Xian 1980 3 Degree GK CM 75E.prj（三度分带法，中央经线东经75度，横坐标前不加带号）Xian 1980 3 Degree GK Zone 25.prj（三度分带法，带号25，横坐标前加带号）Xian 1980 GK CM 75E.prj（六度分带法，中央经线东经75度，横坐标前不加带号）Xian 1980 GK Zone 13.prj（六度分带法，中央经线东经75度，横坐标前加带号）。
自定义坐标系(北京54、西安80、2000坐标系)
网址：
技术支持邮箱： ENVI-IDL@
电话： 010-57632288
技术支持热线： 400-819-2881-5
遥感事业部
1.2 国内坐标系介绍
先了解大地坐标的概念。大地坐标，在地面上建立一系列相连接的三角形，量取一段精确 的距离作为起算边，在这个边的两端点，采用天文观测的方法确定其点位（经度、纬度和方位 角） ，用精密测角仪器测定各三角形的角值，根据起算边的边长和点位，就可以推算出其他各点 的坐标。这样推算出的坐标，称为大地坐标。 我国 1954 年在北京设立了大地坐标原点， 由此计算出来的各大地控制点的坐标， 称为 1954 年北京坐标系。 为了适应大地测量的发展， 我国于 1978 年采用国际大地测量协会推荐的 IAG-75 地球椭球体建立了我国新的大地坐标系， 并在 1986 年宣布在陕西省泾阳县设立了新的大地坐标 原点，由此计算出来的各大地控制点坐标，称为 1980 年大地坐标系。随着社会的进步，国民经 济建设、国防建设和社会发展、科学研究等对国家大地坐标系提出了新的要求，迫切需要采用 原点位于地球质量中心的坐标系统作为国家大地坐标系。2000 国家大地坐标系（China Geodetic Coordinate System 2000，CGCS2000）是全球地心坐标系在我国的具体体现，其原点为包括海洋 和大气的整个地球的质量中心，CGCS2000 是我国当前最新的国家大地坐标系。 我们经常给影像投影时用到的北京 54、西安 80 和 2000 坐标系是投影直角坐标系，如下表 所示为国内坐标系采用的主要参数。从中可以看到我们通常称谓的北京 54 坐标系、西安 80 坐 标系实际上指的是我国的大地基准面。 表1.2 坐标名称 北京 54 西安 80 CGCS2000 投影类型 Gauss Kruger（Transverse Mercator） Gauss Kruger（Transverse Mercator） Gauss Kruger（Transverse Mercator） 北京 54、西安 80 和 2000 坐标系参数列表 椭球体 Krasovsky Xian_1980 CGCS2000 基准面 D_Beijing_1954 D_Xian_1980 D_China_2000
ERDAS IMAGINE遥感影像处理操作流程-2010全国土地利用变更调查监测与核查项目
北京望神州科技有限公司 北京海淀区上地三街9号嘉华大厦F座507室 邮编：100085 联系人：王鹤 Tel：(010)82781747/49‐107 手机：13439429782 邮箱：wanghe@
6
在状态栏可正要求 GCP 总体误差一般平坦区域要小于 1，山区小于 2。如 果误差较大，需要进行修改，删除点或增加新的控制点降低误差。 3、影像重采样 我们采集的点的数量、分布、精度都满足要求后，就可以进行影像的重采样 工作。在工具栏点击影像重采样 图标，打开影像重采样对话框。
RapidEye 原始数据 （5 波段）
每个波段的 NITF数据
DEM
控制资料
逐波段正射纠正 （选用NITF RPC模型）
波段合成
影像镶嵌
DOM成果
图1 使用RapidEye制作DOM的基本流程
北京望神州科技有限公司 北京海淀区上地三街9号嘉华大厦F座507室 邮编：100085 联系人：王鹤 Tel：(010)82781747/49‐107 手机：13439429782 邮箱：wanghe@
二 正射纠正 .......................................................................................................... 2 2.1 2.2 相关概念 ................................................................................................ 2 RapidEye 影像正射纠正 ...................................................................... 3 2.2.1 2.2.2 2.3 使用核心模块 .............................................................................. 3 使用 AutuSync ........................................................................... 9
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
ERDAS中设置西安80坐标系
1.打开ERDAS目录
一般情况下在C:\Leica Geosystems\Geospatial Imaging 9.2\etc中
找到spheroid.tab文件
选择记事本打开
在记事本中输入以下红色内容
"Xian_1980"{
75 6378140 6356755.2882
"D_xian_1980" 0 0 0 0 0 0 0 }
保存并关闭记事本
2.在ERDAS中配置投影参数
打开ERDAS
选择DataPrep
选择Reproject Images（重投影影像）
在Categories处选择自定义
在Custom选项卡中，首先选择投影类型
Gauss Kruger （高斯克鲁格投影）
Spheriod（椭球）中即可选择刚才添加的“Xian_1980”椭球
设置中央经线为东经99°，向东偏移33500000米
设置完成后保存
在Save as中输入Xian_1980_3_Degree_GK_Zone_33 In Category GCS_Xian_1980
双击回车
点击Cancel 退出自定义面板
便可以在Categories中选择刚刚创建的“GCS_Xian_1980”系列，以及33度带的投影了。
3.影像的重投影
选择需要重投影的源文件，并设置导出路径和文件名，点击左下角的OK即可对影像进行重投影
图中为将名为“mz-pinjie.img”的影像重投影到“GCS_Xian_1980”系列下的
“Xian_1980_3_Degree_GK_Zone_33”投影下。
附：（摘自网络）
RDAS中定义西安80坐标系的方法
ERDAS中添加椭球体和基准面：
ERDAS包含了一个能够自定义椭球体、基准面、投影方式的扩展库，通过这个扩展库，可以在ERDAS中添加任何可能存在的投影系统。
基于以上2点，就可以将IAG75椭球的参数添加到ERDAS中，并且应用这个椭球对栅格数据进行投影变换。
下面就来说明一下具体的添加过程。
1 在ERDAS安装目录下的etc/spheroid.tab文件是用来记载椭球体和基准面参数的。
它是一个TXT文本文件，可以用文本编辑器对它进行修改，只要依照它的语法就可以任意添加自定义的椭球体和基准面参数。
基本语法为：
“椭球名称”{
“椭球序号”椭球体长半轴椭球体短半轴
“椭球名称”0 0 0 0 0 0 0
“基准面名称1”dx1 dy1 dz1 rx1 rz1 ds1
“基准面名称2”dx2 dy2 dz2 rx2 rz1 ds2
……….
}
其中：“基准面名称”dx dy dz rx rz ds中，dx、dy、dz是x、y、z3个轴对于WGS84基准点的平移参数，单位为m。
rx、ry、rz是x、y、z、3个轴对于WGS84基准点的旋转参数，单位为rad。
Ds是对于WGS84基准点的比例因子。
在更多的情况下椭球的基准面是基于它本身的。
这时假定椭球的中心点是与没有经过任何平移或旋转的WGS84的基准面相重合，即这时椭球基准面的7个参数均为0，即这时椭球基准面的7个参数均为0。
我国在使用克拉索夫斯基椭球和IAG75椭球时就是用椭球体本身为基准。
在spheroid.tab文件末尾加入如下语句即可，假设spheroid.tab文件中最后一个椭球体序号为73（可以在文件最后一个椭球体中读出序号），则加入：
"IAG 75" {
74 6378140 6356755.2882
"xian 80" 0 0 0 0 0 0 0
}
经过以上的操作IAG75 椭球就会出现在ERDAS的椭球选择列表中。
2. 在Viewer中打开图像数据，Utility--->layer info，在projection info 栏中可以看到目前的数据投影信息还不完整。
点击edit菜单中的change map model，在弹出窗口中将unite参数设为meters，projection，参数设为Tansverse Mercator。
接下来再点击edit菜单中的Add/Change projection，在弹出对话框中将原始投影参数添加进去。
Custom
Projection Type:Transverse Mercator
Spheroid Name:IAG 75
Datum Name:xian80
Scale factor at central meridian:1.000000
Longitude of central meridian:117:00:00.000000000000 E
Latitude of origin of projection:0:00:00.000000000000 N
False easting:39500000.00000000000 meters False northing:0.00000000000000000 meters。