function?[?xy?]?=?BL2XY（?BLaek?）
%?此函數用于進行高斯正算

ee=e/sqrt（1-e^2）;
%?計算中央子午線經度
if?k==6
????L0=6*（fix（L/6）+1）-3;
elseif?k==3
????L0=3*（fix（（L-1.5）/3）+1）;
end
LP=（L-L0）/180*pi;
B=B/180*pi;

%?計算子午圈弧長
ap=1+3/4*e^2+45/64*e^4+175/256*e^6+11025/16384*e^8+43659/65536*e^10;
bp=3/4*e^2+15/16*e^4+525/512*e^6+2205/2048*e^8+72765/65536*e^10;
cp=15/64*e^4+105/256*e^6+2205/4096*e^8+10395/16384*e^10;
dp=35/512*e^6+315/2048*e^8+31185/131072*e^10;
ep=315/16384*e^8+3465/65536*e^10;
fp=639/131072*e^10;
X=a*（1-e^2）*（ap*B-bp/2*sin（2*B）+cp/4*sin（4*B）-dp/6*sin（6*B）+ep/8*sin（8*B）-fp/10*sin（10*B））;

%?計算高斯平面坐標
N=a/sqrt（1-e^2*sin（B）*sin（B））;
eta=ee*cos（B）;
t=tan（B）;
m=cos（B）*LP;
x=X+N*t*（0.5*m^2+1/24*（5-t^2+9*eta^2+4*eta^4）*m^4+1/720*（61-58*t^2+t^4）*m^6）;
y=N*（m+1/6*（1-t^2+eta^2）*m^3

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????901??2016-06-04?16:39??投影與坐標轉換\BL2XY.m
?????文件?????????375??2016-06-22?22:59??投影與坐標轉換\LBH2XYZ.m
?????文件?????????464??2016-06-22?23:05??投影與坐標轉換\main.m
?????文件????????1100??2016-06-22?23:17??投影與坐標轉換\Result.txt
?????文件????????1220??2016-06-22?22:25??投影與坐標轉換\XY2BL.m
?????文件????????1094??2016-06-22?23:01??投影與坐標轉換\XYZ2LBH.m
?????文件???????26112??2016-05-30?18:34??投影與坐標轉換\投影與坐標轉換-程序.v12.suo
?????文件?????????790??2016-06-22?21:40??投影與坐標轉換\原始數據.txt