Thanks for the reply. I had wondered how much better sigma clipping might be. Looks easy enough to implement.
I have wondered if the skewness of the sample pixels could be calculated and that used as a threshold for pixel removal – though I suspect that might be more computationally intensive than sigma clip.
I shall have a play on some data I’m working on.