国家大地坐标系.ppt
国家大地坐标系 Still waters run deep.流静水深流静水深,人静心深人静心深 Where there is life,there is hope。有生命必有希望。有生命必有希望现有测绘成果转换到2000国家大地坐标系技术指南一、一、2000国家大地坐标系的定义国家大地坐标系的定义二、点位坐标转换方法二、点位坐标转换方法三、三、1:2.5-1:25万数据库的转换万数据库的转换四、四、1:1万及万及1:5千基础地理信息数据库的转换千基础地理信息数据库的转换五、相对独立的平面坐标系统与五、相对独立的平面坐标系统与2000国家大地坐国家大地坐标系建立联系的方法标系建立联系的方法六、坐标转换方法六、坐标转换方法七、附录七、附录一、一、2000国家大地坐标系的定义国家大地坐标系的定义国国家家大大地地坐坐标标系系的的定定义义包包括括坐坐标标系系的的原原点点、三三个个坐坐标标轴轴的的指指向向、尺尺度度以以及及地地球球椭椭球球的的4个个基基本参数的定义。本参数的定义。2000国国家家大大地地坐坐标标系系的的原原点点为为包包括括海海洋洋和和大气的整个地球的质量中心;大气的整个地球的质量中心;2000国国家家大大地地坐坐标标系系的的Z轴轴由由原原点点指指向向历历元元2000.0的的地地球球参参考考极极的的方方向向,该该历历元元的的指指向向由由国国际际时时间间局局给给定定的的历历元元为为1984.0的的初初始始指指向向推推算算,定定向向的的时时间间演演化化保保证证相相对对于于地地壳壳不不产产生生残余的全球旋转,残余的全球旋转,X轴轴由由原原点点指指向向格格林林尼尼治治参参考考子子午午线线与与地地球球赤赤道道面面(历历元元2000.0)的交点,)的交点,Y轴与轴与Z轴、轴、X轴构成右手正交坐标系。轴构成右手正交坐标系。v采用广义相对论意义下的尺度。采用广义相对论意义下的尺度。v2000国家大地坐标系采用的地球椭球参数的数值为:国家大地坐标系采用的地球椭球参数的数值为:v长半轴长半轴 a6378137mv扁率扁率 f=1/298.257222101 v地心引力常数地心引力常数 GM3.9860044181014m3s-2 v自转角速度自转角速度 7.292l1510-5rad s-1v其它参数见下表:其它参数见下表:短半径短半径b(m)6356752.31414极曲率半径极曲率半径c(m)6399593.62586第一偏心率第一偏心率e0.0818191910428第一偏心率平方第一偏心率平方e20.00669438002290第二偏心率第二偏心率e0.0820944381519第二偏心率平方第二偏心率平方e20.006739496775481/4子午圈的长度子午圈的长度Q(m)10001965.7293椭球平均半径椭球平均半径R1(m)6371008.77138相同表面积的球半径相同表面积的球半径R2(m)6371007.18092相同体积的球半径相同体积的球半径R3(m)6371000.78997椭球的正常位椭球的正常位U0(m2s-2)62636851.7149动力形状因子动力形状因子J20.001082629832258球谐系数球谐系数J4-0.00000237091126球谐系数球谐系数J60.00000000608347球谐系数球谐系数J8-0.000000000014270.00344978650678赤道正常重力值赤道正常重力值e(伽)伽)9.7803253361两极正常重力值两极正常重力值p(伽)(伽)9.8321849379正常重力平均值正常重力平均值(伽)(伽)9.7976432224纬度纬度45度的正常重力值度的正常重力值45(伽)(伽)9.8061977695采用采用2000国家大地坐标系后仍采用无潮汐系统。国家大地坐标系后仍采用无潮汐系统。返返回回主主菜菜单单二、点位坐标转换方法二、点位坐标转换方法v(一)模型选择(一)模型选择全国及省级范围的坐标转换选择二维七参数转换模型;省全国及省级范围的坐标转换选择二维七参数转换模型;省级以下的坐标转换可选择三维四参数模型或平面四参数模型。级以下的坐标转换可选择三维四参数模型或平面四参数模型。对于相对独立的平面坐标系统与对于相对独立的平面坐标系统与2000国家大地坐标系的联系可国家大地坐标系的联系可采用平面四参数模型或多项式回归模型。坐标转换模型详见本采用平面四参数模型或多项式回归模型。坐标转换模型详见本指南第六部分。指南第六部分。v(二)重合点选取(二)重合点选取坐标重合点可采用在两个坐标系下均有坐标成果的点。但坐标重合点可采用在两个坐标系下均有坐标成果的点。但最终重合点还需根据所确定的转换参数,计算重合点坐标残差,最终重合点还需根据所确定的转换参数,计算重合点坐标残差,根据其残差值的大小来确定,若残差大于根据其残差值的大小来确定,若残差大于3倍中误差则剔除,倍中误差则剔除,重新计算坐标转换参数,直到满足精度要求为止;用于计算转重新计算坐标转换参数,直到满足精度要求为止;用于计算转换参数的重合点数量与转换区域的大小有关,换参数的重合点数量与转换区域的大小有关,但不得少于但不得少于5个。个。v(三)模型参数计算(三)模型参数计算用所确定的重合点坐标,根据坐标转换模型利用最小二乘用所确定的重合点坐标,根据坐标转换模型利用最小二乘法计算模型参数。法计算模型参数。v(四)精度评估与检核(四)精度评估与检核用上述模型进行坐标转换时必须满足相应的精度指标,具用上述模型进行坐标转换时必须满足相应的精度指标,具体精度评估指标及评估方法见附件中相关内容。选择部分重合体精度评估指标及评估方法见附件中相关内容。选择部分重合点作为外部检核点,不参与转换参数计算,用转换参数计算这点作为外部检核点,不参与转换参数计算,用转换参数计算这些点的转换坐标与已知坐标进行比较进行外部检核。些点的转换坐标与已知坐标进行比较进行外部检核。应选定至应选定至少少6个均匀分布的重合点对坐标转换精度进行检核个均匀分布的重合点对坐标转换精度进行检核。v(五)数据库中点位坐标转换模型参数计算(五)数据库中点位坐标转换模型参数计算的区域选取的区域选取对于对于1980西安坐标系下的数据库,采用全国数西安坐标系下的数据库,采用全国数据计算的一套模型参数可满足据计算的一套模型参数可满足1:5万及万及1:25万万比例尺数据库转换的精度要求;采用全国数比例尺数据库转换的精度要求;采用全国数据计算的六个分区的模型参数可满足据计算的六个分区的模型参数可满足1:1万比万比例尺数据库转换的精度要求。对于例尺数据库转换的精度要求。对于1954年北年北京坐标系下的数据库的转换,采用全国数据京坐标系下的数据库的转换,采用全国数据计算的六个分区的模型参数可满足计算的六个分区的模型参数可满足1:5万及万及1:25万比例尺数据库转换的精度要求;按万比例尺数据库转换的精度要求;按(23)进行分区计算模型参数可满足)进行分区计算模型参数可满足1:1万比例尺数据库转换的精度要求。万比例尺数据库转换的精度要求。返返回回主主菜菜单单三、三、1:2.5-1:25万数据库的转换万数据库的转换v(一)按国家基本比例尺地形图分幅组织的数据库(一)按国家基本比例尺地形图分幅组织的数据库按国家基本比例尺地形图分幅组织的图形数据(按国家基本比例尺地形图分幅组织的图形数据(DLG、DEM、DRG),依据以下方案进行转换。),依据以下方案进行转换。v1、1:2.5-1:10万万DLG数据库转换数据库转换(1)1954年北京坐标系下年北京坐标系下1:2.5-1:10万万DLG数据库转换数据库转换a、依据相应比例尺分幅进行区域划分,分两步完成坐标转换。依据相应比例尺分幅进行区域划分,分两步完成坐标转换。首先进行椭球体变换,再利用对应的比例尺图幅区域的首先进行椭球体变换,再利用对应的比例尺图幅区域的X、Y坐标平移量进行坐标平移;坐标平移量进行坐标平移;b、依据依据2000国家大地坐标系下对应的比例尺标准分幅图廓进国家大地坐标系下对应的比例尺标准分幅图廓进行数据裁切,区域边缘图幅中的数据空白区利用相邻图幅数行数据裁切,区域边缘图幅中的数据空白区利用相邻图幅数据进行补充;据进行补充;c、添加添加2000国家大地坐标系下的方里格网层,删除原方里格国家大地坐标系下的方里格网层,删除原方里格网数据层;网数据层;d、完成图廓更改、数据编辑、数据接边、拓扑重建、数据入完成图廓更改、数据编辑、数据接边、拓扑重建、数据入库等数据后处理及建库工作;库等数据后处理及建库工作;b、依据依据2000国家大地坐标系下对应的比例尺标准分幅图廓进国家大地坐标系下对应的比例尺标准分幅图廓进行数据裁切,区域边缘图幅中的数据空白区利用相邻图幅数行数据裁切,区域边缘图幅中的数据空白区利用相邻图幅数据进行补充;据进行补充;c、添加添加2000国家大地坐标系下的方里格网层,删除原方里格国家大地坐标系下的方里格网层,删除原方里格网数据层;网数据层;d、完成图廓更改、数据编辑、数据接边、拓扑重建、数据入完成图廓更改、数据编辑、数据接边、拓扑重建、数据入库等数据后处理及建库工作;库等数据后处理及建库工作;e、图幅换带接边:采用右图(图幅换带接边:采用右图(1954年北京坐标系)年北京坐标系)接左图(接左图(2000国家大地坐标系)时,先进行右图的国家大地坐标系)时,先进行右图的椭球体与换带转换,在左带中利用左图的平移量进椭球体与换带转换,在左带中利用左图的平移量进行右图的坐标平移,完成接边后保存在左带中的右行右图的坐标平移,完成接边后保存在左带中的右图(备份)成果。返回右图取消先前换带接边加入图(备份)成果。返回右图取消先前换带接边加入的平移量,并进行投影变换,最后利用右带自身的的平移量,并进行投影变换,最后利用右带自身的平移量完成平移后,方可与其相邻的右图接边;平移量完成平移后,方可与其相邻的右图接边;f、对基础地理信息数据库元数据相关条目进行更改。对基础地理信息数据库元数据相关条目进行更改。(2)1980西安坐标系下西安坐标系下1:2.5-1:10万万DLG数据库转换数据库转换依据相应比例尺分幅进行区域划分,不考虑椭球体变依据相应比例尺分幅进行区域划分,不考虑椭球体变换,直接利用对应的比例尺图幅区域的换,直接利用对应的比例尺图幅区域的X、Y坐标平坐标平移量进行坐标平移;然后按照移量进行坐标平移;然后按照1954年北京坐标系下年北京坐标系下DLG数据库转换的数据库转换的bf对应步骤进行。对应步骤进行。v2、1:2.5-1:10万万DRG数据库转换数据库转换原数据为原数据为300500dpi的原版印刷地图经的原版印刷地图经扫描纠正生成的扫描纠正生成的RGB栅格数据,无图幅间要素栅格数据,无图幅间要素的接边处理。的接边处理。(1)1954年北京坐标系下年北京坐标系下1:2.5-1:10万万DRG数据库转换数据库转换a、考虑椭球变换及对应图廓角点的考虑椭球变换及对应图廓角点的X、Y坐标平移量,计算坐标平移量,计算1954年北京坐标系分幅图廓年北京坐标系分幅图廓角点在角点在2000国家大地坐标系下的坐标,并修改国家大地坐标系下的坐标,并修改数据头文件中相应的定位坐标;数据头文件中相应的定位坐标;b、在在DRG数据上叠加数据上叠加2000国家大地坐标系下新国家大地坐标系下新的大地控制基础层(图廓及方里格网等),新图廓中的大地控制基础层(图廓及方里格网等),新图廓中数据空白或数据出图区域不做图纹补充和裁减;数据空白或数据出图区域不做图纹补充和裁减;c、在图例中添加在图例中添加2000国家大地坐标系下新的控国家大地坐标系下新的控制基准说明条款;制基准说明条款;d、完成数据合层,并保持完成数据合层,并保持DRG数据的原有分辨数据的原有分辨率;率;e、更改元数据中相关内容,增加更改元数据中相关内容,增加1954年北京坐年北京坐标系标准分幅的图廓四角点在标系标准分幅的图廓四角点在2000国家大地坐标系下国家大地坐标系下坐标,计算坐标,计算2000国家大地坐标系标准分幅的图廓四角国家大地坐标系标准分幅的图廓四角点的坐标。点的坐标。转换后数据为转换后数据为2000国家大地坐标系坐标、国家大地坐标系坐标、1954年年北京坐标系分幅。北京坐标系分幅。v(2)1980西安坐标系下西安坐标系下1:2.5-1:10万万DRG数据库转换数据库转换v获取图幅对应比例尺图幅图廓角点的获取图幅对应比例尺图幅图廓角点的X、Y坐标平移量,根据平移量计算图幅定位坐标,坐标平移量,根据平移量计算图幅定位坐标,修改数据头文件;然后按照修改数据头文件;然后按照1954年北京坐标系年北京坐标系到到2000国家大地坐标系的国家大地坐标系的1:2.5-1:10万万DRG数数据库转换的据库转换的be步骤进行。步骤进行。v转换后数据为转换后数据为2000国家大地坐标系坐标、国家大地坐标系坐标、1980年西安坐标系分幅。年西安坐标系分幅。v3、1:2.5-1:10万万DEM数据库转换数据库转换原数据为原数据为25米分辨率的灰阶(米分辨率的灰阶(256个)栅个)栅格数据,建库数据图幅间接边处理完好。此数格数据,建库数据图幅间接边处理完好。此数据转换可有两种方式:一种是依据据转换可有两种方式:一种是依据2000国家大国家大地坐标系下地坐标系下DLG相关图层数据(等高线、高程相关图层数据(等高线、高程点)重新生成点)重新生成DEM(见(见DEM数据生产规范)数据生产规范),一种是进行一种是进行DEM数据的转换。以下给出数据的转换。以下给出DEM数数据转换方法。据转换方法。v(1)1954年北京坐标系下年北京坐标系下1:2.5-1:10万万DEM数数据库转换据库转换a、按照比例尺对应图幅分块,在需补充内容的按照比例尺对应图幅分块,在需补充内容的邻接边各增加一个相应比例尺图幅;邻接边各增加一个相应比例尺图幅;b、考虑椭球变换及相应的比例尺图幅的考虑椭球变换及相应的比例尺图幅的X、Y坐坐标平移量,求得标平移量,求得X、Y坐标改正值;坐标改正值;c、根据坐标改正值进行图幅坐标平移,同时,根据坐标改正值进行图幅坐标平移,同时,参考像素分辨率确定起算坐标进行数据重采样;参考像素分辨率确定起算坐标进行数据重采样;d、按按2000国家大地坐标系新的图廓及重叠像素国家大地坐标系新的图廓及重叠像素进行图幅裁切,更改数据头文件中定位坐标;进行图幅裁切,更改数据头文件中定位坐标;e、修改元数据相关条目。修改元数据相关条目。v(2)1980西安坐标系下西安坐标系下1:2.5-1:10万万DEM数据库转换数据库转换a、按照相应比例尺对应图幅分块,在需补、按照相应比例尺对应图幅分块,在需补充内容的邻接边各增加一个相应比例尺图幅;充内容的邻接边各增加一个相应比例尺图幅;b、依据相应的比例尺图幅的、依据相应的比例尺图幅的X、Y坐标平坐标平移量,进行图幅坐标平移,并参考像素分辨率移量,进行图幅坐标平移,并参考像素分辨率确定起算坐标完成数据重采样;确定起算坐标完成数据重采样;c、d按按1954年北京坐标系年北京坐标系1:2.5-1:10万万DEM数据库转换的数据库转换的d、e步骤进行。步骤进行。v4、1:25万万DLG数据库转换数据库转换(1)将)将1:25万分幅的平面坐标平移量转换为对万分幅的平面坐标平移量转换为对应的经、纬度平移量或直接获取对应图幅的经、纬度应的经、纬度平移量或直接获取对应图幅的经、纬度平移量;平移量;(2)根据)根据1:25万分幅的经、纬度平移量,完成万分幅的经、纬度平移量,完成1:25万经纬度数据到万经纬度数据到2000国家大地坐标系经纬度数据国家大地坐标系经纬度数据的转换(的转换(1954年北京坐标系需同时考虑椭球体变化和年北京坐标系需同时考虑椭球体变化和平移量);平移量);(3)依据)依据2000国家大地坐标系下对应的国家大地坐标系下对应的1:25万万标准分幅图廓进行数据裁切,区域边缘图幅中的数据标准分幅图廓进行数据裁切,区域边缘图幅中的数据空白区利用相邻图幅数据进行补充;空白区利用相邻图幅数据进行补充;(4)数据后处理,包括:图廓更改、新格网层)数据后处理,包括:图廓更改、新格网层添加、数据编辑、数据接边、拓扑重建、数据入库等;添加、数据编辑、数据接边、拓扑重建、数据入库等;(5)更改元数据文件)更改元数据文件 v5、1:25万万DEM数据库转换数据库转换(1)利用)利用2000国家大地坐标系对应的国家大地坐标系对应的DLG数据层,重新内插生成数据层,重新内插生成DEM;(2)依据新的)依据新的DEM更改元数据文件。更改元数据文件。(二)按其它方式建立的数据库(二)按其它方式建立的数据库v1、按区域建立的图形数据库、按区域建立的图形数据库v按区域(省、地区、流域等)建立的图形数据库按区域(省、地区、流域等)建立的图形数据库(DLG、DEM、DRG),可先分带分块分层完成转),可先分带分块分层完成转换,参照以上相应比例尺基础地理信息数据库的转换换,参照以上相应比例尺基础地理信息数据库的转换方案转换后拼接合成。方案转换后拼接合成。v1:10万万-1:25万数据库,依万数据库,依1:25万数据库转换方案万数据库转换方案逐块进行转换,再整体拼接合成;按非高斯投影方式逐块进行转换,再整体拼接合成;按非高斯投影方式组织的,将原数据经纬网组织的,将原数据经纬网3030或或1515交点作为交点作为坐标转换参考点,计算这些参考点在坐标转换参考点,计算这些参考点在2000国家大地坐国家大地坐标系下的坐标,利用地理信息软件进行图形纠正,完标系下的坐标,利用地理信息软件进行图形纠正,完成数据转换。成数据转换。v2、按线性条带建立的图形数据库、按线性条带建立的图形数据库v按线性条带(境界、河流、交通线、管道线等)建立的图形按线性条带(境界、河流、交通线、管道线等)建立的图形数据库,可依据条带的方向、长短等分段进行,再拼接合成;数据库,可依据条带的方向、长短等分段进行,再拼接合成;也可通过条带中一定密度地物点的两套坐标,通过软件逐点进也可通过条带中一定密度地物点的两套坐标,通过软件逐点进行纠正。具体方法:行纠正。具体方法:v分块纠正:对于分块纠正:对于1:1万分块,按万分块,按1:1万数据转换方案逐块纠正万数据转换方案逐块纠正后接边合成;对于后接边合成;对于1:5万分块,按万分块,按1:2.5-1:10万数据转换方案逐万数据转换方案逐块纠正后接边合成;块纠正后接边合成;v逐点纠正:依据数据精度,建立一定密度(逐点纠正:依据数据精度,建立一定密度(1:1万万100米格网米格网点、点、1:5万万2000米格网点)的坐标转换参考点,计算这些参考米格网点)的坐标转换参考点,计算这些参考点在新坐标系下的坐标,利用地理信息软件完成数据转换。点在新坐标系下的坐标,利用地理信息软件完成数据转换。v3、按无固定分幅分区建立的图形数据库、按无固定分幅分区建立的图形数据库v按无固定分幅分区建立的图形数据,根据坐按无固定分幅分区建立的图形数据,根据坐标系、比例尺及数据主体所在的图幅、数据的标系、比例尺及数据主体所在的图幅、数据的组织方式、产品类型(组织方式、产品类型(DLG、DEM、DRG)等,)等,参照相应比例尺的转换方案,实施数据转换。参照相应比例尺的转换方案,实施数据转换。v4、DOM数据库转换数据库转换v原数据为航空或航天遥感获取的黑白或彩色原数据为航空或航天遥感获取的黑白或彩色影像数据,是连续的灰度(全色)或影像数据,是连续的灰度(全色)或RGB(彩(彩色)栅格数据,分辨率有多种方式(主要包括色)栅格数据,分辨率有多种方式(主要包括用于用于1:5地形图测绘的各种分辨率航空影像,以地形图测绘的各种分辨率航空影像,以及用于专题调查的及用于专题调查的10米、米、15米、米、30米等卫星影米等卫星影像)。影像数据转换可参照下列方式进行。像)。影像数据转换可参照下列方式进行。v对于已按数据库组织方式加工与处理的对于已按数据库组织方式加工与处理的DOM数据,可采用数据,可采用1:2.5-1:10万万DEM的数据转换方的数据转换方法,也可采用计算各景影像有效图边的法,也可采用计算各景影像有效图边的4点在点在2000国家大地坐标系下的坐标来重新定位的方国家大地坐标系下的坐标来重新定位的方式。式。v对于尚未按数据库组织方式加工与处理的对于尚未按数据库组织方式加工与处理的DOM数据,数据,可采用可采用1:2.5-1:10万万DRG的数据转换方法,不再添加的数据转换方法,不再添加新的控制基础信息。新的控制基础信息。v分辨率分辨率5米米-30米的数据,需依据其数据主体所在的米的数据,需依据其数据主体所在的1:25万图幅区域来选用万图幅区域来选用1:25万对应图幅的综合坐标改万对应图幅的综合坐标改正值;对于分辨率在正值;对于分辨率在2米到米到5米间的数据,需依据其数米间的数据,需依据其数据主体所在的据主体所在的1:5万图幅区域来选用万图幅区域来选用1:5万对应图幅的万对应图幅的综合坐标改正值;由此确定各自的综合坐标改正值;由此确定各自的X、Y方向平移像素方向平移像素数对应的坐标值(直接取数对应的坐标值(直接取1:25万或万或1:5万综合坐标改正万综合坐标改正值,或由像素数值,或由像素数像素分辨率求得)。按高斯投影、像素分辨率求得)。按高斯投影、分像对(分景)组织的高分辨率影像数据,参照分像对(分景)组织的高分辨率影像数据,参照1:1万万DOM转换技术方案进行转换。转换技术方案进行转换。返返回回主主菜菜单单四、四、1:1万及万及1:5千基础地理信息数据库的转换千基础地理信息数据库的转换v(一)(一)1:1万及万及1:5千格网点坐标转换改正量计算千格网点坐标转换改正量计算v1、1980西安坐标系坐标转换改正量计算西安坐标系坐标转换改正量计算v1:1万以上大比例尺一般按万以上大比例尺一般按(23)进行分区,并进行分区,并对每个分区向外扩充约对每个分区向外扩充约20,分别解算出各分区的转换,分别解算出各分区的转换参数后,利用确定的转换方法与转换模型分别计算全参数后,利用确定的转换方法与转换模型分别计算全国国1:1万及万及1:5千格网点的千格网点的2000国家大地坐标系坐标国家大地坐标系坐标B2000,L2000,进而求出各点的,进而求出各点的1980西安坐标系与西安坐标系与2000国家大地坐标系的差值国家大地坐标系的差值DB802000,DL802000(B2000-B80,L2000-L80),形成全国),形成全国1:1万及万及1:5千格网点的千格网点的1980西安坐标系与西安坐标系与2000国家大地坐标系的转换改正量国家大地坐标系的转换改正量DB802000,DL802000。v2、1954年北京坐标系坐标转换改正量计算年北京坐标系坐标转换改正量计算v全国全国1954年北京坐标系向年北京坐标系向2000国家大地坐标系转换国家大地坐标系转换改正量计算采用两步法:首先计算改正量计算采用两步法:首先计算1954年北京坐标系年北京坐标系转换向转换向1980西安坐标系转换改正量,其次计算西安坐标系转换改正量,其次计算1980西西安坐标系向安坐标系向2000国家大地坐标系转换改正量,最后将国家大地坐标系转换改正量,最后将两改正量叠加形成两改正量叠加形成1954年北京坐标系向年北京坐标系向2000国家大地国家大地坐标系转换坐标转换改正量。坐标系转换坐标转换改正量。v 1954年北京坐标系向年北京坐标系向1980西安坐标系转换坐标西安坐标系转换坐标改正量计算改正量计算v新旧坐标系统(新旧坐标系统(1954年北京坐标系与年北京坐标系与1980年西安坐年西安坐标系)的转换标系)的转换 大地坐标改正量计算公式:式中式中:分别为分别为IAG-75椭球与克拉索夫斯基椭球椭球与克拉索夫斯基椭球长半径,第一偏心率平方之差。即长半径,第一偏心率平方之差。即则各个点在则各个点在1980西安坐标系中的大地坐标为西安坐标系中的大地坐标为:v根据转换的,采用高斯投影正算公式计算相应的高斯平面坐标。v求取全国1:1万以大比例尺格网点的转换改正量平差改正量的计算平差改正量的计算v1954年北京坐标系所提供的大地点成果没有经过整体平差,年北京坐标系所提供的大地点成果没有经过整体平差,1980西安坐标系提供的大地点成果是经过整体平差的数据,西安坐标系提供的大地点成果是经过整体平差的数据,所以新旧系统转换还要考虑平差改正量的问题。计算平差改所以新旧系统转换还要考虑平差改正量的问题。计算平差改正量比较麻烦,没有一定的数学模式,不同地区,平差改正正量比较麻烦,没有一定的数学模式,不同地区,平差改正量差别很大,在我国中部某些地区,平差改正量在量差别很大,在我国中部某些地区,平差改正量在1米以下,米以下,而在东北地区的某些图幅则在而在东北地区的某些图幅则在10米以上。在实际计算中,在米以上。在实际计算中,在全国均匀地选择一定数量的一、二等大地点,利用它们新全国均匀地选择一定数量的一、二等大地点,利用它们新(1980西安坐标系)旧(西安坐标系)旧(1954年北京坐标系)坐标系的坐年北京坐标系)坐标系的坐标进行多种分析试算并剔除粗差点,然后分别计算它们的坐标进行多种分析试算并剔除粗差点,然后分别计算它们的坐标差值,根据这些差值和它们的大地坐标分别绘制两张平差标差值,根据这些差值和它们的大地坐标分别绘制两张平差改正量分布图(即改正量分布图(即dX,dY分布图),这样在分布图上可以分布图),这样在分布图上可以直接内插出全国直接内插出全国1:1万以大比例尺格网点的平差改正量万以大比例尺格网点的平差改正量DX2,DY2。v根据全国根据全国1:1万以大比例尺格网点的转换改正万以大比例尺格网点的转换改正量量DX1,DY1和平差改正量和平差改正量DX2,DY2按下列按下列公式计算公式计算1954年北京坐标系向年北京坐标系向1980西安坐标西安坐标系转换坐标转换改正量系转换坐标转换改正量DX,DY。将将DX,DY换算成换算成1:1万以大比例尺格网点大地坐标万以大比例尺格网点大地坐标转换改正量转换改正量DB5480,DL5480。v1954年北京坐标系向年北京坐标系向2000国家大地坐标系转换坐国家大地坐标系转换坐标转换改正量计算标转换改正量计算将全国将全国1:1万以大比例尺格网点的万以大比例尺格网点的1954年北京坐标系年北京坐标系向向1980西安坐标系的转换改正量西安坐标系的转换改正量DB5480与与1980西安西安坐标系向坐标系向2000国家大地坐标系的转换改正量国家大地坐标系的转换改正量DB802000叠加,得到全国叠加,得到全国1:1万以大比例尺格网点万以大比例尺格网点1954年北京坐标系向年北京坐标系向2000国家大地坐标系转换的坐国家大地坐标系转换的坐标转换改正量标转换改正量DB542000。即:即:DB542000=DB5480+DB802000 DL542000=DL5480+DL802000(二)1:1万及1:5千DLG数据库转换转换流程如图1所示:图图1 1:1万及万及1:5千基础地理信息数据库转换技术流程千基础地理信息数据库转换技术流程v1、1954年北京坐标系下年北京坐标系下1:1万、万、1:5千千DLG数据库转换数据库转换(1)每个图幅的四个图廓点坐标改正量选用)每个图幅的四个图廓点坐标改正量选用1954年北京坐标年北京坐标系向系向2000国家大地坐标系转换方法计算;国家大地坐标系转换方法计算;(2)图幅内各要素点的坐标改正量根据选用的本图幅的四个)图幅内各要素点的坐标改正量根据选用的本图幅的四个图廓点坐标改正量,按双线性内插等方法计算;图廓点坐标改正量,按双线性内插等方法计算;(3)根据图幅四个图廓点坐标改正量和图幅内各要素点的坐)根据图幅四个图廓点坐标改正量和图幅内各要素点的坐标改正量,计算标改正量,计算2000国家大地坐标系下的图幅四个图廓点坐国家大地坐标系下的图幅四个图廓点坐标和图幅内各要素点的坐标;标和图幅内各要素点的坐标;(4)与周边图幅拼接;)与周边图幅拼接;(5)按照)按照2000国家大地坐标系下对应国家大地坐标系下对应1:1万、万、1:5千标准分幅千标准分幅计算新的公里格网数据,即添加计算新的公里格网数据,即添加2000国家大地坐标系下新的国家大地坐标系下新的公里格网层;公里格网层;(6)完成图廓更改、数据编辑、换带接边、拓扑重建;)完成图廓更改、数据编辑、换带接边、拓扑重建;(7)对空间数据库元数据相关条目进行更改;)对空间数据库元数据相关条目进行更改;(8)数据入库等数据后处理工作。)数据入库等数据后处理工作。v2、1980西安坐标系下西安坐标系下1:1万、万、1:5千千DLG数数据库转换据库转换(1)每个图幅的四个图廓点坐标改正量选用每个图幅的四个图廓点坐标改正量选用1980西安坐标系向西安坐标系向2000国家大地坐标系转换国家大地坐标系转换方法计算;方法计算;(2)()(8)参照)参照1954年北京坐标系到年北京坐标系到2000国国家大地坐标系家大地坐标系1:1万、万、1:5千千DLG数据库转换数据库转换的对应步骤进行。的对应步骤进行。(三)(三)1:1万及万及1:5千千DRG数据库转换数据库转换在保持原分辨率不变的情况下,利用逐格网纠正的方法进行数在保持原分辨率不变的情况下,利用逐格网纠正的方法进行数据转换。据转换。v1、1954年北京坐标系下年北京坐标系下1:1万、万、1:5千千DRG数据库转换数据库转换(1)在)在2000国家大地坐标系下生成图廓坐标及公里格网,逐国家大地坐标系下生成图廓坐标及公里格网,逐公里格网点纠正公里格网点纠正1:1万、万、1:5千千DRG数据;数据;(2)修改元数据相关条目;)修改元数据相关条目;(3)修改相关的图外整饰。)修改相关的图外整饰。v2、1980西安坐标系下西安坐标系下1:1万、万、1:5千千DRG数据库转换数据库转换v参照参照1954年北京坐标系到年北京坐标系到2000国家大地坐标系国家大地坐标系1:1万、万、1:5千千DRG数据库转换的对应步骤进行。数据库转换的对应步骤进行。(四)(四)1:1万及万及1:5千千DEM数据库转换数据库转换v利用DEM生产过程中形成的矢量数据与DEM离散点数据完成数据转换。v1、1954年北京坐标系下年北京坐标系下1:1万、万、1:5千千DEM数据库转换数据库转换v(1)矢量数据与DEM离散点数据转换方法;va、每个图幅的四个图廓点坐标改正量选用1954年北京坐标系向2000国家大地坐标系转换方法计算;vb、图幅内各要素点的坐标改正量根据选用的本图幅的四个图廓点坐标改正量,按双线性内插等方法计算;vc、根据图幅四个图廓点坐标改正量和图幅内各要素点的坐标改正量,计算2000国家大地坐标系下的图幅四个图廓点坐标和图幅内各要素点的坐标;vd、与周边图幅拼接。v(2)构造TIN;v(3)按相关规范或规定内插DEM;v(4)对检查点坐标文件进行点对点坐标转换;v(5)修改元数据条目。v2、1980西安坐标系下西安坐标系下1:1万、万、1:5千千 DEM数据库转数据库转换换v(1)每个图幅的四个图廓点坐标改正量选用1980西安坐标系向2000国家大地坐标系转换方法计算;v(2)(5)参照1954年北京坐标系1:1万、1:5千 DEM数据库转换的对应步骤进行。(五)(五)1:1万及万及1:5千千DOM数据库转换数据库转换v原数据为航空或航天遥感获取的黑白或彩色原数据为航空或航天遥感获取的黑白或彩色影像数据,是连续的灰度或影像数据,是连续的灰度或RGB栅格数据,栅格数据,在转换中应保持原影像分辨率。在转换中应保持原影像分辨率。1、在原、在原DOM上叠加相应坐标系统的内图廓及上叠加相应坐标系统的内图廓及公里格网,在公里格网,在2000国家大地坐标系下生成图国家大地坐标系下生成图廓坐标及公里格网,逐公里格网点纠正廓坐标及公里格网,逐公里格网点纠正1:1万、万、1:5千千DOM数据;数据;2、转换后,删除内图廓及公里格网;、转换后,删除内图廓及公里格网;3、修改元数据相关条目。、修改元数据相关条目。返返回回主主菜菜单单五、相对独立的平面坐标系统与五、相对独立的平面坐标系统与2000国家大地国家大地坐标系建立联系的方法坐标系建立联系的方法v(一)相对独立的平面坐标系统控制点建立联系的方法(一)相对独立的平面坐标系统控制点建立联系的方法v可通过现行国家大地坐标系的平面坐标过渡,利用坐标转换可通过现行国家大地坐标系的平面坐标过渡,利用坐标转换方法将相对独立的平面坐标系统下控制点成果转换到方法将相对独立的平面坐标系统下控制点成果转换到2000国国家大地坐标系下。家大地坐标系下。v选取相对独立的平面坐标系统与选取相对独立的平面坐标系统与2000国家大地坐标系的重合国家大地坐标系的重合点的原则如下:择优选取地方控制网的起算点及高精度控制点的原则如下:择优选取地方控制网的起算点及高精度控制点、周围国家高精度的控制点,大中城市至少选取点、周围国家高精度的控制点,大中城市至少选取5个重合个重合点(城外点(城外4个,市内中心个,市内中心1个);小城市在城市外围至少选取个);小城市在城市外围至少选取4个重合点,重合点要分布均匀,包围城市区域,并在城市个重合点,重合点要分布均匀,包围城市区域,并在城市内部选定至少内部选定至少6个均匀分布的重合点对坐标转换精度进行检个均匀分布的重合点对坐标转换精度进行检核。核。v建立相对独立的平面坐标系统与建立相对独立的平面坐标系统与2000国家大国家大地坐标系联系时,坐标转换模型要同时适用地坐标系联系时,坐标转换模型要同时适用于地方控制点转换和城市数字地图的转换。于地方控制点转换和城市数字地图的转换。一般采用平面四参数转换模型,重合点较多一般采用平面四参数转换模型,重合点较多时可采用多元逐步回归模型。当相对独立的时可采用多元逐步回归模型。当相对独立的平面坐标系统控制点和数字地图均为三维地平面坐标系统控制点和数字地图均为三维地心坐标时,采用心坐标时,采用Bursa七参数转换模型。坐七参数转换模型。坐标转换中误差应小于标转换中误差应小于0.05米米。(二)相对独立的平面坐标系统下数字地形图(二)相对独立的平面坐标系统下数字地形图转换转换v采用点对点转换法完成相对独立的平面坐标系统下采用点对点转换法完成相对独立的平面坐标系统下数字地形图到数字地形图到2000国家大地坐标系的转换,转换后相国家大地坐标系的转换,转换后相邻图幅不存在接边问题。具体步骤如下:邻图幅不存在接边问题。具体步骤如下:v利用控制点的转换模型和参数,对相对独立的平面利用控制点的转换模型和参数,对相对独立的平面坐标系统下数字地形图进行转换,形成坐标系统下数字地形图进行转换,形成2000国家大地国家大地坐标系地形图。坐标系地形图。v根据转换后的图幅四个图廓点在根据转换后的图幅四个图廓点在2000国家大地坐标国家大地坐标系下的坐标,重新划分公里格网线,原公里格网线删系下的坐标,重新划分公里格网线,原公里格网线删除。除。v根据根据2000国家大地坐标系下的图廓坐标,对每幅图国家大地坐标系下的图廓坐标,对每幅图进行裁剪和补充。进行裁剪和补充。返返回回主主菜菜单单六、坐标转换方法六、坐标转换方法v(一)坐标转换模型(一)坐标转换模型v1、二维七参数转换模型、二维七参数转换模型其中:其中:同同一一点点位位在在两两个个坐坐标标系系下下的的纬纬度度差差、经经度度差差,单位为弧度单位为弧度 椭球长半轴差(单位米)、扁率差(无量纲)椭球长半轴差(单位米)、扁率差(无量纲)平移参数,单位为米平移参数,单位为米 旋转参数,单位为弧度旋转参数,单位为弧度尺度参数(无量纲)尺度参数(无量纲)。2、平面四参数转换模型、平面四参数转换模型属于两维坐标转换,对于三维坐标,需将坐标通过高斯投属于两维坐标转换,对于三维坐标,需将坐标通过高斯投影变换得到平面坐标再计算转换参数。影变换得到平面坐标再计算转换参数。平面直角坐标转换模型:平面直角坐标转换模型:其中,其中,x0,y0为平移参数,为平移参数,为旋转参数,为旋