/
Rgelus
/
science_processing
Обзор
Документация
Войти
/
Rgelus
/
science_processing
Код
Запросы
0
Задачи
Вики
Пакеты
0
Релизы
0
CI/CD
Аналитика
Безопасность
master
extraction/main.py
65 строк
1 KB
Rgelus
merger 3 projects
29 июн 2025, 09:53
29 июн 2025, 09:53
1593f54
Код
Авторство
О чём код?
# https://oafa.master.sai.msu.ru/viewer.php?id=103170443 # Откуда взял исходники import matplotlib.pyplot as plt import numpy as np from astropy.visualization import astropy_mpl_style plt.style.use(astropy_mpl_style) from astropy.utils.data import get_pkg_data_filename from astropy.io import fits # image_file = fits.open('fits/1116955.fit') # # image_data = fits.getdata('fits/1116955.fit', ext=0) # # print(image_data.shape) # # # plt.figure() # plt.imshow(np.log10(image_data), vmin = 3, vmax = 3.5, cmap='gray') # # plt.imshow(image_data,vmin = 0, vmax = 8000, cmap='gray') # plt.scatter(5530, 4227, s=3, c='red', marker='o') # plt.colorbar() # plt.show() # image_data = fits.getdata('fits/1116969.fit', ext=0) # # print(image_data.shape) # # # plt.figure() # # plt.imshow(np.log10(image_data), vmin = 3.0, vmax = 3.4, cmap='gray') # # plt.imshow(np.log10(image_data), cmap='gray') # # plt.imshow(np.log10(image_data), cmap='gray') # plt.imshow(image_data,vmin = 1500, vmax = 2000, cmap='gray') # plt.scatter(5507.75, 4227.01, s=3, c='red', marker='o') # plt.colorbar() # plt.show() image_data = fits.getdata('fits/1116985.fit', ext=0) print(image_data.shape) plt.figure() plt.imshow(np.log10(image_data), vmin = 2.5, vmax = 3.5, cmap='gray') # plt.imshow(np.log10(image_data), cmap='gray') # plt.imshow(np.log10(image_data), cmap='gray') # plt.imshow(image_data,vmin = 1500, vmax = 2000, cmap='gray') plt.scatter(5507.75, 4227.01, s=3, c='red', marker='o') plt.colorbar() plt.show()