国土资源数据2000国家大地坐标系统转换技术要求

附件国土资源数据2000国家大地坐标系转换技术要求国土资源部国家测绘地理信息局2017年2月目录一、坐标转换的数据内容...............................................................................................................2二、坐标转换基本要求...................................................................................................................2三、矢量数据的转换.......................................................................................................................3（一）转换工作流程...........................................................................................................4（二）转换方法..................................................................................................................41.管理单元（以县或者单图幅）转换方法............................................................52.空间数据库转换方法...........................................................................................6四、栅格数据转换...........................................................................................................................7（一）分幅转换流程...........................................................................................................7（二）分景数据转换流程...................................................................................................8（三）转换方法..................................................................................................................81.文件形式栅格数据转换方法................................................................................82. 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 分幅栅格数据转换方法................................................................................9五、相对独立的平面坐标系与2000国家大地坐标系建立联系的方法.....................................9（一）相对独立的平面坐标系统控制点建立联系的方法...............................................9（二）相对独立的平面坐标系统下空间图形转换.........................................................11附录A：点位坐标转换方法.........................................................................................................12附录B：坐标转换改正量计算.....................................................................................................16附录C：双线性内插方法.............................................................................................................18附录D：常用坐标转换模型.........................................................................................................19附录E：高斯投影正反算公式.....................................................................................................22附录F：子午线弧长和底点纬度 计算公式 六西格玛计算公式下载结构力学静力计算公式下载重复性计算公式下载六西格玛计算公式下载年假计算公式 .................................................................................23本技术要求规定了国土资源数据内容、转换基本要求、国土资源存量数据及增量数据由1980西安坐标系到2000国家大地坐标系的技术流程、转换方法及转换步骤，相对独立的平面坐标系与2000国家大地坐标系建立的联系方法等内容。一、坐标转换的数据内容全面梳理、合理评估国土资源各项调查、勘界、评价、资源管理等空间数据，根据实际需要，按照“应转尽转”的原则，转换为2000国家大地坐标系。国土资源数据应涵盖实际应用需要的各级各类国土资源空间数据，主要包括遥感影像、土地利用现状、土地利用总体规划、矿产资源总体规划、土地整治规划、农用地分等、基本农田、土地资源批、供、用、补、矿产资源勘查、开发、基础地质、区域地质、地球物理、地球化学等各级各类相关数据。二、坐标转换基本要求坐标转换应遵循以下基本要求：1.1:5万及以小比例尺数据库转换可利用国家测绘地理信息局提供的1:5万1980西安坐标系到2000国家大地坐标系图幅改正量，点位坐标按双线性内插方法（见附录C）进行逐点转换，点位数据及矢量数据也可利用两个坐标系下的重合点作为控制点计算转换参数，使用此参数实现数据转换（见附录A）。栅格数据按本要求中第四部分介绍的方法实施转换。2.1:1万比例尺空间数据的转换采用国家测绘地理信息局提供的1:1万比例尺1980西安坐标系到2000国家大地坐标系图幅改正量通过双线性内插的方法逐点计算改正量。也可采用按（2°×3°）进行分区，逐区计算转换参数，按点位转换方法进行转换。计算模型见附录A。3.1:1万以大比例尺下点位数据按点位坐标的转换方法逐点进行坐标转换和转换精度评定，见附录A精度评定（不包括建立相对独立的平面坐标系的各类数据对应的比例尺）。4.原1980西安坐标系下建立的相对独立的平面坐标系按与1980西安坐标系建立联系的方法建立与2000国家大地坐标系的联系，模型和方法见第五部分。5.1954北京坐标系下的数据，先转换为1980西安坐标系，再转换为2000国家大地坐标系。1954北京坐标系转换为1980西安坐标系的方法参照测绘部门发布的技术方法。三、矢量数据的转换矢量数据的转换，以1:1万比例尺数据为例。（一）转换工作流程根据1980西安坐标系向2000国家大地坐标系转换相应比例尺地形图坐标转换改正量，采用逐要素逐点转换法进行坐标转换或平移方法进行坐标转换，见图1。重新生成2000国家大地坐标系方里格2000国家大地坐标系国土资源数据（X,Y）修改图廓点及方里网标注坐标大地坐标转换为高斯平面坐标坐标转换高斯平面坐标转换到大地坐标2000国家大地坐标系数据库（B,L）1980西安坐标系国土资源数据库（B,L）图幅改正量图1矢量数据转换技术流程图（二）转换方法国土资源数据以空间数据库或管理单元（以县或者单图幅）存放。其存储方式不同，转换到2000国家大地坐标系下方法有所不同。1980西安坐标系国土资源数据（X,Y）1.管理单元（以县或者单图幅）转换方法（1）获取1980西安坐标系各要素的坐标，计算其2000国家大地坐标系下各要素的坐标；（2）将2000国家大地坐标系下各要素的坐标写回原要素；（3）添加2000国家大地坐标系下新的方里格网层及标注，删除原方里格网数据层及方里网标注坐标、图廓标注。该数据成果为2000国家大地坐标系，其有效图范围为原1980西安坐标系范围。计算2000国家大地坐标系坐标流程图见图2，其中1980西安坐标系高斯平面直角坐标转换流程见图(a)，1980西安坐标系大地坐标转换流程见图(b)。图2矢量数据转换流程图2.空间数据库转换方法连接后台数据库，加载空间数据库中每个要素类，读取各要素1980西安坐标系坐标，逐点计算2000国家大地坐标系下各要素坐标，将2000国家大地坐标系下的要素存储到空间数据库中新建的要素类，具体方法如下：（1）新建一个与原要素类结构相同的新要素类；（2）获取各要素1980西安坐标系坐标，逐点计算2000国家大地坐标系下各要素的坐标；读要素高斯平面直角坐标（X、Y）根据大地坐标（B、L）查找其所在格网四个角点向2000国家大地坐标系转换的坐标改正量按双线性内插方法计算该点向2000国家大地坐标系转换的坐标改正量db、dl将高斯平面直角坐标（X、Y）转换到大地坐标（B、L）2000国家大地坐标系下大地坐标B2000=B+dbL2000=L+dl将大地坐标B2000、L2000转换到高斯平面直角坐标X2000、Y2000根据大地坐标（B、L）查找其所在格网四个角点向2000国家大地坐标系转换的坐标改正量按双线性内插方法计算该点向2000国家大地坐标系转换的坐标改正量db、dl读要素大地坐标（B、L）2000国家大地坐标系下大地坐标B2000=B+dbL2000=L+dl（a）高斯平面坐标图（b）大地坐标图（3）将2000国家大地坐标系下要素写入新建要素类；（4）添加2000国家大地坐标系下新的方里格网层，删除原方里格网数据层；（5）更新相关字段属性值。四、栅格数据转换（一）分幅转换流程计算栅格中心点坐标改正量按照像素关系计算移动量（像素数）2000国家大地坐标系分幅栅格数据修改头文件坐标1980西安坐标系标准分幅栅格数据图幅改正量图3栅格数据转换流程图（二）分景数据转换流程图4分景数据转换流程图（三）转换方法栅格数据以文件形式和数据库形式存放，其存储方式不同，转换到2000国家大地坐标系下方法有所不同。1.文件形式栅格数据转换方法（1）获取1980西安坐标系栅格数据坐标范围，可依据区域大小、区域形状、精度要求高低等技术参数均匀内插部分控制点（5公里），计算其2000国家大地坐标系下的坐标；（2）用1980西安坐标系坐标及2000国家大地坐标系计算控制点2000国家大地坐标系坐标用1980西安坐标系坐标及2000国家大地坐标系坐标进行纠正、重采样2000国家大地坐标系分景数据1980西安坐标系分景数据1：10000改正量读1980西安坐标系坐标范围、均匀内插控制点下坐标，通过工具软件进行配准，完成数据的坐标转换；（3）更改数据头文件中定位坐标。2.标准分幅栅格数据转换方法对标准分幅栅格数据读取1980西安坐标系坐标头文件中定位坐标，计算数据中心点坐标“改正量”，按照像素关系计算移动量（像素数），避免图幅之间接边数据重新采样；（1）读取1980西安坐标系坐标头文件中定位坐标，计算数据中心点坐标“改正量”；（2）更改数据头文件中定位坐标。五、相对独立的平面坐标系与2000国家大地坐标系建立联系的方法（一）相对独立的平面坐标系统控制点建立联系的方法利用坐标转换方法将相对独立的平面坐标系统下控制点成果转换到2000国家大地坐标系下。（1）相对独立的平面坐标系与2000国家大地坐标系转换技术流程如图5。两坐标系重合点的 分析 定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析 与确定坐标转换方法与模型的确定独立坐标系成果转换到2000系分析独立坐标系建立方法收集与分析资料图5相对独立的平面坐标系转换技术流程（2）重合点选取原则择优选取地方控制网的起算点及高精度控制点、周围国家高精度的控制点，大中城市至少选取5个重合点（城外4个，市内中心1个）；小城市在城市外围至少选取4个重合点，重合点要分布均匀，包围城市区域，并在城市内部选定至少6个均匀分布的重合点对坐标转换精度进行检核。（3）转换模型确定建立相对独立的平面坐标系统与2000国家大地坐标系联系时，坐标转换模型要同时适用于地方控制点转换和城市数字地图的转换。一般采用平面四参数转换模型，重合点较多时可采用多元逐步回归模型。当相对独立的平面坐标系统控制点和数字地图均为三维地心坐标时，采用Bursa七参数转换模型。坐标转换中误差应小于0.05米。（二）相对独立的平面坐标系统下空间图形转换采用点对点转换法完成相对独立的平面坐标系统下空间数据数字地形图到2000国家大地坐标系的转换，转换后相邻图幅不存在接边问题。具体步骤如下：利用控制点的转换模型和参数，对相对独立的平面坐标系统下数字地形图进行转换，形成2000国家大地坐标系地形图。根据转换后的图幅四个图廓点在2000国家大地坐标系下的坐标，重新划分公里格网线，原公里格网线删除。根据2000国家大地坐标系下的图廓坐标，对每幅图进行裁剪和补充。附录A：点位坐标转换方法1.坐标转换流程图6点位坐标转换流程2.坐标转换步骤（1）重合点选取坐标重合点可采用在两个坐标系下均有坐标成果的点。但最终重合点还需根据所确定的转换参数，计算重合点坐标残差，根据其残差值的大小来确定，若残差大于3倍中误差则剔除，重新计算坐标转换参数，直到满足精度要求为止；用于计算转换参数的重合点数量与转换区域的大小有关，但不得少于5个。（2）模型选择全国及省级范围的坐标转换选择二维七参数转换模型；省级以下的坐标转换可选择平面四参数模型。对于相对独立的平面坐标系统与2000国家大地坐标系的联系可采用平面四参数模型或多项式回归模型。模型选取可参照下 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 1。表1点位坐标转换法模型选取表重合点所属坐标系坐标类型转换模型适用区域范围三维七参数二维七参数大地坐标椭球面多项式拟合全国和省级范围椭球面经纬差≥3区域布尔莎模型全国及省级范围莫洛金斯基模型省级以下范围空间直角坐标三维四参数＜2×2局部区域1980西安坐标系1954年北京坐标系平面坐标二维四参数局部小区域二维四参数局部小区域相对独立的平面坐标系平面坐标平面多项式拟合局部小区域（3）模型参数计算用所确定的重合点坐标，根据坐标转换模型利用最小二乘法计算模型参数。（4）精度检核选择部分重合点作为外部检核点，不参与转换参数计算，用转换参数计算这些点的转换坐标与已知坐标进行比较进行外部检核。应选定至少6个均匀分布的重合点对坐标转换精度进行检核。3.精度评定和评估方法分区转换及数据库转换点位的平均精度应小于图上的0.1mm。具体：1:5万空间数据库坐标转换精度≤5.0m；1:1万空间数据库坐标转换精度≤1.0m；1：2000空间数据库坐标转换精度≤0.2m依据计算坐标转换模型参数的重合点的残差中误差评估坐标转换精度。对于n个点，坐标转换精度估计公式如下：①V（残差）=重合点转换坐标-重合点已知坐标②空间直角坐标X残差中误差[]1XXvvMn=-③空间直角坐标Y残差中误差[]1YYvvMn=-④空间直角坐标Z残差中误差[]1ZZvvMn=-点位中误差222ZYXpMMMM++=⑤平面坐标x残差中误差[]1xxvvMn=-⑥平面坐标y残差中误差[]1yyvvMn=-⑦大地高H残差中误差[]1HHvvMn=-平面点位中误差为22yxpMMM+=附录B：坐标转换改正量计算1:1万格网点坐标转换改正量计算时一般按(2°×3°)进行分区，并对每个分区向外扩充约20′，分别解算出各分区的转换参数后，利用确定的转换方法与转换模型分别计算全国1:1万格网点的2000国家大地坐标系坐标B2000，L2000，进而求出各点的1980西安坐标系与2000国家大地坐标系的差值DB802000，DL802000（B2000-B80，L2000-L80），形成全国1:1万格网点的1980西安坐标系与2000国家大地坐标系的转换改正量DB802000，DL802000。大地坐标改正量计算公式：BBeBeWNaWeMBMzLBMyLBMxdBcossin])sin2(2[1cossinsincossin22222D-+D+D+D-D-=)cossin(cos1LyLxBNdLD-D-=式中:2,eaDD分别为IAG-75椭球与2000国家大地坐标系椭球长半径，第一偏心率平方之差。即Da=a2000-a80,De2=e20002-e802则各个点在2000国家大地坐标系中的大地坐标为：B2000=B80+dBL2000=L80+dL根据转换的B2000,L2000，采用高斯投影正算公式计算相应DX1=X2000-X80DY1=Y2000-Y80的高斯平面坐标X2000,Y2000。求取全国1:1万以大比例尺格网点的转换改正量全国1954年北京坐标系向2000国家大地坐标系转换改正量计算采用两步法：首先计算1954年北京坐标系向1980西安坐标系转换改正量，其次计算1980西安坐标系向2000国家大地坐标系转换改正量，最后将两改正量叠加形成1954年北京坐标系向2000国家大地坐标系转换改正量。附录C：双线性内插方法以图所示标记说明方法，已知待插点(x,y)周围4个邻点的Z值，设待插值Z(x,y)与它们在x、y两个方向上均线性相关，则在x方向(或y方向)内插两次，得到如图中插值Z(m,y)和Z(x,n)，再在y方向(或x方向)内插一次，得到(x,y)点的改正值Z(x,y)。该方法称为双线性内插。图7双线性内插示意图结合图，双线性内插公式:Z(x,m)=Z(i,m)+[Z(j,m)-Z(i,m)]×[(j-x)÷(j-i)]Z(x,n)=Z(i,n)+[Z(j,n)-Z(i,n)]×[(j-x)÷(j-i)]Z(x,y)=Z(x,n)+[Z(x,n)-Z(x,m)]×[(n-y)÷(n-m)]附录D：常用坐标转换模型1.二维七参数转换模型2222sincos""0coscossincossinsincos"""0cossin1sincos0sincos"00(2sin)sincos"1xyzLLXLNBNBYBBLBLBZMMMtgBLtgBLmNLLeBBMNeBeBBMaD-D=D+D--D-+--+--sincos"aBBffDD其中：,BLDD同一点位在两个坐标系下的纬度差、经度差，单位为弧度，,afDD椭球长半轴差（单位米）、扁率差（无量纲），,,XYZDDD平移参数，单位为米，,,xyz旋转参数，单位为弧度，m尺度参数（无量纲）。2.平面四参数转换模型属于两维坐标转换，对于三维坐标，需将坐标通过高斯投影变换得到平面坐标再计算转换参数。平面直角坐标转换模型：其中，x0，y0为平移参数，α为旋转参数，m为尺度参数。x2，y2为2000国家大地坐标系下的平面直角坐标，x1，y1为原坐标系下平面直角坐标。坐标单位为米。3.综合法坐标转换所谓综合法即就是在相似变换（Bursa七参数转换）的基础上，再对空间直角坐标残差进行多项式拟合，系统误差通过多项式系数得到消弱，使统一后的坐标系框架点坐标具有较好的一致性，从而提高坐标转换精度。综合法转换模型及转换方法：利用重合点先用相似变换转换Bursa七参数坐标转换模型式中，3个平移参数TZYXDDD，3个旋转参数000TSSXSSTSSYSSTSSZSSXXZYXXYYZXmYYZZYXZZD-=D+-++D-021021cossin(1)sincosxxxmyyy-=++TZZX和1个尺度参数m。对相似变换后的重合点残差ZYXVVV,,采用多项式拟合式中：B，L单位：弧度；K为拟合阶数；ija为系数，通过最小二乘求解。00KiijjXYZijSSijVVVaBL-===或或附录E：高斯投影正反算公式1.高斯投影正算公式+-+++-++=44442222422222cos)5861(7201cos)495(2415.0coslBttlBtlBNtx-++-++-+=4442224222222cos)5814185(1201cos)1(611coslBNtttlBtlBNy公式中子午线弧长X及公式中相关变量的计算见附录F。2.高斯投影反算公式+++-++--=44222222))(459061(3601))(935(1211)(2ffffffffffffNyttNyttNyyMtBB+++++++-=422242222))(8624285(1201))(21(611)(cosfffffffffffNytttNytNyBl式中f、ft分别为按fB值计算的相应量，fB及公式中相关变量的计算见附件F。附录F：子午线弧长和底点纬度计算公式1.常用量定义a为椭球长半轴，1980西安坐标系为6378140mb为椭球短半轴f为椭球扁率，1980西安坐标系为1/298.25721eab-=e─第一偏心率e─第二偏心率22abea-=222eff=-22'abeb-=abfa-=Be222cos'=tgBt=BeV22cos'1+=221+=VBeW22sin1-=B为纬度，单位弧度M─子午圈曲率半径233(1)aecMWV-==N─卯酉圈曲率半径acWV==2.子午线弧长X设有子午线上两点p1和p2，p1在赤道上p2的纬度为B，p1、p2间的子午线弧长X计算公式)12sin10sin8sin6sin4sin2sin)(1(2BGBFBEBDBCBBarcBAeaX+-+-+--=式中2acb=1212101210812108612108641210864212108642838860810015242880900913107206938388608990995242883465131072315104857655055262144103954096105307235838860814864856553610395163842205102410525615524288297297131072727654096220510245253215831048576693693655364365916384110252561756445431eGeeFeeeEeeeeDeeeeeCeeeeeeBeeeeeeA=+=++=+++=++++=+++++=++++++=3.底点纬度fB迭代公式)()(,)1(120iiiiBFBFXBBAeaXB-+=-=+直到iiBB--1小于某一个指定数值，即可停止迭代。式中]12sin10sin8sin6sin4sin2sin)[1()(2BGBFBEBDBCBBarcBAeaBF+-+-+--=]12cos1210cos108cos86cos64cos42cos2)[1()(2BGBFBEBDBCBBAeaBF+-+-+--=