Transformasi Fourier Untuk Peningkatan Kualitas Citra Erna Zuni Astuti
Abstract : Fourier Transform is an important trasformasi on signal processing, especially in image processing. Fourier transform is tools to convert from spatial domain to frequency domain, which is a function of time domains (spasial) into the frequency domain. In general, Fourier transformation is divided into two, spatial transformation and domain transformation. Spatial transformation is to change intensity of pixel (brightness, contrast, negation, thresholding), while domain transformation is a process of change of image from one domain to another domain. In this way, transformed first digital image with the Fourier transformation, followed by manipulation of the results of the Fourier transformation.Once manipulation is complete, do inverse Fourier transformation to get the image back. Fourier transformation changing the spatial domain into the frequency domain in the lowpass filter and highpass filter. Lowpass filter to work at low frequency, the greater the radius of the filter, the more low frequencies are passed so that the image of the original image mendekatai results, conversely the smaller the radius of the circle on the filter, the less amount of low frequency is passed so that the image becomes blur of cint original. highpass filter while working at high frequency, the greater the radius of the filter, the more high frequencies are passed so that the image looks darker not even seen the original image. Meanwhile, if the radius of the filter is smaller, then the high frequency passed the less, it produces a dark image but it still looks the image sketched the original image. Keywords : Transformasi Fourier, Domain Frekuensi, Highpass Filter, Lowpass Filter
PENDAHULUAN Transformasi Fourier adalah sub pokok bahasan dalam mata kuliah yang sering dipelajari di Fakultas MIPA, maupun Fakultas Teknik, Di Fakultas Ilmu Komputer transformasi ini diterapkan pada materi pengolahan citra. Aplikasi yang memanfaatkan transformasi fourier yaitu pengolahan sinyal digital dan pengolahan citra digital. Sebagai aplikasi transformasi fourier dalam pengolahan citra digital, sebuah citra memegang peranan sangat penting sebagai bentuk informasi visual. Citra mempunyai karakteristik yang kaya dengan informasi. Meskipun sebuah citra kaya informasi, namun seringkali citra mengalami penurunan mutu (degradasi), misalnya
Erna Zuni Astuti adalah Dosen Fakultas Ilmu Komputer UDINUS Semarang 53
54
Techno.Com, Vol. 9 No. 1, Februari 2010
mengandung cacat atau derau (noise), warnanya terlalu kontras, kurang tajam, kabur (bluring) dan sebagainya. Citra yang demikian ini menjadi sulit diinterprestasikan karena informasi yang disampaikan oleh citra tersebut menjadi berkurang. Agar citra mudah diinterprestasikan, maka citra tersebut perlu dimanipulasi menjadi citra lain dengan cara ditingkatkan kualitasnya sesuai kebutuhan masing – masing pengguna. Dalam proses pengolahan citra, transformasi fourier dapat digunakan untuk perbaikan citra (image restoration) atau peningkatan kualitas citra (image enhancement). Perbaikan citra diartikan sebagai proses untuk mengolah citra digital yang didapat agar lebih mendekati bentuk citra aslinya, atau sering disebut sebagai proses mendapatkan kembali citra asli dari suatu citra yang telah mengalami proses degradasi. Sedangkan peningkatan kualitas citra adalah suatu proses untuk mengubah sebuah citra baru sesuai dengan kebutuhan yang diinginkan melalui berbagai cara. Pada transformasi spasial yang diubah adalah intensitas piksel (brightness, kontras, negasi, thresholding) atau posisi piksel (rotasi, translasi, scaling, shear, dan lain- lain). Transformasi jenis ini relatif mudah diimplementasikan dan banyak aplikasi yang dapat melakukannya (Paint, ACDSee, dan lain-lain). Sedangkan transformasi domain yaitu proses perubahan citra dari suatu domain ke domain lainnya, sebagai contoh domain spasial ke domain frekuensi. Transformasi yang mengubah dari domain spasial ke domain frekuensi disebut transformasi fourier.
TINJAUAN PUSTAKA Pengolahan Citra Definisi citra (image) menurut kamus Webster adalah suatu representasi, kemiripan atau imitasi dari suatu objek atau benda. Secara harfiah, citra adalah gambar pada bidang dua dimensi (dwimatra). Ditinjau dari sudut pandang matematis, merupakan fungsi menerus (continue) dari intensitas cahaya pada bidang dua dimensi. Citra yang terlihat merupakan cahaya yang direfleksikan dari objek, sumber cahaya menerangi objek, objek memantulkan kembali sebagian dari berkas cahaya tersebut dan pantulan cahaya ditangkap oleh alat–alat optik, misalnya mata manusia, kamera, scanner, sensor satelit dan sebagainya, sehingga bayangan objek yang disebut citra tersebut terekam. Citra sebagai keluaran dari suatu sistem perekaman data dapat bersifat Digital, artinya dapat langsung disimpan pada media penyimpan magnetic.Citra yang dimaksudkan dalam penulisan ini adalah citra diam (still image). Citra diam adalah citra tunggal yang tidak bergerak. Untuk selanjutnya citra diam akan disebut citra saja. Perbaikan Citra Jenis operasi ini bertujuan untuk memperbaiki citra dengan cara memanipulasi perameter-parameter citra. Perbaikan citra diartikan sebagai proses untuk mengolah citra digital yang didapat agar lebih mendekati bentuk citra aslinya, atau sering disebut sebagai proses mendapatkan kembali citra asli dari suatu citra yang telah mengalami proses degradasi.
Transformasi Fourier (Erna)
55
Peningkatan Kualitas Citra (Image Enhancement) Proses peningkatan kualitas citra atau image enhancement berasal dari enhancement yang artinya mempertinggi atau meningkatkan kualitas citra dengan metode–metode tertentu. Peningkatan kualitas citra adalah suatu proses untuk mengubah sebuah citra menjadi citra baru sesuai dengan kebutuhan yang diinginkan melalui berbagai cara, bisa berupa fungsi transformasi, operasi matematis, pemfilteran, dan lain-lain. Tujuan utama dari peningkatan kualitas citra adalah untuk memproses citra sehingga citra yang dihasilkan lebih baik dari pada citra aslinya Pemugaran Citra (image Restoration) Operasi ini bertujuan menghilangkan atau meminimalkan cacat pada citra. Tujuanya hampir sama dengan operasi perbaikan citra, bedanya pada pemugaran citra penyebab degradasi gambar diketahui.
Tranformasi Fourier Transformasi Fourier 2D Transformasi Fourier 2D dari suatu fungsi f(x,y) didefinisikan sebagai berikut
F (u, v)
f ( x, y). exp[ j 2 .(u.x v. y)]dxdy
5-2(a)
Invers dari transformasi tersebut adalah
f ( x, y )
F (u, v). exp[ j 2 .(u.x v. y)]dudv
5-2(b)
Discrete Fourier Transform (DFT) 1D Transformasi Fourier bersifat kontinu, oleh karena itu sulit untuk dilakukan komputasi karena ada operasi integral dan sifat kontinu itu sendiri. Untuk kebutuhan pengolahan citra, fungsi yang akan ditransformasi haruslah fungsi diskrit. Transformasi Fourier Diskrit atau disebut dengan Discrete Fourier Transform (DFT) adalah Transformasi Fourier yang dikenakan pada fungsi diskrit, dan hasil transformasinya juga merupakan fungsi diskrit. DFT 1D didefinisikan sebagai N 1
F (u ) f ( x). exp[ j x 0
Invers dari transformasi tersebut adalah
2 .u.x ] N
6-3(a)
Techno.Com, Vol. 9 No. 1, Februari 2010
56
f ( x)
1 N 1 2 .u.x F (u ). exp[ j ] N x 0 N
Dimana u = 0,1,2,..........N−1
dan
x = 0,1,2,..........N−1
Contoh 6.3 Diketahui f(x) dalam bentuk diskrit berikut f(x) 4 2
0
1
2
3
4
5
6
7
x
Fungsi tersebut mempunyai N = 8 (dari 0,1,2,.....7). DFT dari f(x) adalah sebagai berikut 7
F (u ) f ( x). exp[ j x 0
2 .u.x ] 8
Untuk u = 0, 7
F (0) f ( x). exp[0] x 0 7
F (0) f ( x) x 0
F (0) f (0) f (1) f (2) f (3) f (4) f (5) f (6) f (7) F (0) 2 4 4 2 2 4 4 2 F (0) 24 Untuk u = 1, 7
F (1) f ( x). exp[ j x 0
2 .x ] 8
6-3(b)
Transformasi Fourier (Erna)
57
2 3 F (1) f (0). exp[0] f (1). exp[ j ] f (2). exp[ j ] f (3). exp[ j ] 4 4 4 f (4). exp[ j
4 5 6 7 ] f (5). exp[ j ] f (6). exp[ j ] f (7). exp[ j ] 4 4 4 4
3 3 F (1) 2.1 4.cos( ) j sin( ) 4.cos( ) j sin( ) 2.cos( ) j sin( ) 4 4 2 2 4 4 5 5 2.cos( ) j sin( ) 4.cos( ) j sin( ) 4 4 3 3 7 7 4.cos( ) j sin( ) 2.cos( ) j sin( ) 2 2 4 4
F (1) 2 2 2 2 2 j 4 j 2 2 2 2 j 2 2 2 2 2 j 4 j 2 2 2 2 j F(1) = 0
Untuk u = 2, 7
F (2) f ( x). exp[ j x 0
2 2.x ] 8
7
.x
x 0
2
F (2) f ( x). exp[ j
]
2 3 F (2) f (0). exp[0] f (1). exp[ j ] f (2). exp[ j ] f (3). exp[ j ] 2 2 2 f (4). exp[ j
4 5 6 7 ] f (5). exp[ j ] f (6). exp[ j ] f (7). exp[ j ] 2 2 2 2
F(2) = − (4+4j)
Untuk u = 3, 7
F (3) f ( x). exp[ j x 0
7
F (3) f ( x). exp[ j x 0
2 3.x ] 8
3 .x ] 4
3 6 9 ] f (2). exp[ j ] f (3). exp[ j ] 4 4 4 12 15 18 21 f (4). exp[ j ] f (5). exp[ j ] f (6). exp[ j ] f (7). exp[ j ] 4 4 4 4
F (3) f (0). exp[0] f (1). exp[ j
Techno.Com, Vol. 9 No. 1, Februari 2010
58
F (3) 0 Untuk u = 4, 7
F (4) f ( x). exp[ jx] x 0
F (4) f (0). exp[0] f (1). exp[ j ] f (2). exp[ j 2 ] f (3). exp[ j3 ] f (4). exp[ j 4 ] f (5). exp[ j5 ] f (6). exp[ j 6 ] f (7). exp[ j 7 ] F (4) 0 Untuk u = 5, 7
F (5) f ( x). exp[ j x 0 7
F (5) f ( x). exp[ j x 0
2 .5.x ] 8 5 .x ] 4
F (5) f (0). exp[0] f (1). exp[ j f (4). exp[ j
5 10 15 ] f (2). exp[ j ] f (3). exp[ j ] 4 4 4
20 25 30 35 ] f (5). exp[ j ] f (6). exp[ j ] f (7). exp[ j ] 4 4 4 4
F (5) 0 Untuk u = 6, 7
F (6) f ( x). exp[ j x 0 7
F (6) f ( x). exp[ j x 0
2 .6.x ] 8 3 .x ] 2
F (6) f (0). exp[0] f (1). exp[ j f (4). exp[ j
3 6 9 ] f (2). exp[ j ] f (3). exp[ j ] 2 2 2
12 15 18 21 ] f (5). exp[ j ] f (6). exp[ j ] f (7). exp[ j ] 2 2 2 2
F (6) (4 4 j )
Transformasi Fourier (Erna)
59
Untuk u = 7, 7
F (7) f ( x). exp[ j x 0 7
F (7) f ( x). exp[ j x 0
2 .7.x ] 8 7 .x ] 4
F (7) f (0). exp[0] f (1). exp[ j f (4). exp[ j
7 14 21 ] f (2). exp[ j ] f (3). exp[ j ] 4 4 4
28 35 42 49 ] f (5). exp[ j ] f (6). exp[ j ] f (7). exp[ j ] 4 4 4 4
F (7) 0
Hasil DFT fungsi f(t) di atas, setelah unsur real dan kompleksnya dipisahkan adalah k
Real {F(k)}
Im {F(k)}
0
24
0
1
0
0
2
−4
− 4j
3
0
0
4
0
0
5
0
0
6
−4
− 4j
7
0
0
Discrete Fourier Transform (DFT) 2D Transformasi Fourier diskrit dua dimensi dari sebuah fungsi diskrit f(x,y) dinyatakan segabai berikut
F (u, v)
x M 1 y N 1
x 0
Inversnya adalah
y 0
ux vy f ( x, y). exp 2j dxdy M N
6-4(a)
Techno.Com, Vol. 9 No. 1, Februari 2010
60
f ( x, y)
1 MN
M 1 N 1
u 0 v 0
ux
vy
F (u, v). exp 2j M N dxdy
Dimana u = 1,2,....M−1; x = 1,2,....M−1; v = 1,2,....N−1
6-4(b)
y = 1,2,....N−1.
Misalkan R(u,v) dan I(u,v) adalah komponen real dan imajiner dari F(u,v), spektrum fourier didefinisikan sebagai
F (u, v) R 2 (u, v) I 2 (u, v)
1/ 2
6-5
Dan sudut fase transformasi definisikan sebagai
I (u, v) (u, v) tan 1 R(u, v)
6-6
Sehingga F(u,v) dapat ditulis sebagai
F (u, v) F (u, v) e i (u ,v )
6-7
Spektrum daya didefinisikan sebagai
P(u, v) F (u, v) R 2 (u, v) I 2 (u, v) 2
6-8
Transformasi Fourier untuk Peningkatan kualitas Citra Karena citra digital merupakan besaran diskrit dua dimensi (2D), maka untuk analisis citra hanya dibutuhkan transformasi fourier diskrit 2D. Untuk menganalisis citra pada domain frekuensi, hasil transformasi fourier dapat ditampilkan sebagai citra, dimana intensitasnya sebanding dengan besarnya |F(u,v)| atau spektrum fourier. Ideal Lowpass Filter Ideal Lowpass Filter adalah filter yang menghilangkan semua komponen frekuensi tinggi dari transformasi fourier pada jarak yang lebih besar daripada jarak spesifik D0 (dari pusat transformasi). Filter semacam ini dinamakan Ideal Lowpass filter (ILPF).
1 if H (u, v) 0 if
D(u, v) D0 D(u, v) D0
D(u,v) adalah jarak dari titik (u,v) ke titik pusat transformasi.
6-10
Transformasi Fourier (Erna)
61
Butterworth Lowpass Filters Fungsi dari Butterworth LowPass Filter (BLPF) dengan order n, dan dengan cutoff frequency pada jarak D0 dari titik asal adalah:
H (u, v)
1 2n 1 D(u, v) / D0
Gaussian Lowpass Filters Fungsi Gaussian LowPass Filters (GLPF) dinyatakan oleh persamaan berikut
H (u, v) e D
2
( u ,v ) / 2 2
Bila σ = D0 , persamaan nya akan menjadi , H (u, v) e
D 2 ( u ,v ) / 2 Do2
dimana D0 adalah cutoff frequency. Bila D(u,v)=D0, filter akan turun hingga 0,607 dari nilai maximumnya. Aplikasi yang menerapkan lowpass filtering, di antaranya adalah: pengenalan karakter, printing and publishing industry, dan processing satellite and aerial images. Lowpass filtering juga dapat dimanfaatkan untuk memperbaiki citra agar lebih indah dipandang, misalnya: pada citra wajah seseorang, cacat atau kelemahan pada wajahnya dapat “ditutupi” dengan cara diblurr (disamarkan). Lowpass filtering membantu memudahkan analisa dengan merata-rata fitur yang penting agar jumlah fitur yang perlu diperhatikan menjadi lebih kecil. Ideal Highpass Filters Ideal Highpass Filter (IHPF) merupakan kebalikan dari ILPF. IHPF memberikan nilai 0 untuk semua frequency di dalam lingkaran radius D0 ketika dilewati, tanpa pengurangan, semua frequency di luar lingkaran diset menjadi 1. Butterworth Highpass Filters Butterworth Highpass Filter (BHPF) mempunyai persamaan sebagai berikut
H (u, v) 1 H BLPF (u, v)
Gaussian Highpass Filters Fungsi Gaussian Highpass Filter (GHPF) dengan cutoff frequency terletak pada jarak D0 dari origin, adalah:
Techno.Com, Vol. 9 No. 1, Februari 2010
62
H (u, v) 1 H GLPF (u, v) Pada proses Inverse, citra asli akan ditransformasi dengan transformasi fourier untuk mengubah domain spasial ke dalam domain frekuensi. Sehingga hasilnya akan dikalikan dengan citra lowpass filter 2D atau highpass filter 2D. Selanjutnya dilakukan proses inverse untuk mengembalikan citra dari domain frekuensi ke dalam domain spasial.
Rancangan proses Inverse ini dapat digambarkan dalam FlowChart berikut: Peningkatan_Kualitascitr a File VISUALISASI TRANSFORMASI FOURIER Domain Spasial
DomainFrekuen si
Buka Gambar
Bersih Layar
Masukkan jari2 lingkaran
Lowpass Filter
Filter 2D
Reset Highpass Filter
ILPF
ILPF
BLP F
BLP F
GLP F
GLP F
Inverse
Keluar
Filter 3D
Gambar 2: Rancangan Output Cita Asli, Lowpass Filter 2D / Highpass Filter 2D, dan Lowpass Filter 3 D / Highpass Filter 3D.
HASIL DAN PEMBAHASAN Visualisasi Transformasi Fourier untuk jari-jari = 50, 40, 30, 20 dan 10
Transformasi Fourier (Erna)
63
Dilakukan ujicoba untuk citra dengan jari-jari 50, 40, 30, 20 dan 10, dengan salah satu Lowpass Filter dan salah satu Highpass Filter. Lowpass Filter terdiri Ideal Lowpass Filter, Butterworth Lowpass Filter dan Gaussian Lowpass Filter. Sedangkan Highpass Filter terdiri dari Ideal Highpass Filter, Butterworth Highpass Filter, Gaussian Highpass Filter. Citra yang diujicobakan adalah citra celsi.jpg. Hasil dari percobaan ini bisa dilihat di dalam ringkasan hasil akhir dari semua percobaan. Citra celsi.jpg dengan Ideal Lowpass Filter Jari-
Citra
ILPF
ILPF
Transformasi
Transformasi
jari
Asli
2D
3D
Fourier 2D
Fourier 3D
Citra Hasil
50
40
30
20
10
Gambar3: Hasil Pengujian Transformasi fourier pada Ideal Lowpass Filter citra celsi.jpg
Techno.Com, Vol. 9 No. 1, Februari 2010
64
Jari-
Citra Asli
jari
IHPF 2 D
IHPF 3 D
Transformas
Transformas
i Fourier 2D
i Fourier 3D
Citra Hasil
50
40
30
20
10
Gambar 4: Hasil Pengujian Transformasi fourier pada Ideal Highpass Filter citra celsi.jpg
KESIMPULAN 1. Transformasi fourier adalah transformasi yang mengubah dari domain spasial ke domain frekuensi. 2. Peningkatan kualitas citra adalah suatu proses untuk mengubah sebuah citra menjadi citra baru sesuai dengan kebutuhan yang diinginkan melalui berbagai cara. 3. Transformasi fourier dapat memproses citra dalam domain frekuensi, hasil transformasi fourier, filter ideal, butterworth filter, gaussian filter baik lowpass filter maupun highpass filter.
Transformasi Fourier (Erna)
65
4. Dengan transformasi fourier pada lowpass filter jika jari-jari pada filter semakin besar maka frekuensi rendah yang diloloskan semakin banyak sehingga menghasilkan citra hasil yang hampir mendekati citra asli. Sebaliknya jika pada lowpass filter jari-jari pada filter semakin kecil maka frekuensi rendah yang diloloskan semakin sedikit sehingga menghasilkan citra hasil yang blur atau menjadi tidak jelas dan berbeda dengan citra asli.
DAFTAR PUSTAKA [1]. [2]. [3]. [4]. [5]. [6]. [7].
Murni, Aniati. 1992. Pengantar Pengolahan Citra. Jakarta: PT Elex Media Komputindo Firmansyah A. 2003. Dasar - dasar Pemrograman Matlab. Indonesia: IlmuKomputer. Paulus, Erick S.Si, M.Kom, dan Nataliani, Yessica, S.Si., M.Kom. 2007. Cepat Mahir GUI Matlab.Yogyakarta : Andi Munir, Rinaldi. 2004. Pengolahan Citra Digital dengan Pendekatan Algoritmik. Bandung: Informatika Ch.Wijaya, Marvin dan Prijono, Agus. 2007. Pengolahan Citra Digital Bandung : http://www.pdf-search-engine.com/transformasi-fourier-pdf.html, diakses tanggal 10 Januari 2009 http://elearning.gunadarma.ac.id/docmodul/matematika_lanjut/bab10-deret_fourier.pdf, diakses tanggal 12 Februari 2009