From 92271302c179b17bd08d400bd13a752973d6effd Mon Sep 17 00:00:00 2001
From: stnava <stnava@gmail.com>
Date: Mon, 12 Aug 2024 10:15:04 -0400
Subject: [PATCH] BUG: wrong type for idx

---
 R/makePowersPointsImage.R | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/R/makePowersPointsImage.R b/R/makePowersPointsImage.R
index 7c5e1c41..ecfb2918 100644
--- a/R/makePowersPointsImage.R
+++ b/R/makePowersPointsImage.R
@@ -31,7 +31,7 @@ makePowersPointsImage <- function(pts, mask, radius = 5) {
       for (j in seq(-n[2], n[2], by = 0.5)) {
         for (k in seq(-n[3], n[3], by = 0.5)) {
           local <- idx + c(i, j, k)
-          localpt <- antsTransformIndexToPhysicalPoint(mask, local)
+          localpt <- antsTransformIndexToPhysicalPoint(mask, as.numeric(local))
           dist <- sqrt(sum((localpt - pt) * (localpt - pt)))
           inImage <- (prod(idx <= dim(mask)) == 1) && (length(which(idx < 1)) == 0)
           if ((dist <= rad) && (inImage == TRUE)) {
@@ -83,7 +83,7 @@ makePointsImage <- function(pts, mask, radius = 5) {
         if (dim == 3) {
           for (k in seq(-n[3], n[3], by = 0.5)) {
             local <- idx + c(i, j, k)
-            localpt <- antsTransformIndexToPhysicalPoint(mask, local)
+            localpt <- antsTransformIndexToPhysicalPoint(mask, as.numeric(local))
             dist <- sqrt(sum((localpt - pt) * (localpt - pt)))
             inImage <- all(local >= 1 & local <= dim(powersLabels))
             if ((dist <= rad) && (inImage == TRUE)) {
@@ -95,7 +95,7 @@ makePointsImage <- function(pts, mask, radius = 5) {
         } # if dim == 3
         if (dim == 2) {
           local <- idx + c(i, j)
-          localpt <- antsTransformIndexToPhysicalPoint(mask, local)
+          localpt <- antsTransformIndexToPhysicalPoint(mask, as.numeric(local))
           dist <- sqrt(sum((localpt - pt) * (localpt - pt)))
           inImage <- (prod(local <= dim(mask)) == 1) && (length(which(local < 1)) == 0)
           if ((dist <= rad) && (inImage == TRUE)) {