诸君安! geodist:计算测地距离,即最短曲线的长度,是在沿着地球表面的两个点之间的距离 语法: geodist lat1 lon1 lat2 lon2 [if] [in] , generate(newvar) [options] 注: 纬度的范围在-90°至90°之间,经度的范围在-180°到180°之间。北纬0-90°与东经0-180°为正,南纬0-90°与西经0-180°为负 坐标可以用变量,标量或数字坐标字符串来表示。如果坐标是标量,或数字坐标字符串,则无须添加generate(newvar) 这里,geodist具有如下几个选项option: miles:用miles来测度距离 ellipsoid(a,f):可参考椭圆体中两个定义常数(a,f),a即为长半轴,f即为短半轴,长半轴是地心到赤道的球半径长度,短半轴是地心到两级的球半径长度。这些是以米为单位的半长 轴的长度(即到赤道的半径)和扁平率。具体的计算方法可参考: http://www.movable-type.co.uk/scripts/latlong-vincenty.html maxiter(#):表示执行计算椭圆距离时的最大迭代次数。默认值为25。请注意实际上所有情况都需要少于10次迭代 sphere:计算大圆距离。geodist默认地球近似为球体,其平均半径为6371公里。具体的计算方法可参考: http://www.movable-type.co.uk/scripts/latlong.html radius:在计算大圆距离时,自定义半径为多少(以千米为单位)的球体 那么我们该如何测算距离呢?首先,我们得安装geodist命令: ssc install geodist,replace 然后,这里我们选择了一组经纬度的数据分别用不同的option测量一下,会得到不同的结果: geodist 42.270873 -83.726329 45.545448 -73.639075 以英里来表示距离: geodist 42.270873 -83.726329 45.545448 -73.639075,miles 计算大圆距离: geodist 42.270873 -83.726329 45.545448 -73.639075, sphere 自定义半径为6399.6km来测度距离: geodist 42.270873 -83.726329 45.545448 -73.639075,radius(6399.6) 定义椭球的长半轴为6378136.49米,短半轴为6356755.00米: geodist 42.270873 -83.726329 45.545448 -73.639075, ellipsoid(6378136.49, 6356755.00) 而在我们的公众号在2016年9月24号发表了一文《爬虫俱乐部周末送大礼—chinagcode提取中文地址的经纬度》,我们可以通过中文地址来提取经纬度,那么,我们是不是就可以利用geodist命令来将经纬度转化成距离呢? 那么,我们就举个下面这个简单的例子来说明: 我们可以分别输入两所大学的所在具体地点,其中具体的地点中得包括省、市、区、地址,因为有的大学同时存在多个分校与校区,所以我们就必须得精确其地址: clear set more off cap mkdir 'D:/distance' cd D:/distance input str100 fulladdr1 str100 fulladdr2 '湖北省武汉市洪山区中南财经政法大学南湖校区' '湖北省武汉市洪山区武汉理工大学南湖校区' end list 我们将百度地图API的秘钥输入选项baidukey()中,得到我们需要的经纬度: chinagcode, baidukey(40GD1PPONTXxKdrktzhGtC0NUyXsHP5c) /// fulladdress(fulladdr1) lat(lat1) long(lon1) chinagcode, baidukey(40GD1PPONTXxKdrktzhGtC0NUyXsHP5c) /// fulladdress(fulladdr2) lat(lat2) long(lon2) list 最后,我们生成一个新的变量distance来测度其距离: geodist lat1 lon1 lat2 lon2,gen(distance) list fulladdr1 fulladdr2 distance save 'distance.dta',replace 接下来是空气质量报告 全国空气质量如下 全国空气质量较佳 天津稍有逊色 以上就是今天给大家分享的内容了,说得好就赏个铜板呗!有钱的捧个钱场,有人的捧个人场~,点赞打赏随您心意,么么哒~ 应广大粉丝要求,爬虫俱乐部的推文公众号打赏功能可以开发票啦,累计打赏超过1000元我们即可给您开具发票,发票类别为“咨询费”。用心做事,只为做您更贴心的小爬虫。第一批发票已经寄到各位小主的手中,大家快来给小爬虫打赏呀~ 文字编辑:张 欢 技术总编:刘贝贝 |
geodist——测测我们相距多远
特别重申:本篇文档资料为 “好网角收藏夹” 注册用户(收藏家)上传共享,仅供参考之用,请谨慎辨别,不代表本站任何观点。
好网角收藏夹为网友提供资料整理云存储服务,仅提供信息存储共享平台。
如发现不良信息删除、涉嫌侵权,请 点击这里举报 ,或发送邮件到:dongye2016qq.com。
如发现不良信息删除、涉嫌侵权,请 点击这里举报 ,或发送邮件到:dongye2016qq.com。