全国省会城市经纬度.doc

全国省会城市经纬度2010-04-25 18:19北京E11628N3954上海E12129N3114天津E11711N3909重庆E10632N2932哈尔滨E12641N4545长春E12519N4352沈阳E12324N4150呼和浩特E11148N4049石家庄E11428N3802太原E11234N3752济南E117N3638郑州E11342N3448西安E10854N3416兰州E10349N3603银川E10616N3820西宁E10145N3638乌鲁木齐E8736N4348合肥E11718N3151南京E11850N3202杭州E12009N3014长沙E113N2811南昌E11552N2841武汉E11421N3037成都E10405N3039贵阳E10642N2635福州E11918N2605台北E12131N2503广州E11315N2308海口E11020N2002南宁E10820N2248昆明E10241N25拉萨E9008 N2939香港E11410N2218澳门E11335N2214台北E12131N2503把地球当作一个正常的球体（其实它是椭球）来说，球面两点之间的距离计算并不复杂，运用球坐标很容易就能计算出两点之间的弧长。当然这都是高中的知识，我 和你一样，也没有那个耐心来将其推导，所以我就利用google map的经纬度到距离计算的js脚本，将球面弧长的公式给还原出来（估计这个公式是经过部分修正的），还原出来的公式如下：公式解释如下： 公式中经纬度均用弧度表示，角度到弧度的转化应该是很简单的了吧，若不会，依然请参考这个这个经纬度算距离的工具；Lat1 Lung1 表示A点经纬度，Lat2 Lung2 表示B点经纬度；a=Lat1 Lat2 为两点纬度之差 b=Lung1 -Lung2 为两点经度之差；6378.137为地球半径，单位为公里；matlab代码：adress=北京;上海;天津;重庆;哈尔滨;长春;沈阳;呼和浩特;石家庄;太原;济南;郑州;西安;兰州;银川;西宁;乌鲁木齐;合肥;南京;杭州;长沙;南昌;武汉;成都;贵阳;福州;台北;广州;海口;南宁;昆明;拉萨;香港;香港;澳门;台北;%经纬度，第一列是东经E，第二列是北纬NEN=116.28,39.54;121.29,31.14;117.11,39.09;106.32,29.32;126.41,45.45;125.19,43.52;123.24,41.5;111.48,40.49;114.28,38.02;112.34,37.52;117,36.38;113.42,34.48;108.54,34.16;103.49,36.03;106.16,38.2;101.45,36.38;87.36,43.48;117.18,31.51;118.5,32.02;120.09,30.14;113,28.11;115.52,28.41;114.21,30.37;104.05,30.39;106.42,26.35;119.18,26.05;121.31,25.03;113.15,23.08;110.2,20.02;108.2,22.48;102.41,25;90.08,29.39;114.1,22.18;114.1,22.18;113.35,22.14;121.31,25.03;%计算任意两点间的直线距离ENNew(:,1)=(floor(EN(:,1)+(EN(:,1)-floor(EN(:,1)/60)*pi/180;ENNew(:,2)=(floor(EN(:,2)+(EN(:,2)-floor(EN(:,2)/60)*pi/180;%计算任意两点的距离，使用公式$2*6378.137*arcsinsqrtsin2(a)+cos(Lat2)*cos(Lat2)*sin2(b/2)$%Lat1 Lung1 表示A点经纬度，Lat2 Lung2 表示B点经纬度；%a=Lat1 Lat2 为两点纬度之差 b=Lung1 -Lung2 为两点经度之差；%Lat1 Lung1 表示A点经纬度，Lat2 Lung2 表示B点经纬度；SizeEN=size(EN,1);Dis=zeros(SizeEN,SizeEN);for i=1:SizeENfor j=1:SizeENWeiDuCha=EN(i,2)-EN(j,2);JingDuCha=EN(i,1)-EN(j,1)