单点经纬度转换坐标

MATLAB—经纬度坐标转换为平面坐标；Mercator投影变换。这个如果出现不合理转换请查看我的另一篇（正常是可以的）单点经纬度转换坐标，那个需要积分多一点（后台审核积分）。博文中含有这几个坐标转换的文档解释

function [XY]=transformation（latlonB0L0）
% lat为纬度，lon为经度（单位为弧度）
% B0：坐标原点的纬度，L0：坐标原点的经度，单位为弧度（（rad）;
% X：水平直角坐标， Y：纵向直角坐标，单位为米（m）;
B=lat*pi/180;L=lon*pi/180;%弧度转换为度数
a=6378137;% a：地球椭球体长半轴，单位为米（m）;
b=6356752.3