-
Notifications
You must be signed in to change notification settings - Fork 0
/
program_specialty_radar_charts_2012-2019.py
158 lines (133 loc) · 6.03 KB
/
program_specialty_radar_charts_2012-2019.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
# Author: Eni Awowale
# Date first written: June 20, 2019
# Date last updated: August 14, 2019
# Summary: Create Radar Charts of University Geography Departments Program Specialities
"""
Purpose:
Using a Radar chart (spider or star chart) visualize university program specialties.
Each year since 2012 the AAG (American Association of Geographers)
publishes a Guide to Geography Programs in the Americas, Program Specialties section.
The data used is from Program Specialties published from 2012-2019, excluding 2013,
(the Guide was not published in 2013).
"""
import numpy as np
import os
import sys
import ast
import matplotlib.pyplot as plt
from matplotlib import rcParams
from matplotlib.patches import Circle, RegularPolygon
from matplotlib.path import Path
from matplotlib.projections.polar import PolarAxes
from matplotlib.projections import register_projection
from matplotlib.spines import Spine
from matplotlib.transforms import Affine2D
#Set font
rcParams['font.family'] = 'DejaVu Sans'
folder = r'C:\Users\oawowale\Documents\GitHub\affinities-of-geog-departments'
os.chdir(folder)
def radar_factory(num_vars, frame='circle'):
"""Create a radar chart with `num_vars` axes.
This function creates a RadarAxes projection and registers it.
Parameters
----------
num_vars : int
Number of variables for radar chart.
frame : {'circle' | 'polygon'}
Shape of frame surrounding axes.
"""
# calculate evenly-spaced axis angles
theta = np.linspace(0, 2*np.pi, num_vars, endpoint=False)
class RadarAxes(PolarAxes):
name = 'radar'
# use 1 line segment to connect specified points
RESOLUTION = 1
def __init__(self, *args, **kwargs):
super().__init__(*args, **kwargs)
# rotate plot such that the first axis is at the top
self.set_theta_zero_location('N')
def fill(self, *args, closed=True, **kwargs):
"""Override fill so that line is closed by default"""
return super().fill(closed=closed, *args, **kwargs)
def plot(self, *args, **kwargs):
"""Override plot so that line is closed by default"""
lines = super().plot(*args, **kwargs)
for line in lines:
self._close_line(line)
def _close_line(self, line):
x, y = line.get_data()
# FIXME: markers at x[0], y[0] get doubled-up
if x[0] != x[-1]:
x = np.concatenate((x, [x[0]]))
y = np.concatenate((y, [y[0]]))
line.set_data(x, y)
def set_varlabels(self, labels):
self.set_thetagrids(np.degrees(theta), labels)
def _gen_axes_patch(self):
# The Axes patch must be centered at (0.5, 0.5) and of radius 0.5
# in axes coordinates.
if frame == 'circle':
return Circle((0.5, 0.5), 0.5)
elif frame == 'polygon':
return RegularPolygon((0.5, 0.5), num_vars,
radius=.5, edgecolor="k")
else:
raise ValueError("unknown value for 'frame': %s" % frame)
def _gen_axes_spines(self):
if frame == 'circle':
return super()._gen_axes_spines()
elif frame == 'polygon':
# spine_type must be 'left'/'right'/'top'/'bottom'/'circle'.
spine = Spine(axes=self,
spine_type='circle',
path=Path.unit_regular_polygon(num_vars))
# unit_regular_polygon gives a polygon of radius 1 centered at
# (0, 0) but we want a polygon of radius 0.5 centered at (0.5,
# 0.5) in axes coordinates.
spine.set_transform(Affine2D().scale(.5).translate(.5, .5)
+ self.transAxes)
return {'polar': spine}
else:
raise ValueError("unknown value for 'frame': %s" % frame)
register_projection(RadarAxes)
return theta
#University program specialties data
data_for_reading = open('final_radar_chart_data.txt').readline()
data_list = ast.literal_eval(data_for_reading)
'''Actual List:
['Urban and Economic Geography', 'Physical Geography', 'Methods', 'Human-Environmental Interactions', 'Geospatial Technologies', 'Human Geography']
'''
if __name__ == '__main__':
N = 6
theta = radar_factory(N, frame='polygon')
for data_index in range(1, len(data_list)):
data = data_list[data_index]
spoke_labels = data_list[0]
#Formatting tick names
spoke_labels[0] = '\n\nUrban and Economic\nGeography'
spoke_labels[1] = 'Physical\nGeography'
spoke_labels[2] = 'Methods'
spoke_labels[3] = 'Human-Environmental\nInteractions'
spoke_labels[4] = 'Geospatial\nTechnologies'
spoke_labels[5] = 'Human\nGeography'
ax = plt.subplot(111, projection='radar')
#Tick parameters
ax.tick_params(axis='x', which='major', grid_alpha=0, direction='out', labelsize='small', pad=12)
#Colors and line thickness
linewidths = [1, 1.5, 2, 2.5, 3, 3.5, 4][::-1]
colors = ['b', 'r', 'g', 'm', 'y', 'c', '#fc9403'][::-1]
ax.set_rgrids([0.2, 0.4, 0.6, 0.8], size='small')
ax.set_title(data[0].title(), weight='bold', size='medium', position=(0.5, 1.1),
horizontalalignment='center', verticalalignment='center')
for d, color, widths in zip(data[1], colors, linewidths):
ax.plot(theta, d, color=color, linewidth=widths, linestyle=':')
ax.fill(theta, d, facecolor='#02113A', alpha=0.25)
ax.set_varlabels(spoke_labels)
#Add legend
labels = ('2012', '2014', '2015', '2016', '2017', '2018', '2019')
legend = ax.legend(labels, loc=(1.12, 0),
labelspacing=0.1, fontsize='small')
#Save current figure
plt.savefig(r'radar_charts\radar_chart_{0}.png'.format(data[0]))
#Clear plot for new data
plt.clf()