There's a new cool package out there: ridge_map
. Let's see what it can do!
%matplotlib inline
from ridge_map import RidgeMap
RidgeMap().plot_map()
/Users/kappamaki/anaconda/lib/python3.6/site-packages/skimage/util/dtype.py:135: UserWarning: Possible precision loss when converting from float64 to uint8 .format(dtypeobj_in, dtypeobj_out))
<matplotlib.axes._subplots.AxesSubplot at 0x10b1ac198>
Let's use http://bboxfinder.com as suggested in the README.
import matplotlib.pyplot as plt
rm = RidgeMap((6.996918,43.628123,7.358437,43.734639))
values = rm.get_elevation_data(num_lines=100)
values[values == 0] = 'nan'
values=rm.preprocess(
values=values,
lake_flatness=0,
water_ntile=0,
vertical_ratio=340)
rm.plot_map(values=values,
label='Cagnes-sur-Mer',
label_y=0.1,
label_x=0.55,
label_size=50,
linewidth=.5,
kind='gradient',
line_color='red')
/Users/kappamaki/anaconda/lib/python3.6/site-packages/skimage/util/dtype.py:135: UserWarning: Possible precision loss when converting from float64 to uint8 .format(dtypeobj_in, dtypeobj_out))
<matplotlib.axes._subplots.AxesSubplot at 0x11f762d30>
rm = RidgeMap((5.286140,43.420759,5.725594,43.604262))
values = rm.get_elevation_data(num_lines=100)
values=rm.preprocess(
values=values,
lake_flatness=0,
water_ntile=0,
vertical_ratio=200)
rm.plot_map(values=values,
label='Aix-en-Provence',
label_y=0.25,
label_x=0.55,
label_size=50,
linewidth=.5,
kind='gradient',
line_color='green')
/Users/kappamaki/anaconda/lib/python3.6/site-packages/skimage/util/dtype.py:135: UserWarning: Possible precision loss when converting from float64 to uint8 .format(dtypeobj_in, dtypeobj_out))
<matplotlib.axes._subplots.AxesSubplot at 0x120681f60>
rm = RidgeMap((1.856689,48.714525,2.796021,49.039668))
values = rm.get_elevation_data(num_lines=100)
values=rm.preprocess(
values=values,
lake_flatness=2,
water_ntile=10,
vertical_ratio=100)
rm.plot_map(values=values,
label='Paris',
label_y=0.25,
label_x=0.55,
label_size=50,
linewidth=.5,
kind='gradient',
line_color='black')
/Users/kappamaki/anaconda/lib/python3.6/site-packages/skimage/util/dtype.py:135: UserWarning: Possible precision loss when converting from float64 to uint8 .format(dtypeobj_in, dtypeobj_out))
<matplotlib.axes._subplots.AxesSubplot at 0x11f6b10b8>
What about a better zoom?
rm = RidgeMap((2.132378,48.793973,2.571831,48.960384))
values = rm.get_elevation_data(num_lines=150)
values=rm.preprocess(
values=values,
lake_flatness=2,
water_ntile=10,
vertical_ratio=150)
rm.plot_map(values=values,
label='Paris',
label_y=0.25,
label_x=0.55,
label_size=50,
linewidth=.5,
kind='gradient',
line_color='black')
/Users/kappamaki/anaconda/lib/python3.6/site-packages/skimage/util/dtype.py:135: UserWarning: Possible precision loss when converting from float64 to uint8 .format(dtypeobj_in, dtypeobj_out))
<matplotlib.axes._subplots.AxesSubplot at 0x117f21128>
rm = RidgeMap((-5.548096,46.430285,-0.417480,48.936935))
values = rm.get_elevation_data(num_lines=150)
values[values == 0] = 'nan'
values=rm.preprocess(
values=values,
lake_flatness=2,
water_ntile=10,
vertical_ratio=100)
rm.plot_map(values=values,
label='Bretagne',
label_y=0.25,
label_x=0.05,
label_size=100,
linewidth=1,
kind='gradient',
line_color='black')
/Users/kappamaki/anaconda/lib/python3.6/site-packages/skimage/util/dtype.py:135: UserWarning: Possible precision loss when converting from float64 to uint8 .format(dtypeobj_in, dtypeobj_out))
<matplotlib.axes._subplots.AxesSubplot at 0x117db1198>
rm = RidgeMap((7.932129,41.310824,9.931641,43.141078))
values = rm.get_elevation_data(num_lines=200)
values[values == 0] = 'nan'
values=rm.preprocess(
values=values,
lake_flatness=2,
water_ntile=5,
vertical_ratio=100)
rm.plot_map(values=values,
label='La Corse',
label_y=0.85,
label_x=0.05,
label_size=100,
linewidth=1,
line_color=plt.get_cmap('ocean'),
kind='elevation')
/Users/kappamaki/anaconda/lib/python3.6/site-packages/skimage/util/dtype.py:135: UserWarning: Possible precision loss when converting from float64 to uint8 .format(dtypeobj_in, dtypeobj_out))
<matplotlib.axes._subplots.AxesSubplot at 0x11c6943c8>
/Users/kappamaki/anaconda/lib/python3.6/site-packages/matplotlib/colors.py:512: RuntimeWarning: invalid value encountered in less xa[xa < 0] = -1
rm = RidgeMap((1.669922,44.629573,3.002014,45.176229))
values = rm.get_elevation_data(num_lines=100)
values=rm.preprocess(
values=values,
lake_flatness=1,
water_ntile=5,
vertical_ratio=100)
rm.plot_map(values=values,
label='Aurillac',
label_y=0.10,
label_x=0.72,
label_size=50,
linewidth=1,
kind='elevation',
line_color=plt.get_cmap('magma'))
/Users/kappamaki/anaconda/lib/python3.6/site-packages/skimage/util/dtype.py:135: UserWarning: Possible precision loss when converting from float64 to uint8 .format(dtypeobj_in, dtypeobj_out))
<matplotlib.axes._subplots.AxesSubplot at 0x11f7f1940>
/Users/kappamaki/anaconda/lib/python3.6/site-packages/matplotlib/colors.py:512: RuntimeWarning: invalid value encountered in less xa[xa < 0] = -1
from skimage.filters import gaussian
rm = RidgeMap((9.757490,50.557888,9.867353,50.598058))
values = rm.get_elevation_data(num_lines=75)
values=rm.preprocess(
values=values,
lake_flatness=0,
water_ntile=0,
vertical_ratio=200)
rm.plot_map(values=values,
label='Niederbieber',
label_y=0.08,
label_x=0.60,
label_size=50,
linewidth=1,
kind='gradient',
line_color='black')
/Users/kappamaki/anaconda/lib/python3.6/site-packages/skimage/util/dtype.py:135: UserWarning: Possible precision loss when converting from float64 to uint8 .format(dtypeobj_in, dtypeobj_out))
<matplotlib.axes._subplots.AxesSubplot at 0x138a27cc0>
rm = RidgeMap((9.757490,50.557888,9.867353,50.598058))
values = rm.get_elevation_data(num_lines=75)
values = gaussian(values, sigma=(0, 1))
values=rm.preprocess(
values=values,
lake_flatness=0,
water_ntile=0,
vertical_ratio=200)
rm.plot_map(values=values,
label='Niederbieber',
label_y=0.08,
label_x=0.60,
label_size=50,
linewidth=1,
kind='gradient',
line_color='black')
/Users/kappamaki/anaconda/lib/python3.6/site-packages/skimage/util/dtype.py:135: UserWarning: Possible precision loss when converting from float64 to uint8 .format(dtypeobj_in, dtypeobj_out))
<matplotlib.axes._subplots.AxesSubplot at 0x13720a8d0>
rm = RidgeMap((9.835510,50.521394,9.945374,50.561596))
values = rm.get_elevation_data(num_lines=50)
values = gaussian(values, sigma=(0, .75))
values=rm.preprocess(
values=values,
lake_flatness=0,
water_ntile=0,
vertical_ratio=100)
rm.plot_map(values=values,
label='die Milseburg',
label_y=0.08,
label_x=0.60,
label_size=50,
linewidth=1,
line_color=plt.get_cmap('inferno'),
kind='elevation')
/Users/kappamaki/anaconda/lib/python3.6/site-packages/skimage/util/dtype.py:135: UserWarning: Possible precision loss when converting from float64 to uint8 .format(dtypeobj_in, dtypeobj_out))
<matplotlib.axes._subplots.AxesSubplot at 0x1377c2a90>
rm = RidgeMap((9.461975,50.456630,9.901428,50.617449))
values = rm.get_elevation_data(num_lines=100)
values=rm.preprocess(
values=values,
lake_flatness=0,
water_ntile=0,
vertical_ratio=200)
rm.plot_map(values=values,
label='Fulda',
label_y=0.78,
label_x=0.10,
label_size=100,
linewidth=1,
line_color=plt.get_cmap('hot'),
kind='elevation')
/Users/kappamaki/anaconda/lib/python3.6/site-packages/skimage/util/dtype.py:135: UserWarning: Possible precision loss when converting from float64 to uint8 .format(dtypeobj_in, dtypeobj_out))
<matplotlib.axes._subplots.AxesSubplot at 0x12e80a5f8>