Я новичок в использовании R. У меня есть 2 изображения FITS, и я хочу использовать одно в качестве маски для другого. На практике я хочу сохранить значения изображения B в формате FITS выше определенного значения изображения A в формате FITS (и поместить остальные пиксели изображения B в NaN).
До сих пор я мог читать изображение FITS и тренировать маску. Я хотел бы знать, есть ли способ сделать это, избегая использования двойной петли. Спасибо
### READ FITS File ################
Mono <- readFITS("Mono.fits")
Vr <- readFITS("Vr.fits")
### Put the FITS image in a matrix #####
TempMono <- as.matrix(Mono$imDat)
TempVr <- as.matrix(Vr$imDat)
TempMono[TempMono < 40.0] <- 0.0
TempMono[TempMono > 40.0] <- 1.0