はじめに
skimage.filtersのthreshold_otsu、morphologyのopening, closingなどによって、画像中の対象物をラベリングする方法について説明する。
コード
![](https://sabopy.com/wp/wp-content/uploads/2019/09/plot_label-800x600.jpg)
解説
モジュールのインポート
画像データの読み込み
昨年蒔いたアストロフィツム・ストロンギロゴナムの幼苗の画像を用いる。
![](https://sabopy.com/wp/wp-content/uploads/2019/09/strongi-2.jpg)
大津の2値化
![](https://sabopy.com/wp/wp-content/uploads/2019/09/plot_label_threshold_otsu-800x600.jpg)
大津の2値化により単一のしきい値を求め、しきい値処理を行うことでバイナリ画像を生成する。
クロージングによる穴埋め
![](https://sabopy.com/wp/wp-content/uploads/2019/09/plot_label_closing-800x600.jpg)
膨張処理後に収縮処理を行うクロージングによって、白い領域中の黒点を除去することができる。
オープニングによる小さい領域の除去
![](https://sabopy.com/wp/wp-content/uploads/2019/09/plot_label_opening-800x600.jpg)
収縮処理後に膨張処理を行うオープニングによって、小さな白い領域を除去することができる。
画像の端に接している領域の除去
![](https://sabopy.com/wp/wp-content/uploads/2019/09/plot_label_clear_border-800x600.jpg)
clear_borderにより画像の端に接している白い領域を除去する。
ラベリングし、ラベルごとに色を変えて表示
![](https://sabopy.com/wp/wp-content/uploads/2019/09/plot_label_label2rgb-800x600.jpg)
label(バイナリ画像)で連結している部分ごとにラベル付けを行う。label_imageはバイナリ画像と同じ形状で整数の配列となる。
label2rgbにより、ラベリングされた部分に色をつけることができる。
ラベリングした領域を四角で囲む
![](https://sabopy.com/wp/wp-content/uploads/2019/09/plot_label_rectangle-800x600.jpg)
regionprops(label_image)によってラベリングした領域に関する多くの情報を得ることができる。region.area によってラベリングした領域の面積が得られる。ここでは、100ピクセル以上の面積をもつものに対して、region.bboxでラベル画像を囲む最小の四角の座標を取得し、緑線でRectangleを描いた。
regionpropsで得られる情報
areaで面積、labelでラベル番号、centroidで重心の値を得ることができる。その他については下記が詳しい。
skimage.measure — skimage 0.24.0 documentation
ラベリングした画像をそれぞれ表示
![](https://sabopy.com/wp/wp-content/uploads/2019/09/plot_label_cactus-800x600.jpg)
面積が100以上のものに対してregion.bboxでラベル画像を囲む最小の四角の座標を取得し、その座標で切り取りを行い、切り取った画像のリストを作成した。
そのリストをImagegridを用いて表示した。タイトルには切り取った画像のラベル番号を記している。
参考
Label image regions — skimage 0.24.0 documentation
skimage.filters — skimage 0.24.0 documentation
skimage.segmentation — skimage 0.24.0 documentation
skimage.measure — skimage 0.24.0 documentation
skimage.measure — skimage 0.24.0 documentation
![](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn.qiita.com%2Fassets%2Fpublic%2Farticle-ogp-background-412672c5f0600ab9a64263b751f1bc81.png?ixlib=rb-4.0.0&w=1200&mark64=aHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZ3PTk3MiZoPTM3OCZ0eHQ9bWF0cGxvdGxpYiVFMyU4MSVBN2NvbG9yJTIwY3ljbGUlRTMlODElQUVOJUU3JTk1JUFBJUU3JTlCJUFFJUUzJTgxJUFFJUU4JTg5JUIyJUUzJTgyJTkyJUU2JThDJTg3JUU1JUFFJTlBJUUzJTgxJTk5JUUzJTgyJThCJUUzJTgxJTg0JUUzJTgxJThGJUUzJTgxJUE0JUUzJTgxJThCJUUzJTgxJUFFJUU2JTk2JUI5JUU2JUIzJTk1JnR4dC1hbGlnbj1sZWZ0JTJDdG9wJnR4dC1jb2xvcj0lMjMyMTIxMjEmdHh0LWZvbnQ9SGlyYWdpbm8lMjBTYW5zJTIwVzYmdHh0LXNpemU9NTYmcz0wZTc1YTI4MmQ1NzQ4YmRkZTFlNmRmZDYyYjUwNGQwNw&mark-x=142&mark-y=57&blend64=aHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZoPTc2Jnc9NzcwJnR4dD0lNDBza290YXJvJnR4dC1jb2xvcj0lMjMyMTIxMjEmdHh0LWZvbnQ9SGlyYWdpbm8lMjBTYW5zJTIwVzYmdHh0LXNpemU9MzYmdHh0LWFsaWduPWxlZnQlMkN0b3Amcz0xODIzNjQwMGY3OWM2YTFhYjhiY2NkMDEwYWU3YzdhNQ&blend-x=142&blend-y=486&blend-mode=normal&s=34c2577b95e3ccd9e7f986ac6c7e1d81)
matplotlibでcolor cycleのN番目の色を指定するいくつかの方法 - Qiita
実際の活用場面に即したタイトルにすると__「同じ図に複数データをプロットする際にcolor cycleから同じ色を指定する方法」とでもなりますが、少し一般性を持たせるためにどういう作業をするかに注目…
コメント