Metode Otsu

Dalam penglihatan komputer dan pengolahan citra digital, **metode Otsu** () dipakai untuk melakukan pengambangan (*thresholding*) citra otomatis. Metode ini dinamai dari .

Sederhananya, algoritme ini mengembalikan nilai ambang intensitas tunggal yang membagi piksel-piksel menjadi dua kelas, yaitu latar depan dan latar belakang. Nilai ambang ini ditentukan dengan meminimalkan ragam intensitas dalam kelas atau memaksimalkan ragam intensitas antarkelas. Perluasan untuk pengambangan banyak tingkat telah dijelaskan dalam karya tulis asli dan implementasi yang efisien telah diusulkan.

## Metode Otsu

Algoritma ini secara menyeluruh mencari nilai ambang yang meminimalkan ragam dalam kelas yang didefinisikan sebagai jumlah berbobot ragam kedua kelas.

: \sigma^2_w(t) = \omega_0(t)\sigma^2_0(t) + \omega_1(t)\sigma^2_1(t)

Bobot \omega_{0} dan \omega_{1} adalah peluang kedua kelas dipisahkan oleh nilai ambang t. Nilai \sigma^2_{0} dan \sigma^2_{1} adalah ragam kedua kelas.

Peluang kelas \omega_{0, 1}(t) dihitung dari L tempat dalam histogram.

: \begin{align}
\omega_0(t) &= \sum_{i = 0}^{t - 1} p(i) \\ [4pt]
\omega_1(t) &= \sum_{i = t}^{L - 1} p(i)
\end{align}

Untuk dua kelas, meminimalkan ragam dalam kelas sama dengan memaksimalkan ragam antarkelas:

: \begin{align}
\sigma^2_b(t) &= \sigma^2 - \sigma^2_w(t) \\
&= \omega_0(\mu_0 - \mu_T)^2 + \omega_1(\mu_1 - \mu_T)^2 \\
&= \omega_0(t) \omega_1(t) \left[\mu_0(t) - \mu_1(t)\right]^2
\end{align}

dengan \omega adalah peluang kelas dan \mu adalah rata-rata kelas. Hubungan nilai \mu_0(t), \mu_1(t), dan \mu_T adalah sebagai berikut.

: \begin{align}
\mu_0(t) &= \frac{\sum_{i = 0}^{t - 1} i p(i)}{\omega_0(t)} \\ [4pt]
\mu_1(t) &= \frac{\sum_{i = t}^{L - 1} i p(i)}{\omega_1(t)} \\
\mu_T &= \sum_{i = 0}^{L - 1} i p(i)
\end{align}

Hubungan berikut dapat dibuktikan dengan mudah.

: \begin{align}
\omega_0 \mu_0 + \omega_1 \mu_1 &= \mu_T \\
\omega_0 + \omega_1 &= 1
\end{align}

Peluang dan rata-rata kelas dapat dihitung secara iteratif. Cara ini menghasilkan algoritma yang efektif.

### Algoritma
# Hitung histogram dan peluang-peluang pada tiap tingkat intensitas.
# Hitung nilai awal \omega_i(0) dan \mu_i(0).
# Hitung semua kemungkinan nilai ambang t = 1, \ldots, intensitas maksimal.
## Perbarui \omega_i dan \mu_i.
## Hitung \sigma^2_b(t).
# Nilai ambang yang diminta adalah yang menghasilkan \sigma^2_b(t) maksimal.

### Implementasi Octave atau MATLAB
`cacahHistogram` adalah histogram 256 elemen dari citra berderajat keabuan (8 bit). `level` adalah nilai ambang citra (`double`).


```matlab
function ambang = otsu(cacahHistogram)
	total = sum(cacahHistogram); % banyak piksel dalam citra
	tingkat = size(cacahHistogram, 1);
	jumlahB = 0;
	wB = 0;
	maks = 0.0;
	jumlah1 = dot(0:(tingkat - 1), cacahHistogram);
	for (ii = 1:tingkat)
		wF = total - wB;
		if (wB > 0 && wF > 0)
			mF = (jumlah1 - jumlahB) / wF;
			val = wB * wF * ((jumlahB / wB) - mF) * ((jumlahB / wB) - mF);
			if (val >= maks)
				ambang = ii;
				maks = val;
			end;
		end;
		wB = wB + cacahHistogram(ii);
		jumlahB = jumlahB + (ii-1) * cacahHistogram(ii);
	end;
end;
```


Matlab memiliki fungsi `graythresh()` dan `multithresh()` dalam *Image Processing Toolbox* yang menggunakan metode Otsu dan metode multi-Otsu. Octave memiliki fungsi `graythresh()` dalam paket *Image Processing* yang secara bawaan menggunakan metode Otsu. Untuk menghitung dengan metode Otsu (khusus), Octave memiliki fungsi `otsuthresh()`.

## Kelemahan
Metode Otsu cukup baik jika histogram dianggap memiliki persebaran dua puncak serta memiliki lembah yang curam dan dalam di antara dua puncak. Namun, bila luas objek (latar depan) berukuran kecil dibandingkan luas latar belakang, histogram tidak lagi memiliki sifat dua puncak. Jika ragam objek dan ragam latar belakang cukup besar dibanding selisih rata-rata atau jika citra sangat rusak akibat derau aditif, kecuraman lembah pada histogram menurun. Akibatnya, nilai ambang yang dihasilkan menyebabkan kesalahan segmentasi.

Hasil empiris menunjukkan bahwa kinerja pengambangan global yang dipakai untuk segmentasi citra (termasuk metode Otsu) dibatasi oleh ukuran objek yang kecil, selisih rata-rata antara latar depan dan belakang yang kecil, ragam latar depan dan latar belakang masing-masing yang besar, derau yang banyak, dan lain-lain.

## Pengembangan

Berbagai perluasan telah dikembangkan untuk mengatasi batasan-batasan metode Otsu. Salah satunya adalah **metode Otsu dua dimensi** yang bekerja lebih baik untuk citra berderau. Dalam cara ini, intensitas tiap piksel dibandingkan dengan rata-rata intensitas tetangganya untuk memperbaiki hasil segmentasi.

## Lihat pula
* Pengambangan (pengolahan citra)
* Pengambangan histogram seimbang

## Referensi

## Pranala luar
* Implementasi metode Otsu sebagai pengaya GIMP dengan Script-Fu
*  Catatan kuliah tentang pengambangan – membahas metode Otsu
* Pengaya ImageJ dengan metode Otsu
*  Penjelasan lengkap tentang metode Otsu dengan contoh dan implementasi dalam Java
* Implementasi metode Otsu dalam ITK
*  Pengambangan Otsu dalam C# – implementasi dalam C# beserta penjelasan
* Metode Otsu dengan MATLAB
* Pengambangan Otsu dengan scikit-image dalam Python