Wrote new python script to compute averages from a set of images to r…

[Gab-D-G]

…eplace the currently slow median implementation.

[gdevenyi]

One more request, please add the effcient trimean as an estimator: https://en.wikipedia.org/wiki/Trimean

R code:

efficient_trimean <- function(X){
  # Boosted effciency (sensitivity) trimean, average of 20th, 50th and 80th percentiles
  a <- quantile(X, c(0.2,0.5,0.8))
  result <- mean(a)
  return(result)
}

[Gab-D-G]

is this alright?

elif opts.method == 'efficient_trimmed_mean':
    # computes the average from the 20th,50th and 80th percentiles
    average = np.quantile(concat_array, (0.2,0.5,0.8),axis=0).mean(axis=0)

[gdevenyi]

Yes, except its not a trimmed mean, its a "trimean". Please include the wikipedia link in the comment.