Plot PHI Blos overlay on EUI maps#

This example demonstrates how to overlay relevant information from line-of-sight magnetic field: (L2 Blos) data from the Solar Orbiter PHI instrument on maps from the Solar Orbiter EUI instrument.

import sunpy_soar
from sunpy.net import Fido, attrs as a
from astropy.time import Time
import sunpy.map
import sunpy.visualization.colormaps
import matplotlib.pyplot as plt
from astropy.io import fits
import numpy as np

from sunpy.coordinates import SphericalScreen
from astropy.coordinates import SkyCoord
import astropy.units as u
from sunpy.coordinates import propagate_with_solar_surface

Searching for PHI-HRT Blos Data and co-temporal EUI-HRI Data#

We first search for Solar Orbiter PHI-HRT (High Resolution Telescope) Blos data in a given time range. The search results will return metadata about available files. no PHI-FDT available during this time range

t_start_hrt = Time('2024-10-15T18:00:00', format='isot', scale='utc')
t_end_hrt = Time('2024-10-15T18:00:35', format='isot', scale='utc')

search_results_phi_hrt = Fido.search(a.Instrument('PHI'), a.Time(t_start_hrt.value, t_end_hrt.value), (a.soar.Product('phi-hrt-blos')))
print(search_results_phi_hrt)

search_results_eui_hri = Fido.search(a.Instrument('EUI'), a.Time(t_start_hrt.value, t_end_hrt.value), (a.soar.Product('eui-hrieuv174-image')), a.Level('L2'))
print(search_results_eui_hri)

#FSI for context
t_start_fsi = Time('2024-10-15T18:00:00', format='isot', scale='utc')
t_end_fsi = Time('2024-10-15T18:02:00', format='isot', scale='utc')

search_results_eui_fsi = Fido.search(a.Instrument('EUI'), a.Time(t_start_fsi.value, t_end_fsi.value), (a.soar.Product('eui-fsi174-image')))
print(search_results_eui_fsi)

Results from 1 Provider:

1 Results from the SOARClient:

Instrument Data product Level        Start time               End time        Filesize          SOOP Name           Detector Wavelength
                                                                               Mbyte
---------- ------------ ----- ----------------------- ----------------------- -------- ---------------------------- -------- ----------
       PHI phi-hrt-blos    L2 2024-10-15 18:00:03.834 2024-10-15 18:01:26.596   13.084 L_BOTH_HRES_HCAD_Major-Flare      HRT   6173.341


Results from 1 Provider:

2 Results from the SOARClient:

Instrument     Data product    Level        Start time               End time        Filesize          SOOP Name           Detector Wavelength
                                                                                      Mbyte
---------- ------------------- ----- ----------------------- ----------------------- -------- ---------------------------- -------- ----------
       EUI eui-hrieuv174-image    L2 2024-10-15 18:00:12.242 2024-10-15 18:00:14.242    5.057 L_BOTH_HRES_HCAD_Major-Flare  HRI_EUV      174.0
       EUI eui-hrieuv174-image    L2 2024-10-15 18:00:28.242 2024-10-15 18:00:30.242     5.06 L_BOTH_HRES_HCAD_Major-Flare  HRI_EUV      174.0


Results from 1 Provider:

2 Results from the SOARClient:

Instrument   Data product   Level        Start time               End time        Filesize SOOP Name Detector Wavelength
                                                                                   Mbyte
---------- ---------------- ----- ----------------------- ----------------------- -------- --------- -------- ----------
       EUI eui-fsi174-image    L1 2024-10-15 18:00:55.242 2024-10-15 18:01:05.242    4.308      none      FSI      174.0
       EUI eui-fsi174-image    L2 2024-10-15 18:00:55.242 2024-10-15 18:01:05.242    4.939      none      FSI      174.0

sr_eui_fsi = search_results_eui_fsi[0,0]
eui_fsi_file = Fido.fetch(sr_eui_fsi)
eui_fsi_map = sunpy.map.Map(eui_fsi_file[0])

sr_eui_hri = search_results_eui_hri[0,0]
eui_hri_file = Fido.fetch(sr_eui_hri)
eui_hri_map = sunpy.map.Map(eui_hri_file[0])

sr_phi_hrt = search_results_phi_hrt[0,0]
blos_file = Fido.fetch(sr_phi_hrt)
blos_map = sunpy.map.Map(blos_file[0])

