(2017-07-02 银河统计)

地理编码(Geocoding)指的是将统计资料或是地址信息建立空间坐标关系的过程。地址经过地理编码后,便可在地图中通过经纬度显示各地址的位置(正向地理编码:地址 - 坐标)。反之,也可以通过经纬度解析出对应地址(逆向地理编码:坐标 - 地址)。地理编码也称为地址解析。

一、标准地址规范

当我们书写地址是,应该从大到小分级书写,地址等级可参考如下分级:

  • 第一级省
  • 第二级市
  • 第三级区、县
  • 第四级镇、乡、街道
  • 第五级路、巷、行政村
  • 第六级小区、自然村
  • 第七级门牌、村组
  • 第八级楼号、办公楼房(单位)
  • 第九级单元号、房号

如,黑龙江省哈尔滨市道里区通达街138号。如果城市较大,可以省略省份。另外,地理编码一般只解析到门牌号,没必要书写楼层或房间号等细节。

二、地址解析

本文根据百度地图开发平台介绍地址和经纬度间的相互转换。

1、在网页中引用百度地图API

<!DOCTYPE html>
<html>
<head>
   <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
   <script type="text/javascript" src="http://api.map.baidu.com/api?v=2.0&ak=您的密钥">
   </script>
   <title>地址解析</title>
</head>
<body></body>
</html>

