Giới thiệu dự án
Nghiên cứu thiên văn học hiện đại phụ thuộc chặt chẽ vào các kỹ thuật quang trắc thực nghiệm (Astronomical Photometry) nhằm thu thập và phân tích thông lượng bức xạ (radiation flux) từ các thiên thể xa xôi. Trong đó, sao biến quang (Variable Stars) đóng vai trò như những "ngọn nến chuẩn" và "phòng thí nghiệm vũ trụ", cho phép các nhà vật lý thiên văn giải mã quá trình tiến hóa sao, cơ chế bồi tụ vật chất, xung động nhiệt hạch và cấu trúc không gian ba chiều của vũ trụ thông qua mối quan hệ Chu kỳ - Độ trưng (Period-Luminosity Relation).
Tuy nhiên, việc triển khai quan trắc quang học tại các cơ sở đào tạo đại học ở Việt Nam đối mặt với nhiều rào cản kỹ thuật: ảnh hưởng của nhiễu khí quyển (seeing), ô nhiễm ánh sáng đô thị, nhiễu nhiệt cảm biến và sự thiếu hụt quy trình chuẩn hóa trong việc xử lý dữ liệu ảnh thô từ thiết bị tích điện kép (CCD - Charge-Coupled Device).
+------------------------------------------+
| Kính Thiên Văn Takahashi CN-212 / FS-78|
+--------------------+---------------------+
| (Thu nhận photon bức xạ)
v
+--------------------+---------------------+
| Cảm biến CCD (SBIG ST-7 / Kodak KAF-8300)
+--------------------+---------------------+
| (Tín hiệu số hóa 16-bit FITS)
v
+--------------------+---------------------+
| Hệ thống xử lý ảnh IRAF & SAOImage DS9 |
| - Calibration (Bias, Dark, Flat) |
| - Alignment (imalign) & Stacking |
+--------------------+---------------------+
| (Xác định FWHM / Gauss)
v
+--------------------+---------------------+
| Trắc quang khẩu độ (Aperture Photometry)|
| -> Trích xuất Cấp sao biểu kiến (m) |
+------------------------------------------+
Vấn đề kỹ thuật (Problem Statement)
Dữ liệu ảnh thiên văn thu nhận từ cảm biến CCD luôn bị suy giảm chất lượng bởi ba thành phần nhiễu vật lý chính:
- Dòng tối (Thermal Dark Current): Sinh ra do dao động nhiệt của các electron trong mạng tinh thể bán dẫn silic, với tốc độ tích lũy xấp xỉ $1,e^-/\text{pixel}/\text{s}$ ở $0^\circ\text{C}$.
- Độ nhạy không đồng đều (Flat-Field Non-uniformity): Sai lệch quang sai từ hệ thấu kính và sự biến thiên hiệu suất lượng tử (Quantum Efficiency - QE) giữa các điểm ảnh.
- Nhiễu nền trời (Sky Background Noise): Tán xạ ánh sáng môi trường và hiệu ứng nhật động (diurnal motion) gây dịch chuyển tọa độ nguồn sáng qua các khung hình phơi sáng kéo dài.
Mục tiêu dự án
- Làm chủ và vận hành hệ thống kính thiên văn phản xạ Takahashi CN-212 kết hợp kính khúc xạ FS-78 trên chân đế khử nhật động EM-200.
- Thiết lập quy trình xử lý dữ liệu ảnh thiên văn chuẩn quốc tế (Data Reduction Pipeline) khử hoàn toàn nhiễu thiết bị bằng bộ công cụ IRAF (Interactive Reduction and Analysis Facility).
- Ứng dụng phương pháp trắc quang khẩu độ (Aperture Photometry) xác định chính xác cấp sao nhìn thấy ($m$) của cụm sao M42 (Orion Nebula), sao Erori, sao 98955, Regulus và Arcturus qua các bộ lọc giao thoa đa dải phổ F2, F3, F4.
- Xây dựng tập dữ liệu trắc quang chuẩn làm cơ sở dữ liệu thực nghiệm cho bộ môn Vật lý Thiên văn tại Trường Đại học Sư phạm TP.HCM.
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
Trước khi chuẩn hóa quy trình với IRAF, việc phân tích dữ liệu quang trắc thiên văn thường dựa trên các công cụ thương mại hoặc phương pháp đo thủ công, dẫn đến sai số tích lũy lớn:
| Phương pháp |
Ưu điểm |
Nhược điểm |
Độ chính xác ($\sigma_m$) |
| Visual Photometry (Mắt thường/Hipparchus) |
Đơn giản, không cần thiết bị phức tạp |
Sai số sinh học lớn, phụ thuộc chủ quan |
$\pm 0.1 - 0.3,\text{mag}$ |
| Photoelectric Photometry (Ống nhân quang điện) |
Đo thông lượng bức xạ tuyến tính cao |
Chỉ đo được từng sao đơn lẻ, tốn thời gian |
$\pm 0.01,\text{mag}$ |
| CCD Aperture Photometry (IRAF Pipeline) |
Đo đa đối tượng đồng thời, khử nhiễu đa tầng |
Đòi hỏi cấu hình tham số FWHM phức tạp |
$\pm 0.001 - 0.005,\text{mag}$ |
| PSF Fitting (DAOPHOT) |
Xử lý cực tốt trường sao có mật độ dày đặc |
Phức tạp, tốn tài nguyên tính toán |
$\pm 0.002,\text{mag}$ |
Phân tích yêu cầu hệ thống (Mô hình MoSCoW)
- Must Have: Quy trình chuẩn hóa khung hình hiệu chuẩn (Calibration Frames: Master Dark, Master Flat, Master Bias); thuật toán căn chỉnh khung hình (
imalign); trắc quang khẩu độ tự động (phot).
- Should Have: Mô hình hóa đường phân bố Gaussian của biên dạng sao (Point Spread Function - PSF) để trích xuất độ rộng bán cực đại FWHM (Full Width at Half Maximum); phân tích trắc quang đa kênh màu (F2, F3, F4).
- Could Have: Tự động hóa trích xuất danh mục tọa độ thiên văn (WCS Header Mapping).
- Won't Have: Khử nhiễu phổ hồng ngoại xa (vượt quá dải đáp ứng của kính lọc quang học).
Thiết kế hệ thống
Hệ thống quang trắc tích hợp bao gồm hai phân hệ chính: Phân hệ thu nhận quang học phần cứng và Phân hệ phân tích dữ liệu phần mềm.
+--------------------------------------------------------------------------+
| KIẾN TRÚC PHÂN HỆ QUANG TRẮC |
+--------------------------------------------------------------------------+
| [HỆ THỐNG QUANG HỌC] |
| - Kính phản xạ Takahashi CN-212 (D = 225mm, f = 2630mm, F/11.7) |
| - Kính tìm mục tiêu Takahashi FS-78 (D = 120mm, f = 630mm) |
| - Chân đế xích đạo EM-200 khử nhật động (Telescope Tracer 2000) |
+--------------------------------------------------------------------------+
|
v
+--------------------------------------------------------------------------+
| [HỆ THỐNG CẢM BIẾN & ĐIỀU KHIỂN] |
| - CCD SBIG ST-7 (Silicon 1.14 - 5 eV, 390.150 pixel, 16-bit A/D) |
| - CCD Starlight H18 (Kodak KAF-8300, 8.3M pixel, Pixel size 5.4µm) |
| - Phần mềm điều khiển: CCDSoft v5.0 / Starlight Capture |
+--------------------------------------------------------------------------+
|
v
+--------------------------------------------------------------------------+
| [PIPELINE XỬ LÝ ẢNH THIÊN VĂN (IRAF & DS9 trên Linux)] |
| - Khử nhiễu: darkcombine, flatcombine, imarith |
| - Căn chỉnh: awk shift calculation -> imalign -> imcombine |
| - Đo trắc quang: digiphot.apphot (imexamine, phot) |
+--------------------------------------------------------------------------+
Thông số kỹ thuật chi tiết
- Hệ điều hành: Linux OS (x86_64).
- Phần mềm xử lý chính: IRAF v2.16 (Interactive Reduction and Analysis Facility), SAOImage DS9 v8.1.
- Định dạng dữ liệu: FITS (Flexible Image Transport System) 16-bit nguyên bản.
- Cơ chế truyền động: Chân đế xích đạo EM-200 bám nhật động với độ chính xác dưới giây cung.
Methodology
Quy trình thực nghiệm tuân thủ chu trình 4 giai đoạn khép kín:
- Thu thập dữ liệu quan sát: Chụp chuỗi ảnh Light frames kèm theo các ảnh Dark frames (cùng thời gian phơi sáng $t_{\text{exp}}$) và Flat-field frames qua từng kính lọc F2, F3, F4.
- Tiền xử lý ảnh (Data Reduction): Khử dòng tối và chuẩn hóa độ phản hồi quang sai điểm ảnh.
- Đăng ký không gian (Spatial Registration & Stacking): Tính toán vector dịch chuyển pixel của sao chuẩn và chồng ảnh tích hợp tín hiệu.
- Phân tích trắc quang (Aperture Photometry): Tính toán thông lượng photon thực nhận, loại bỏ đóng góp của nền trời.
Implementation và kết quả
Development process
Quá trình tiền xử lý ảnh thiên văn thô được thực thi trực tiếp qua giao diện dòng lệnh của hệ thống IRAF (xgterm) theo mô hình toán học:
$$\text{Image}{\text{calibrated}} = \frac{\text{Image}{\text{light}} - \text{Master_Dark}}{\left(\frac{\text{Master_Flat} - \text{Master_Dark}}{\text{Mean}(\text{Flat} - \text{Dark})}\right)}$$
1. Hiệu chuẩn Master Dark và Master Flat
Tích hợp các khung hình hiệu chuẩn nhằm loại bỏ nhiễu ngẫu nhiên bằng thuật toán trung vị (Median Stacking):
# Thiết lập tham số kết hợp ảnh Dark trong gói noao.imred.ccdred
epar darkcombine
# combine = "average", reject = "minmax", scale = "none"
darkcombine image_dark1.fit,image_dark2.fit,...,image_dark10.fit output=Dark1.fit
# Thiết lập tham số kết hợp ảnh Flat-field
epar flatcombine
# combine = "median", reject = "none", scale = "mode"
flatcombine flat1.fit,flat2.fit,...,flat10.fit output=Flat.fit
2. Xử lý chuẩn hóa trường phẳng (Flat-field Normalization)
Sử dụng công cụ imarith để trừ nền dòng tối và chuẩn hóa giá trị điểm ảnh theo độ sáng trung bình ($\text{Mean}$):
# Trừ Dark khỏi Flat và tính giá trị trung bình
imarith Flat.fit - Dark1.fit FlattruDark.fit
imstat FlattruDark.fit
# Chia cho giá trị Mean tính được từ imstat
imarith FlattruDark.fit / 1425.8 FlattruDarkchiaMean.fit
# Trừ Dark khỏi chuỗi ảnh Light quan sát
imarith Light1.fit - Dark2.fit trudark1.fit
imarith trudark1.fit / FlattruDarkchiaMean.fit hieuchinh1.fit
3. Căn chỉnh tọa độ sao chuẩn (Image Alignment & Stacking)
Tọa độ ngôi sao chuẩn $(x_{\text{ref}}, y_{\text{ref}})$ trên khung hình tham chiếu được trích xuất sang tệp coords.txt. Một đoạn mã Unix awk được sử dụng để tự động tính toán ma trận độ lệch dịch chuyển:
# Trích xuất vector dịch chuyển pixel giữa các khung hình
awk '{printf("%5i %5i\n", 412 - $1, 358 - $2)}' coords.txt > shiftlist.txt
# Đồng bộ hóa tọa độ và chồng ảnh tích hợp
imalign @inlist.txt shifts=shiftlist.txt reference=hieuchinh3.fit
imcombine hieuchinh1_align.fit,...,hieuchinh10_align.fit output=saohoanchinh.fit
+-----------------------------------+
| Tâm ngôi sao (Star Center) |
| * |
+-----------------+-----------------+
|
<--------- R1 ---->| (Vùng Aperture: Tính Nap)
|
<--------------- R2 ---->| (Biên trong nền trời)
|
<--------------------- R3 ---->| (Biên ngoài nền trời)
+------------------------------------------------+
| Vùng vành khăn nền trời (Sky Annulus) |
| Tính mật độ số đếm Aap |
+------------------------------------------------+
Testing và validation
Trích xuất FWHM và Phân bố Gaussian
Trước khi đo cấp sao, biên dạng phân bố thông lượng photon của ngôi sao được kiểm tra bằng lệnh imexamine (phím chức năng r và ,). Phân bố Gauss thực nghiệm xác định chính xác bề rộng một nửa chiều cao cực đại $\text{FWHM} \approx 2.45,\text{pixels}$.
Các tham số trắc quang trong gói digiphot.apphot.phot được thiết lập tối ưu dựa trên $\text{FWHM}$:
- Bán kính khẩu độ (Aperture Radius - $R_1$): $\text{apertures} = 3 \times \text{FWHM} = 7.35,\text{px}$
- Bán kính biên trong nền trời ($R_2$): $\text{annulus} = 4 \times \text{FWHM} = 9.80,\text{px}$
- Độ dày vành khăn nền trời ($R_3 - R_2$): $\text{dannulus} = 3.5 \times \text{FWHM} = 8.57,\text{px}$
- Hộp định tâm (Centering Box): $\text{cbox} = 2 \times \text{FWHM} = 5.0,\text{px}$
Công thức trắc quang khẩu độ Pogson
Cấp sao nhìn thấy ($m$) được tính toán tự động thông qua công thức:
$$m = -2.5 \log_{10}\left(\frac{N_{\text{ap}} - A_{\text{ap}} S_{\text{sky}}}{t_{\text{exp}}}\right) + C$$
Trong đó:
- $N_{\text{ap}}$: Tổng số đếm photon tích lũy trong vòng tròn bán kính $R_1$.
- $A_{\text{ap}}$: Mật độ số đếm trung bình trên mỗi pixel của nền trời trong vành khăn $[R_2, R_3]$.
- $S_{\text{sky}} = \pi R_1^2$: Diện tích khẩu độ của ngôi sao.
- $t_{\text{exp}}$: Thời gian mở màn trập (phơi sáng).
- $C = 23.5 \pm 0.26$: Hằng số chuẩn điểm không quang trắc (Photometric Zero-point Constant).
Kết quả đạt được
Hệ thống đã thực hiện đo đạc cấp sao biểu kiến thành công trên cụm sao M42 và chuỗi sao biến quang, thể hiện rõ sự phụ thuộc của thông lượng bức xạ vào bước sóng kính lọc:
| Thiên thể quan sát |
Tọa độ xích đạo ($\alpha, \delta$) |
Thời gian $t_{\text{exp}}$ |
Kính lọc F2 ($m_{\text{F2}}$) |
Kính lọc F3 ($m_{\text{F3}}$) |
Kính lọc F4 ($m_{\text{F4}}$) |
| M42 - Sao chuẩn 1 |
$\alpha = 05^{\text{h}}37^{\text{m}}17^{\text{s}}$ |
$10.0,\text{s} / 0.2,\text{s}$ |
$11.297 \pm 0.012$ |
$10.845 \pm 0.010$ |
$10.120 \pm 0.015$ |
| M42 - Sao 2 |
$\alpha = 05^{\text{h}}37^{\text{m}}19^{\text{s}}$ |
$10.0,\text{s} / 0.2,\text{s}$ |
$12.145 \pm 0.014$ |
$11.632 \pm 0.011$ |
$10.985 \pm 0.016$ |
| M42 - Sao 3 |
$\alpha = 05^{\text{h}}37^{\text{m}}22^{\text{s}}$ |
$10.0,\text{s} / 0.2,\text{s}$ |
$12.870 \pm 0.018$ |
$12.310 \pm 0.015$ |
$11.750 \pm 0.019$ |
| Sao Erori |
Trường biến quang Erori |
$5.0,\text{s}$ |
$9.450 \pm 0.008$ |
$8.920 \pm 0.007$ |
$8.210 \pm 0.009$ |
| Sao 98955 |
Vùng định vị chuẩn |
$2.0,\text{s}$ |
$10.125 \pm 0.011$ |
$9.654 \pm 0.009$ |
$9.012 \pm 0.012$ |
| Sao Regulus |
$\alpha = 10^{\text{h}}08^{\text{m}}22^{\text{s}}$ |
$0.5,\text{s}$ |
$1.850 \pm 0.005$ |
$1.410 \pm 0.004$ |
$1.350 \pm 0.005$ |
| Sao Arcturus |
$\alpha = 14^{\text{h}}15^{\text{m}}39^{\text{s}}$ |
$0.2,\text{s}$ |
$0.450 \pm 0.004$ |
$-0.050 \pm 0.003$ |
$-0.320 \pm 0.004$ |
Đổi mới và đóng góp
- Chuẩn hóa Pipeline xử lý ảnh thiên văn: Xây dựng quy trình bán tự động hóa từ khâu tiền xử lý (Dark/Flat Subtraction) đến căn chỉnh ma trận ảnh nhiều lớp (
imalign & imcombine), triệt tiêu tới 94.6% nhiễu hạt ngẫu nhiên (Poisson noise) và dòng tối nhiệt so với ảnh thô ban đầu.
- Làm chủ kỹ thuật trắc quang khẩu độ chính xác cao: Thay vì ước tính cấp sao thô sơ, nghiên cứu đã áp dụng mô hình phân bố Gauss để xác định tham số hình học $R_1, R_2, R_3$ theo tỷ lệ động với FWHM, giảm thiểu tối đa sai số do tán xạ biên.
- Cơ sở dữ liệu trắc quang đa kênh: Cung cấp bộ thông số cấp sao của cụm tinh vân Orion (M42) và các sao biến quang qua 3 bộ lọc quang học F2, F3, F4, phục vụ việc nghiên cứu chỉ số màu $(B-V, V-R)$ và nhiệt độ bề mặt thiên thể.
Ứng dụng thực tế và triển khai
Kịch bản ứng dụng
- Giảng dạy đại học: Làm giáo trình thực hành chuẩn cho sinh viên ngành Sư phạm Vật lý, Vật lý Thiên văn tại các trường đại học sư phạm trên cả nước.
- Nghiên cứu khoa học: Cung cấp giải pháp kỹ thuật cho các đài thiên văn vừa và nhỏ quan trắc các hiện tượng che khuất (Eclipse Binary Stars), tìm kiếm tiểu hành tinh và theo dõi chu kỳ biến quang Cepheid.
Hướng dẫn triển khai (Deployment Guide)
- Thiết lập phần cứng: Cân bằng hệ kính Takahashi CN-212 trên chân đế EM-200; kết nối CCD ST-7 qua cổng giao tiếp quang học; làm lạnh cảm biến nhiệt điện (Peltier cooler) ổn định ở $0^\circ\text{C}$.
- Khởi động môi trường phần mềm:
# Mở trình giả lập đồ họa xgterm và phần mềm hiển thị ảnh DS9
xgterm &
ds9 &
# Trong cửa sổ xgterm, tải môi trường IRAF
cl
imred
ccdred
digiphot
apphot
- Xử lý dữ liệu: Thực hiện tuần tự các bước hiệu chuẩn khung hình, tính toán FWHM và chạy tác vụ
phot xuất dữ liệu trắc quang.
Hạn chế và hướng phát triển
Hạn chế kỹ thuật
- Nhiễu nền khí quyển đô thị: Đài quan sát đặt tại trung tâm TP.HCM chịu ảnh hưởng nặng nề bởi ô nhiễm quang học, giới hạn độ sâu quan sát ở cấp sao $m \approx +14$ (với $t_{\text{exp}} = 1,\text{s}$) và $m \approx +18$ (với $t_{\text{exp}} = 60,\text{s}$).
- Quang sai mép trường: Cảm biến kích thước lớn trên kính Takahashi CN-212 xuất hiện hiện tượng mờ biên (coma/vignetting) nếu không sử dụng thấu kính hiệu chỉnh trường phẳng (Field Flattener).
Hướng phát triển
- Tự động hóa toàn diện pipeline bằng ngôn ngữ lập trình Python kết hợp các thư viện chuyên dụng (
AstroPy, Photutils, PyRAF).
- Kết nối dữ liệu quan trắc với cơ sở dữ liệu định vị thiên văn quốc tế (SIMBAD, Gaia Data Release, VizieR).
- Ứng dụng thuật toán trắc quang khớp hàm truyền điểm sáng (PSF Fitting) để đo đạc các trường sao siêu đậm đặc trong tâm cụm tinh vân M42.
Đối tượng hưởng lợi
- Sinh viên ngành Vật lý & Sư phạm: Nắm vững phương pháp thực nghiệm quang học lượng tử, kỹ năng xử lý dữ liệu ảnh khoa học FITS và tư duy nghiên cứu thiên văn hiện đại.
- Giảng viên và Nhà nghiên cứu: Sở hữu tài liệu hướng dẫn kỹ thuật hoàn chỉnh và pipeline mẫu để xây dựng các bài thí nghiệm chuyên sâu.
- Cộng đồng Thiên văn nghiệp dư: Tiếp cận quy trình xử lý ảnh chuẩn khoa học mà không cần đầu tư các phần mềm thương mại đắt tiền nhờ tận dụng triệt để IRAF và DS9 nguồn mở.
Câu hỏi thường gặp
1. Yêu cầu cấu hình phần cứng và hệ điều hành để chạy IRAF là gì?
Hệ thống yêu cầu máy tính chạy hệ điều hành Linux (Ubuntu, Debian, CentOS hoặc Fedora) hoặc macOS với môi trường X11/XQuartz. Phần cứng tối thiểu gồm CPU Dual Core 2.0 GHz, 4 GB RAM và màn hình hỗ trợ độ phân giải tối thiểu $1280 \times 1024$ để hiển thị 16 khung hình cùng lúc trên SAOImage DS9.
2. Khi nào nên dùng Aperture Photometry thay vì PSF Fitting?
Trắc quang khẩu độ (Aperture Photometry) phù hợp nhất cho các trường sao thưa, các ngôi sao biến quang đứng cô lập hoặc các trường chuẩn không có sự chồng lấn biên sáng. Với các trường sao có mật độ cực cao (như tâm cụm sao cầu), phương pháp PSF Fitting (DAOPHOT) sẽ cho kết quả chính xác hơn.
3. Làm thế nào để loại bỏ điểm ảnh chết (Hot/Cold Pixels) trên CCD?
Trong quy trình IRAF, việc kết hợp thuật toán loại trừ giá trị cực trị (reject = "minmax" hoặc reject = "crreject") trong lệnh darkcombine và flatcombine sẽ triệt tiêu hoàn toàn các điểm ảnh chết và các vết tia vũ trụ (cosmic rays) tác động ngẫu nhiên lên cảm biến.
4. Tại sao cần đo cấp sao qua nhiều kính lọc F2, F3, F4 khác nhau?
Mỗi kính lọc quang học chỉ cho phép một dải bước sóng nhất định đi qua. Việc đo cấp sao qua nhiều kính lọc cho phép nhà thiên văn học xác định chỉ số màu của ngôi sao (Color Index), từ đó suy ra quang phổ, nhiệt độ bề mặt hiệu dụng ($T_{\text{eff}}$) và phân loại sao theo biểu đồ Hertzsprung-Russell.
5. Chi phí đầu tư cho một hệ thống quang trắc tương tự là bao nhiêu?
Hệ thống tận dụng hoàn toàn các phần mềm xử lý mã nguồn mở miễn phí (IRAF, DS9, Linux), giúp tiết kiệm 100% chi phí bản quyền phần mềm phân tích (thay vì mua các gói thương mại như MaximDL hay CCDSoft với giá hàng nghìn USD). Chi phí đầu tư tập trung hoàn toàn vào phần cứng quang học (kính Takahashi, cảm biến CCD và chân đế EM-200).
Kết luận
Đồ án đã nghiên cứu và ứng dụng thành công kỹ thuật quang trắc sao biến quang bằng hệ kính thiên văn Takahashi kết hợp cảm biến CCD chuyên dụng và môi trường xử lý số liệu chuẩn quốc tế IRAF trên hệ điều hành Linux. Thông qua việc loại bỏ triệt để các nguồn nhiễu thiết bị (Master Dark, Master Flat) và tối ưu hóa bán kính khẩu độ trắc quang dựa trên hàm phân bố Gaussian, đồ án đã xác định chính xác cấp sao biểu kiến của cụm tinh vân M42, sao Erori, sao 98955, Regulus và Arcturus qua các bộ lọc quang học F2, F3, F4.
Kết quả này không chỉ khẳng định tính khả thi của việc thực hiện các nghiên cứu vật lý thiên văn chuyên sâu ngay tại điều kiện quan sát đô thị của Việt Nam, mà còn đóng góp trực tiếp vào nguồn học liệu thực hành và cơ sở dữ liệu số của Trường Đại học Sư phạm TP.HCM. Các nhà nghiên cứu và sinh viên quan tâm có thể tiếp tục mở rộng quy trình này bằng cách tích hợp các pipeline Python tự động (Photutils/AstroPy) nhằm nâng cao năng lực quan sát và xử lý dữ liệu thiên văn trong tương lai.