Example maps¶

This notebook shows all the built-in example maps in maria.

[1]:

import maria

crab_nebula = maria.map.get("maps/M1.h5")
crab_nebula.plot(slices="all")
print(crab_nebula)

ProjectionMap:
  data(1, 3, 1205, 1187):
    min: -6.076e+06
    max: 9.216e+05
    units: mJy beam^-1
    quantity: spectral_flux_density_per_beam
  stokes(1):
    components: ['I']
  nu(3):
    values: [149.8962 214.1375 272.5386] GHz
  eta(1205):
    height: 20.07’
    res: -1”
  xi(1187):
    width: 19.77’
    res: 1”
  frame: ra/dec
  center:
    ra: 05ʰ34ᵐ31.95ˢ
    dec: 22°00’52.16”
  beam(maj, min, psi): ragged
  memory: 34.33 MB

../../../_images/usage_inputs_maps_example-maps_1_1.png

[2]:

CO_spectral_line = maria.map.get("maps/12CO(2-1).fits")
CO_spectral_line.plot(slices={"v": [9, 10, 11]})
print(CO_spectral_line)

2026-06-02 15:58:14.169 INFO: Fetching https://github.com/thomaswmorris/maria-data/raw/master/maps/12CO(2-1).fits
Downloading: 100%|██████████| 14.1M/14.1M [00:00<00:00, 141MB/s]

ProjectionMap:
  data(1, 20, 357, 493):
    min: -5.769e+00
    max: 5.732e+00
    units: K_RJ
    quantity: rayleigh_jeans_temperature
  nu(1):
    values: [230.538] GHz
  v(20):
    min: -4.463 km/s
    max: 1.714 km/s
    res: -325.1 m/s
  eta(357):
    height: 1.618°
    res: -16.36”
  xi(493):
    width: 2.236°
    res: 16.36”
  frame: galactic
  center:
    glon: 124.92499999999998
    glat: 26.26333
  beam(maj, min, psi): (34.49”, 34.49”, 0 rad)
  memory: 14.08 MB

../../../_images/usage_inputs_maps_example-maps_2_2.png

[3]:

time_evolving_sun = maria.map.get("maps/sun.h5")
time_evolving_sun.plot(slices={"t": [0, -1]})
print(time_evolving_sun)

time_evolving_sun_change = time_evolving_sun[:, :, -1]
time_evolving_sun_change.data -= time_evolving_sun[:, :, 0].data
time_evolving_sun_change.plot()

2026-06-02 15:58:17.219 INFO: Fetching https://github.com/thomaswmorris/maria-data/raw/master/maps/sun.h5
Downloading: 100%|██████████| 1.67M/1.67M [00:00<00:00, 61.0MB/s]

ProjectionMap:
  data(1, 1, 16, 256, 256):
    min: 0.000e+00
    max: 1.339e+04
    units: K_RJ
    quantity: rayleigh_jeans_temperature
  stokes(1):
    components: ['I']
  nu(1):
    values: [950.] GHz
  t(16):
    min: 1970-01-01 00:00:00.000+00:00
    max: 1970-01-01 00:03:00.000+00:00
    res: 12 s
  eta(256):
    height: 40.81’
    res: -9.603”
  xi(256):
    width: 40.81’
    res: 9.603”
  frame: ra/dec
  center:
    ra: 00ʰ00ᵐ0.00ˢ
    dec: 00°00’0.00”
  beam(maj, min, psi): (0 rad, 0 rad, 0 rad)
  memory: 4.194 MB

../../../_images/usage_inputs_maps_example-maps_3_2.png

../../../_images/usage_inputs_maps_example-maps_3_3.png

[4]:

whirlpool_galaxy = maria.map.get("maps/M51HA.fits")
whirlpool_galaxy.plot(slices="all")
print(whirlpool_galaxy)

2026-06-02 15:58:19.506 INFO: Fetching https://github.com/thomaswmorris/maria-data/raw/master/maps/M51HA.fits
Downloading: 100%|██████████| 2.10M/2.10M [00:00<00:00, 96.2MB/s]

ProjectionMap:
  data(1, 512, 512):
    min: 0.000e+00
    max: 1.361e+01
    units: Jy beam^-1
    quantity: spectral_flux_density_per_beam
  nu(1):
    values: [100.] GHz
  eta(512):
    height: 4.258’
    res: -500 marcsec
  xi(512):
    width: 4.258’
    res: 500 marcsec
  frame: ra/dec
  center:
    ra: 13ʰ29ᵐ5.31ˢ
    dec: 47°11’43.00”
  beam(maj, min, psi): (3.36”, 3.181”, -89.2°)
  memory: 1.049 MB

