画像分割 問題とは、何らかの均質性基準に基づいて画像を複数の領域に分割する問題である。本稿では主に、最小カット または最大カット による グラフ分割 を適用した、画像分割に対するグラフ理論的アプローチについて述べる。分割に基づくオブジェクト分類は、画像分割に適用された スペクトルクラスタリング の特殊なケースと見なすことができる。
画像セグメンテーションの応用 画像圧縮 画像を均質な成分に分割し、各成分に最適な圧縮アルゴリズムを使用して圧縮率を向上させる。 医学的診断 MRI画像から癌領域を自動的にセグメンテーションする。 マッピングと測定 交通機関 輸送ネットワークを分割することで、均質な交通状態を特徴とする領域を特定することが可能になります。[ 1 ]
正規化カットを用いたセグメンテーション
正規化されたカット G = ( V , E , w ) を重み付きグラフとする。A {\displaystyle A} そしてB {\displaystyle B} 2つの頂点部分集合とする。
させて:
w ( A 、 B ) = ∑ 私 ∈ A 、 j ∈ B w 私 j {\displaystyle w(A,B)=\sum \limits _{i\in A,j\in B}w_{ij}} カット ( A 、 B ) = w ( A 、 B ) w ( A 、 V ) + w ( A 、 B ) w ( B 、 V ) {\displaystyle \operatorname {ncut} (A,B)={\frac {w(A,B)}{w(A,V)}}+{\frac {w(A,B)}{w(B,V)}}} ナソック ( A 、 B ) = w ( A 、 A ) w ( A 、 V ) + w ( B 、 B ) w ( B 、 V ) {\displaystyle \operatorname {nassoc} (A,B)={\frac {w(A,A)}{w(A,V)}}+{\frac {w(B,B)}{w(B,V)}}} 正規化カットアプローチでは、[ 2 ] 任意のカットに対して( S 、 S ¯ ) {\displaystyle (S,{\overline {S}})} でG {\displaystyle G} 、カット ( S 、 S ¯ ) {\displaystyle \operatorname {ncut} (S,{\overline {S}})} 異なる部分間の類似性を測定し、ナソック ( S 、 S ¯ ) {\displaystyle \operatorname {nassoc} (S,{\overline {S}})} 同じ部分内の頂点の全体的な類似性を測定します。
以来カット ( S 、 S ¯ ) = 2 − ナソック ( S 、 S ¯ ) {\displaystyle \operatorname {ncut} (S,{\overline {S}})=2-\operatorname {nassoc} (S,{\overline {S}})} カット( S * 、 S ¯ * ) {\displaystyle (S^{*},{\overline {S}}^{*})} 最小限に抑えるカット ( S 、 S ¯ ) {\displaystyle \operatorname {ncut} (S,{\overline {S}})} また最大限に活用するナソック ( S 、 S ¯ ) {\displaystyle \operatorname {nassoc} (S,{\overline {S}})} 。
カットを計算する( S * 、 S ¯ * ) {\displaystyle (S^{*},{\overline {S}}^{*})} 最小限に抑えるカット ( S 、 S ¯ ) {\displaystyle \operatorname {ncut} (S,{\overline {S}})} これはNP困難 問題である。しかし、多項式時間でカットを見つけることができる。( S 、 S ¯ ) {\displaystyle (S,{\overline {S}})} 正規化された小さな重量カット ( S 、 S ¯ ) {\displaystyle \operatorname {ncut} (S,{\overline {S}})} スペクトル技術 を用いて。
計算複雑性 すべての固有ベクトルに対する標準的な固有値問題を解くには(例えばQRアルゴリズムを使用して)、 O ( n 3 ) {\displaystyle O(n^{3})} 時間。これは、次のような画像セグメンテーション アプリケーションには非現実的です。n {\displaystyle n} は画像内のピクセル数です。
非カットアルゴリズムでは、2番目に小さい一般化固有値に対応する固有ベクトルのみが使用されるため、対応する固有値問題の解を行列フリー方式、つまり、 ランチョスアルゴリズム のように行列Wを明示的に操作したり計算したりすることなく実行すれば、効率を劇的に向上させることができます。行列フリー方式では 、各反復で与えられたベクトルに対して行列とベクトルの積を実行する関数のみが必要です。画像セグメンテーションの場合、行列Wは通常疎行列であり、非ゼロ要素が多数存在します。O ( n ) {\displaystyle O(n)} なので、このような行列ベクトル積はO ( n ) {\displaystyle O(n)} 時間。
高解像度画像の場合、第2固有値はしばしば条件が悪く、 ランチョス法 などの反復固有値ソルバーの収束が遅くなります。前処理は 、例えば行列フリーLOBPCG 法において、収束を加速する重要な技術です。最適に前処理された行列フリー法を使用して固有ベクトルを計算するには、O ( n ) {\displaystyle O(n)} 最適な複雑さである時間、固有ベクトルはn {\displaystyle n} コンポーネント。
OBJカット OBJ CUT [ 7 ] は、オブジェクトを自動的にセグメント化する効率的な方法です。OBJ CUT メソッドは汎用的な方法であるため、任意のオブジェクトカテゴリモデルに適用できます。既知のオブジェクトカテゴリのインスタンス (たとえば牛) を含む画像 D が与えられると、OBJ CUT アルゴリズムはオブジェクトのセグメント化を計算します。つまり、ラベルのセットm を推論します。
m をバイナリラベルの集合とし、Θ {\displaystyle \Theta } 形状パラメータである(Θ {\displaystyle \Theta } は、階層化された画像構造 (LPS)モデルのラベルに対する形状事前分布です。エネルギー関数E ( m 、 Θ ) {\displaystyle E(m,\Theta )} は以下のように定義される。
E ( m 、 Θ ) = ∑ ϕ x ( D | m x ) + ϕ x ( m x | Θ ) + ∑ Ψ x y ( m x 、 m y ) + ϕ ( D | m x 、 m y ) {\displaystyle E(m,\Theta )=\sum \phi _{x}(D|m_{x})+\phi _{x}(m_{x}|\Theta )+\sum \Psi _{xy}(m_{x},m_{y})+\phi (D|m_{x},m_{y})} (1)用語ϕ x ( D | m x ) + ϕ x ( m x | Θ ) {\displaystyle \phi _{x}(D|m_{x})+\phi _{x}(m_{x}|\Theta )} は単項項と呼ばれ、Ψ x y ( m x 、 m y ) + ϕ ( D | m x 、 m y ) {\displaystyle \Psi _{xy}(m_{x},m_{y})+\phi (D|m_{x},m_{y})} はペアワイズ項と呼ばれます。単項項は尤度で構成されます。ϕ x ( D | m x ) {\displaystyle \phi _{x}(D|m_{x})} 色と単項ポテンシャルに基づくϕ x ( m x | Θ ) {\displaystyle \phi _{x}(m_{x}|\Theta )} からの距離に基づいてΘ {\displaystyle \Theta } ペアワイズ項は事前情報から構成される。Ψ x y ( m x 、 m y ) {\displaystyle \Psi _{xy}(m_{x},m_{y})} そして対照的な用語ϕ ( D | m x 、 m y ) {\displaystyle \phi (D|m_{x},m_{y})} 。
最高のラベルm * {\displaystyle m^{*}} 最小限に抑える∑ 私 w 私 E ( m 、 Θ 私 ) {\displaystyle \sum \limits _{i}w_{i}E(m,\Theta _{i})} 、 どこw 私 {\displaystyle w_{i}} はパラメータの重みですΘ 私 {\displaystyle \Theta _{i}} 。
m * = 引数 ミニ m ∑ 私 w 私 E ( m 、 Θ 私 ) {\displaystyle m^{*}=\arg \min \limits _{m}\sum \limits _{i}w_{i}E(m,\Theta _{i})} (2)
参考文献 ↑ Lopez, Clélia; Leclercq, Ludovic; Krishnakumari, Panchamy; Chiabaut, Nicolas; Van Lint, Hans (2017年10月25日). 「3D速度マップによる都市交通渋滞パターンの日々の規則性の解明」 . Scientific Reports . 7 (14029): 14029. Bibcode : 2017NatSR...714029L . doi : 10.1038/ s41598-017-14237-8 . PMC 5656590. PMID 29070859 . ↑ Jianbo Shi およびJitendra Malik (1997): "Normalized Cuts and Image Segmentation", IEEE Conference on Computer Vision and Pattern Recognition, pp 731–737 ↑ 「スペクトルクラスタリング — scikit-learn ドキュメント」 。 ↑ Knyazev, Andrew V. (2003). Boley; Dhillon; Ghosh; Kogan (編). スペクトル画像分割とグラフ二分法のための最新の前処理付き固有値ソルバー . 大規模データセットのクラスタリング; 第 3 回 IEEE 国際データマイニング会議 (ICDM 2003) フロリダ州メルボルン: IEEE コンピュータ協会. pp. 59–62 . ↑ Knyazev, Andrew V. (2006). マルチスケールスペクトル画像セグメンテーション 画像セグメンテーションにおけるグラフラプラシアンの固有値を計算するマルチスケール前処理 。Fast Manifold Learning Workshop、WM Williamburg、VA。doi : 10.13140/RG.2.2.35280.02565 。 ↑ Knyazev, Andrew V. (2006). マルチスケールスペクトルグラフ分割と画像セグメンテーション 。スタンフォード大学とYahoo!リサーチ主催の「現代の大規模データセットのためのアルゴリズム」ワークショップ。 ↑ MP Kumar、PHS Torr、A. Zisserman。「Obj cut」。IEEE Conference on Computer Vision and Pattern Recognition 論文集 、サンディエゴ、18~25ページ、2005年。 ↑ E. Borenstein、S. Ullman:クラス固有のトップダウンセグメンテーション。第7回ヨーロッパコンピュータビジョン会議議事録、デンマーク、コペンハーゲン、109~124ページ、2002年。 ↑ Z. Tu、X. Chen、AL Yuille、SC Zhu:画像解析:セグメンテーション、検出、認識の統合。カテゴリレベルの物体認識に向けて 2006: 545–576 ↑ B. Leibe、A. Leonardis、B. Schiele:複合オブジェクト分類とセグメンテーションのための暗黙的形状モデル。カテゴリレベルオブジェクト認識に向けて 2006: 508–524 ↑ J. Winn、N. Joijic。Locus : 教師なしセグメンテーションによるオブジェクトクラスの学習。IEEE 国際コンピュータビジョン会議議事録、北京、2005 年。 ↑ JM Winn、J. Shotton:部分的に遮蔽された物体を認識およびセグメント化するためのレイアウト一貫性ランダムフィールド。CVPR (1) 2006: 37–44