经纬度坐标与高斯坐标的转换代码.docx

} } / 功能说明： 将绝对高斯坐标 (y,x) 转换成绝对的地理坐标 (wd,jd) X_3 = x / - 3 ; + * X_3 - * pow(X_3,2) - * pow(X_3,3) + * pow(X_3,4) + * pow(X_3,5) - * pow(X_3,6) ; tf = tan(Bf0*PI/180); jd_hd = jd / * PI / ; // 将以秒为单位的经度转换成弧度 wd_hd = wd / * PI / ; // 将以秒为单位的纬度转换成弧度 // 如果不设中央经线(缺省参数 : -1000 )，则计算中央经线， // 否则，使用传入的中央经线，不再使用带号和带宽参数 //L = (DH - * DH_width ; // 计算中央经线的经度 if (LP == -1000) { L = (DH - * DH_width ; // 计算中央经线的经度 } else { L = LP ; } } l0 = jd / - L ; // 计算经差 tsin = sin(wd_hd); // 计算 sinB tcos = cos(wd_hd); // 计算 cosB // 计算克拉索夫斯基椭球中子午弧长 X X = / * wd - * tsin + * pow(tsin,3) + * pow(tsin,5) + * pow(tsin,7) ) * tcos; et2 = b_e2 * pow(tcos,2) ; // et2 = (e ** 2) * (cosB ** 2) N = b_c / sqrt( 1 + et2 ) ; // N = C / sqrt(1 + et2) t = tan(wd_hd); // t=tgB m = PI/180 * l0 * tcos; // m = cosB * PI/180 * l0 *x = X + N * t * ( * pow(m,2) + - pow(t,2) + * et2 + 4 * pow(et2,2)) * pow(m,4)/ + - * pow(t,2) + pow(t,4)) * pow(m,6) / ) ; *y = N * ( m + ( - pow(t,2) + et2 ) * pow(m,3) / + ( - * pow(t,2) + pow(t,4) + * et2 - * et2 * pow(t,2) ) * pow(m,5) / ); 3 度带与 6 度带 1．我国采用 6 度分带和 3 度分带： 1：万及1 : 5万的地形图采用6度分带投影，即经差为6度，从零度子午线开始，自西向东 每个经差6度为一投影带，全球共分 60个带，用1，2，3，4，5，……表示.即东经 0?6 度为第一带，其中央经线的经度为东经 3度，东经6?12度为第二带，其中央经线的经度为 9度。 1 : 1万的地形图采用3度分带，从东经度的经线开始，每隔 3度为一带，用1, 2, 3,…… 表示，全球共划分 120个投影带，即东经? 度为第 1带，其中央经线的经度为东经 3度， 东经?度为第 2带，其中央经线的经度为东经 6度.我省位于东经 113度－东经 120度之间， 跨第 38、39、40共计 3个带，其中东经度以西为第 38带，其中央经线为东经 114度；东经? 度为 39带，其中央经线为东经 117 度；东经度以东到山海关为 40带，其中央经线为东经 120 度。 地形图上公里网横坐标前 2位就是带号，例如： 1： 5万地形图上的横坐标为，其中 20即为 带号， 345486 为横坐标值。 当地中央经线经度的计算 六度带中央经线经度的计算： 当地中央经线经度=6°x当地带号一3° ,例如：地形图上的 横坐标为20345,其所处的六度带的中央经线经度为： 6°X 20-3°= 117° （适用于1 : 2.5 万和 1： 5万地形图）。 三度带中央经线经度的计算：中央经线经度= 3°^当地带号（适用于 1 : 1万地形图）。 3、 如何计算当地的中央子午线 当地中央子午线决定于当地的直角坐标系统，首先确定您的直角坐标系统是 3 度带还是 6 度带投影公式推算 : 6度带中央子午线计算公式：当地经度 /6=N ;中央子午线L=6 * N （带号） 中央子午线 L=6*N - 3当没有除尽， 中央子午线 L=6*N - 3 3 度带中央子午线计算公式： 当地经度 /3=N ；中央子午线 L=3 X N 我国的经度范围西起 73 °东至 135°，可分成 六度带^一个（13号带一23号带），各带中央经线依次为 （75°、81 °、……123°、129°、 135°）； 三度带二十二 个（24号带一45号带）。各带中央经线依次为（72°、7