datatrans.m

MATLAB语言下的(有源码):Mike Zero系列软件的前后处理 GE中岸线与MIKE岸线相互转换 Mapinfo与MIKE中数据相互转换.(对学习MIKE与学习MATLAB编程人皆有帮助,且有说

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;