datatrans.m

来自「MATLAB语言下的(有源码):Mike Zero系列软件的前后处理 GE中岸线」· M 代码 · 共 19 行

function datatrans(X0,Y0)
Data=load('data.txt','-ascii');
X=Data(:,1);
Y=Data(:,2);
Z=data(:,3);
N=length(X);
K_Line=(Y0(1)-Y0(2))./(X0(1)-X0(2));
AngleLine=atan(K_Line); %弧度
for c1=1:N
    xData=X(c1);
    yData=Y(c1);
    K_Data=(Y0(1)-yData)./(X0(1)-xData);
    Angle_Data=atan(K_Data); %弧度
    Angle_Use=Angle_Data-AngleLine; %弧度
    Distance_Data=sqrt((Y0(1)-yData).^2+(X0(1)-xData).^2);
    Result(c1,1)=Distance_Data.*cos(Angle_Use);
end
Result(:,2)=Z;
save result.txt Result -ascii;