[論文レビュー] Gradient-Domain Processing for Large EM Image Stacks

Michael Kazhdan, Randal Burns+4 more|arXiv (Cornell University)|Sep 30, 2013

Advanced Optical Imaging Technologies参考文献 18被引用数 3

ひとこと要約

本稿では、大規模な電子線画像スタックにおけるスライス間の不連続性を除去しつつスライス内の詳細を保持する二段階の勾配ドメイン法を提案する。本手法はスライス軸に沿った異方的拡散を実行し、その後各スライスでスクリーネッド・ポアソン再構成を適用する。さらに、12コアと45GBのRAMを用いて60時間未満でテラボクセルデータセットを処理できる、新規の外部メモリ対応マルチグリッドソルバにより高速化されている。

ABSTRACT

We propose a new gradient-domain technique for processing registered EM image stacks to remove the inter-image discontinuities while preserving intra-image detail. To this end, we process the image stack by first performing anisotropic diffusion to smooth the data along the slice axis and then solving a screened-Poisson equation within each slice to re-introduce the detail. The final image stack is both continuous across the slice axis (facilitating the tracking of information between slices) and maintains sharp details within each slice (supporting automatic feature detection). To support this editing, we describe the implementation of the first multigrid solver designed for efficient gradient domain processing of large, out-of-core, voxel grids.

研究の動機と目的

• 照明の変動、カメラパラメータの違い、切断に起因するアーティファクトによって生じる大規模なEM画像スタックにおけるスライス間の不連続性を解消すること。
• 処理中に個々のスライス内の高周波数の解剖学的詳細を保持すること。
• テラボクセル規模の外部メモリ対応3次元ボクセルデータに対する効率的でスケーラブルな勾配ドメイン処理を可能にすること。
• 大規模で外部メモリ対応のボクセルグリッドに特化したマルチグリッドソルバを開発すること。
• 2次元の勾配ドメイン画像処理を、接続体解析や神経科学分野の応用を想定して3次元ボリュームデータに拡張すること。

提案手法

• 本手法はまず、z軸（スライス軸）に沿って異方的拡散を適用し、スライス間の不連続性を平滑化しつつ、xおよびy方向のスライス内勾配を保持する。
• 異方的拡散を、出力の勾配と目標勾配場との差の重み付きL2ノルムの最小化として定式化し、z成分の重みを低く設定する。
• その後、各スライス内で独立にスクリーネッド・ポアソン方程式を解き、元のデータからの高周波数の詳細を再導入しながらスライス間の連続性を維持する。
• z成分を抑制した重み付き勾配場を用い、スクリーネッド・ポアソン方程式は平滑化されたデータへの忠実度と元の勾配とのバランスをとる。
• スクリーネッド・ポアソン方程式に由来する大規模な線形方程式系を、テラボクセル規模のデータセット上で効率的に解けるように、新規の並列処理・外部メモリ対応マルチグリッドソルバを実装する。
• メモリ内処理とディスクベース処理の両方をサポートし、適応的Vサイクル戦略と大規模ボクセルグリッドに最適化されたメモリレイアウトを採用する。

実験結果

リサーチクエスチョン

• RQ1 勾配ドメイン処理により、スライス内解剖学的詳細のぼかしを伴わずに、大規模なEM画像スタックにおけるスライス間不連続性を効果的に除去できるか？
• RQ2 スライス軸に沿った異方的拡散をどのように適用することで、スライス間アーティファクトを低減しつつ、平面内勾配を保持できるか？
• RQ3 勾配ドメイン処理の文脈において、外部メモリ対応の大規模3次元ボクセルグリッドに対するマルチグリッドソルバの性能とスケーラビリティはいかがなものか？
• RQ4 異方的拡散とスクリーネッド・ポアソン再構成の組み合わせにより、EMデータにおいてスライス間連続性とスライス内シャープネスの両立が可能か？
• RQ5 定数、ハイブリッド、線形の各ソルバの選択が、大規模3次元ボリューム処理における収束性と実行時間に与える影響は？

主な発見

• 提案された二段階手法は、元データ、拡散後データ、再構成後データの視覚的比較により、スライス間不連続性を効果的に除去するとともに、微細な解剖学的詳細を保持していることが実証された。
• 異方的拡散段階では、ADおよびSPフェーズの両方で残差ノルムが約4.7×10⁻⁴に低下し、ハイブリッドソルバが速度と収束性のバランスが最良であった。
• マルチグリッドソルバは、1850×21496×25792のEMスタックを12コアと45GBのRAMを用いて60時間未満で処理を完了し、テラボクセルデータセットへのスケーラビリティを示した。
• ハイブリッドマルチグリッドソルバは、10⁵秒あたり0.53の残差減衰率を達成し、実行時間の観点からも定数および線形バージョンを上回る効率性を示した。
• 本手法は、1920×1080×742の動画データに対しても非線形フィルタリングに応用可能であり、小さな勾配を抑制することでノイズを低減した。1チャンネルあたりのピークメモリ使用量はわずか185MBであった。
• 実装は、外部メモリ対応3次元データにおける一般化された勾配ドメインフィルタリング、例えばHDR圧縮、シャープニング、非フォトリアリスティックレンダリングをサポートしている。