はじめに
画像の輪郭をsobelフィルタにより抽出し、抽出した輪郭ごとにラベリングする方法について、都道府県の境界が表示された日本地図のイラストを例として、説明する。
コード
![](https://sabopy.com/wp/wp-content/uploads/2019/12/plot_label_label2rgb.jpg)
![](https://sabopy.com/wp/wp-content/uploads/2019/12/plot_label_japan_2.jpg)
解説
モジュールのインポート
画像の読み込み
用いる画像は下記サイトから取得した。
![](https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhyrjymJd2UQmx5oKVED5CZHKtsdN_X9CPt2xkD1fTW1FqiPDG1RhHPvAX4iB_Pxzsle80c2Q-aE59XytLEUaQ6cVdzuVf-Y369nvVvJy_owqSSkpOvV79_WpykKJUsAYgwyNfgsJLG1gp5/s400/nihonchizu_area.png)
以下の画像を用いる。
![](https://sabopy.com/wp/wp-content/uploads/2019/12/nihonchizu_area.jpg)
img_cはRGB画像であり、それをrgb2grayでグレースケール化したものをimageとする。
sobelフィルタの適用
sobelフィルタで画像の輪郭を抽出する。sobelフィルタを適用した画像は以下のようになる。
![](https://sabopy.com/wp/wp-content/uploads/2019/12/plot_sobel.jpg)
sobelフィルタを適用した画像の2値化
値として0.04より小さい部分をTrueとすると以下のような画像が得られる。
![](https://sabopy.com/wp/wp-content/uploads/2019/12/plot_sobel_binali.jpg)
ノイズの除去
opening処理を2値化像に適用することでノイズを除去する。ここでは2回openingした。結果として以下の画像が得られる。
![](https://sabopy.com/wp/wp-content/uploads/2019/12/plot_opening.jpg)
画像のラベリング
label(openimg)でさっきの画像をラベリングする。
label2rgbでラベルごとにcmapsの色を適用して表示する。結果は以下のようになる。
![](https://sabopy.com/wp/wp-content/uploads/2019/12/plot_label_label2rgb.jpg)
重心とラベル番号の表示
regionprops(label_image)のcentroidでラベル領域の重心が得られるので、その部分に赤丸とラベル番号を表示すると下のようになる。
![](https://sabopy.com/wp/wp-content/uploads/2019/12/plot_label_center.jpg)
ラベル領域を切り取って表示
regionprops(label_image).bboxでラベル領域を囲む四角形の座標が得られるので、その部分で画像を切り取り、一覧を表示したのが下の図となる。
この場合、他の領域も周りに混ざってくるのでわかり難い図となっている。
![](https://sabopy.com/wp/wp-content/uploads/2019/12/plot_label_japan_1.jpg)
ラベリングした部分のみを切り取る
im_label = (licrop==region.label)でラベル番号が一致する部分のみをTrueとしたbool配列を得る。
この配列をもとのRGB画像と掛け合わせる。bool配列はTrueが1、Falseが0で計算されるのでTrueの部分以外は0となり黒くなり、ラベリングした部分のみが残る。例を以下に示す。
![](https://sabopy.com/wp/wp-content/uploads/2019/12/mawari_jokyo.jpg)
各都道府県を表示
すべてのラベリングした領域に上記の処理を行うことで各都道府県を分離できる。
![](https://sabopy.com/wp/wp-content/uploads/2019/12/plot_label_japan_2.jpg)
参考
![](https://sabopy.com/wp/wp-content/uploads/2021/01/r_s_filter_obesa-1-160x90.jpg)
![](https://sabopy.com/wp/wp-content/uploads/2019/09/plot_label-160x90.jpg)
![](https://sabopy.com/wp/wp-content/uploads/2019/11/plotregionprops-160x90.png)
コメント