本文介绍如何使用R语言中的EBImage库进行图像处理，包括图像显示、灰度转换及快速傅里叶变换(FFT)等操作。特别地，文中详细展示了FFT的Shift变换过程，以及如何获取FFT变换后的幅度和相位图像。
setwd(".../Project/R/Workspace/Task1")
library("EBImage" , lib.loc="~/R/win-library/3.2")
# Image
img
display(img, title='Image')
# Grayscaled
img_gray
# FFT
img_ff
###################################################
###################################################
# FFT SHIFT
fftshift
rows
cols
swap_up_down
rows_half
return(rbind(img_ff[((rows_half+1):rows), (1:cols)], img_ff[(1:rows_half), (1:cols)]))
}
swap_left_right
cols_half
return(cbind(img_ff[1:rows, ((cols_half+1):cols)], img_ff[1:rows, 1:cols_half]))
}
if (dim == -1) {
img_ff
return(swap_left_right(img_ff))
}
else if (dim == 1) {
return(swap_up_down(img_ff))
}
else if (dim == 2) {
return(swap_left_right(img_ff))
}
else {
stop("Invalid dimension parameter")
}
}
ifftshift
rows
cols
swap_up_down
rows_half
return(rbind(img_ff[((rows_half+1):rows), (1:cols)], img_ff[(1:rows_half), (1:cols)]))
}
swap_left_right
cols_half
return(cbind(img_ff[1:rows, ((cols_half+1):cols)], img_ff[1:rows, 1:cols_half]))
}
if (dim == -1) {
img_ff
return(swap_up_down(img_ff))
}
else if (dim == 1) {
return(swap_up_down(img_ff))
}
else if (dim == 2) {
return(swap_left_right(img_ff))
}
else {
stop("Invalid dimension parameter")
}
}
###################################################
###################################################
# FFT SHIFT
# Magnitude and Phase
magntd
phase
img_fftsh
display(log(img_fftsh),title="FFT")
扫一扫
5886
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
