我需要读第三coloumn和matlab文本文件的所有行。请给你的建议。
1 Ansicht(30天)的
Altere Kommentare anzeigen
xyz是文本文件包含100行3列。大地坐标的XYZ坐标的转换代码是下面讨论:
j = 1;
而(< 100)
X = fopen (“xyz.txt (j, 1) ');
Y = fopen (“xyz.txt (j, 2) ');
Z = fopen (“xyz.txt (j, 3) ');
%估算曲率半径;(WGS84模型被认为是)
一个= 6378137;%计;%轴;
f = 1/298.257223563;%椭球压扁;
b = *(行进);%定义长半轴;
%估计辅助值;
P =√X ^ 2 + Y ^ 2);
θ=每股(Z * / P * b);
e =√((^ 2)——(b ^ 2)) / a ^ 2);%第一地球的偏心;
e2 =√(^ 2 - b ^ 2) / b ^ 2);%地球第二偏心;
%纬度的初始值;
ppi =量化((Z) * (1 + e2 ^ 2)), P);
N = /√(单电子^ 2 * sin (ppi) ^ 2);% '垂直;
h = (P / cos (ppi)) - N;%的高度;
ppi =量化(Z) * (1 + e2 ^ 2 * (N / (N + h))), P);
λ=(量化(Y, X)) * 180 /π;
ppi = ppi * 180 /π;
%保存llnew_rad.txt
%显示结果
disp (“Longitutde是:”);
disp(λ);
disp(“纬度是:”);
disp (ppi);
disp(“高度(米):”);
disp (h);
流(“j的价值:% d \ n”, j);
j = + 1;
结束