zbzh.m

这是用于遥感数据处理中的相对定向及绝对定向及坐标的七参数转换的程序。具有一定代表性

[FileName,PathName]=uigetfile('*.txt','打开文件');
if FileName==0
    return
end
fn=fopen([PathName,FileName]);
%------------------------------------------------------模型点的转换---------------------------------------------

IJz=fscanf(fn,'%f',[2,26]);
IJz=IJz';
[c,r]=size(IJz);
XY=[[linspace(1,1,c)]' IJz]*NEI;
XYZz=[XY [linspace(-f,-f,c)]']';
IJy=fscanf(fn,'%f',[2,26]);
IJy=IJy';
XY=[[linspace(1,1,c)]' IJy]*NEI;
XYZy=[XY [linspace(-f,-f,c)]']';
Bx=0.09216*(1-0.6)*10000;
R1(1,1)=cos(phy1)*cos(kapa1)-sin(phy1)*sin(0)*sin(kapa1);
R1(1,2)=-cos(phy1)*sin(kapa1)-sin(phy1)*sin(0)*cos(kapa1);
R1(1,3)=-sin(phy1)*cos(0);
R1(2,1)=cos(0)*sin(kapa1);
R1(2,2)=cos(0)*cos(kapa1);
R1(2,3)=-sin(0);
R1(3,1)=sin(phy1)*cos(kapa1)+cos(phy1)*sin(0)*sin(kapa1);
R1(3,2)=-sin(phy1)*sin(kapa1)+cos(phy1)*sin(0)*cos(kapa1);
R1(3,3)=cos(phy1)*cos(0);

R2(1,1)=cos(phy2)*cos(kapa2)-sin(phy2)*sin(om2)*sin(kapa2);
R2(1,2)=-cos(phy2)*sin(kapa2)-sin(phy2)*sin(om2)*cos(kapa2);
R2(1,3)=-sin(phy2)*cos(om2);
R2(2,1)=cos(om2)*sin(kapa2);
R2(2,2)=cos(om2)*cos(kapa2);
R2(2,3)=-sin(om2);
R2(3,1)=sin(phy2)*cos(kapa2)+cos(phy2)*sin(om2)*sin(kapa2);
R2(3,2)=-sin(phy2)*sin(kapa2)+cos(phy2)*sin(om2)*cos(kapa2);
R2(3,3)=cos(phy2)*cos(om2);

XYZ1=R1*XYZz;
XYZ2=R2*XYZy;

for i=1:c
    N1=Bx*XYZ2(3,i)/(XYZ1(1,i)*XYZ2(3,i)-XYZ2(1,i)*XYZ1(3,i));
    N2=Bx*XYZ1(3,i)/(XYZ1(1,i)*XYZ2(3,i)-XYZ2(1,i)*XYZ1(3,i));
    XYZa(1,i)=N1*XYZ1(1,i);
    XYZa(2,i)=0.5*(N1*XYZ1(2,i)+N2*XYZ2(2,i));
    XYZa(3,i)=N1*XYZ1(3,i)+1200;
end
dXYZ=[dx;dy;dz];
for i=1:c-1
    dXYZ=[dXYZ [dx;dy;dz]];
end
R(1,1)=cos(ey)*cos(ez);
R(1,2)=cos(ex)*sin(ez)+sin(ex)*sin(ey)*cos(ez);
R(1,3)=sin(ex)*sin(ez)-cos(ex)*sin(ey)*cos(ez);
R(2,1)=-cos(ey)*sin(ez);
R(2,2)=cos(ex)*cos(ez)-sin(ex)*sin(ey)*sin(ez);
R(2,3)=sin(ex)*cos(ez)+cos(ex)*sin(ey)*sin(ez);
R(3,1)=sin(ey);
R(3,2)=-sin(ex)*cos(ey);
R(3,3)=cos(ex)*cos(ey);

XYZ=dXYZ+(1+dk)*R*XYZa;
XYZ=XYZ';

fud=fopen(strcat(PathName,'坐标转换结果.txt'),'wt');
fprintf(fud,'坐标转换结果:\n   IL     JL     IR     JR             X              Y\n');
for i=1:size(XYZ)
fprintf(fud,'%6.0f %6.0f %6.0f %6.0f       %10.6f    %10.6f\n',[IJz(i,:) IJy(i,:) XYZ(i,2) XYZ(i,1)]);
end

fclose all;
open(strcat(PathName,'坐标转换结果.txt'));