antsAverageImages error

[dorianps]

It doesn't seem the variable "template" is needed for anything.

average = antsAverageImages(trn.imgCon, normalize = T)
Error: object 'template' not found

function (imageList, normalize = FALSE)
{
avg <- imageList[[1]] * 0
for (i in imageList) {
if (normalize) {
i <- i/mean(i)
}
template <- template + i
}
avg <- avg/length(imageList)
return(avg)
}
<environment: namespace:ANTsR>

[dorianps]

Corrected:

antsAverageImages = function (imageList, normalize = FALSE)
{
avg <- imageList[[1]] * 0
for (i in imageList) {
if (normalize) {
i <- i/mean(i)
}
avg <- avg + i
}
avg <- avg/length(imageList)
return(avg)
}

[bkandel]

If the help file for this function had a reproducible example, this bug would have been caught automatically by Travis.

[stnava]

agreed - the example should be added

brian

On Tue, Mar 24, 2015 at 10:22 AM, bkandel notifications@github.com wrote:

If the help file for this function had a reproducible example, this bug
would have been caught automatically by Travis.

—
Reply to this email directly or view it on GitHub
#35 (comment).

[dorianps]

Normalization step was missing in fix above:

antsAverageImages = function (imageList, normalize = FALSE)
{
avg <- imageList[[1]] * 0
for (i in imageList) {
if (normalize) {
i <- i/mean(i)
}
avg <- avg + i
}
avg <- avg/length(imageList)
if (normalize) avg = avg %>% iMath('Normalize')
return(avg)
}

[bkandel]

The fix I pushed took care of this: 3e97beb

[dorianps]

Good, cause mine above is not good enough, needed another fix.