Hoi Han,
ik heb nog even naar je stukje Pascal code gekeken. Je hoeft een negatieve error niet anders af te handelen dan een positieve. De correcte workflow voor het testen op outliers en vinden van de uiteindelijke focuspositie is als volgt:
1) Stel de kappa-waarde in (2 of 3) en maak een inschatting van parameters a, b en Focus
2) Bepaal het verschil tussen ieder datapunt en de op a, b en Focus gebaseerde hyperbool (dit zijn de errors)
3) Bepaal de wortel van de kwadratensom van deze errors
4) [START ITERATIELOOP] Bepaal de hyperbolische best fit op basis van het minimaliseren van de kwadratensom door a, b en Focus aan te passen en stappen 2) en 3) te herhalen, dit resulteert in nieuwe waarden voor a, b en Focus
5) Bepaal de mediaan van de errors: MEDIAN(Error)
6) Bepaal voor ieder datapunt ABS(Error-MEDIAN(Error))
7) Bepaal de MAD: 1.4826*MEDIAN([alle gevonden waarden in 6)])
8) Bepaal voor ieder datapunt ABS(Error-Median(Error))/MAD
9) Toets voor ieder datapunt of de uitkomst van 8) de kappa-waarde overschrijdt: indien dit het geval is dan vervalt het datapunt
10) Herhaal vanaf stap 4) totdat er geen verwerpingen meer plaatsvinden.
11) Presenteer de gevonden waarden voor a en Focus (de waarde van b is minder interessant voor de gebruiker).
Nicolàs
https://www.dehilster.info/astronomy
In the observatory: Mount: 10Micron GM3000HPS, OTAs mounted: SW Esprit 80ED & Esprit 150ED, Celestron C11 XLT EdgeHD, Lunt LS80THA single stack, GTT60 (60mm aperture Galilean Type Telescope), Cameras: ZWO ASI1600MM Cool (2x), ASI174MM, ASI290MM & MC, QHYCCD QHY163M, OTAs on the ground: SW Explorer 300PDS, Bresser Messier 130/650 & 90/500.