将经纬度数组转换为xy坐标

如果您有一个经纬度数组，可以使用MATLAB的矢量化运算，将其转换为xy坐标数组。具体的代码如下：

% 定义中央子午线经度和地球半径
lon0 = 0;
R = 6371000; % 地球半径，单位为米

% 将经纬度转换为弧度制
lon = deg2rad(lon + 180);
lat = deg2rad(lat + 90);

% 计算墨卡托投影下的x和y坐标
x = R * (lon - lon0);
y = R * log(tan(lat/2 + pi/4));

% 对x和y坐标进行缩放和平移，使得它们都落在指定的范围内
x = (x - min(x)) / (max(x) - min(x)) * (xmax - xmin) + xmin;
y = (y - min(y)) / (max(y) - min(y)) * (ymax - ymin) + ymin;

% 将x和y坐标组成一个二维数组
xy = [x(:) y(:)];

其中，`lon`和`lat`分别表示经度和纬度，`xmin`、`xmax`、`ymin`和`ymax`分别表示x和y坐标的范围。最后，将x和y坐标组成一个二维数组`xy`，其中第一列是x坐标，第二列是y坐标。需要注意的是，如果您需要在绘图中使用这个数组，可以使用`scatter`函数或者`plot`函数进行绘制，例如：

scatter(xy(:,1), xy(:,2));

或者

plot(xy(:,1), xy(:,2));

这样就可以将经纬度数组转换为xy坐标数组，并且在图中绘制出来。