| 基于MATLAB的坐标转换系统的设计与实现 |
|
|
|
网上，需对旧坐标加以平移、旋转和尺度因子改正，以保持形状不变。
已知新旧坐标系的坐标转换方程为：
式中a、b、c、d为所求的未知量，即平差参数。由于平差过程中可能出现病态矩阵的问题，故平差前应使新旧坐标同时减去各自平均值进行中心化处理。用处理后的数据列出误差方程，计算四参数。
2.3 高程坐标转换模型
在本系统中，采用多项式拟合模型进行GPS水准的高程异常拟合。假设是点i的高程异常（似大地水准面至地球椭球面的高度），(xi,yi)为点的平面坐标，视为无误差，并认为是关于平面坐标的函数，即可取拟合函数为：
在的条件下，求出参数b0,b1,...,bn。
3 系统设计流程和程序的主要代码
3.1 系统设计流程
3.2 平面四参数部分求解代码
x1=pt(1,1);y1=pt(1,2);x2=pt(2,1);y2=pt(2,2);
xx1=ps(1,1);yy1=ps(1,2);xx2=ps(2,1);yy2=ps(2,2);
cscz=inv([1 0 x1 -y1;0 1 y1 x1;1 0 x2 -y2;0 1 y2 x2])*[xx1;yy1;xx2;yy2];
bb=[];ll=[];
for i= 1:a
bbi=[1 0 pt(i,1) -pt(i,2);0 1 pt(i,2) pt(i,1)];
li=[ps(i,1);ps(i,2)]-bbi*cscz;
bb=[bb;bbi];ll=[ll;li];
end
csgz=inv(bb'*bb)*bb'*ll;
pmcs=cscz+csgz;
3.3 高程拟合部分求解代码
x1=x.*x;
y1=y.*y;
y2=x.*y;
n=[ones(b,1),x',y',x1',y2',y1'];
[xs,bint,r,rint,s]=regress(z',n);
gccs=xs(6)*y.*y+xs(5)上一页 [1] [2] [3] [4] 下一页 |
|
|
| |
|
上一个论文： 西瓜常见生理病害与防治
下一个论文： 基于J2EE在线考试系统的设计与实现 |
|
用户登录 
新用户注册 