Tengo un archivo de ajustes del galaxy m87 producido por SCUBA 2 y me pregunto cómo obtengo los flujos a 450 y 850 micras. la imagen de salida da magnitudes de píxeles en mJy/beam donde quiero encontrar el flujo eh para una parte determinada de la galaxia. Llegué a la cima de los valores de todos los píxeles dentro de una apertura, pero no estoy seguro de adónde ir desde allí para obtener el flujo, ya que hay más de un píxel dentro del FWHM. Digamos que tengo una apertura de radio de 1,5 píxeles que da una suma de 6985,04885203 mJy/haz con un tamaño de píxel de 3,99999999 segundos de arco y un tamaño de haz de 13 segundos de arco. ¿Cómo obtengo el flujo para esa apertura en esa longitud de onda?
En esencia, lea el valor de flujo máximo y, si es una fuente puntual, será el mismo que su flujo integrado.
para profundizar en esto:
En esta pregunta, su fuente parece ser comparable a su haz. así que lo trataría como una fuente puntual. - tenga en cuenta que es posible que también desee considerar que muchas personas aplican un "filtro de coincidencia" a sus datos si se trata de una emisión de fuente puntual débil - esto suaviza el mapa - aumenta el tamaño efectivo del haz dentro del mapa y también puede haber alguna pérdida de flujo ( ver por ejemplo: http://adsabs.harvard.edu/cgi-bin/bib_query?arXiv:1707.00990 )
Verá que esto (mJy/beam to mJy) también se discute aquí: http://www.eaobservatory.org/jcmt/help/
Para una fuente puntual real, un valor de brillo máximo informado en unidades de mJy es lo mismo que un valor de brillo máximo informado en mJy/haz.
Pero, ¿qué ocurre si tenemos un mapa en mJy/haz y queremos obtener un valor de intensidad integrado, un valor de flujo total? Primero sumamos una cantidad de píxeles y ahora queremos obtener nuestras unidades correctas de mJy/beam a mJy...
Flujo total = flujo sumado sobre un número de píxeles/(número de píxeles en un haz)
Entonces tus unidades son:
[mJy*píxeles/haz] / [píxeles/haz] = [mJy].
Ahora sabemos que para una Gaussiana:
Área del haz = 2 × π × σ² [arcsec]
donde el σ del haz gaussiano se puede calcular a partir de los valores JCMT FWHM a 850 y 450 micrones (recuerde que los componentes del haz se proporcionan en el artículo de Dempsey de 2013).
FWHM = 2 σ √(2 ln 2) [seg de arco]
Entonces podemos usar el FWHM para obtener σ para calcular el área del haz e informar el área del haz en términos de píxeles:
número de píxeles en un haz = Área del haz [arcsec] / (longitud de píxel)²
Finalmente, debo agregar que también querrá hacer una corrección de apertura: consulte el documento de calibración de Dempsey et al 2013 si está utilizando una apertura de un tamaño diferente a 60 ".
Carlos Witthoft