1 / 31

Pengolahan dalam Domain Frekuensi dan Restorasi Citra

Pengolahan dalam Domain Frekuensi dan Restorasi Citra. Pengolahan Citra Digital Materi 5. Eko Prasetyo Teknik Informatika Universitas Muhamamdiyah Gresik 2011. Domain Spasial vs Domain Frekuensi. Domain Spasial Konsep koordinat baris dan kolom Pemrosesan pixel-by-pixel

vea
Télécharger la présentation

Pengolahan dalam Domain Frekuensi dan Restorasi Citra

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Pengolahandalam Domain FrekuensidanRestorasi Citra Pengolahan Citra Digital Materi 5 Eko Prasetyo TeknikInformatika UniversitasMuhamamdiyah Gresik 2011

  2. Domain Spasialvs Domain Frekuensi Domain Spasial • Konsepkoordinatbarisdankolom • Pemrosesan pixel-by-pixel • Komputasi lama (terutamacitradenganukuranspasialtinggi) Domain Frekuensi • Konsepfrekuensi, perubahanintensitaspikselkepiksel (frekuensirendahdantinggi) • Pemrosesanberdasarkanpemilihanfrekuensi yang akandifilteratautidak • Komputasirelatifcepat (terutamacitradenganukuranspasialtinggi)

  3. KonsepFrekuensidalamcitra • Sembarangsinyalspasialmempunyairepresentasifrekuensi • Maknafrekuensidalamcitra: • Komponenfrekuensitinggidikaitkandenganperubahanpikselkepikselsecaracepatsepanjangcitra. Misal: teks, tekstur, dsb. • Komponenfrekuensitinggidikaitkandenganfiturberskalabesarpadacitra. Misal: daerahdenganintensitaskonstan, ataupiksel yang jumlahnyamendominasidalamseluruhdaerahcitra.

  4. Transformasi Fourier • Fungsiperiodikdapatdinyatakansebagaijumlah sinus dan/ataucosinusdarperbedaanfrekuensisetiapperkaliannyadengankoefisien yang berbeda

  5. Transformasi Fourier • Fungsi yang tidakperiodiktetapidengandaerahkurva yang terbatasdapatdinyatakansebagai integral sinus dan/ataucosinusdikalikandenganfungsibobot. • Transformasi Fourier 1 dimensi: • Transformasi Fourier 2 dimensi:

  6. Transformasi Fourier Diskrit • Karenacitraadalahgelombangdiskrit, makafungsi f(x), x=0,1,…,M-1, untuksatudimensikitamendapatkan: • Formula Euler: • Sehinggadidapatkan: • Untuku = 0,…,M-1 • f(x) adalahnilaiintensitassetiappiksel • Nilai u adalahkomponendalam domain frekuensi • Setiap F(u) adalahnilaifrekuensidalamtransformasi

  7. Transformasi Fourier Diskrit 2-D • Untukcitra 2 dimensi, DFT yang digunakan: • Untuk u=0,…,M-1 and v=0,…,N-1 daniDFTdidefinisikan: • Karenanilai FT adalahbilangankompleks, kadang-kadangkitanyatakan F(u) dalamkoordinat polar: • Dimanajarakatauspektrumdinyatakandengan: • Sudutfasednyatakanoleh:

  8. Transformasi Fourier Diskrit 2-D • Untuk u=0, v=0, didapatkan: • Samadengan rata-rata nilaiintensitas. • Lokasiinijugaadalahtitik origin pada domain frekuensi.

  9. MendapatkanSpektrum Fourier Citra >> f = imread('rice.tif'); >> f = im2double(f); >> F = fft2(f); >> figure, imshow(F); >> F2 = log(1+abs(F)); >> figure, imshow(F2,[ ]); >> Fs = fftshift(F2); >> figure, imshow(Fs,[ ]); >> f2 = ifft2(F); Citra asli Spektrumasli Frekuensitinggi Frekuensirendah Origin Spektrumsetelahdi-enhance dengan log Setelahdigeser (memusatkan origin)

  10. Filter dalam Domain Frekuensi Dasaruntuk filter linear dalam domain spasialdanfrekuensiadalahteorikonvolusi, yang dapatdituliskandengan: • Pemfilterandalam domain spasialberisikonvolusicitra f(x,y) mask filter h(x,y). • Sepertihalnyateorikonvolusi, jugabisamendapatkanhasil yang samadalam domain frekuensidenganperkalianantara F(u,v) dengan H(u,v), transformasi Fourier filter spasial. • Dasarnya, idedalampemfilteran domain frekuensiadalahuntukmemilihfungsi transfer filter yang memodifikasi F(u,v) dengancaratertentu.

  11. Langkahpemfilteran • Tentukan parameter padding menggunakanfungsipaddedsize: PQ = paddedsize(size(f)); • Tentukantransformasi Fourier dengan padding: F = fft2(f, PQ(1), PQ(2)); • Keluarkanfungsi filter H berukuran PQ(1) x PQ(2) menggunakanmetode yang akandibahaspadababini. Filter harusdalam format tidakterpusat. Jikamasihterpusatmakaharusdigeserdenganfungsifftshiftsebelummenggunakan filter. • Kalikantransformasidengan filter: G = H.*F; • Tentukanbagian real dariinvers FFT dari G: g = real(ifft2(G)); • Potongpersegipanjangpadabagiankiriataspadaukuranasli: g = g(1:size(f, 1), 1:size(f, 2));

  12. Teknik Filter dalam Domain Frekuensi Filter Penghalusan (Smoothing) • Ideal Lowpass Filter (ILPF) • Butterworth Lowpass Filter (BLPF) • Gaussian Lowpass Filter (GLPF) Filter Penajaman (Sharpening) • Ideal Highpass Filter (IHPF) • Butterworth Highpass Filter (BHPF) • Gaussian Highpass Filter (GHPF)

  13. Filter Penghalusan (Smoothing) • Smoothing (blurring) dicapaidalam domain frekuensidenganpelemahanfrekuensitinggi; yang disebutdenganlowpass filter. • Ideal Lowpass Filter (ILPF) • Filter lowpass 2-D yang melewatkantanpapelemahansemuafrekuensidalamlingkaran radius D0dari origin danmeng-“cut off” semuafrekuensidiluarlingkarandisebutIdeal Lowpass Filter (ILPF) yang ditentukanolehfungsi: • dimana D0adalahkonstantapositifdan D(u,v) adalahjarakantaratitik (u,v) dalam domain frekuensidanpusatpersegipanjangfrekuensi, maka: • D(u,v) = [(u – P/2)2 + (v – Q/2)2]1/2

  14. Filter Penghalusan (Smoothing) Butterworth Lowpass Filter • Fungsi transfer filter lowpass Butterworth (BLPF) dari order n, dandengan cutoff frekuensipadajarak D0dari origin, didefinisikansebagai: • dimana D(u,v) dinyatakanolehpersamaansebelumnya.

  15. Filter Penghalusan (Smoothing) • Gaussian Lowpass Filter • Bentuk Gaussian Lowpass Filters (GLPF) dalamduadimensididefinisikandengan: • D(u,v) adalahjarakdaripusatpersegipanjangfrekuensi

  16. Ideal Lowpass Filter x = Spektrumasli ILPF, D0 = 20 Spektrumhasil Citra hasil x = ILPF, D0 = 60

  17. Butterworth Lowpass Filter x = Spektrumasli BLPF, D0 = 20, sig = 2 Spektrumhasil Citra hasil x = BLPF, D0 = 60, sig = 2

  18. Gaussian Lowpass Filter x = Spektrumasli GLPF, D0 = 20 Spektrumhasil Citra hasil x = GLPF, D0 = 60

  19. %SOURCE CODE ILPF f = imread(‘rice.png'); f = im2double(f); [M, N] = size(f); F = fft2(f); D0 = 20; %Ideal Li = lpfilter('ideal', M, N, D0); fli = dftfilt(f,Li); figure, imshow(fftshift(log(1+abs(fft2(f)))),[ ]); F2 = fftshift(log(1+abs(fft2(fli)))); figure, imshow(F2,[ ]); figure, imshow(fftshift(Li)); figure, imshow(fli); %SOURCE CODE BLPF f = imread(‘rice.png'); f = im2double(f); [M, N] = size(f); F = fft2(f); D0 = 20; sig = 2; %Ideal Lb = lpfilter(‘btw', M, N, D0, sig); flb = dftfilt(f,Lb); figure, imshow(fftshift(log(1+abs(fft2(f)))),[ ]); F2 = fftshift(log(1+abs(fft2(flb)))); figure, imshow(F2,[ ]); figure, imshow(fftshift(Lb)); figure, imshow(flb);

  20. %SOURCE CODE GLPF f = imread(‘rice.png'); f = im2double(f); [M, N] = size(f); F = fft2(f); D0 = 20; %Ideal Lg = lpfilter(‘gaussian', M, N, D0); flg = dftfilt(f,Lg); figure, imshow(fftshift(log(1+abs(fft2(f)))),[ ]); F2 = fftshift(log(1+abs(fft2(flg)))); figure, imshow(F2,[ ]); figure, imshow(fftshift(Lg)); figure, imshow(flg);

  21. Mau difilter HIGHPASS ?

  22. Ideal Highpass Filter x = Spektrumasli IHPF, D0 = 20 Spektrumhasil Citra hasil x = IHPF, D0 = 80

  23. Butterworth Highpass Filter x = Spektrumasli BHPF, D0 = 20, sig = 5 Spektrumhasil Citra hasil x = BHPF, D0 = 80, sig = 5

  24. Gaussian Highpass Filter x = Spektrumasli GHPF, D0 = 20 Spektrumhasil Citra hasil x = GHPF, D0 = 80

  25. %SOURCE CODE IHPF f = imread('m.png'); f = im2double(f); [M, N] = size(f); F = fft2(f); D0 = 20; %Ideal Li = hpfilter('ideal', M, N, D0); fli = dftfilt(f,Li); figure, imshow(fftshift(log(1+abs(fft2(f)))),[ ]); F2 = fftshift(log(1+abs(fft2(fli)))); figure, imshow(F2,[ ]); figure, imshow(fftshift(Li)); figure, imshow(fli); %SOURCE CODE BHPF f = imread('m.png'); f = im2double(f); [M, N] = size(f); F = fft2(f); D0 = 20; sig = 2; %Ideal Lb = hpfilter(‘btw', M, N, D0, sig); flb = dftfilt(f,Lb); figure, imshow(fftshift(log(1+abs(fft2(f)))),[ ]); F2 = fftshift(log(1+abs(fft2(flb)))); figure, imshow(F2,[ ]); figure, imshow(fftshift(Lb)); figure, imshow(flb);

  26. %SOURCE CODE GHPF f = imread('m.png'); f = im2double(f); [M, N] = size(f); F = fft2(f); D0 = 20; %Ideal Lg = hpfilter(‘gaussian', M, N, D0); flg = dftfilt(f,Lg); figure, imshow(fftshift(log(1+abs(fft2(f)))),[ ]); F2 = fftshift(log(1+abs(fft2(flg)))); figure, imshow(F2,[ ]); figure, imshow(fftshift(Lg)); figure, imshow(flg);

  27. Selective Filtering • Filter yang dibahassebelumnyaberoperasipadasemuabagianpersegipanjangfrekuensi. • Adaaplikasi yang memprosespada band frekuensitertentuatau region kecildaripersegipanjangfrekuensi. • Filter dalamkategoriinidisebutdenganbandrejectataubandpass filterdannotch filter. • Jenis filter inimudahdibuatmenggunakankonsep filter sebelumnya. • Formula untuk Ideal Bandreject Filter: • Formula untuk Butterworth Bandreject Filter: • Formula untuk Gaussian Bandreject Filter:

  28. Bandreject Perspective plot bandreject Ideal Butterworth Gaussian Tampilancitra

  29. Penguranganefekkorandengan filter Bandreject f = imread('car.tif'); F = fft2(f); %i adalah citra abu-abu newspaper figure, imshow(fftshift(log(1+abs(F))), [ ]); H=bandreject('ideal',size(F, 1), size(F,2), 50, 5); figure, imshow(fftshift(H), [ ]); g = H.*F; figure, imshow(fftshift(log(1+abs(g))), [ ]); G = real(ifft2(g)); G = (G-min(min(G)))./(max(max(G))-min(min(G))); %menormalisasi menjadi 0-1 figure, imshow(G); Citra asli Setelahdifilter

  30. Penguranganefekkorandengan filter Bandreject x = Spektrumasli Filter ideal bandreject, D0 = 50, W = 5 Spektrumhasil

  31. To Be Continued … Materi 6 – Morfologi Citra ANY QUESTION ?

More Related