পাইথনে SunPy ব্যবহার করে সৌর ছবি প্লট করা

SunPy কী?

SunPy হলো সৌর фізиكس ও ডেটা বিশ্লেষণের জন্য পাইথন লাইব্রেরি। এতে সৌর অভ{Y}derabad (SDO, SOHO ইত্যাদি) থেকে ডেটা লোড, ভিজ্যুয়ালাইজ, অ্যানালাইসিস করার টুলস রয়েছে।

স্থাপন

pip install sunpy

উদাহরণ: AIA ১৭১ Å ছবি প্লট


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

# স্যাম্পল AIA 171 Å ছবি লোড
aia_map = sunpy.map.Map(sunpy.data.sample.AIA_171_IMAGE)

# প্লট
fig = plt.figure(figsize=(10, 8))
ax = plt.subplot(111, projection=aia_map)
aia_map.plot(axes=ax, cmap='sdoaia171')  # SDO AIA 171 কালার ম্যাপ
aia_map.draw_limb()      # সৌর লিম্ব (किनারা)
aia_map.draw_grid()      # গ্রিড (দ্রাঘিমাংশ/অক্ষাংশ)
plt.colorbar(ax=ax, label='Intensity (DN/s)')
plt.title('SDO/AIA 171 Å')
plt.show()

মূল ফাংশনসমূহ

sunpy.map.Map() — FITS/JP2 ফাইল থেকে Map অবজেক্ট বানায়।
.plot() — ম্যাটপ্লটলিবে প্লট করে।
.draw_limb() — সূর্যের কিনারা আঁকে।
.draw_grid() — হেলিওগ্রাফিক গ্রিড আঁকে।
.submap() — নির্দিষ্ট অঞ্চল ক্রপ করা।
.rotate() — রোটেশন কোরেকশন।

বাস্তব ডেটা ডাউনলোড (Fido)


from sunpy.net import Fido, attrs as a

# SDO/AIA 171 Å, নির্দিষ্ট সময়ের ডেটা সার্চ
results = Fido.search(
    a.Time('2023-01-01 00:00', '2023-01-01 01:00'),
    a.Instrument('AIA'),
    a.Wavelength(171 * a.u.Angstrom)
)
print(results)

# ডাউনলোড
files = Fido.fetch(results)
solar_map = sunpy.map.Map(files[0])

প্রয়োজনীয় প্যাকেজ

pip install sunpy matplotlib astropy