Radial Profile

TODO: move this entirely to a colab notebook, and only keep some key figures on this page. — also add more code to get radial averages (maybe even with errors), not just the coordinates.

Here's a sample plot of output ring coordinates around Newton crater, where each ring is a different color and we dramatically undersample the points on each ring for clarity (20 rings, 10km between points on each ring):

The code to generate this plot is as follows:

 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
import numpy as np
import matplotlib.pyplot as plt
import redplanet as rp

## get crater data
crater = rp.Craters.get(name='Newton', as_dict=True)[0]

## get concentric ring coordinates
ring_radius_km__per_ring, ring_coords__per_ring = rp.analysis.radial_profile.get_concentric_ring_coords(
    lon = crater['lon'],
    lat = crater['lat'],
    radius_km = (crater['diam'] / 2) * 2.5,
    num_rings = 20,
    dist_btwn_points_km = 10,
)

## plot hillshade (base)
r = 10
lons = np.linspace(crater['lon']-r, crater['lon']+r, 1000)
lats = np.linspace(crater['lat']-r, crater['lat']+r, 1000)
rp.Crust.topo.load('DEM_463m')

fig, ax = rp.plot(
    lons, lats,
    dat = rp.Crust.topo.get(lons, lats),
    figsize = (7,7),
    hillshade = True,
    cbar = False,
    alpha_dat = 0,
    alpha_hs = 0.5,
    show = False,
)

## plot concentric rings
colors__per_ring = plt.cm.jet(
    np.linspace(0, 1, len(ring_radius_km__per_ring))
)

for i in range(len(ring_radius_km__per_ring)):
    ax.scatter(
        ring_coords__per_ring[i][:,0],
        ring_coords__per_ring[i][:,1],
        color = colors__per_ring[i],
        s = 1,
        alpha = 0.5,
    )

ax.set_title(f'Concentric rings around {crater["name"]} crater')

plt.show()