python地理处理包——geopy使用之地理编码与反地理编码

2023年1月19日 202次阅读 来源: whgiser

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

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

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

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

OpenStreetMapNominatim,ESRIArcGIS, GoogleGeocodingAPI(V3), BaiduMaps, BingMapsAPI,Yahoo!PlaceFinder, Yandex, IGNFrance, GeoNames, NaviData,OpenMapQuest, What3Words, OpenCage, SmartyStreets, geocoder.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
    原文作者：whgiser
    原文地址: https://www.cnblogs.com/giserliu/p/4982187.html
    本文转自网络文章，转载此文章仅为分享知识，如有侵权，请联系博主进行删除。
点赞