经纬度坐标转换为平面坐标的matlab实现

% function B=matrix_transformation(A,B0,L0) % 调用格式：[X,Y]=transformation(lat,lon,B0,L0) % 输入：A为经纬度坐标矩阵（第一列为纬度，第二列为经度），B0：坐标原点的纬度，L0：坐标原点的经度，单位为弧度((rad); % 输出：B为平面坐标系的坐标矩阵，%第一列为水平坐标，第二列为垂直坐标 B0=0;L0=120; A=[0,120; 30,100; -30,100; -25,150];%A第一列为纬度，第二列为经度； B=zeros(c,2);%对应坐标矩阵 k=length(A(:,1)); c=1; for i=1:k if c<=k lat=A(c,1);lon=A(c,2); [X,Y]=transformation(lat,lon,B0,L0);%调用单点经纬度坐标转换函数 B(c,1)=X; B(c,2)=Y; c=c+1; end end B