WGS84和BJ54坐标转换源程序

最新推荐文章于 2017-07-31 16:03:15 发布

nj66666666

最新推荐文章于 2017-07-31 16:03:15 发布

阅读量4.2k

点赞数 2

分类专栏： ArcGIS desktop

本文链接：https://blog.csdn.net/nj198624/article/details/8160415

版权

该博客介绍了如何使用7参数进行WGS84和BJ54坐标系之间的转换，包括平移、旋转和缩放操作。提供了一个名为CoordTrans7Param的类，包含了设置转换参数、坐标转换和误差计算的方法。

摘要由CSDN通过智能技术生成

两个坐标系转换一般需要平移，旋转，缩放共七参数。
Y=（1＋k)*M(x,y,z)*X+dX;
X,Y为3*1矩阵,M(x,y,z)为3*3的旋转矩阵.

public class CoordTrans7Param
{
public double[,] values=new double[7,1];
//{ {dx},{dy},{dz},{rx},{ry},{rz},{k}};
//public double dx,dy,dz,rx,ry,rz,k;
public void Set4Param(double dx,double dy,double dz,double k)
{
this.dx=dx;
this.dy=dy;
this.dz=dz;
this.k=k;
this.rx=this.ry=this.rz=0;
}
public void SetRotationParamRad(double rx,double ry,double rz)
{
this.rx=rx;
this.ry=ry;
this.rz=rz;
}
public void SetRotationParamMM(double rx,double ry,double rz)
{
SetRotationParamRad(rx*Math.PI/648000,ry*Math.PI/648000,rz*Math.PI/648000);
}
private double[,] GetMx()
{
double [,] Mx=new double[,]
{ {1,0,0},
{0,Math.Cos(rx),Math.Sin(rx)},
{0,-Math.Sin(rx),Math.Cos(rx)}};
return Mx;
}
private double[,] GetMy()
{
double [,] My=new double[,]
{ {Math.Cos(ry),0,-Math.Sin(ry)},
{0,1,0},
{Math.Sin(ry),0,Math.Cos(ry)}};
return My;
}
private double[,] GetMz()
{
double [,] Mz=new double[,]
{ {Math.Cos(rz),Math.Sin(rz),0},
{-Math.Sin(rz),Math.Cos(rz),0},
{0,0,1}};
return Mz;
}

private double[,] GetM() //M=Mx*My*Mz? or M=Mz*My*Mx?
{
double [,] M=new double[3,3];

MatrixTool.Multi(GetMz(),GetMy(),ref M);
MatrixTool.Multi(M,GetMx(),ref M);
return M;
}
private double[,] GetMdx()
{
double[,] mt = { { 0, 0, 0 },
{ 0, -Math.Sin(rx), Math.Cos(rx) },
{ 0, -Math.Cos(rx), -Math.Sin(rx) }};

double[,] m=new dou