-
Notifications
You must be signed in to change notification settings - Fork 53
Open
Description
This code is really cool! Here's Bozeman, with love from the MSU physics department:
rm = RidgeMap((-111.373901,45.502497,-110.916595,45.908167))
msu = (-45.46202, -41.665885)
msu_coords = ((msu[0] - rm.longs[0])/(rm.longs[1] - rm.longs[0]),(msu[1] - rm.lats[0])/(rm.lats[1] - rm.lats[0]))
values = rm.get_elevation_data(num_lines=150)
newmap = colors.LinearSegmentedColormap.from_list('trunc({n},{a:.2f},{b:.2f})'.format(n='Blues', a=0.7, b=1), plt.get_cmap('Blues')(np.linspace(0.7, 1, 100)))
values=rm.preprocess(
values=values,
lake_flatness=1,
water_ntile=5,
vertical_ratio=180)
ridges = rm.plot_map(values=values,
label='Bozeman',
label_y=0.9,
label_x=0.05,
label_size=40,
line_color=newmap,
kind = 'elevation',
linewidth=1)
ridges.scatter(msu_coords[0], msu_coords[1], marker='$\heartsuit$',
color='red',
s=80,
zorder=len(values)+10)
plt.show()```
Metadata
Metadata
Assignees
Labels
No labels