GALFIT 3용 고급 표면 밝기 분석 도구 EllipSect

초록

EllipSect는 Python으로 구현된 GALFIT 3 출력 데이터를 바로 활용해 1차원 표면 밝기 프로파일, 비모수 반경(Petrosian, Kron 등) 및 총 유효 반경 등을 자동으로 계산하고, 고품질 그래프와 CSV 형태의 결과물을 제공하는 사용자 친화적인 분석 툴이다.

상세 분석

EllipSect는 GALFIT 3이 생성하는 FITS 데이터 큐브(원본 이미지, 모델 이미지, 잔차 이미지)와 파라미터 파일(galfit.nn)을 입력으로 받아, 사용자가 별도 변환 과정을 거칠 필요 없이 바로 분석을 수행한다. 핵심 알고리즘은 mgefit 라이브러리의 sectors_photometry 함수를 차용해 은하 중심을 기준으로 일정 각도 간격의 섹터를 정의하고, 각 섹터 내에서 지정된 최소 카운트 수준까지 반경을 확장하면서 평균 표면 밝기를 측정한다. 이때 GALFIT 3에서 제공된 마지막 컴포넌트(보통 지수형 디스크)의 타원율과 위치각을 기본값으로 사용하지만, 사용자는 명령행 옵션으로 모든 파라미터를 재정의할 수 있다.

표면 밝기 계산식은
μ_I = zp_t + 5·log(PlateScale) + 2.5·log(ExpTime) – 2.5·log(C_0) + 0.1
을 따르며, 여기서 zp_t는 매직 제로포인트, PlateScale는 픽셀당 각도, ExpTime은 노출 시간, C_0는 카운트‑픽셀 변환 계수이다. 이 식은 무한 aperture 보정(0.1 mag)까지 포함한다.

EllipSect는 GALFIT 3이 제공하지 않는 비모수 측정값을 추가로 산출한다. 주요 파라미터는 다음과 같다.

• 총 유효 반경 (r_te): 다중 컴포넌트 모델의 전체 광량을 합산한 뒤, 절반 광량에 해당하는 반경을 이분법으로 찾는다. 이는 단일 Sérsic 모델에서 r_e를 사용하는 기존 관행을 크게 개선한다.
• Petrosian, Kron 반경: 각각 Petrosian η와 Kron R_K 값을 이용해 정의된 반경을 자동 계산한다.
• Bumpiness (B): 관측 이미지와 모델 이미지 간의 표면 밝기 차이를 2 r_e 영역 내에서 RMS 형태로 정량화한다. 이는 Blakeslee et al. (2006)의 형태분류에 활용될 수 있다.
• Tidal parameter (T_gal): I(x,y)와 모델 M(x,y) 간의 차이를 정규화해 은하의 외부 교란 정도를 평가한다.
• AIC/BIC: 모델 복잡도와 χ²_ν를 결합해 최적 모델을 선택하도록 돕는다. 특히 BIC는 해상도 요소(N_res) 기반 변형도 제공한다.

스카이 배경 추정은 두 가지 방법을 제공한다. Sky gradient method는 Barden et al. (2012)를 따라 원형 링을 확장하며 평균 스카이 값을 구하고, 처음 두 번 양의 기울기가 나타날 때까지 진행한다. Random box method는 무작위 박스를 이미지 전역에 배치해 통계적 평균을 산출한다. 두 방법 모두 마스크 파일을 적용해 불필요한 픽셀을 배제한다.

출력은 (1) 주요 1D 표면 밝기 프로파일 그래프(주축 및 여러 방위각), (2) 원본‑모델‑잔차 PNG 합성 이미지, (3) CSV 형태의 정량적 파라미터 파일이다. 사용자는 –comp 옵션으로 개별 컴포넌트의 프로파일을 별도 표시하거나, 외부 플롯 라이브러리와 연동해 맞춤형 시각화를 구현할 수 있다. 또한 NED API를 호출해 은하의 거리 모듈러스, 은하 외부 소멸, Kpc 변환 계수를 자동으로 가져와 절대광도와 광도 등을 계산한다(단, k‑correction은 사용자가 직접 제공해야 함).

전체적으로 EllipSect는 GALFIT 3의 강력한 2D 피팅 능력을 보완해, 후처리 단계에서 반복적인 포맷 변환과 수동 계산을 없애고, 비모수 반경 및 통계적 적합도 지표를 일관되게 제공함으로써 대규모 은하 표면 밝기 연구와 파이프라인 자동화에 큰 가치를 제공한다.