GEUS-Glaciology-and-Climate/pypromice

Latest coordinates shifted by 200m at DY2 compared to station location

BaptisteVandecrux opened this issue · 3 comments

red dots are monthly averaged transmitted positions
image

During the bufr processing, a linear regression is fitted to the available transmitted coordinates. Unfortunately, this regression and the OLS minimization is to sensitive to outliers. E.g. these outliers:
image

... lead to this linear regression:
image

Switching to a RANSAC regressor leads to more robust estimation:
image

More testing on the way before I make this update.

Examples where RANSAC don't work:
image

image

image

another problem at KAN_M where only coordinates from SBD messages are available:
image

For the moment I am switching off the gap-filling with SBD coordinates and using SBD average when no coordinates are available.

This is now fixed:
billede