本文作者：宁刘莹，中南财经政法大学金融学院
文字编辑：赵冰洁
技术总编：余术玲
导读
在之前的推文《轻松获取通行数据——cntraveltime命令介绍》里，我们介绍了用cntraveltime命令来获取两地之间的交通时间和交通路程，那么如果我们从地理的角度出发想要得到两地之间的直线距离该怎么做呢？根据中学的数学物理知识，在已知地球半径、两个地点的经纬度的情况下，列一个计算公式是可以计算出在地球表面两地之间的直线距离的，但列式子未免过于麻烦，这个由懒人驱动的世界就开发了一个命令—geodist—来帮助我们实现这个功能。给我几列经纬度数据，还你一份直线距离。
一、命令介绍
该命令的语法如下：
如果纬度/经度坐标是数值变量：
geodist lat1 lon1 lat2 lon2 [if] [in] , generate(new_dist_var) [options]
该命令有这样几个选项：
（1）miles：指定以英里为单位报告距离；
（2）ellipsoid(a,f)：自定义椭球体参数（a，f），a为一个椭球体的长半轴，f为短半轴，长半轴是地心到赤道的球半径长度，短半轴是地心到两级的球半径长度。这些是以米为单位的半长轴的长度(即到赤道的半径）和扁平率。具体的计算方法可参考: http://www.movable-type.co.uk/scripts/latlong-vincenty.html ；
（3）maxiter(#)：表示执行计算椭圆距离时的最大迭代次数，默认值为25。请注意实际上所有情况都需要少于10次迭代；
（4）sphere：计算球体的大圆距离。默认地球近似为球体，其平均半径为6371公里。具体的计算方法可参考: http://www.movable-type.co.uk/scripts/latlong.html ；
（5）radius：在计算大圆距离时，指定球体的半径（以千米为单位）。
二、举例介绍
下面我们通过一个具体的例子来说明。首先要安装这个命令：
ssc install geodist
结合爬虫俱乐部之前开发的 cngcode 命令，我们先来得到中南财经政法大学和武汉大学的经纬度。首先输入两所学校的详细地址，包括省、市、区、校名：
clear all
input str100 fulladdr1 str100 fulladdr2
"湖北省武汉市洪山区中南财经政法大学南湖校区" "湖北省武汉市湖北省武汉市武昌区武汉大学"
end
list
image
接下来使用 cngcode 得到两个地点的经纬度：
cngcode, baidukey(your secret key) fulladdress(fulladdr1) lat(lat1) long(long1)
cngcode, baidukey(your secret key) fulladdress(fulladdr2) lat(lat2) long(long2)
list
image
然后就可以使用 geodist 用得到的经纬度获取两地之间的距离：
geodist lat1 long1 lat2 long2, gen(dist)
list
image
看到距离是7.2568235公里。与地图软件结果基本一致：
image
接下来我们用 cntraveltime 命令来看一下交通距离：
cntraveltime, baidukey(your secret key) startlat(lat1) startlng(long1) endlat(lat2) endlng(long2) mode("car")
list
image
得到结果，以驾车为交通方式的路程为11103米，与苹果手机自带的地图软件对比也是基本一致的：
image
以上就是本文对geodist命令的介绍了，感兴趣的一起动手试试吧！
网友评论