| 查看: 1132 | 回复: 7 | |||
[交流]
【求助】坐标从6°带转为3°带的坐标
|
目前要转一批6°分带的坐标数据为3°分带的坐标，不知道有没有什么公式，不涉及比例尺的改变，谢谢各位大侠~ |
回复此楼» 猜你喜欢
河南大学纳米科学与材料工程学院2025年科研人员招聘启事
已经有0人回复
-
江苏科技大学2025年招聘公告
已经有4人回复
-
地球化学论文润色/翻译怎么收费?
已经有55人回复
-
中国科学院华南植物园2025年各类人才招聘宣公告
已经有3人回复
-
ACS期刊投稿快要正式接受了，副主编提出引用已发表图片或数据引用的书面许可，求助！
已经有11人回复
-
香港中文大学（深圳)2025年招聘简章
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
经纬度数据如何转换为北京六分带投影坐标系
已经有3人回复- 求助Fe2O3原子坐标问题，谢谢~
已经有3人回复
- 3度带坐标转为1.5度带坐标怎么搞？
已经有4人回复
- 【分享】“北京54坐标系”转“西安80坐标系”的转换方法和步骤
已经有11人回复
- 【分享】[09-24]绿色精品软件下载[清清心整理]【已搜索无重复】
已经有3人回复
- 【求助】关于In2O3的氧空位坐标的问题
已经有9人回复
- 【转帖】GPS的坐标系统及坐标系转换-转贴
已经有6人回复
Ky_20060924
银虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 407.7
- 散金: 95
- 红花: 4
- 帖子: 325
- 在线: 63.2小时
- 虫号: 1094983
- 注册: 2010-09-10
- 专业: 微量元素地球化学
gisml(金币+3):多谢！我看看是否可以~ 2010-10-22 13:27:01
2楼2010-10-22 13:07:00
mjq240301
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 4276.6
- 红花: 2
- 帖子: 166
- 在线: 32.6小时
- 虫号: 540976
- 注册: 2008-04-07
- 专业: 工程地质学
3楼2010-10-22 14:50:33
Ky_20060924
银虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 407.7
- 散金: 95
- 红花: 4
- 帖子: 325
- 在线: 63.2小时
- 虫号: 1094983
- 注册: 2010-09-10
- 专业: 微量元素地球化学
4楼2010-10-22 15:35:24
mjq240301
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 4276.6
- 红花: 2
- 帖子: 166
- 在线: 32.6小时
- 虫号: 540976
- 注册: 2008-04-07
- 专业: 工程地质学
5楼2010-10-22 19:29:30
mjq240301
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 4276.6
- 红花: 2
- 帖子: 166
- 在线: 32.6小时
- 虫号: 540976
- 注册: 2008-04-07
- 专业: 工程地质学
6楼2010-10-22 19:32:45
Ky_20060924
银虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 407.7
- 散金: 95
- 红花: 4
- 帖子: 325
- 在线: 63.2小时
- 虫号: 1094983
- 注册: 2010-09-10
- 专业: 微量元素地球化学
|
这里有公式 http://htsoft.net.cn/html/y2009/ ... ormation-class.html 我记得当时大一的时候上测量课，所有的计算都是手算的，偶尔借助一下casio计算器！ |
7楼2010-10-23 07:29:26
Ky_20060924
银虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 407.7
- 散金: 95
- 红花: 4
- 帖子: 325
- 在线: 63.2小时
- 虫号: 1094983
- 注册: 2010-09-10
- 专业: 微量元素地球化学
gisml(金币+5):已经解决~ 十分感谢您的回复 2010-10-25 09:24:09
|
或者这里还有 3度带6度带 北京54坐标系3度带的坐标系，54北京坐标但是是6度带。 1∶2.5万及1∶5万的地形图采用6度分带投影，即经差为6度，从零度子午线开始，自西向东每个经差6度为一投影带，全球共分60个带，用1，2，3，4，5，……表示．即东经0～6度为第一带，其中央经线的经度为东经3度，东经6～12度为第二带，其中央经线的经度为9度……。 1∶1万的地形图采用3度分带，从东经1.5度的经线开始，每隔3度为一带，用1，2，3，……表示，全球共划分120个投影带，即东经1.5～ 4.5度为第1带，其中央经线的经度为东经3度，东经4.5～7.5度为第2带，其中央经线的经度为东经6度． 如何计算当地的中央子午线？ 当地中央子午线决定于当地的直角坐标系统，首先确定您的直角坐标系统是3度带还是6度带投影 ----------------------------------------------------------------------------------------------------- 公式推算: 6度带中央子午线计算公式：当地经度/6=N；中央子午线L=6 * N （带号） 当没有除尽，N有余数时， 中央子午线L=6*N - 3 ---------------------------------------------------------------------------------------------------- 3度带中央子午线计算公式： 当地经度/3=N；中央子午线L=3 X N -------------------------------------------------------------------------------------------------------- 我国的经度范围西起 73°东至135°，可分成 六度带十一个（13号带—23号带），各带中央经线依次为（75°、81°、……123°、129°、135°）； 三度带二十二 个（24号带—45号带）。各带中央经线依次为（72°、75°、……132°、135°）； 六度带可用于中小比例尺（如 1：250000）测图，三度带可用于大比例尺（如 1：10000）测图，城建坐标多采用三度带的高斯投影 --------------------------------------------------------------------------------------------------------- 如何判断投影坐标是3度带坐标还是6度带坐标 如(4231898,21655933)其中21即为带号，同样所定义的东伪偏移值也需要加上带号，如21带的东伪偏移值为21500000米。 假如你的工作区经度在120度至126度范围，则该坐标系为6度带坐标系，该带的中央经度为123度。 如(2949320,36353822)其中36即为带号，已知该地点位于贵阳市附近，而从地图上我们看到贵阳大概的经度是东经108度左右，因此可以36*3=108，所以该坐标系为3度带坐标系，该带的中央经度为108度。而不可能为6度带：36*6=216。 高斯投影是正形投影的一种，同一坐标系中的高斯投影换带计算公式是根据正形投影原理推导出的两个高斯坐标系间的显函数式。在同一大地坐标系中(例如1954北京坐标系或1980西安坐标系)，如果两个高斯坐标系只是主子午线的经度不同，那么显函数式前的系数可以根据坐标系使用的椭球元素和主子午线经度唯一确定。但如果两个高斯坐标系除了主子午线的经度不同以外，还存在其他线性系，则将线性变换公式代入换带计算的显函数式中，仍然可以得到严密的坐标变换公式。此时显函数式前的系数等价于使用两个坐标系主子午线的经度和线性变换参数联合求解得到的，可以唯一确定。 //6度带宽 54北京坐标系 //高斯投影由大地坐标(Unit:Metres)反算经纬度(Unit  D) void GaussProjInvCal(double X, double Y, double *longitude, double *latitude) { int ProjNo; int ZoneWide; ////带宽 double longitude1,latitude1, longitude0,latitude0, X0,Y0, xval,yval; double e1,e2,f,a, ee, NN, T,C, M, D,R,u,fai, iPI; iPI = 0.0174532925199433; ////3.1415926535898/180.0; a = 6378245.0; f = 1.0/298.3; //54年北京坐标系参数 ////a=6378140.0; f=1/298.257; //80年西安坐标系参数 ZoneWide = 6; ////6度带宽 ProjNo = (int)(X/1000000L) ; //查找带号 longitude0 = (ProjNo-1) * ZoneWide + ZoneWide / 2; longitude0 = longitude0 * iPI ; //中央经线 X0 = ProjNo*1000000L+500000L; Y0 = 0; xval = X-X0; yval = Y-Y0; //带内大地坐标 e2 = 2*f-f*f; e1 = (1.0-sqrt(1-e2))/(1.0+sqrt(1-e2)); ee = e2/(1-e2); M = yval; u = M/(a*(1-e2/4-3*e2*e2/64-5*e2*e2*e2/256)); fai = u+(3*e1/2-27*e1*e1*e1/32)*sin(2*u)+(21*e1*e1/16-55*e1*e1*e1*e1/32)*sin( 4*u)+(151*e1*e1*e1/96)*sin(6*u)+(1097*e1*e1*e1*e1/512)*sin(8*u); C = ee*cos(fai)*cos(fai); T = tan(fai)*tan(fai); NN = a/sqrt(1.0-e2*sin(fai)*sin(fai)); R = a*(1-e2)/sqrt((1-e2*sin(fai)*sin(fai))*(1-e2*sin(fai)*sin(fai))*(1-e2*sin (fai)*sin(fai))); D = xval/NN; //计算经度(Longitude) 纬度(Latitude) longitude1 = longitude0+(D-(1+2*T+C)*D*D*D/6+(5-2*C+28*T-3*C*C+8*ee+24*T*T)*D *D*D*D*D/120)/cos(fai); latitude1 = fai -(NN*tan(fai)/R)*(D*D/2-(5+3*T+10*C-4*C*C-9*ee)*D*D*D*D/24 +(61+90*T+298*C+45*T*T-256*ee-3*C*C)*D*D*D*D*D*D/720); //转换为度 DD *longitude = longitude1 / iPI; *latitude = latitude1 / iPI; } //高斯投影由经纬度(UnitD)正算平面坐标(含带号，Unit:Metres) void GaussProjCal(double longitude, double latitude, double *X, double *Y) { int ProjNo=0; int ZoneWide; ////带宽 double longitude1,latitude1, longitude0,latitude0, X0,Y0, xval,yval; double a,f, e2,ee, NN, T,C,A, M, iPI; iPI = 0.0174532925199433; ////3.1415926535898/180.0; ZoneWide = 6; ////6度带宽 a=6378245.0; f=1.0/298.3; //54年北京坐标系参数 ////a=6378140.0; f=1/298.257; //80年西安坐标系参数 ProjNo = (int)(longitude / ZoneWide) ; longitude0 = ProjNo * ZoneWide + ZoneWide / 2; longitude0 = longitude0 * iPI ; latitude0=0; longitude1 = longitude * iPI ; //经度转换为弧度 latitude1 = latitude * iPI ; //纬度转换为弧度 e2=2*f-f*f; ee=e2*(1.0-e2); NN=a/sqrt(1.0-e2*sin(latitude1)*sin(latitude1)); T=tan(latitude1)*tan(latitude1); C=ee*cos(latitude1)*cos(latitude1); A=(longitude1-longitude0)*cos(latitude1); M=a*((1-e2/4-3*e2*e2/64-5*e2*e2*e2/256)*latitude1-(3*e2/8+3*e2*e2/32+45*e2*e2 *e2/1024)*sin(2*latitude1)+(15*e2*e2/256+45*e2*e2*e2/1024)*sin(4*latitude1)-(35*e2*e2*e2/3072)*sin(6*latitude1)); xval = NN*(A+(1-T+C)*A*A*A/6+(5-18*T+T*T+72*C-58*ee)*A*A*A*A*A/120); yval = M+NN*tan(latitude1)*(A*A/2+(5-T+9*C+4*C*C)*A*A*A*A/24 +(61-58*T+T*T+600*C-330*ee)*A*A*A*A*A*A/720); X0 = 1000000L*(ProjNo+1)+500000L; Y0 = 0; xval = xval+X0; yval = yval+Y0; *X = xval; *Y = yval; } NN卯酉圈曲率半径，测量学里面用N表示 M为子午线弧长，测量学里用大X表示 fai为底点纬度，由子午弧长反算公式得到，测量学里用Bf表示 R为底点所对的曲率半径，测量学里用Nf表示。估计有了以上公式，不用软件，自己在EXCEL中也可以算了！！！ 总结，使用软件的前提要懂得基础理论，“殊途同归”，不过怎么样也能算出来！ |
8楼2010-10-23 07:32:39