实验二，图像变换

gaohaikuo1

内容一：对一维信号[0，1，0，1]作傅里叶变换（套用公式）
学习matlab矩阵运算的方法
内容二：（1）建立一个400*400的全黑图像（RGB）;（2）图像的最左上角的20*20调整为红色。
I = zeros(400,400,3);
I = uint8(I);
subplot(1,2,1);imshow(I);
I(1:20,1:20,1) = 255;
subplot(1,2,2);imshow(I);
内容三：（1）建立一个400*400的全黑图像（RGB）;（2）把中间宽为81的竖条改为红色（3）图像灰度化
内容四：（1）建立一个400*400的全黑图像（灰度图）;（2）把中间宽为81的竖条改为白色（3）对图像快速傅里叶变换。
内容五：（1）建立一个400*400的全黑图像（灰度图）;（2）把中间宽为81高81的方块改为白色（3）对图像快速傅里叶变换。
内容六：对上述频率图做反傅里叶变换；
内容七：把原始图片换成lena.jpg的灰度图，重复内容五和六。
内容八：在lena.jpg的频谱图上添加白点，观察反变换之后的图像

gaohaikuo1

I = imread('lena.jpg');
I = rgb2gray(I);subplot(2,3,1);imshow(I);title('original image');
F = fft2(I);subplot(2,3,2);imshow(log(1+abs(F)),[]);title('frequency image');
F_c  = fftshift(F);
subplot(2,3,3);imshow(log(1+abs(F_c)),[]);title('centered frequency image');
F_c(255-30:255-20,254-30:254-20) = 1000000;
subplot(2,3,4);imshow(log(1+abs(F_c)),[]);title('centered frequency image');
F1 = ifftshift(F_c);
subplot(2,3,5);imshow(log(1+abs(F1)),[]);title('new frequency image');
I1 = uint8(ifft2(F1));
subplot(2,3,6);imshow(I1);title('new image');