以物流行业为例，分析了 postgresql 与 greenplum 在地理位置信息处理，最佳路径算法，机器学习等方面的物流行业应用方法。其中提到了地址转换成坐标的问题，更专业些的名词应该是“地理编码”，即知道一个地址，如北京市海淀区上地十街10号，怎么样可以获取到对应的经纬度位置信息(40,116)，或者反过来。
地理编码概念
很多地图相关的厂商都提供了相关的api，我们可以直接利用这些api得到这些信息。比如百度的geocoding api。
geocoding api是一类接口，用于提供从地址到经纬度坐标或者从经纬度坐标到地址的转换服务，用户可以使用c# 、c++、java等开发语言发送请求且接收json、xml的返回数据。geocoding api包括地址解析和逆地址解析功能：
安装
pip install geocoder
地理编码
import geocoder
g = geocoder.google(“1403 washington ave, new orleans, la 70130″)
g = geocoder.arcgis(u”北京市海淀区上地十街10号”)
g.latlng
输出为
[29.9287839, -90.08421849999999]
也可以查看完整的geojson
g.geojson
输出为
{‘bbox’: [-90.0855674802915,
29.9274349197085,
-90.0828695197085,
29.9301328802915],
‘geometry’: {‘coordinates’: [-90.08421849999999, 29.9287839],
‘type’: ‘point’},
‘properties’: {‘accuracy’: u’rooftop’,
‘address’: u’1403 washington ave, new orleans, la 70130, usa’,
‘bbox’: [-90.0855674802915,
29.9274349197085,
-90.0828695197085,
29.9301328802915],
‘city’: u’new orleans’,
‘confidence’: 9,
‘country’: u’us’,
‘county’: u’orleans parish’,
‘encoding’: ‘utf-8’,
‘housenumber’: u’1403′,
‘lat’: 29.9287839,
‘lng’: -90.08421849999999,
‘location’: ‘1403 washington ave, new orleans, la 70130’,
‘neighborhood’: u’garden district’,
‘ok’: true,
‘place’: u’chijgyfhwc2liiyrysoneaxauiw’,
‘postal’: u’70130′,
‘provider’: ‘google’,
‘quality’: u’street_address’,
‘state’: u’la’,
‘status’: ‘ok’,
‘status_code’: 200,
‘street’: u’washington ave’},
‘type’: ‘feature’}
直接用google尝试查询中文地址时失败
g = geocoder.google(u”北京市海淀区上地十街10号”)
g.ok
输出为
false
用百度应该没问题，不过我没有申请相应的key。切换到arcgis，能够成功编码
g = geocoder.arcgis(u”北京市海淀区上地十街10号”)
g.latlng
输出为
[40.050934, 116.30079]
逆地理编码
g = geocoder.google([29.9287839, -90.08421849999999], method=’reverse’)
print g.address
print g.city
print g.state
print g.country
输出为
1403 washington ave, new orleans, la 70115, usa
new orleans
la
us
换成中国的地址
g = geocoder.google([40.050934, 116.30079], method=’reverse’)
print g.address
print g.city
print g.state
print g.country
输出为
bai du da sha, haidian qu, beijing shi, china, 100193
beijing
beijing shi
cn
用arcgis的服务试试
g = geocoder.arcgis([40.050934, 116.30079], method=’reverse’)
print g.address
print g.city
print g.state
print g.country
输出为
none
北京市
北京市
chn
google转换成的是英文，但地址比较全。arcgis虽然是中文，但是详细的地址居然输出为了none，这有个x用。
其他
geocoder 的功能不止于此，它还可以查ip（包括自己的）。
g = geocoder.ip(‘199.7.157.0’)
print g.latlng
print g.city
g = geocoder.ip(‘me’)
print g.latlng
print g.city
输出为
[43.6934, -79.4857]
toronto
[51.05, 13.75]
dresden
查询一个城市的空间包围盒
g = geocoder.arcgis(u”山东”)
g.bbox
输出为
{‘northeast’: [38.976997, 121.976998], ‘southwest’: [33.022997, 116.022998]}
小结
空间信息既可以利用行政区划、自然地理区域等文本信息描述，也可以用坐标系统、数字（邮编等）来标识。利用地理编码技术，可以将空间信息的地理定位要素与相应的文本信息关联起来。本文主要介绍了geocoder地理编码这一小工具，可以方便快捷的利用地图等相关厂商提供的地理编码服务，将文字描述的位置转换成地图上的经纬度，或者通过地图上的某个位置坐标获得相应的位置信息文字描述。