geopy使用详解

geopy使用详解

  由于专业需要,经常接触一些地理处理的工具包,文档都是英文的,自己看的同时将其翻译一下,一方面自己学习的同时有个记录,要是能同时给一起的学习的童鞋们一些帮助,想想也是极好的。以下的文档内容主要翻译自官方文档,水平有限,错误的地方希望大家指出。

  处理地理数据经常会涉及到地理编码的问题。地理编码指的是将地理信息转化成坐标关系的过程。分为正向和反向的编码。正向的是指将地址信息转换为坐标点,比如:武汉市武汉大学-->(114.3594147, 30.5401222);反向地理编码就是将地理坐标转换为具体的地址,是一个与前面相反的过程。

  基于python的地理编码库geopy 是用于地理编码的常用工具,使用它 可获取多种地图服务的坐标。目前Python2和Python3下都支持。Python开发者可以使用geopy很容易的获取全球的某个街道地址,城市, 国家和地块的地理坐标,它是通过第三方的地理编码器和数据源来解析的。

geopy可以使用的地理定位服务如下:

OpenStreetMapNominatim,ESRIArcGISGoogleGeocodingAPI(V3)BaiduMapsBingMapsAPI,Yahoo!PlaceFinderYandexIGNFranceGeoNamesNaviData,OpenMapQuestWhat3WordsOpenCageSmartyStreetsgeocoder.us, and GeocodeFarm等. 这些丰富的地理编码器在geopy.geocoders模块下面,它提供各个地理定位服务的API类。每个地理编码器至少定义了一个geocode方法, 用于将字符串转换为地理位置;也定义了一个reverse方法,用于将地理坐标转换为具体地址。每个地理编码器需要接受认证和设置才能使用它的服务,例 如:在初始化时,需要一个API key或者locale。

  geopy通过了CPython 2.7, CPython 3.2, CPython 3.4, PyPy, and PyPy3下的测试。

安装:

可以通过pip或者easy_install方式安装:

pip install geopy

easy_install geopy

或者通过从PyPI下载wheel或源文件安装

下面将通过实例演示一下geopy的使用

地理编码

复制代码

>>> from geopy.geocoders import Nominatim>>> geolocator = Nominatim()>>> location = geolocator.geocode("175 5th Avenue NYC")>>> print(location.address)
Flatiron Building, 175, 5th Avenue, Flatiron, New York, NYC, New York, ...>>> print((location.latitude, location.longitude))
(40.7410861, -73.9896297241625)>>> print(location.raw)
{'place_id': '9167009604', 'type': 'attraction', ...}

复制代码

反地理编码

复制代码

>>> from geopy.geocoders import Nominatim>>> geolocator = Nominatim()>>> location = geolocator.reverse("52.509669, 13.376294")>>> print(location.address)
Potsdamer Platz, Mitte, Berlin, 10117, Deutschland, European Union>>> print((location.latitude, location.longitude))
(52.5094982, 13.3765983)>>> print(location.raw)
{'place_id': '654513', 'osm_type': 'node', ...}

复制代码

计算距离

geopy能使用经纬度距离公式(Vincenty distance) 或球面距离(great-circle distance)公式在两点间计算测地距离。在geopy中用的经纬度距离是默认的方式,类为geopy.distance.distance,计算距离为其属性(e.g., miles, meters, etc)。

给出计算经纬度距离的例子如下:

>>> from geopy.distance import vincenty>>> newport_ri = (41.49008, -71.312796)>>> cleveland_oh = (41.499498, -81.695391)>>> print(vincenty(newport_ri, cleveland_oh).miles)538.3904451566326

使用球面距离:

>>> from geopy.distance import great_circle
>>> newport_ri = (41.49008, -71.312796)
>>> cleveland_oh = (41.499498, -81.695391)
>>> print(great_circle(newport_ri, cleveland_oh).miles)
537.1485284062816


转载于:https://my.oschina.net/u/2306127/blog/601152

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值