【问题标题】:Find and Replace number to NA in array R在数组R中查找并将数字替换为NA
【发布时间】:2021-08-25 21:57:53
【问题描述】:

我有一些脑图像的 NifTi 文件,其中有很多我想用 NA 替换的零,但我不知道该怎么做。我在is.na() 函数的描述中读到:“通用函数 is.na

library(RNifti)

in.path <- "R_Things/Voxel/"
out.path <- "R_Things/Outputs/"

ids <- c('100','101','102','103','104')

for (i in ids) {

a <- readNifti(paste0(in.path, i, ".nii.gz"))

is.na(a) <- 0


writeNifti(image=a, file=paste0(out.path, i, "_with_NAs.nii.gz"))

}

任何关于我可以做些什么不同的想法都会非常有帮助,谢谢!

【问题讨论】:

标签: r arrays for-loop nifti


【解决方案1】:

这里有一些可以用NA替换0的选项

> replace(a, a == 0, NA)
[1] NA  1  2 NA  3 -1

> a * NA^(a == 0)
[1] NA  1  2 NA  3 -1

虚拟数据

> a <- c(0, 1, 2, 0, 3, -1)

【讨论】:

  • a * NA^(a == 0) 很迷人,因为它可以看到指数中的真值。
  • 感谢@ThomasIsCoding,但信不信由你。我想这可能与数据在数组中的事实有关?
【解决方案2】:

感谢@Baroque,这是我的最终代码

library(RNifti)

in.path <- "R_Things/Voxel/"
out.path <- "R_Things/Outputs/"

ids <- c('100','101','102','103','104')

for (i in ids) {

a <- readNifti(paste0(in.path, i, ".nii.gz"))

a[a==0] <- NA

writeNifti(image=a, file=paste0(out.path, i, "_with_NAs.nii.gz"))

}

【讨论】:

    猜你喜欢
    • 1970-01-01
    • 2021-07-02
    • 1970-01-01
    • 2018-09-11
    • 1970-01-01
    • 2020-10-26
    • 1970-01-01
    • 1970-01-01
    • 2019-05-27
    相关资源
    最近更新 更多