-
Notifications
You must be signed in to change notification settings - Fork 12
/
makeLawrenceMap.py
271 lines (203 loc) · 10.1 KB
/
makeLawrenceMap.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
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
""" ***************************************************************************
# * File Description: *
# * Makes a street map of Lawrence, Kansas where the streets are colored by *
# * their length. *
# * *
# * The contents of this script are: *
# * 1. Importing Libraries *
# * 2. Version Check *
# * 3. Get Data *
# * 4. Unpack Data *
# * 5. Assign Each Segment a Color Based on its Length *
# * 6. Assign Each Segment a Width Based on its type *
# * 7. Make Map *
# * 8. Helper Functions: Add Border to the Map *
# * 9. Adding Border to Map *
# * 10. Adding Text to Map *
# * *
# * --------------------------------------------------------------------------*
# * AUTHORS(S): Frank Ceballos <[email protected]> *
# * --------------------------------------------------------------------------*
# * DATE CREATED: September 30, 2019 *
# * --------------------------------------------------------------------------*
# * NOTES: None *
# * ************************************************************************"""
###############################################################################
# 1. Importing Libraries #
###############################################################################
import networkx as nx
import osmnx as ox
import requests
import matplotlib.cm as cm
import matplotlib.colors as colors
from matplotlib.lines import Line2D
# To add text and a border to the map
import PIL # You dont really need this import
from PIL import Image, ImageOps, ImageColor, ImageFont, ImageDraw
# Remove if you're working in a Jupyter Notebook
%matplotlib auto
###############################################################################
# 2. Version Check #
###############################################################################
print(f"The NetworkX package is version {nx.__version__}")
print(f"The OSMNX package is version {ox.__version__}")
print(f"The Request package is version {requests.__version__}")
print(f"The PIL package is version {PIL.__version__}")
###############################################################################
# 3. Get Data #
###############################################################################
# Define city/cities
places = ["Lawrence, Kansas, USA"]
# Get data for places
G = ox.graph_from_place(places, network_type="all", simplify = True)
###############################################################################
# 4. Unpack Data #
###############################################################################
u = []
v = []
key = []
data = []
for uu, vv, kkey, ddata in G.edges(keys=True, data=True):
u.append(uu)
v.append(vv)
key.append(kkey)
data.append(ddata)
###############################################################################
# 5. Assign Each Segment a Color Based on its Length #
###############################################################################
# List to store colors
roadColors = []
# The length is in meters
for item in data:
if "length" in item.keys():
if item["length"] <= 100:
color = "#d40a47"
elif item["length"] > 100 and item["length"] <= 200:
color = "#e78119"
elif item["length"] > 200 and item["length"] <= 400:
color = "#30bab0"
elif item["length"] > 400 and item["length"] <= 800:
color = "#bbbbbb"
else:
color = "w"
roadColors.append(color)
###############################################################################
# 6. Assign Each Segment a Width Based on its type #
###############################################################################
# List to store linewidths
roadWidths = []
for item in data:
if "footway" in item["highway"]:
linewidth = 1
else:
linewidth = 2.5
roadWidths.append(linewidth)
###############################################################################
# 7. Make Map #
###############################################################################
# Center of map
latitude = 38.960323
longitude = -95.263223
# Bbox sides
north = latitude + 0.035
south = latitude - 0.035
east = longitude + 0.05
west = longitude - 0.05
# Make Map
fig, ax = ox.plot_graph(G, node_size=0, bbox = (north, south, east, west), margin = 0,
fig_height=40, fig_width=40, dpi = 300, bgcolor = "#061529",
save = False, edge_color=roadColors,
edge_linewidth=roadWidths, edge_alpha=1)
# Text and marker size
markersize = 16
fontsize = 16
# Add legend
legend_elements = [Line2D([0], [0], marker='s', color="#061529", label= 'Length < 100 m',
markerfacecolor="#d40a47", markersize=markersize),
Line2D([0], [0], marker='s', color="#061529", label= 'Length between 100-200 m',
markerfacecolor="#e78119", markersize=markersize),
Line2D([0], [0], marker='s', color="#061529", label= 'Length between 200-400 m',
markerfacecolor="#30bab0", markersize=markersize),
Line2D([0], [0], marker='s', color="#061529", label= 'Length between 400-800 m',
markerfacecolor="#bbbbbb", markersize=markersize),
Line2D([0], [0], marker='s', color="#061529", label= 'Length > 800 m',
markerfacecolor="w", markersize=markersize)]
l = ax.legend(handles=legend_elements, bbox_to_anchor=(0.0, 0.0), frameon=True, ncol=1,
facecolor = '#061529', framealpha = 0.9,
loc='lower left', fontsize = fontsize, prop={'family':"Georgia", 'size':fontsize})
# Legend font color
for text in l.get_texts():
text.set_color("w")
# Save figure
fig.savefig("Lawrence.png", dpi=300, bbox_inches='tight', format="png", facecolor=fig.get_facecolor(), transparent=True)
###############################################################################
# 8. Helper Functions: Add Border to the Map #
###############################################################################
# Get color
def _color(color, mode):
color = ImageColor.getcolor(color, mode)
return color
# Expand image
def expand(image, fill = '#e0474c', bottom = 50, left = None, right = None, top = None):
"""
Expands image
Parameters
----------
image: The image to expand.
bottom, left, right, top: Border width, in pixels.
param fill: Pixel fill value (a color value). Default is 0 (black).
return: An image.
"""
if left == None:
left = 0
if right == None:
right = 0
if top == None:
top = 0
width = left + image.size[0] + right
height = top + image.size[1] + bottom
out = Image.new(image.mode, (width, height), _color(fill, image.mode))
out.paste(image, (left, top))
return out
# Add border
def add_border(input_image, output_image, fill = '#e0474c', bottom = 50, left = None, right = None, top = None):
""" Adds border to image and saves it.
Parameters
----------
input_image: str,
String object for the image you want to load. This is the name of the file you want to read.
output_image: str,
String object for the output image name. This is the name of the file you want to export.
fill: str,
Hex code for border color. Default is set to reddish.
bottom, left, right, top: int,
Integer object specifying the border with in pixels.
"""
if left == None:
left = 0
if right == None:
right = 0
if top == None:
top = 0
img = Image.open(input_image)
bimg = expand(img, bottom = bottom, left = left, right = right, top = top, fill= '#e0474c')
bimg.save(output_image)
###############################################################################
# 9. Adding Border to Map #
###############################################################################
# Input image
in_img = 'Lawrence.png'
# Output Image
add_border(in_img, output_image='Lawrence.png', fill = '#e0474c', bottom = 800)
###############################################################################
# 10. Adding Text to Map #
###############################################################################
# Open Image
img = Image.open("Lawrence.png")
draw = ImageDraw.Draw(img)
# Get font from working directory. Visit https://www.wfonts.com/search?kwd=pmingliu to download fonts
font = ImageFont.truetype("PMINGLIU.ttf", 650)
# Add font: position, text, color, font
draw.text((1650,8500),"LAWRENCE, KANSAS", (255,255,255), font=font)
# Save image
img.save('Lawrence.png')