Astropy: FWHM négatif dans Moffat1D

Créé le 5 juin 2019  ·  3Commentaires  ·  Source: astropy/astropy

Je fais une demande pour obtenir la meilleure mise au point pour un instrument. J'utilise un modèle Moffat1D pour ajuster le profil et pour une raison quelconque, le gamma est négatif (ce qui n'a pas d'importance et je suppose que c'est correct) mais quand j'obtiens le FWHM est également négatif, ce qui, à ma connaissance, ne fait pas sens dans les images astronomiques (ou autre). Il me semble qu'il y a une erreur concernant la formule, voir ce commentaire : https://github.com/astropy/astropy/pull/6027#issuecomment -300380649

de l'avant-dernière à la dernière ligne
CodeCogsEqn
et donc (devrait être)
CodeCogsEqn_fwhm

et en utilisant la notation dans la documentation d'astropie, le FWHM devrait être :

CodeCogsEqn_astropy_fwhm

J'espère que quelqu'un a une meilleure idée de cela, je vais contourner ce problème en faisant toujours une valeur absolue de gamma .

Je peux faire un PR si c'est vraiment quelque chose qui en vaut la peine.

tracé et modèle de référence

moffat_fitted

Le modèle

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

Tous les 3 commentaires

@ysBach , qu'en pensez-vous ?

@pllim Désolé de ne pas avoir répondu pendant quelques jours.
Je pense que @simontorres a tout à fait raison.

En raison de la mise au carré, le signe du paramètre "noyau" ( gamma de l'astropie = alpha de l'IRAF = R du travail original de Moffat) n'a pas d'importance dans le calcul du profil de Moffat, mais c'est le cas pour le calcul FWHM actuel...

@simontorres , la pull request est la bienvenue. Veuillez vous référer à http://docs.astropy.org/en/stable/development/workflow/development_workflow.html . Merci!

Cette page vous a été utile?
0 / 5 - 0 notes