-
Notifications
You must be signed in to change notification settings - Fork 25
/
get_grid.py
98 lines (77 loc) · 3.59 KB
/
get_grid.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
from get_distance import get_distance
import transCoordinateSystem
import settings
#TODO
require_distance = 1 #要求的地图网格间隔距离(千米为单位) -- 百度地图
print(settings.city_bound_point_A)
#武汉
#四个点分别代表矩阵的四个角的坐标
#point_A = [114.286652,30.642638] #左上角点，x619
#point_D = [114.239273,30.580588] #左下角点，农场十一队
#point_B = [114.462433,30.574677] #右上角点，木妙
#point_C = [114.418488,30.479746] #右下角点，牛场右下角
point_A = settings.city_bound_point_A #左上角点，x619
point_D = settings.city_bound_point_D #左下角点，农场十一队
point_B = settings.city_bound_point_B #右上角点，木妙
point_C = settings.city_bound_point_C #右下角点，牛场右下角
#AB间的距离
distance_AB = get_distance(point_A[1], point_A[0], point_B[1], point_B[0])
print(distance_AB) #13.381562710218205
#AD间的距离
distance_AD = get_distance(point_A[1], point_A[0], point_D[1], point_D[0])
print(distance_AD) #7.735240077970707
#AB边上平分点的个数
num_split_AB = int(distance_AB / require_distance if distance_AB % require_distance == 0 else distance_AB // require_distance + 1) - 1
print('AB边上平分点的个数', str(num_split_AB))
#AD边上评分的个数
num_split_AD = int(distance_AD / require_distance if distance_AD % require_distance == 0 else distance_AD // require_distance + 1) - 1
print('AD边上评分的个数', str(num_split_AD))
#计算AB线段上各个点的坐标
points_coordinate_AB = []
for i in range(num_split_AB + 1):
x = point_A[0] + (point_B[0] - point_A[0]) / (num_split_AB + 1) * i
y = point_A[1] + (point_B[1] - point_A[1]) / (num_split_AB + 1) * i
points_coordinate_AB.append([x, y])
points_coordinate_AB.append(point_B)#加上B点的坐标
print('AB线段上各个平分点的坐标集合：', points_coordinate_AB)
#计算DC线段上各个点的坐标
points_coordinate_DC = []
for i in range(num_split_AB + 1):
x = point_D[0] + (point_C[0] - point_D[0]) / (num_split_AB + 1) * i
y = point_D[1] + (point_C[1] - point_D[1]) / (num_split_AB + 1) * i
points_coordinate_DC.append([x, y])
points_coordinate_DC.append(point_C)#加上B点的坐标
print('DC线段上各个平分点的坐标集合：', points_coordinate_DC)
#计算中间网格交汇点坐标
points_coordinate_result = []
for i in range(num_split_AD + 2):
for j in range(num_split_AB + 2):
x = points_coordinate_AB[j][0] + (points_coordinate_DC[j][0] - points_coordinate_AB[j][0]) / (num_split_AD + 1) * (i + 1)
y = points_coordinate_AB[j][1] + (points_coordinate_DC[j][1] - points_coordinate_AB[j][1]) / (num_split_AD + 1) * (i + 1)
points_coordinate_result.append([x, y])
points_coordinate_result.extend(points_coordinate_AB)
points_coordinate_result.extend(points_coordinate_DC)
print('======================================')
print(points_coordinate_result)
'''
latA = 116.371067
lonA = 39.959193
latB = 116.483032
lonB = 39.973018
d2 = get_distance(lonA, latA, lonB, latB)
print(d2) #9.664978035037414
'''
def create_yichuxing_data(points_coordinate_result):
'''生成宜出行爬取的数据各式'''
text = 'OBJECTID,x,y\n'
for i, point in enumerate(points_coordinate_result):
lng = points_coordinate_result[i][0]
lat = points_coordinate_result[i][1]
res = transCoordinateSystem.bd09_to_wgs84(lng, lat)
lng = res[0]
lat = res[1]
text += str(i) + ',' + str(lng) + ',' + str(lat) + '\n'
with open('data.txt', 'w') as f:
f.write(text.strip())
print('写入成功')
create_yichuxing_data(points_coordinate_result)