空间直角坐标直接求解大地纬度_第1页
空间直角坐标直接求解大地纬度_第2页
空间直角坐标直接求解大地纬度_第3页
空间直角坐标直接求解大地纬度_第4页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、空间直角坐标直接求解大地纬度1 计算公式依次使用下面的公式进行复数运算：2 VC+代码#include <math.h>#include <complex>/*对一个复数 x 开 n 次方*/std:complex<double> sqrtn(const std:complex<double>&x,double n) double r = _hypot(x.real(),x.imag(); /模 double a = atan2(x.imag(),x.real(); /辐角 n = 1.0 / n; r = pow(r,n); a *=

2、n; return std:complex<double>(r * cos(a),r * sin(a);/*空间直角坐标变换为大地坐标a in 椭球长半轴f in 扁率倒数，不能为零X,Y,Z in 空间直角坐标B,L,H in 大地纬度，大地经度，大地高。纬度和经度的单位为弧度*/void XYZtoBLH(double a,double f,double X,double Y,double Z,double&B,double&L,double&H) doubleb=a * (1.0 - 1.0 / f); doublex=_hypot(X,Y); dou

3、bley=Z; L=atan2(Y,X); f=2.0 / (a + b) * (a - b); doubleax=a * x * f; doubleby=b * y * f; std:complex<double>b4(-ax,by); std:complex<double>d4( ax,by); std:complex<double>P= (b4 * d4 + 4.0) / -3.0; std:complex<double>Q= (d4 + b4) * (d4 - b4) / 2.0; std:complex<double>D=s

4、qrtn(Q * Q - P * P * P,2.0); std:complex<double>M(-0.5,-0.86602540378443864676372317075294); M=M * sqrtn(Q + D,3.0); std:complex<double>y3=M + P / M; std:complex<double>m=sqrtn(b4 * b4 + 4.0 * y3,2.0); std:complex<double>n= (b4 * y3 - 2.0 * d4) / m; std:complex<double>z1; D=sqrtn(b4 + m) * (b4 + m) - 8.0 * (y3 + n),2.0); z1=(-(b4 + m) + D) * 0.25; B=atan2(a * z1.imag(),b * z1.real(); ax=a * z1.imag(); by=b * z1.r

人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论