单个点经纬度转换坐标

MATLAB—经纬度坐标转换为平面坐标；Mercator投影变换。（这个程序考虑第一偏心率、第二偏心率）原创度高所以积分多。博文中含有这几个坐标转换的文档解释https://blog.csdn.net/weixin_43428682/article/details/87889753。
批量点转换的链接https://download.csdn.net/download/weixin_43428682/11155590

function [XY]=transformation（latlonB0L0）
% 调用格式：[XY]=transformation（latlonB0L0）
% 输入： lat为纬度，lon为经度（单位为弧度），B0：坐标原点的纬度，L0：坐标原点的经度，单位为弧度（（rad）;
% 输出：X为水平直角坐标， Y为纵向直角坐标，单位为米（m）;
B=lat*pi/180;L=lon*pi/180;
B0=B0*pi/180;L0=L0*pi/1