如題所說，直接上程式。驗證自己做一下，結果應該是對的。誠不我欺。
（注意：程式名與函式名保持一致！！！）
%地理經緯度座標轉換為utm座標
size_shuzu=size(lat_shuru);
for i=1:size_shuzu(2)
%輸入經緯度
% lat=29.819206;
% lon=116.133243;
lat=lat_shuru(i);
lon=lon_shuru(i);
%% %%unit：km
a=6378.137;
e=0.0818192;
k0=0.9996;
e0=500;
n0=0;
zonenum=fix(lon/6)+31;
lamda0=(zonenum-1)*6-180+3;%degree
lamda0=lamda0*pi/180;%radian
phi=lat*pi/180;
lamda=lon*pi/180;%radian
v=1/sqrt(1-e^2*(sin(phi)^2));
a=(lamda-lamda0)*cos(phi);
t=tan(phi)*tan(phi);
c=e^2*cos(phi)*cos(phi)/(1-e^2);
s=(1-e^2/4-3*e^4/64-5*e^6/256)*phi-(3*e^2/8+3*e^4/32+45*e^6/1024)*sin(2*phi)+(15*e^4/256+45*e^6/1024)*sin(4*phi)-35*e^6/3072*sin(6*phi);
utme=e0+k0*a*v*(a+(1-t+c)*a^3/6+(5-18*t+t^2)*a^5/120);
utmn=n0+k0*a*(s+v*tan(phi)*(a^2/2+(5-t+9*c+4*c^2)*a^4/24+(61-58*t+t^2)*a^6/720));
shuchu(1,i)=utme;
shuchu(2,i)=utmn;
end
經緯度座標系轉東北天 經緯度座標系轉換
網際網路地圖的座標系現狀 地球座標 wgs84 國際標準，從 gps 裝置中取出的資料的座標系 國際地圖提供商使用的座標系 火星座標 gcj 02 也叫國測局座標系 中國標準，從國行移動裝置中定位獲取的座標資料使用這個座標系 國家規定 國內出版的各種地圖系統 包括電子形式 必須至少採用gcj 02對...
經緯度座標轉換
經緯度與xy座標轉換工具類 public class latlngxyzconverthelper xy轉經緯度 需要轉換的x座標 需要轉換的y座標 地圖級別 轉換後的經度 轉換後的緯度 public static void tiletolatlng double tile x,double til...
不同座標系經緯度轉換演算法（JS）
墨卡托投影座標轉經緯度 export function mercatortolonlat mercator var x mercator.x 20037508.34 180 var y mercator.y 20037508.34 180 y 180 math.pi 2 math.atan math...