Files Downloaded:   0%|          | 0/1 [00:00<?, ?file/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:   0%|          | 0.00/4.31M [00:00<?, ?B/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:   0%|          | 1.02k/4.31M [00:00<12:59, 5.52kB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:   1%|▏         | 63.9k/4.31M [00:00<00:15, 274kB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:   3%|▎         | 120k/4.31M [00:00<00:11, 380kB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:   5%|▍         | 207k/4.31M [00:00<00:07, 545kB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:   6%|▋         | 276k/4.31M [00:00<00:06, 586kB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:   9%|▉         | 391k/4.31M [00:00<00:05, 755kB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  12%|█▏        | 521k/4.31M [00:00<00:04, 916kB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  15%|█▌        | 651k/4.31M [00:00<00:03, 1.02MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  18%|█▊        | 782k/4.31M [00:01<00:03, 1.10MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  21%|██        | 911k/4.31M [00:01<00:02, 1.15MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  24%|██▍       | 1.04M/4.31M [00:01<00:02, 1.18MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  27%|██▋       | 1.17M/4.31M [00:01<00:02, 1.20MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  30%|███       | 1.30M/4.31M [00:01<00:02, 1.22MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  33%|███▎      | 1.43M/4.31M [00:01<00:02, 1.23MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  36%|███▌      | 1.56M/4.31M [00:01<00:02, 1.24MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  39%|███▉      | 1.69M/4.31M [00:01<00:02, 1.25MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  42%|████▏     | 1.82M/4.31M [00:01<00:01, 1.25MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  45%|████▌     | 1.95M/4.31M [00:01<00:01, 1.26MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  48%|████▊     | 2.08M/4.31M [00:02<00:01, 1.25MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  51%|█████▏    | 2.21M/4.31M [00:02<00:01, 1.26MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  54%|█████▍    | 2.34M/4.31M [00:02<00:01, 1.26MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  57%|█████▋    | 2.47M/4.31M [00:02<00:01, 1.26MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  60%|██████    | 2.60M/4.31M [00:02<00:01, 1.26MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  63%|██████▎   | 2.73M/4.31M [00:02<00:01, 1.26MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  66%|██████▋   | 2.86M/4.31M [00:02<00:01, 1.26MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  69%|██████▉   | 2.99M/4.31M [00:02<00:01, 1.26MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  72%|███████▏  | 3.12M/4.31M [00:02<00:00, 1.26MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  75%|███████▌  | 3.25M/4.31M [00:02<00:00, 1.26MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  78%|███████▊  | 3.38M/4.31M [00:03<00:00, 1.26MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  82%|████████▏ | 3.51M/4.31M [00:03<00:00, 1.26MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  85%|████████▍ | 3.64M/4.31M [00:03<00:00, 1.26MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  88%|████████▊ | 3.77M/4.31M [00:03<00:00, 1.26MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  91%|█████████ | 3.90M/4.31M [00:03<00:00, 1.26MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  94%|█████████▎| 4.03M/4.31M [00:03<00:00, 1.26MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits:  97%|█████████▋| 4.16M/4.31M [00:03<00:00, 1.26MB/s]

solo_L1_eui-fsi174-image_20241015T180055242_V01.fits: 100%|█████████▉| 4.29M/4.31M [00:03<00:00, 1.26MB/s]

Files Downloaded: 100%|██████████| 1/1 [00:04<00:00,  4.10s/file]
Files Downloaded: 100%|██████████| 1/1 [00:04<00:00,  4.10s/file]

Files Downloaded:   0%|          | 0/1 [00:00<?, ?file/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:   0%|          | 0.00/5.06M [00:00<?, ?B/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:   0%|          | 1.02k/5.06M [00:00<16:04, 5.24kB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:   1%|▏         | 64.0k/5.06M [00:00<00:19, 259kB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:   2%|▏         | 121k/5.06M [00:00<00:13, 358kB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:   5%|▍         | 250k/5.06M [00:00<00:07, 650kB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:   8%|▊         | 381k/5.06M [00:00<00:05, 833kB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  10%|█         | 511k/5.06M [00:00<00:04, 948kB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  13%|█▎        | 641k/5.06M [00:00<00:04, 1.02MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  15%|█▌        | 771k/5.06M [00:00<00:04, 1.07MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  18%|█▊        | 901k/5.06M [00:01<00:03, 1.10MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  20%|██        | 1.03M/5.06M [00:01<00:03, 1.13MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  23%|██▎       | 1.16M/5.06M [00:01<00:03, 1.14MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  26%|██▌       | 1.29M/5.06M [00:01<00:03, 1.15MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  28%|██▊       | 1.42M/5.06M [00:01<00:03, 1.16MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  31%|███       | 1.55M/5.06M [00:01<00:03, 1.17MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  33%|███▎      | 1.68M/5.06M [00:01<00:02, 1.17MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  36%|███▌      | 1.81M/5.06M [00:01<00:02, 1.17MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  38%|███▊      | 1.94M/5.06M [00:01<00:02, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  41%|████      | 2.07M/5.06M [00:02<00:02, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  44%|████▎     | 2.20M/5.06M [00:02<00:02, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  46%|████▌     | 2.33M/5.06M [00:02<00:02, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  49%|████▊     | 2.46M/5.06M [00:02<00:02, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  51%|█████     | 2.59M/5.06M [00:02<00:02, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  54%|█████▍    | 2.72M/5.06M [00:02<00:01, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  56%|█████▋    | 2.85M/5.06M [00:02<00:01, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  59%|█████▉    | 2.98M/5.06M [00:02<00:01, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  61%|██████▏   | 3.11M/5.06M [00:02<00:01, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  64%|██████▍   | 3.24M/5.06M [00:03<00:01, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  67%|██████▋   | 3.37M/5.06M [00:03<00:01, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  69%|██████▉   | 3.50M/5.06M [00:03<00:01, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  72%|███████▏  | 3.63M/5.06M [00:03<00:01, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  74%|███████▍  | 3.76M/5.06M [00:03<00:01, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  77%|███████▋  | 3.89M/5.06M [00:03<00:00, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  79%|███████▉  | 4.02M/5.06M [00:03<00:00, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  82%|████████▏ | 4.15M/5.06M [00:03<00:00, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  85%|████████▍ | 4.28M/5.06M [00:03<00:00, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  87%|████████▋ | 4.41M/5.06M [00:04<00:00, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  90%|████████▉ | 4.54M/5.06M [00:04<00:00, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  92%|█████████▏| 4.67M/5.06M [00:04<00:00, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  95%|█████████▍| 4.80M/5.06M [00:04<00:00, 1.18MB/s]

solo_L2_eui-hrieuv174-image_20241015T180012242_V01.fits:  97%|█████████▋| 4.93M/5.06M [00:04<00:00, 1.18MB/s]

Files Downloaded: 100%|██████████| 1/1 [00:04<00:00,  4.82s/file]
Files Downloaded: 100%|██████████| 1/1 [00:04<00:00,  4.82s/file]

Files Downloaded:   0%|          | 0/1 [00:00<?, ?file/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:   0%|          | 0.00/13.1M [00:00<?, ?B/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:   0%|          | 1.02k/13.1M [00:00<40:22, 5.40kB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:   0%|          | 64.0k/13.1M [00:00<00:48, 271kB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:   1%|          | 121k/13.1M [00:00<00:34, 377kB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:   2%|▏         | 251k/13.1M [00:00<00:18, 691kB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:   3%|▎         | 381k/13.1M [00:00<00:14, 884kB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:   4%|▍         | 512k/13.1M [00:00<00:12, 1.01MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:   5%|▍         | 642k/13.1M [00:00<00:11, 1.09MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:   6%|▌         | 772k/13.1M [00:00<00:10, 1.14MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:   7%|▋         | 902k/13.1M [00:01<00:10, 1.18MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:   8%|▊         | 1.03M/13.1M [00:01<00:10, 1.20MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:   9%|▉         | 1.16M/13.1M [00:01<00:09, 1.22MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  10%|▉         | 1.29M/13.1M [00:01<00:09, 1.23MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  11%|█         | 1.42M/13.1M [00:01<00:09, 1.24MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  12%|█▏        | 1.55M/13.1M [00:01<00:09, 1.24MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  13%|█▎        | 1.68M/13.1M [00:01<00:09, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  14%|█▍        | 1.81M/13.1M [00:01<00:08, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  15%|█▍        | 1.94M/13.1M [00:01<00:08, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  16%|█▌        | 2.07M/13.1M [00:01<00:08, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  17%|█▋        | 2.20M/13.1M [00:02<00:08, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  18%|█▊        | 2.33M/13.1M [00:02<00:08, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  19%|█▉        | 2.46M/13.1M [00:02<00:08, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  20%|█▉        | 2.59M/13.1M [00:02<00:08, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  21%|██        | 2.72M/13.1M [00:02<00:08, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  22%|██▏       | 2.85M/13.1M [00:02<00:08, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  23%|██▎       | 2.98M/13.1M [00:02<00:08, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  24%|██▍       | 3.11M/13.1M [00:02<00:07, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  25%|██▍       | 3.24M/13.1M [00:02<00:07, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  26%|██▌       | 3.37M/13.1M [00:02<00:07, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  27%|██▋       | 3.50M/13.1M [00:03<00:07, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  28%|██▊       | 3.63M/13.1M [00:03<00:07, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  29%|██▊       | 3.76M/13.1M [00:03<00:07, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  30%|██▉       | 3.89M/13.1M [00:03<00:07, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  31%|███       | 4.02M/13.1M [00:03<00:07, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  32%|███▏      | 4.15M/13.1M [00:03<00:07, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  33%|███▎      | 4.28M/13.1M [00:03<00:07, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  34%|███▎      | 4.41M/13.1M [00:03<00:06, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  35%|███▍      | 4.54M/13.1M [00:03<00:06, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  36%|███▌      | 4.67M/13.1M [00:04<00:06, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  37%|███▋      | 4.80M/13.1M [00:04<00:06, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  38%|███▊      | 4.93M/13.1M [00:04<00:06, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  39%|███▊      | 5.06M/13.1M [00:04<00:06, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  40%|███▉      | 5.19M/13.1M [00:04<00:06, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  41%|████      | 5.32M/13.1M [00:04<00:06, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  42%|████▏     | 5.45M/13.1M [00:04<00:06, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  43%|████▎     | 5.58M/13.1M [00:04<00:05, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  44%|████▎     | 5.71M/13.1M [00:04<00:05, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  45%|████▍     | 5.84M/13.1M [00:04<00:05, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  46%|████▌     | 5.97M/13.1M [00:05<00:05, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  47%|████▋     | 6.10M/13.1M [00:05<00:05, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  48%|████▊     | 6.23M/13.1M [00:05<00:05, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  49%|████▊     | 6.36M/13.1M [00:05<00:05, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  50%|████▉     | 6.49M/13.1M [00:05<00:05, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  51%|█████     | 6.62M/13.1M [00:05<00:05, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  52%|█████▏    | 6.75M/13.1M [00:05<00:05, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  53%|█████▎    | 6.88M/13.1M [00:05<00:04, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  54%|█████▎    | 7.01M/13.1M [00:05<00:04, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  55%|█████▍    | 7.14M/13.1M [00:05<00:04, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  56%|█████▌    | 7.27M/13.1M [00:06<00:04, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  57%|█████▋    | 7.40M/13.1M [00:06<00:04, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  58%|█████▊    | 7.53M/13.1M [00:06<00:04, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  59%|█████▊    | 7.66M/13.1M [00:06<00:04, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  60%|█████▉    | 7.79M/13.1M [00:06<00:04, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  61%|██████    | 7.92M/13.1M [00:06<00:04, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  62%|██████▏   | 8.05M/13.1M [00:06<00:03, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  63%|██████▎   | 8.18M/13.1M [00:06<00:03, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  64%|██████▎   | 8.31M/13.1M [00:06<00:03, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  65%|██████▍   | 8.44M/13.1M [00:07<00:03, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  66%|██████▌   | 8.57M/13.1M [00:07<00:03, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  66%|██████▋   | 8.70M/13.1M [00:07<00:03, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  67%|██████▋   | 8.83M/13.1M [00:07<00:03, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  68%|██████▊   | 8.96M/13.1M [00:07<00:03, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  69%|██████▉   | 9.09M/13.1M [00:07<00:03, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  70%|███████   | 9.22M/13.1M [00:07<00:03, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  71%|███████▏  | 9.35M/13.1M [00:07<00:02, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  72%|███████▏  | 9.48M/13.1M [00:07<00:02, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  73%|███████▎  | 9.61M/13.1M [00:07<00:02, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  74%|███████▍  | 9.74M/13.1M [00:08<00:02, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  75%|███████▌  | 9.87M/13.1M [00:08<00:02, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  76%|███████▋  | 10.0M/13.1M [00:08<00:02, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  77%|███████▋  | 10.1M/13.1M [00:08<00:02, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  78%|███████▊  | 10.3M/13.1M [00:08<00:02, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  79%|███████▉  | 10.4M/13.1M [00:08<00:02, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  80%|████████  | 10.5M/13.1M [00:08<00:02, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  81%|████████▏ | 10.6M/13.1M [00:08<00:01, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  82%|████████▏ | 10.8M/13.1M [00:08<00:01, 1.25MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  83%|████████▎ | 10.9M/13.1M [00:08<00:01, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  84%|████████▍ | 11.0M/13.1M [00:09<00:01, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  85%|████████▌ | 11.2M/13.1M [00:09<00:01, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  86%|████████▋ | 11.3M/13.1M [00:09<00:01, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  87%|████████▋ | 11.4M/13.1M [00:09<00:01, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  88%|████████▊ | 11.6M/13.1M [00:09<00:01, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  89%|████████▉ | 11.7M/13.1M [00:09<00:01, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  90%|█████████ | 11.8M/13.1M [00:09<00:01, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  91%|█████████▏| 11.9M/13.1M [00:09<00:00, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  92%|█████████▏| 12.1M/13.1M [00:09<00:00, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  93%|█████████▎| 12.2M/13.1M [00:10<00:00, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  94%|█████████▍| 12.3M/13.1M [00:10<00:00, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  95%|█████████▌| 12.5M/13.1M [00:10<00:00, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  96%|█████████▋| 12.6M/13.1M [00:10<00:00, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  97%|█████████▋| 12.7M/13.1M [00:10<00:00, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  98%|█████████▊| 12.9M/13.1M [00:10<00:00, 1.26MB/s]

solo_L2_phi-hrt-blos_20241015T180003_V01_0450150701.fits:  99%|█████████▉| 13.0M/13.1M [00:10<00:00, 1.26MB/s]

Files Downloaded: 100%|██████████| 1/1 [00:10<00:00, 10.95s/file]
Files Downloaded: 100%|██████████| 1/1 [00:10<00:00, 10.95s/file]

# Check HRT WCS is corrected (using history and or drawing limb)
history = blos_map.fits_header['HISTORY']
for line in history:
    if 'WCS' in line:
        print(line)

/home/docs/checkouts/readthedocs.org/user_builds/solar-orbiter-docs/envs/latest/lib/python3.12/site-packages/sunpy/io/_fits.py:246: SunpyMetadataWarning: The meta key  is not valid ascii, dropping from the FITS header
  warn_metadata(f'The meta key {k} is not valid ascii, dropping from the FITS header')
WCS updated by HRT pipeline using FDT, S/W version: 1.0.1. Check parent
file for old WCS.

# Crop EUI FSI Map
center = SkyCoord(-2500*u.arcsec, -2500*u.arcsec, frame=eui_fsi_map.coordinate_frame)
width = 5000*u.arcsec
height = 5000*u.arcsec

eui_fsi_map = eui_fsi_map.submap(center, width=width, height=height)

Plot Full Sun context from EUI-FSI and show EUI-HRI and PHI-HRT FoV#

fig, ax = plt.subplots(figsize=(10,10), subplot_kw={'projection': eui_fsi_map})

eui_fsi_map.plot(axes=ax)

with SphericalScreen(eui_hri_map.observer_coordinate, only_off_disk=False):
    eui_hri_map.draw_extent(axes=ax, color='k', lw=1, ls='dotted', label='HRI')

with SphericalScreen(blos_map.observer_coordinate, only_off_disk=False):
    blos_map.draw_extent(axes=ax, color='red', lw=1, ls='dotted', label='HRT')

eui_hri_map.draw_extent(axes=ax, color='k', lw=1)

ax.legend(loc='upper right')

eui_fsi_map.draw_limb()

# Draw the HRT limb, if it is close to the FSI one, then the HRT WCS will be ok.
blos_map.draw_limb(color='red')

ax.grid(False)
plt.show()

$FSI $174 \; \mathrm{\mathring{A}}$ 2024-10-15 18:00:55$

Plot PHI-HRT Blos contours on top of EUI-HRI#

The EUI-HRI and PHI-HRT computed limb does not quite line up with the limb as seen in the EUI-HRI image. It may be that the WCS in both are not perfect.

#propagage_with_solar_surface corrects for solar rotation if time between datasets
with propagate_with_solar_surface():
    fig = plt.figure(figsize=(8,8), layout='tight')
    ax = fig.add_subplot(projection=eui_hri_map)
    eui_hri_map.plot(annotate=True, axes=ax)
    xlim = ax.get_xlim()
    ylim = ax.get_ylim()

    blos_map.draw_limb(color='red')
    eui_hri_map.draw_limb(color='k')

    for contour in blos_map.find_contours(-200*u.Gauss):
        ax.plot_coord(contour,color='k', alpha=.75)
    for contour in blos_map.find_contours(200*u.Gauss):
        ax.plot_coord(contour,color='w',alpha=.75)

    ax.set_xlim(xlim); ax.set_ylim(ylim);

    plt.show()

$HRI_EUV $174 \; \mathrm{\mathring{A}}$ 2024-10-15 18:00:12$

Total running time of the script: (1 minutes 12.988 seconds)

Gallery generated by Sphinx-Gallery