将一个城市的地图按照每1平方公里进行拆分为若干个区域
将一个城市的地图按照每1平方公里进行拆分为若干个区域（地图使用高德）
如图：
核心代码（Python 3.5）：
map_zoning.py
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
import json
import math
import os
import time
from amap import Map
class MapZoning:
def __init__(self, city, distance=None, map_service=None):
self.city = city
self.distance = distance if distance else 1 # 单位km，默认0.5km
self.earth_radius = 40075.04 / (2 * math.pi)
self.lat_unit = self.get_lat_unit(self.distance)
self.map_service = map_service if map_service else Map(os.environ.get('MAP_KEY'))
self.location_validated = {} # key：[经度,维度] value: 1验证成功，2超出范围
def get_lat_unit(self, distance):
"""
获取一定距离的维度单位
:param distance:
:return:
"""
return distance * 360 / (2 * math.pi) / self.earth_radius
def get_lng_unit(self, lat, distance):
"""
获取一定距离的经度单位
:param lat:
:param distance:
:return:
"""
return distance * 360 / (2 * math.pi) / math.cos(math.radians(lat)) / self.earth_radius
def _get_area_points(self, quadrant, location):
"""
获取区域的点
:param int quadrant:
:param list location:
:return list : 矩形四个点，顺序：左上、右上、右下、左下
"""
if quadrant == 1:
left_lng = float(location[0]) # 经度
bottom_lat = float(location[1]) # 维度
top_lat = round(bottom_lat + self.lat_unit, 6)
right_lng = round(left_lng + self.get_lng_unit(float(top_lat), self.distance), 6)
elif quadrant == 2:
right_lng = float(location[0])
bottom_lat = float(location[1])
top_lat = round(bottom_lat + self.get_lat_unit(self.distance), 6)
left_lng = round(right_lng - self.get_lng_unit(top_lat, self.distance), 6)
elif quadrant == 3:
right_lng = float(location[0])
top_lat = float(location[1])
left_lng = round(right_lng - self.get_lng_unit(top_lat, self.distance), 6)
bottom_lat = round(top_lat - self.get_lat_unit(self.distance), 6)
elif quadrant == 4:
left_lng = float(location[0])
top_lat = float(location[1])
right_lng = round(left_lng + self.get_lng_unit(float(top_lat), self.distance), 6)
bottom_lat = round(top_lat - self.lat_unit, 6)
else:
return None
return [[left_lng, top_lat], [right_lng, top_lat], [right_lng, bottom_lat], [left_lng, bottom_lat]]
def get_area_points_x(self, quadrant, res_points, points, i_temp):
if i_temp > 60:
return
if quadrant == 1:
location = points[2]
elif quadrant == 2:
location = points[3]
elif quadrant == 3:
location = points[0]
elif quadrant == 4:
location = points[1]
else:
return
tmp_pos = self._get_area_points(quadrant, location)
if self.validate_points(tmp_pos):
res_points.append(tmp_pos)
i_temp += 1
self.get_area_points_x(quadrant, res_points, tmp_pos, i_temp)
def get_area_points_y(self, quadrant, res_points, points, i_temp):
if i_temp > 60:
return
if quadrant == 1:
location = points[0]
elif quadrant == 2:
location = points[1]
elif quadrant == 3:
location = points[2]
elif quadrant == 4:
location = points[3]
else:
return
tmp_pos = self._get_area_points(quadrant, location)
if self.validate_points(tmp_pos):
res_points.append(tmp_pos)
i_temp += 1
self.get_area_points_x(quadrant, res_points, tmp_pos, 1)
self.get_area_points_y(quadrant, res_points, tmp_pos, i_temp)
def get_all_area_points(self, quadrant, res_points, location):
tmp_pos = self._get_area_points(quadrant, location)
res_points.append(tmp_pos)
self.get_area_points_x(quadrant, res_points, tmp_pos, 1)
self.get_area_points_y(quadrant, res_points, tmp_pos, 1)
def validate_points(self, points):
b_result = False
for location in points:
if location[0] < 113.74181 or location[1] < 22.438768 or location[1] > 22.887492 or location[0] > 114.649706:
continue
v_key = ','.join([str(x) for x in location])
v_value = self.location_validated.get(v_key)
if not v_value:
res = self.map_service.get_get_address(location)
address_component = res.get('regeocode').get('addressComponent') if res and res.get('regeocode') else {}
if not address_component or address_component.get('citycode') != '0755' or address_component.get(
'seaArea'):
self.location_validated[v_key] = 2
else:
b_result = True
self.location_validated[v_key] = 1
elif v_value == 1:
b_result = True
return b_result
def zoning(self):
print(time.time())
res_points = []
location = [114.057868, 22.543099] # self.map_service.get_geo_code(self.city).split(',')
self.get_all_area_points(1, res_points, location)
self.get_all_area_points(2, res_points, location)
self.get_all_area_points(3, res_points, location)
self.get_all_area_points(4, res_points, location)
with open('tmp/points.js', 'w', encoding='utf-8') as _file:
_file.write('var points = ' + json.dumps(res_points) + ';')
print(time.time())
# points = [['113.980092', '22.475719'], ['113.984953', '22.475719'], ['113.984953', '22.471227'],
# ['113.980092', '22.471227']]
# print(self.validate_points(points))
# print(self.map_service.get_get_address(['113.980092', '22.475719']))
return res_points
amap.py
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
import logging
import requests
class Map:
def __init__(self, key):
if not key:
raise Exception('初始map服务错误无效Key')
self.key = key
def get_geo_code(self, address):
"""
获取地理位置编码
:param address:
:return:
"""
api_geocode = 'http://restapi.amap.com/v3/geocode/geo'
try:
response = requests.get(api_geocode, {'key': self.key,
'address': address,
'output': 'JSON'})
if response.status_code != 200:
return None
res = response.json()
if not res:
return None
geocodes = res.get('geocodes')
if not geocodes:
return None
location = geocodes[0].get('location')
if not location:
return None
return location
except Exception as ex:
logging.error('获取地址位置编码错误：' + str(ex))
return None
def get_get_address(self, location):
api = 'http://restapi.amap.com/v3/geocode/regeo'
try:
response = requests.get(api, {'key': self.key,
'location': ','.join([str(x) for x in location]),
'output': 'JSON'})
if response.status_code != 200:
return None
res = response.json()
if not res:
return None
return res
except Exception as ex:
logging.error('获取地址位置错误：' + str(ex))
return None
app.py
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
import logging
from os.path import join, dirname
from dotenv import load_dotenv
from map_zoning import MapZoning
dotenv_path = join(dirname(__file__), '.env')
load_dotenv(dotenv_path)
if __name__ == '__main__':
try:
MapZoning('深圳市').zoning()
except Exception as ex:
logging.error(ex)
.env
# application APP_NAME="split_map" APP_RUNTIME="dev" APP_LOG_LEVEL="INFO" APP_LOG_HANDLER="console" # map MAP_KEY="amap_web_api_key"
