Astropy: FWHM negativo en Moffat1D

Creado en 5 jun. 2019  ·  3Comentarios  ·  Fuente: astropy/astropy

Estoy haciendo una solicitud para obtener el mejor enfoque para un instrumento. Estoy usando un modelo Moffat1D para ajustar el perfil y, por alguna razón, la gamma es negativa (lo cual no importa y supongo que es correcto) pero cuando obtengo el FWHM también es negativo, lo que hasta donde yo sé no hace sentido en imágenes astronómicas (o cualquiera). Me parece que hay un error con respecto a la fórmula, vea este comentario: https://github.com/astropy/astropy/pull/6027#issuecomment -300380649

desde la penúltima hasta la última línea
CodeCogsEqn
y así (debería ser)
CodeCogsEqn_fwhm

y usando la notación en la documentación de Astropy, el FWHM debería ser:

CodeCogsEqn_astropy_fwhm

Con suerte, alguien tiene una mejor idea de esto, lo solucionaré haciendo siempre un valor absoluto de gamma .

Puedo hacer un PR si realmente vale la pena hacer esto.

diagrama y modelo de referencia

moffat_fitted

El modelo

Model: Moffat1D
Inputs: ('x',)
Outputs: ('y',)
Model set size: 1
Parameters:
      amplitude       x_0          gamma         alpha    
    ------------- ------------ -------------- ------------
    248.392233095 2065.9401953 -90.3411844413 1161.3273732
Bug modeling

Todos 3 comentarios

@ysBach , ¿qué opinas?

@pllim Perdón por no haber respondido durante unos días.
Creo que @simontorres tiene toda la razón.

Debido a la cuadratura, el signo del parámetro "núcleo" ( gamma de astropía = alpha de IRAF = R del trabajo original de Moffat) no importa en el cálculo del perfil de Moffat, pero lo hace para el cálculo actual de FWHM ...

@simontorres , la solicitud de extracción es bienvenida. Consulte http://docs.astropy.org/en/stable/development/workflow/development_workflow.html . ¡Gracias!

¿Fue útil esta página
0 / 5 - 0 calificaciones