我是使用R的新手。我有2张FITS图像,我想将其中一个用作另一个的蒙版。在实践中,我想使FITS图像B的像素值保持在FITS图像A的某个值之上(并将图像B的其余px放入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