../../../_images/usage_inputs_maps_example-maps_4_2.png

[5]:

radio_galaxy = maria.map.get("maps/radio_galaxy_3C_288.fits")
radio_galaxy.plot(slices="all")
print(radio_galaxy)

2026-06-02 15:58:20.921 INFO: Fetching https://github.com/thomaswmorris/maria-data/raw/master/maps/radio_galaxy_3C_288.fits
Downloading: 100%|██████████| 268k/268k [00:00<00:00, 19.9MB/s]

ProjectionMap:
  data(1, 1, 257, 256):
    min: -1.833e+00
    max: 1.008e+02
    units: uJy beam^-1
    quantity: spectral_flux_density_per_beam
  nu(1):
    values: [4.8851] GHz
  z(1):
    values: [0.246]
  eta(257):
    height: 38.4”
    res: -150 marcsec
  xi(256):
    width: 38.25”
    res: 150 marcsec
  frame: ra/dec
  center:
    ra: 13ʰ36ᵐ38.57ˢ
    dec: 39°06’21.64”
  beam(maj, min, psi): (600 marcsec, 600 marcsec, 0 rad)
  memory: 263.2 kB

../../../_images/usage_inputs_maps_example-maps_5_2.png

[6]:

tarantula_nebula = maria.map.get("maps/30dor.fits")
tarantula_nebula.plot(slices="all")
print(tarantula_nebula)

ProjectionMap:
  data(251, 251):
    min: 2.573e+00
    max: 6.945e+02
    units: MJy sr^-1
    quantity: spectral_radiance
  eta(251):
    height: 9’
    res: -2.16”
  xi(251):
    width: 9’
    res: 2.16”
  frame: ra/dec
  center:
    ra: 05ʰ38ᵐ44.48ˢ
    dec: -69°06’12.81”
  beam(maj, min, psi): (0 rad, 0 rad, 0 rad)
  memory: 252 kB

../../../_images/usage_inputs_maps_example-maps_6_1.png

[7]:

orion_molecular_cloud = maria.map.get("maps/orion_A.h5")
orion_molecular_cloud.plot(slices="all")
print(orion_molecular_cloud)

2026-06-02 15:58:22.741 INFO: Fetching https://github.com/thomaswmorris/maria-data/raw/master/maps/orion_A.h5
Downloading: 100%|██████████| 1.01M/1.01M [00:00<00:00, 42.9MB/s]
2026-06-02 15:58:22.790 WARNING: 'x_res' is deprecated and will be removed in a future version (use 'xi_res' instead)
2026-06-02 15:58:22.790 WARNING: 'y_res' is deprecated and will be removed in a future version (use 'eta_res' instead)

ProjectionMap:
  data(780, 502):
    min: -1.175e+00
    max: 1.328e+02
    units: MJy sr^-1
    quantity: spectral_radiance
  eta(780):
    height: 5.193°
    res: -24”
  xi(502):
    width: 3.34°
    res: 24”
  frame: ra/dec
  center:
    ra: 05ʰ35ᵐ14.20ˢ
    dec: -5°22’21.50”
  beam(maj, min, psi): (0 rad, 0 rad, 0 rad)
  memory: 3.132 MB

../../../_images/usage_inputs_maps_example-maps_7_2.png

[8]:

monoceros_molecular_cloud = maria.map.get("maps/monoceros_R2.h5")
monoceros_molecular_cloud.plot(slices="all")
print(monoceros_molecular_cloud)

2026-06-02 15:58:24.323 INFO: Fetching https://github.com/thomaswmorris/maria-data/raw/master/maps/monoceros_R2.h5
Downloading: 100%|██████████| 4.57M/4.57M [00:00<00:00, 125MB/s]

ProjectionMap:
  data(1, 3, 595, 633):
    min: -4.915e+05
    max: 6.616e+04
    units: mJy beam^-1
    quantity: spectral_flux_density_per_beam
  stokes(1):
    components: ['I']
  nu(3):
    values: [149.8962 214.1375 272.5386] GHz
  eta(595):
    height: 9.9’
    res: -1”
  xi(633):
    width: 10.53’
    res: 1”
  frame: ra/dec
  center:
    ra: 06ʰ07ᵐ46.30ˢ
    dec: -6°23’9.00”
  beam(maj, min, psi): ragged
  memory: 9.039 MB

