När det gäller FWHM så kan man tex ta de pixlar som ligger horisontellt utmed stjärnans mittposition, men då använder man ju inte all information i bilden. Andra produkter som tex Mira Pro gör en radial profile + gaussian curve fitting för att estimera FWHM. Jag använder också alla pixlar, men har hittills fuskat med curve fitting. Jag kör en mindre exakt, men snabb approximation istället.
I grova drag så gör mitt program följande:
- Läs in råbild (TIFF, CR2 etc)
- Konvertera bilden till greyscale
- Gaussian blur -> Threshold bild där gränsen är median + bias
- Skapa lista över stjärnor i threshold bilden
- Sortera bort överexponerade stjärnor
- Sortera bort stjärnor som ligger för nära bildkanten (kommer orsaka problem senare)
- Sätt estimerat centrum för varje stjärna utifrån maxvärde i blur bild (pixel noggrannhet)
- Utvidga stjärnans radie en pixel i taget, avbryt då median värdet för cirkeln inte minskar nämnvärt längre
- Beräkna exakt mittpunkt för stjärnan utifrån föregående yta med "center of mass" algoritm (sub pixel)
- Beräkna FWHM i pixels och i arcsec per stjärna (se bild nedan)
- FWHM för bilden = Median värdet av alla FWHM
![glad :)](./images/smilies/icon_e_smile.gif)
Hur som helst så stämmer kurvornas utseende överlag så det är jag nöjd med.