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
y así (debería ser)
y usando la notación en la documentación de Astropy, el FWHM debería ser:
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.
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
@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!