Фотометрия / Апертурная фотометрия Как я могу найти ошибку на фотометрии? - PullRequest
Новая
       38

Фотометрия / Апертурная фотометрия Как я могу найти ошибку на фотометрии?

1 голос
/ 21 февраля 2020

Я выполняю апертурную фотометрию с использованием фотоутил. Я следовал инструкциям на странице с фотографиями, вычислял фон из круглого кольца вокруг источника и вычитал его. Теперь мне нужно вычислить погрешность фотометрии, но задача aperture_photometry не дает ошибок апертуры и фоновой ошибки, пока я не укажу ошибку ввода. Какую ошибку ввода я должен дать заданию? Для получения апертурной фотометрии я использовал следующий код: 

    import numpy as np
    import matplotlib.pyplot as plt
    from astropy.io import fits
    from matplotlib.colors import LogNorm
    from astropy.stats import sigma_clipped_stats
    from photutils import aperture_photometry
    from photutils import CircularAnnulus
    from photutils import CircularAperture

    image = fits.open('cut_F555_03.fits')
    hdu_list = image
    image_data = hdu_list[0].data

    import pandas as pd
    sources = pd.read_csv('daofind.tab', delimiter = '\s+', header = 0)


    loc = np.array([sources['xcentroid'], sources['ycentroid']])

    positions = np.transpose(loc)

    apertures = CircularAperture(positions, r = 8.12)
    rawflux = aperture_photometry(image_data, apertures)

    annulus_apertures = CircularAnnulus(positions, r_in = 8.12, r_out = 18.12)

    annulus_masks = annulus_apertures.to_mask(method = 'center')

    bkg_median = []
    for mask in annulus_masks:
        annulus_data = mask.multiply(image_data)
        annulus_data_1d = annulus_data[mask.data > 0]
        _, median_sigclip, _ = sigma_clipped_stats(annulus_data_1d)
        bkg_median.append(median_sigclip)
   bkg_median = np.array(bkg_median)

   rawflux['annulus_median'] = bkg_median / annulus_apertures.area

   rawflux['aper_bkg'] = bkg_median * apertures.area

   rawflux['final_phot'] = rawflux['aperture_sum'] - rawflux['aper_bkg']

Любая идея о том, как получить ошибки, будет очень полезна. С наилучшими пожеланиями. Sara

1 Ответ

0 голосов
/ 22 февраля 2020

Добро пожаловать в Сару!

Это хороший вопрос, и я несколько лет задавал себе вопрос go, когда работал с photutils. В есть раздел *1004* Оценка ошибок , в котором говорится:

Если и только если ключевое слово error введено в aperture_photometry () , возвращаемая таблица будет содержать столбец 'aperture_sum_err' в дополнение к 'aperture_sum'. 'aperture_sum_err 'обеспечивает распространенную неопределенность, связанную с 'aperture_sum'.

Это предполагает, что вы ранее получили общую ошибку и сохранили ее в массиве (через который вы проходите параметр error). Прочитайте документы, это поможет вам понять несколько вещей.

Я открыл вопрос в photutils чуть более двух лет go относительно этой топи c , Пользователь предложил сценарий для получения ошибок фотометрической апертуры так же, как это делает IRAF. К сожалению, я, видимо, не выполнил тестирование предоставленного сценария , поэтому я думаю, что теперь это ваше дело.

К сожалению, Nº 2. IRAF находится на смертном одре и большую часть ресурсов можно использовать чтобы узнать о том, как он выполнял некоторые задачи (например, апертурную фотометрию), теперь нет. Существует группа Facebook для (профессиональных) астрономов, которые могут оказаться полезными для такого рода вопросов.

...