Kode Contrast Stretching di Matlab

Buatlah filter Untuk memperbaiki citra bean pada gambar (a) dengan menggunakan fungsi T sebagaimana terlihat pada gambar (b). Dimana m1, m2 dan m3 adalah gradient. Temukan nilai m1,m2 dan m3 yang paling optimal sehingga menghasilkan gambar (c).

Jawaban :

Code Matlab:

%% PROGRAM CONTRAST STRETCHING %%

%% Hapus layar dan variabel %%

clear all;

close all;

clc;

%% Membaca gambar %%

a=imread('C:\dip\week3\figure.tif');

a=double(a);

s=size(a);

%% Mendefinisikan titik-titik %%

p1=[0,0];

p2=[97,30];

p3=[127,225];

p4=[255,255];

%% Mencari gradien garis tiap titik %%

m1=(p1(1,2)-p2(1,2))/(p1(1,1)-p2(1,1));

m2=(p2(1,2)-p3(1,2))/(p2(1,1)-p3(1,1));

m3=(p3(1,2)-p4(1,2))/(p3(1,1)-p4(1,1));

%% Mencari nilai constanta %%

c1=p1(1,2)-m1*p1(1,1);

c2=p2(1,2)-m2*p2(1,1);

c3=p3(1,2)-m3*p3(1,1);

%% Fungsi transformasi %%

t=[];

for x=0:255

 if(x<=p2(1,1))

 t=[t (m1*x+c1)];

 end

 if(x>p2(1,1) && x<=p3(1,1))

 t=[t (m2*x+c2)];

 end

 if(x>p3(1,1) && x<=p4(1,1))

 t=[t (m3*x+c3)];

 end

end

%% Mendapatkan output gambar %%

for n=1:s(1,1)

 for m=1:s(1,2)

 ot(n,m)=t(a(n,m)+1);

 end

end

%% Menampilkan grafik intensitas %%

plot(t)

grid on;

xlabel('Intensity in input image')

ylabel('Intensity in output image')

title('Piece-wise linear transformation : Contrast stretching function')

%% Menampilkan gambar perbandingan %%

figure()

subplot(1,2,1)

imshow(a/255)

title('Original Image')

subplot(1,2,2)

imshow(ot./255)

title('Contrast stretching')

Outputnya:

Nilai yang optimal adalah

• m1 = 0.3093
• m2 = 6.5
• m3 = 0.2344