使用前登陆百度地图开放平台开发者注册页(http://lbsyun.baidu.com/apiconsole/key),获得自己的密钥。

2、百度地址解析

 地址解析

经度  纬度  墨卡托X  墨卡托Y 

注:在地址栏中输入国内任何地址,可获得该地址经纬度和墨卡托坐标

代码样例

<!DOCTYPE html>
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
<title>地址解析</title>
<script type="text/javascript" src="http://api.map.baidu.com/api?v=2.0&ak=您的密钥"></script>
</head>
<body>
<center>地址 <input  value="北京市海淀区上地10街" style="width:396px;"> 
<button onclick="setJWD(document.getElementById('myAddress').value)">地址解析</button>
<div  style="width:500px; height:300px; border:#aaaaaa solid 1px;"></div>
经度 <input  style="width:80px;"> 纬度 <input  style="width:80px;"> 
墨卡托X <input  style="width:80px;"> 墨卡托Y <input  style="width:80px;"> 
</center>
</body>
</html>

<script type="text/javascript" charset="UTF-8">

var map = new BMap.Map("myMap");
setJWD(document.getElementById("myAddress").value);

function setJWD(addr) {
var point = new BMap.Point(116.331398,39.897445);
map.centerAndZoom(point,12);
var myGeo = new BMap.Geocoder(); // 创建地址解析器实例
var mpObj = new BMap.MercatorProjection(); // 墨卡托坐标对象
myGeo.getPoint(addr, function(point){
if (point) {
    	map.centerAndZoom(point, 16); //设置地图中心点和显示级别
	map.addOverlay(new BMap.Marker(point)); //加载中心点图层
            document.getElementById("myJD").value=point.lng; //获得经度
            document.getElementById("myWD").value=point.lat;  //获得纬度
            var mpoint = mpObj.lngLatToPoint(new BMap.Point(point.lng, point.lat)); //经纬度转换为墨卡托坐标
            document.getElementById("myMX").value=mpoint.x; //获得墨卡托横坐标
            document.getElementById("myMY").value=mpoint.y; //获得墨卡托纵坐标                      
}}, "");
}
</script>

3、百度批量地址解析

批量地址格式文本(地址加引号并用逗号分隔)

"合肥市包河区金寨路1号(金寨路与望江西路交叉口)","合肥市庐阳区凤台路209号(凤台路与蒙城北路交叉口)","合肥市蜀山区金寨路217号(近安医附院公交车站)","合肥市蜀山区梅山路10号(近安徽饭店) ","合肥市蜀山区 长丰南路159号铜锣湾广场312室","合肥市寿春路93号钱柜星乐町KTV(逍遥津公园对面)","合肥市庐阳区长江中路177号","合肥市新站区胜利路89"

注:按格式批量添加地址,然后运行“批量地址解析”按钮

批量地址解析

批量地址经纬度坐标

代码样例

<!DOCTYPE html>
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf8" />
<meta name="viewport" content="initial-scale=1.0, user-scalable=no" />
<title>批量地址</title>
<style type="text/css">
	body, html{width: 100%;height: 100%;margin:0;font-family:"微软雅黑";}
	#l-map{height:300px;width:100%;}
	#r-result{width:100%; font-size:14px;line-height:20px;}
</style>
<script type="text/javascript" src="http://api.map.baidu.com/api?v=2.0&ak=您的密钥"></script>
</head>
<body>
<div id="l-map"></div>
<div >
	<input type="button" value="批量地址解析" onclick="bdGEO()" />
	<div ></div>
</div>
</body>
</html>
<script type="text/javascript">
// 百度地图API功能
var map = new BMap.Map("l-map");
map.centerAndZoom(new BMap.Point(117.269945,31.86713), 13);
map.enableScrollWheelZoom(true);
var index = 0;
var myGeo = new BMap.Geocoder();
var adds = [
	"包河区金寨路1号(金寨路与望江西路交叉口)",
	"庐阳区凤台路209号(凤台路与蒙城北路交叉口)",
	"蜀山区金寨路217号(近安医附院公交车站)",
	"蜀山区梅山路10号(近安徽饭店) ",
	"蜀山区 长丰南路159号铜锣湾广场312室",
	"合肥市寿春路93号钱柜星乐町KTV(逍遥津公园对面)",
	"庐阳区长江中路177号",
	"新站区胜利路89"
];
function bdGEO(){
	var add = adds[index];
	geocodeSearch(add);
	index++;
        }
function geocodeSearch(add){
	if(index < adds.length){
		setTimeout(window.bdGEO,400);
	} 
	myGeo.getPoint(add, function(point){
		if (point) {
			document.getElementById("result").innerHTML +=  index + "、" + add + ":" + point.lng + "," + point.lat + "</br>";
			var address = new BMap.Point(point.lng, point.lat);
			addMarker(address,new BMap.Label(index+":"+add,{offset:new BMap.Size(20,-10)}));
		}
	}, "合肥市");
}
// 编写自定义函数,创建标注
function addMarker(point,label){
	var marker = new BMap.Marker(point);
	map.addOverlay(marker);
	marker.setLabel(label);
}
</script>

注:代码摘自百度API示例文档

4、百度逆地址解析

 纬度  逆地址解析

省份  城市  区镇  街道  牌号 

地址 

注:可替换经纬度(如,经度126.616759、纬度45.74989),点击“逆地址解析”按钮获得明码地址

代码样例

<!DOCTYPE html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
<meta name="viewport" content="initial-scale=1.0, user-scalable=no" />
<style type="text/css">
	body, html{width: 100%;height: 100%;margin:0;font-family:"微软雅黑";font-size:14px;}
	#allmap {width:100%;height:500px;}
</style>
<script type="text/javascript" src="http://api.map.baidu.com/api?v=2.0&ak=您的密钥"></script>
<title>逆地址解析</title>
</head>
<body>
<div id="allmap"></div>
<p>点击地图展示详细地址</p>
</body>
</html>
<script type="text/javascript">
// 百度地图API功能
var map = new BMap.Map("allmap");
var point = new BMap.Point(116.331398,39.897445);
map.centerAndZoom(point,12);
var geoc = new BMap.Geocoder();    

map.addEventListener("click", function(e){        
	var pt = e.point;
	geoc.getLocation(pt, function(rs){
		var addComp = rs.addressComponents;
		alert(addComp.province + ", " + addComp.city + ", " + addComp.district + ", " + addComp.street + ", " + addComp.streetNumber);
	});        
});
</script>

注:代码摘自百度API示例文档

5、百度批量逆地址解析

[[116.307852,40.057031],[116.313082,40.047674],[116.328749,40.026922],[116.347571,39.988698],[116.316163,39.997753],[116.345867,39.998333],[116.403472,39.999411],[116.307901,40.05901]]

注:数组按格式批量添加经纬度,然后运行“批量反地址解析”按钮



代码样例

<!DOCTYPE html>
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=gb2312" />
<meta name="viewport" content="initial-scale=1.0, user-scalable=no" />
<title>批量反地址解析+商圈</title>
<style type="text/css">
	body, html{width: 100%;height: 100%;margin:0;font-family:"微软雅黑";}
	#l-map{height:300px;width:100%;}
	#r-result{width:100%; font-size:14px;line-height:20px;}
</style>
<script type="text/javascript" src="http://api.map.baidu.com/api?v=2.0&ak=您的密钥"></script>
</head>
<body>
<div id="l-map"></div>
<div >
	<input type="button" value="批量反地址解析+商圈" onclick="bdGEO(0)" />
	<div ></div>
</div>
</body>
</html>
<script type="text/javascript">
// 百度地图API功能
var map = new BMap.Map("l-map");
map.centerAndZoom(new BMap.Point(116.328749,40.026922), 13);
map.enableScrollWheelZoom(true);
var index = 0;
var myGeo = new BMap.Geocoder();
var adds = [
	new BMap.Point(116.307852,40.057031),
	new BMap.Point(116.313082,40.047674),
	new BMap.Point(116.328749,40.026922),
	new BMap.Point(116.347571,39.988698),
	new BMap.Point(116.316163,39.997753),
	new BMap.Point(116.345867,39.998333),
	new BMap.Point(116.403472,39.999411),
	new BMap.Point(116.307901,40.05901)
];
for(var i = 0; i<adds.length; i++){
	var marker = new BMap.Marker(adds[i]);
	map.addOverlay(marker);
	marker.setLabel(new BMap.Label("我是商圈:"+(i+1),{offset:new BMap.Size(20,-10)}));
}
function bdGEO(){	
	var pt = adds[index];
	geocodeSearch(pt);
	index++;
}
function geocodeSearch(pt){
	if(index < adds.length-1){
		setTimeout(window.bdGEO,400);
	} 
	myGeo.getLocation(pt, function(rs){
		var addComp = rs.addressComponents;
		document.getElementById("result").innerHTML += index + ". " +adds[index-1].lng + "," + adds[index-1].lat + ":"  + "商圈(" + rs.business + ")  结构化数据(" + addComp.province + ", " + addComp.city + ", " + addComp.district + ", " + addComp.street + ", " + addComp.streetNumber + ")<br/><br/>";
	});
}
</script>

注:代码摘自百度API示例文档

地址正、逆解析精度取决于系统数据库数据采集密度和更新速度。国内地址服务较好的为百度和高德地图。在实践中,地址解析、特别是大数据批量地址解析较为常用。凡是和位置有关的自然和社会现象,通过地址解析获得数量化指标,从而为统计量化分析、特别是空间统计分析奠定坚实基础。

本文提供的各种地址解析工具仅用于个人学习或教学研究用途,不支持大批量商业用途地址解析(由于百度地图key的限制)。

地理编码 python 地理编码数据_地址解析


©哈尔滨商业大学 银河统计工作室


银河统计工作室成员由在校统计、计算机部分师生和企业数据数据分析师组成,维护和开发银河统计网和银河统计博客(技术文档)。专注于数据挖掘技术研究和运用,探索统计学、应用数学和IT技术有机结合,尝试大数据条件下新型统计学教学模式。