Я работаю с данными радиообследований. Мои изображения представляют собой файлы .fits с BUNIT Jy/beam. Я хочу рассчитать чистый поток от непрерывного источника (не точечного источника). Когда я просто сложил все значения интенсивности и перевел в Jy, используя параметры луча, мой профессор сказал, что этого недостаточно. Он сказал: «Вы должны учитывать шум, принимая уровень 3 или 4 сигма в зависимости от распределения шума». Распределение шума (в местах, отличных от источника) почти гауссово, как я могу видеть с помощью специального средства просмотра.
Под «взять уровень 3 сигма» он имеет в виду, что я должен подогнать гауссиан к шуму, найти сигму, а затем выбрать пиксели со значениями интенсивности больше, чем 3 * сигма? Если я сделаю это, очевидно, я потеряю много пикселей. Кроме того, если это правильно, должен ли я вычесть среднее значение интенсивности «шумовых» пикселей из интенсивности источника?
Среднеквадратичное значение шума указано в заголовке с предупреждением о том, что оно сильно варьируется в зависимости от некоторых источников. Связано ли это среднеквадратичное значение шума с сигмой, которую я вычислил ранее?
Может быть, слишком поздно, но, вероятно, он имеет в виду учитывать только пиксели с интенсивностью, превышающей 3 * сигма. Если нет, вы будете считать случайный шум континуумом. Очевидно, вы потеряете много пикселей, но эти пиксели вообще не несут никакой информации.
Профессор также может иметь в виду вычесть потоки фоновой области для вычитания шума. Затем, вычисляя статистику для вычитаемого значения потока и ошибки измерения, скажем, A +/- dA. Затем, применяя критерий 3-сигма, определить, обнаружен ли источник радиоизлучения в этом месте, т. е. если A - 3*dA > 0, источник есть, в противном случае источника нет.
Если источника нет, вместо этого вы захотите указать значение для верхнего предела 3 сигма, например, A + 3*dA.
Также обратите внимание, что для слабого источника регулярная пуассоновская или гауссовская статистика не являются хорошими репрезентативными. В рентгене мы используем статистику наличных денег для случаев с небольшим числом случаев. Я думаю, что та же статистика может быть применима к радионаблюдениям.
StephenG - Помощь Украине