Note
Click here to download the full example code
Downloading and plotting an HMI magnetogram¶
How to download an HMI magnetogram data with Fido and make a plot.
import astropy.units as u
import matplotlib.pyplot as plt
import sunpy.map
from sunpy.net import Fido, attrs as a
To download the required data, we use
Fido, a downloader client.
First define the search variables, a timerange, the instrument,
the observation type,
and a cadence of images spaced every 720 seconds.
result = Fido.search(a.Time('2011/11/09 17:40:00', '2011/11/09 17:55:00'),
a.Instrument('hmi'),
a.Sample(720*u.s),
a.vso.Physobs('LOS_magnetic_field'))
Now we can see what results we obtained from our search. Notice we have two files.
print(result)
Out:
Results from 1 Provider:
2 Results from the VSOClient:
Start Time [1] End Time [1] Source ... Type Wavelength [2]
... Angstrom
str19 str19 str3 ... str8 float64
------------------- ------------------- ------ ... -------- ----------------
2011-11-09 17:40:13 2011-11-09 17:40:14 SDO ... FULLDISK 6173.0 .. 6174.0
2011-11-09 17:52:13 2011-11-09 17:52:14 SDO ... FULLDISK 6173.0 .. 6174.0
The following shows how to download the results. If we don’t provide a path it will download the file into the sunpy data directory. The output provides the path of the downloaded files. The result can be from several data clients, so we have to index the client first client and then index the file.
# Slice the first record returned by the first client.
downloaded_file = Fido.fetch(result[0, 0])
print(downloaded_file)
Out:
Files Downloaded: 0%| | 0/1 [00:00<?, ?file/s]
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 0%| | 0.00/15.7M [00:00<?, ?B/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 0%| | 42.4k/15.7M [00:00<00:36, 424kB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 0%| | 67.2k/15.7M [00:00<00:44, 349kB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 1%| | 116k/15.7M [00:00<00:40, 381kB/s] [A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 1%|1 | 170k/15.7M [00:00<00:37, 419kB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 2%|1 | 239k/15.7M [00:00<00:32, 473kB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 2%|2 | 330k/15.7M [00:00<00:27, 553kB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 3%|2 | 423k/15.7M [00:00<00:24, 629kB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 3%|3 | 534k/15.7M [00:00<00:21, 722kB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 4%|4 | 660k/15.7M [00:00<00:18, 827kB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 5%|5 | 803k/15.7M [00:01<00:15, 945kB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 6%|6 | 963k/15.7M [00:01<00:13, 1.08MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 7%|7 | 1.15M/15.7M [00:01<00:11, 1.23MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 9%|8 | 1.35M/15.7M [00:01<00:10, 1.39MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 10%|# | 1.58M/15.7M [00:01<00:08, 1.58MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 12%|#1 | 1.84M/15.7M [00:01<00:07, 1.78MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 14%|#3 | 2.13M/15.7M [00:01<00:06, 2.02MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 16%|#5 | 2.45M/15.7M [00:01<00:05, 2.26MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 18%|#7 | 2.80M/15.7M [00:01<00:05, 2.52MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 20%|## | 3.17M/15.7M [00:01<00:04, 2.79MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 23%|##2 | 3.58M/15.7M [00:02<00:03, 3.07MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 26%|##5 | 4.03M/15.7M [00:02<00:03, 3.39MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 29%|##8 | 4.53M/15.7M [00:02<00:02, 3.74MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 32%|###2 | 5.07M/15.7M [00:02<00:02, 4.12MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 36%|###6 | 5.68M/15.7M [00:02<00:02, 4.56MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 41%|#### | 6.36M/15.7M [00:02<00:01, 5.05MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 45%|####5 | 7.09M/15.7M [00:02<00:01, 5.55MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 50%|##### | 7.87M/15.7M [00:02<00:01, 6.06MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 56%|#####5 | 8.71M/15.7M [00:02<00:01, 6.61MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 61%|######1 | 9.62M/15.7M [00:02<00:00, 7.17MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 68%|######7 | 10.6M/15.7M [00:03<00:00, 7.82MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 75%|#######4 | 11.7M/15.7M [00:03<00:00, 8.52MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 82%|########2 | 12.9M/15.7M [00:03<00:00, 9.29MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 90%|######### | 14.2M/15.7M [00:03<00:00, 10.1MB/s][A
hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits: 99%|#########9| 15.6M/15.7M [00:03<00:00, 11.0MB/s][A
[A
Files Downloaded: 100%|##########| 1/1 [00:14<00:00, 14.17s/file]
Files Downloaded: 100%|##########| 1/1 [00:14<00:00, 14.17s/file]
['/home/docs/sunpy/data/hmi_m_45s_2011_11_09_17_41_15_tai_magnetogram.fits']
Now load it into a map and plot it
hmi_map = sunpy.map.Map(downloaded_file[0])
fig = plt.figure()
hmi_map.plot()
plt.show()
Total running time of the script: ( 0 minutes 22.648 seconds)