../../../_images/usage_inputs_maps_example-maps_8_2.png

[9]:

protoplanetary_disk = maria.map.get("maps/protoplanetary_disk.fits")
protoplanetary_disk.plot(slices="all")
print(protoplanetary_disk)

2026-06-02 15:58:27.235 INFO: Fetching https://github.com/thomaswmorris/maria-data/raw/master/maps/protoplanetary_disk.fits
Downloading: 100%|██████████| 271k/271k [00:00<00:00, 23.0MB/s]

ProjectionMap:
  data(1, 256, 256):
    min: -8.685e-05
    max: 5.391e-04
    units: Jy beam^-1
    quantity: spectral_flux_density_per_beam
  nu(1):
    values: [220.] GHz
  eta(256):
    height: 1.594”
    res: -6.25 marcsec
  xi(256):
    width: 1.594”
    res: 6.25 marcsec
  frame: ra/dec
  center:
    ra: 00ʰ00ᵐ0.00ˢ
    dec: -23°00’0.00”
  beam(maj, min, psi): (22.17 marcsec, 20.72 marcsec, 92.79°)
  memory: 262.1 kB

../../../_images/usage_inputs_maps_example-maps_9_2.png

[10]:

polarized_quasar = maria.map.get("maps/quasar_3C_286.h5")
polarized_quasar.plot(slices="all")
print(polarized_quasar)

2026-06-02 15:58:29.129 INFO: Fetching https://github.com/thomaswmorris/maria-data/raw/master/maps/quasar_3C_286.h5
Downloading: 100%|██████████| 7.56M/7.56M [00:00<00:00, 143MB/s]

ProjectionMap:
  data(3, 2, 1, 460, 472):
    min: -7.912e+03
    max: 9.685e+03
    units: mJy beam^-1
    quantity: spectral_flux_density_per_beam
  stokes(3):
    components: ['I' 'Q' 'U']
  nu(2):
    values: [214.1375 272.5386] GHz
  z(1):
    values: [0.8493]
  eta(460):
    height: 7.65’
    res: -1000 marcsec
  xi(472):
    width: 7.85’
    res: 1000 marcsec
  frame: ra/dec
  center:
    ra: 13ʰ31ᵐ17.80ˢ
    dec: 30°28’48.13”
  beam(maj, min, psi): ragged
  memory: 10.42 MB

../../../_images/usage_inputs_maps_example-maps_10_2.png

[11]:

circinus = maria.map.get("maps/circinus_galaxy.h5")
circinus.plot(slices=dict(v=[5, 6, 7]))
print(circinus)

2026-06-02 15:58:33.964 INFO: Fetching https://github.com/thomaswmorris/maria-data/raw/master/maps/circinus_galaxy.h5
Downloading: 100%|██████████| 2.24M/2.24M [00:00<00:00, 70.4MB/s]

ProjectionMap:
  data(1, 12, 239, 236):
    min: -7.981e-02
    max: 1.346e+00
    units: K_RJ
    quantity: rayleigh_jeans_temperature
  nu(1):
    values: [345.796] GHz
  v(12):
    min: 64.0 km/s
    max: 614.0 km/s
    res: -50 km/s
  eta(239):
    height: 11.9’
    res: -3”
  xi(236):
    width: 11.75’
    res: 3”
  frame: ra/dec
  center:
    ra: 14ʰ13ᵐ28.85ˢ
    dec: -65°20’21.91”
  beam(maj, min, psi): ragged
  memory: 2.707 MB

../../../_images/usage_inputs_maps_example-maps_11_2.png

[12]:

dust = maria.map.get("maps/dust.fits")
dust.plot(slices="all")
print(dust)

2026-06-02 15:58:36.257 INFO: Fetching https://github.com/thomaswmorris/maria-data/raw/master/maps/dust.fits
Downloading: 100%|██████████| 268k/268k [00:00<00:00, 20.3MB/s]

ProjectionMap:
  data(1, 256, 256):
    min: 1.074e-06
    max: 1.218e-02
    units: K_RJ
    quantity: rayleigh_jeans_temperature
  nu(1):
    values: [344.589] GHz
  eta(256):
    height: 6.044’
    res: -1.422”
  xi(256):
    width: 6.044’
    res: 1.422”
  frame: ra/dec
  center:
    ra: 19ʰ00ᵐ0.00ˢ
    dec: -40°00’0.00”
  beam(maj, min, psi): (1.422”, 1.422”, 0 rad)
  memory: 262.1 kB

