I have some gridded data that i am Z transforming, and an inherent part of the data are low frequency but very high/low values. This is absolutely expected in the context of the data. (It is quite normally distributed)
When i calculate the traditional Z-Score (x - mean/sd) and the Modified Z-Score (using
spatialEco::raster.Zscore in R), the latter is removing outliers.
The resulting surfaces differ quite markedly; Z-scores for normal method are ca. 80% higher than for modified Z. This pushes many values above 2.
Is this purely about whether i want to keep outliers?