Vẽ hình ảnh mặt trời với SunPy trong Python

SunPy là thư viện Python cho vật lý mặt trời, hỗ trợ đọc/hiển thị dữ liệu từ các quan sát điểm (SDO/AIA, SOHO, STEREO...).

Cài đặt

pip install sunpy

Vẽ ảnh AIA 171Å mẫu

import sunpy.map
import matplotlib.pyplot as plt
import sunpy.data.sample

# Load sample data (AIA 171Å)
aia_map = sunpy.map.Map(sunpy.data.sample.AIA_171_IMAGE)

fig = plt.figure()
ax = plt.subplot(111, projection=aia_map)
aia_map.plot()
aia_map.draw_limb()      # Vẽ đường viền mặt trời
aia_map.draw_grid()      # Vẽ lưới tọa độ
plt.colorbar()           # Thanh màu cường độ
plt.show()

Kết quả

Tính năng chính SunPy

• Đọc nhiều format: FITS, CDF, JP2, ASDF...
• Tọa độ hệ mặt trời (Heliographic, Helioprojective...).
• Tìm kiếm/tải dữ liệu từ VSO, JSOC, Fido.
• Tích hợp Matplotlib để vẽ.
• Phân tích flare, CME, active region.

Ví dụ tải dữ liệu thực

from sunpy.net import Fido, attrs as a
results = Fido.search(a.Time('2023-01-01', '2023-01-02'),
                      a.Instrument('AIA'),
                      a.Wavelength(171*u.angstrom))
files = Fido.fetch(results)
solar_map = sunpy.map.Map(files[0])