Astropy: Negative FWHM in Moffat1D

Erstellt am 5. Juni 2019  ·  3Kommentare  ·  Quelle: astropy/astropy

Ich stelle eine Bewerbung um den besten Fokus für ein Instrument zu bekommen. Ich verwende ein Moffat1D-Modell, um das Profil anzupassen, und aus irgendeinem Grund ist Gamma negativ (was egal ist und ich gehe davon aus, dass es richtig ist), aber wenn ich das FWHM bekomme, ist es auch negativ, was meines Wissens nicht ausreicht Sinn in astronomischen Bildern (oder anderen). Mir scheint, dass es einen Fehler in der Formel gibt, siehe diesen Kommentar: https://github.com/astropy/astropy/pull/6027#issuecomment -300380649

von der vorletzten bis zur letzten Zeile
CodeCogsEqn
und so (sollte)
CodeCogsEqn_fwhm

und unter Verwendung der Notation in der Dokumentation von Astropy sollte das FWHM sein:

CodeCogsEqn_astropy_fwhm

Hoffentlich hat jemand einen besseren Einblick darin, ich werde dies umgehen, indem ich gamma immer absoluten Wert mache.

Ich kann eine PR machen, wenn sich das wirklich lohnt.

Referenzplot und Modell

moffat_fitted

Das Model

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

Alle 3 Kommentare

@ysBach , was denkst du?

@pllim Entschuldigung, dass ich einige Tage nicht geantwortet habe.
Ich denke @simontorres hat absolut recht.

Aufgrund der Quadrierung spielt das Vorzeichen des "Kern"-Parameters ( gamma von Astropie = alpha von IRAF = R von Moffats Originalarbeit) bei der Moffat-Profilberechnung keine Rolle, aber für die aktuelle FWHM-Berechnung...

@simontorres , Pull-Request ist willkommen. Siehe http://docs.astropy.org/en/stable/development/workflow/development_workflow.html . Dankeschön!

War diese Seite hilfreich?
0 / 5 - 0 Bewertungen