디지털신호처리및설계(DSP)_MATLAB 프로젝트 - Final Report of Design Project

소개글

디지털신호처리및설계(DSP)_MATLAB 프로젝트 - Final Report of Design Project에 대한 보고서 자료입니다.

목차

Final Report of Design Project

1. Team name and participants
2. Objective
3. Contents of design project
4. Results and discussion
5. Conclusions & Postscripts

본문내용

bilinear : -6.5000e-06 + 1.5831e-05i ~ 0.0294 - 0.0001i
bicubic : -7.5897e-06 - 1.3852e-05i ~ 0.0255 - 0.0023i
6) 주어진 불완전한 레나 영상(2장)을 각자의 방식으로 복원하여 가장 원영상과 흡사하게 (즉, PSNR 수치가 높게) 복원하시오.
gaussian을 이용한 필터링. noise가 많이 제거되어 더욱 선명하게 복원 되었다.
Medan Filter 도 제법 많이 복원 되었지만, Low Pass Filter 가 더 깔끔하게 noise가 필터링 되었다.
5. Conclusions & Postscripts
지금까지 실습수업 때 배운 지식을 가지고 이번 프로젝트에 응용해 봄으로써 MATLAB 프로그램을 이용한 영상처리 기법을 배울 수 있었다. 우리 조는 프로그래밍 언어를 이번 기회에 처음 접했기 때문에 많은 애로사항이 있었다. 도서관에서 관련 서적을 참고하고, 인터넷에서 관련 포스팅 글과 논문을 찾아가며 노력 해 봤지만 결국, 4번과 5번을 완성하지 못해 아쉬움이 남는다. 지금에야 와서 생각 난건데, 모르는 부분을 조교님께 찾아가서 물어봤더라면.. 아쉬움이 많이 남는다. 그래도 결과적으로, MATLAB 프로그램에 더욱 친숙해 질수 있는 기회였고, 영상처리 기술에 더 많은 관심을 갖게 된 계기가 되었다.
<코드>
1)
clear all;
close all;
clc;
image = imread('pic1.bmp');
R=double(image(:,:,1));
G=double(image(:,:,2));
B=double(image(:,:,3));
figure(1);
imshow(image);
gray_image=(rgb2gray(image));
figure(2);
imshow(gray_image);
2)
clear all;
close all;
clc;
image= imread('pic2.jpg');
R =double(image(:,:,1));
G =double(image(:,:,2));
B =double(image(:,:,3));
figure(1);
imshow(image);
figure(2);
x = fspecial('gaussian');
y = imfilter(uint8(R),x)
subplot(131); imshow(y);
title('R')
y1 = imfilter(uint8(G),x)
subplot(132); imshow(y1);
title('G')
z = imfilter(uint8(B),x)
subplot(133); imshow(z);
title('B')
3)
clear all;
close all;
clc;
image = imread('pic1.bmp');
gray=(rgb2gray(image));
f=fspecial('laplacian');
wf=filter2(f,gray);
f1=fspecial('sobel');
wg=filter2(f1,gray);
figure(1)
imshow(image);
figure(2)
imshow(wf);
figure(3)
imshow(wg);
4)
close all;
clear all;
clc;
image = (fopen('Lenna_sample1.raw'));
data = fread(image,[256,256]);
y = data';
rms=(y/(256*256))^(1/2);
y1 = [];
for idx_x1 = 1:128
for idx_y2 = 1:128
y1(idx_x1,idx_y2) = y(idx_x1*2, idx_y2*2);
end
end
rms1=(y1/(idx_x1*idx_y2))^(1/2);
z = imresize(y1,[256,256],'bilinear');
z1 = imresize(y1,[256,256],'bicubic');
rms2=(z/(256*256))^(1/2);
rms3=(z1/(256*256))^(1/2);
show_img(y);
show_img(y1);
show_img(z);
show_img(z1);
6)
clear all;
close all;
clc;
image = (fopen('Lenna_sample2.raw'));
data = fread(image,[256,256]);
BW = data'
BW = BW/255;
[m, n] = size(BW);
figure
subplot(1,3,1);
imshow(BW);
title('Original');
x = zeros(m ,n);
for r = 2:m-1,
for c = 2:n-1,
K = BW((r-1):(r+1), (c-1:c+1));
newv = median(K(:));
x(r, c) = newv;
end
end
subplot(1,3,2);
imshow(x);
title('Median Filter');
k = (1/9)*[1 1 1; 1 1 1; 1 1 1];
Y = conv2(BW, k);
subplot(1,3,3);
imshow(Y);
title('Low Pass Filter');
6_1)
clear all;
close all;
clc;
image = (fopen('Lenna_sample1.raw'));
data = fread(image,[256,256]);
x = data'
[R,C] = size(x);
new_high_boost_one(R,C) = zeros;
x = double(x);
for i=1:R-2
for j=1:C-2
new_high_boost_one(i+1,j+1) = x(i,j)*(-1.5/9) +...
x(i,j+1)*(-1.5/9) + x(i,j+1)*(-1.5/9) + ...
x(i+1,j)*(-1.5/9) + x(i+1,j+1)*(13.5/9) +...
x(i+1,j+1)*(-1.5/9) + ...
x(i+2,j)*(-1.5/9) + x(i+2,j+1)*(-1.5/9) +...
x(i+2,j+1)*(-1.5/9);
end;
end;
figure(1);
subplot(121); imshow(mat2gray(x));
subplot(122); imshow(mat2gray(new_high_boost_one));

키워드

디지털신호,   신호및시스템,   DSP,   MATLAB,   영상추출,   rgb,   gausian,   PSNR,   interpolation,   matlab 프로젝트