../../../_images/usage_inputs_maps_example-maps_12_2.png

[13]:

for i in [1, 2, 3]:

    sz_cluster = maria.map.get(f"maps/cluster{i}.fits")
    sz_cluster.plot(slices="all")
    print(sz_cluster)

2026-06-02 15:58:38.302 INFO: Fetching https://github.com/thomaswmorris/maria-data/raw/master/maps/cluster1.fits
Downloading: 100%|██████████| 4.20M/4.20M [00:00<00:00, 125MB/s]

ProjectionMap:
  data(1024, 1024):
    min: -1.591e-05
    max: -1.117e-08
    units: compton_y
    quantity: compton_y
  eta(1024):
    height: 59.94’
    res: -3.516”
  xi(1024):
    width: 59.94’
    res: 3.516”
  frame: ra/dec
  center:
    ra: 17ʰ20ᵐ0.00ˢ
    dec: -10°00’0.00”
  beam(maj, min, psi): (0 rad, 0 rad, 0 rad)
  memory: 4.194 MB
ProjectionMap:
  data(1024, 1024):
    min: -3.845e-06
    max: -5.005e-09
    units: compton_y
    quantity: compton_y
  eta(1024):
    height: 59.94’
    res: -3.516”
  xi(1024):
    width: 59.94’
    res: 3.516”
  frame: ra/dec
  center:
    ra: 17ʰ20ᵐ0.00ˢ
    dec: -10°00’0.00”
  beam(maj, min, psi): (0 rad, 0 rad, 0 rad)
  memory: 4.194 MB

2026-06-02 15:58:39.738 INFO: Fetching https://github.com/thomaswmorris/maria-data/raw/master/maps/cluster3.fits
Downloading: 100%|██████████| 4.20M/4.20M [00:00<00:00, 120MB/s]

ProjectionMap:
  data(1024, 1024):
    min: -1.155e-05
    max: -2.919e-09
    units: compton_y
    quantity: compton_y
  eta(1024):
    height: 59.94’
    res: -3.516”
  xi(1024):
    width: 59.94’
    res: 3.516”
  frame: ra/dec
  center:
    ra: 17ʰ20ᵐ0.00ˢ
    dec: -10°00’0.00”
  beam(maj, min, psi): (0 rad, 0 rad, 0 rad)
  memory: 4.194 MB

../../../_images/usage_inputs_maps_example-maps_13_4.png

../../../_images/usage_inputs_maps_example-maps_13_5.png

../../../_images/usage_inputs_maps_example-maps_13_6.png

[14]:

einstein = maria.map.get("maps/einstein.h5")
einstein.plot(slices="all")
print(einstein)

ProjectionMap:
  data(4, 1, 685, 685):
    min: -1.000e-02
    max: 2.540e-01
    units: K_RJ
    quantity: rayleigh_jeans_temperature
  stokes(4):
    components: ['I' 'Q' 'U' 'V']
  nu(1):
    values: [90.] GHz
  eta(685):
    height: 60’
    res: -5.263”
  xi(685):
    width: 60’
    res: 5.263”
  frame: ra/dec
  center:
    ra: 00ʰ00ᵐ0.00ˢ
    dec: -23°00’0.00”
  beam(maj, min, psi): (0 rad, 0 rad, 0 rad)
  memory: 7.508 MB

../../../_images/usage_inputs_maps_example-maps_14_1.png

[15]:

mary = maria.map.get("maps/maria.h5")
mary.plot(slices="all")
print(mary)

2026-06-02 15:58:46.889 INFO: Fetching https://github.com/thomaswmorris/maria-data/raw/master/maps/maria.h5
Downloading: 100%|██████████| 411k/411k [00:00<00:00, 33.8MB/s]

ProjectionMap:
  data(4, 1, 500, 500):
    min: -1.000e-02
    max: 2.515e-01
    units: K_RJ
    quantity: rayleigh_jeans_temperature
  stokes(4):
    components: ['I' 'Q' 'U' 'V']
  nu(1):
    values: [90.] GHz
  eta(500):
    height: 59.88’
    res: -7.2”
  xi(500):
    width: 59.88’
    res: 7.2”
  frame: ra/dec
  center:
    ra: 00ʰ00ᵐ0.00ˢ
    dec: -23°00’0.00”
  beam(maj, min, psi): (0 rad, 0 rad, 0 rad)
  memory: 4 MB

../../../_images/usage_inputs_maps_example-maps_15_2.png