FRETImageProcess2.RegionGrowing関数
FRETImageProcess2.SortNLSNESCAAX関数の内部で、
細胞全体の領域を求めるのに使われる関数。
ほかの機会でこの関数を呼び出すことはあまり無い。
NESまたはCAAXのタグがついたFRETプローブ発現細胞に関して、
YFPまたはCFPの画像と、核のバイナリイメージから
領域拡張法アルゴリズムを用いて細胞全体の領域を求め、
バイナリイメージを出力する。
(ソースコードには前バージョンのコードと注釈が
一部コメントアウトされて残っている)
この関数での拡張条件は、seed pointに隣接した座標のintensityが、
seed pointのintensityの0.9倍を超えていることである。
解像度によって最適な条件は変わる。
GrownBinaryImage = RegionGrowing(image, SeedBinaryImage)
引数
image (2次元行列またはTiffオブジェクト)
wavelengthがFRET(YFP)の画像データを格納した数値行列またはTiffオブジェクト。
CFP画像はFRET画像より暗いことが多いのでFRET画像の方が好ましい。
また、核の位置にずれがないよう、
segmentationに用いた核の画像とできるだけ近いtimepointの画像がよい。
数値のデータ型は何でもよい。
小数値を含む場合はround関数により四捨五入される。
SeedBinaryImage (logical型の2次元行列)
一細胞の核の領域を示したバイナリイメージ。
FRETImageProcess2.Segmentation 関数の戻り値の
SegmentationMatrixから以下のようなコードで作成する。
% -------------------------------------------
nucleus_region = SegmentationMatrix == i;
nucleus_region = imdilate(nucleus_region, strel('disk',1));
% -------------------------------------------
iは細胞の番号であるが、核の境界線は0の数値が割り当てられているので、
imdilate関数で1ピクセル分拡張する。
戻り値
GrownBinaryImage(logical型の2次元行列)
SeedBinaryImageの核を持つ
細胞全体を示したバイナリイメージ