Displaying 1 result from an estimated 1 matches for "mediansarr".
Did you mean:
meansarr
2006 Aug 09
2
Speeding indexing and sub-sectioning of 3d array
...(k - radius):(k + radius)] )
Putting that code in 3 for loops
calculateKMedian <- function( arr, radius){
for( i in (radius + 1):(dim(arr)[1] - radius - 1) ){
for( j in (radius + 1):(dim(arr)[2] - radius - 1) )
for( k in (radius + 1):(dim(arr)[3] - radius - 1) ){
mediansArr <- mean( arr[ (i - radius):(i + radius),
(j - radius):(j + radius),
(k - radius):(k + radius)] )
}
}
return(mediansArr)
}
Results in a very slow run:
> system.time( calculateKMedian( a, 3))
[1] 423.468 0.096 423.63...