การพล็อตภาพสุริยะ (Solar Image) ด้วย Python และ SunPy

Python มีแพ็กเกจ SunPy สำหรับสร้างภาพสุริยะ ในแพ็กเกจนี้มีไฟล์ข้อมูลสุริยะต่างๆ เช่นข้อมูลโปรตอน/อิเล็กตรอนจากสังเกตการณ์สุริยะและห้องปฏิบัติการสุริยะ

การติดตั้ง

pip install sunpy

ตัวอย่าง: พล็อตภาพ AIA

AIA = Atmospheric Imaging Assembly อุปกรณ์บนดาวเทียม SDO (Solar Dynamics Observatory)

ใช้ sunpy.Map() สร้างแผนที่จากผลิตภัณฑ์ข้อมูลที่รองรับ

โค้ดตัวอย่าง

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

# โหลดภาพตัวอย่าง AIA 171 Angstrom
my_aia = sunpy.map.Map(sunpy.data.sample.AIA_171_IMAGE)

fig = plt.figure()
ax = plt.subplot(111, projection=my_aia)

my_aia.plot()
my_aia.draw_limb()      # วาดขอบดวงอาทิตย์
my_aia.draw_grid()      # วาดกริดพิกัด
plt.colorbar()          # แสดง colorbar
plt.show()

ผลลัพธ์

ฟังก์ชันเพิ่มเติมที่น่าสนใจ:

• my_aia.draw_limb() - วาดขอบดวงอาทิตย์
• my_aia.draw_grid() - วาดกริดละติจูด/ลองจิจูด
• my_aia.peek() - แสดงภาพพร้อม UI โต้ตอบ
• sunpy.map.Map('file.fits') - โหลดไฟล์ FITS ของตัวเอง

SunPy เป็นไลบรารีที่ทรงพลังสำหรับการวิเคราะห์ข้อมูลสุริยะ รองรับหลายโทรทัศน์และภารกิจอวกาศ