対象の物体の輪郭を半分だけ検出したい

Question

opencvのcv2.findContours()を利用し、マスク画像を通して対象の輪郭線を抽出することには成功したのですが、対象の輪郭線を半分だけ検出したい場合はどのような方法があるのでしょうか？
なかなか実践している人がいらっしゃらないので質問させていただきます。
イメージでは下図のようなものです。

Comments

[fana]

全体的に話が不明瞭です．
対象の輪郭線を「全部」抽出できるのであれば，単に好きなように「半分だけ」を用いればよいのではないでしょうか．

Answer 1

の輪郭線を半分だけ

の意味合いが不明瞭ですが，「輪郭線の全長の半分」という話であれば，

1. 輪郭の座標を（普通に）全部検出し，
2. 検出した座標データを半分だけ使う（不要なら「半分捨てる」）

とすれば良いのではないでしょうか．

半分というのが「どこから」「どこまで」の範囲なのかは知りませんが，とにかく必要な箇所を使って不要な箇所を使わなければよいでしょう．
「単純に，findContours が返してきた結果（輪郭の座標群）の後半を捨てる」としたコードを示します．（C++ですが）

int main()
{
	cv::Mat SrcImg = cv::imread( "Regions.png", cv::IMREAD_GRAYSCALE );	//画像読込
	if( SrcImg.empty() )return 0;	//読込失敗チェック
	cv::threshold( SrcImg, SrcImg, 128, 255, cv::THRESH_BINARY );	//２値化
	
	//ふつーにcontour検出
	std::vector< std::vector< cv::Point > > contours;
	cv::findContours( SrcImg, contours, cv::RETR_EXTERNAL, cv::CHAIN_APPROX_NONE );

	//表示用画像バッファを用意
	cv::Mat ShowImg;
	cv::cvtColor( SrcImg, ShowImg, cv::COLOR_GRAY2BGR );
	ShowImg *= 0.25;

	//contour検出結果をまずはふつーに赤で描画
	cv::drawContours( ShowImg, contours, -1, cv::Scalar(0,0,255), 2 );

	//輪郭データを半分だけにして，シアンで描画
	for( auto &cont : contours ){	cont.resize( cont.size() / 2 );	}	//※各輪郭のデータを半分捨てる
	cv::polylines( ShowImg, contours, false, cv::Scalar(255,255,0), 2 );

	//
	cv::imshow( "ShowImg", ShowImg );
	if( cv::waitKey() == 's' ){	cv::imwrite( "Regions_result.png", ShowImg );	}
	return 0;
}

処理結果例：

• 用いた入力画像

• 処理結果画像

（輪郭の半分を赤，もう半分をシアンで描画している）

Answer 2

あってますか？わからないことは英語で探すと大体見つかります。

import cv2
import numpy as np

img_gray= cv2.cvtColor(cv2.imread("./tri.png"),cv2.COLOR_BGR2GRAY)


ret, thresh = cv2.threshold(img_gray,127,255,0)
cv2.imshow("thresh",thresh)

cnts, _ = cv2.findContours(thresh,cv2.RETR_TREE,cv2.CHAIN_APPROX_SIMPLE)

img_mask = np.zeros_like(img_gray)
img_mask = cv2.drawContours(img_mask, cnts, 1, 255,5)

cv2.imshow("img_mask",img_mask)

out = img_mask.copy()
out_half = img_mask.copy()
# https://stackoverflow.com/questions/28759253/how-to-crop-the-internal-area-of-a-contour
(y, x) = np.where(img_mask == 255)
(topy, topx) = (np.min(y), np.min(x))
(bottomy, bottomx) = (np.max(y), np.max(x))
out = out[topy:bottomy+1, topx:bottomx+1]

out_half = out_half[topy:bottomy+1, int(topx/2):int(bottomx/2)+1]
cv2.imshow("out",out)


cv2.imshow("out_half",out_half)
cv2.imwrite("out_half.png",out_half)

cv2.waitKey(0)

Comments

[fana]

> あってますか？

と問うならば，せめて
「輪郭線を半分だけ検出したい」とかいうぼんやりとした質問文を
「どのような結果を目的とする話なのだと捉えたのか」に関しては述べねばならないのではないでしょうか．

（処理内容に関してはコードを解読しろ，というスタンスだとしても）

[fana]

単にROIを左右に真っ二つにしただけに見えるのですが，
領域形状が　△ ではなく ⊿ みたいな場合でも，「輪郭線を半分」と言える結果が得られるものでしょうか？
（他にも　ヨ　みたいな形のときだとか）

[退会済みユーザー]

> 単にROIを左右に真っ二つにしただけに見えるのですが，
領域形状が　△ ではなく ⊿ みたいな場合でも，「輪郭線を半分」と言える結果が得られるものでしょうか？
> （他にも　ヨ　みたいな形のときだとか）
処理内容はまさに真っ二つにしただけで、「一番左と一番右の座標を調べて、絵に掲載された通り真っ二つに分断した処理」になっています。輪郭線の長さを見て云々等はしていませんので、"ヨ"であれば"≡"とよ"ﾖ"のようになります。

質問した方の絵をなるべく再現する形でプログラムに起こし、質問した方が期待する（と私が思う）結果ができますよ、と計算結果を図示しました。

私の理解では、質問した方の説明文をと絵を見て、質問文の"イメージ"から「△の輪郭の左半分（なのか右半分）を検出したい」というように読み取りました。「全体の輪郭線の長さを見て半分なのか面積を見て半分」なのか、「他の何をもって半分」なのか、もしかしたらいろいろあるかもしれませんが、私にはそこまで深い意図をこの絵から私は読み取れませんでした。

と、ここまでで―－パッと見て左右に半分以上に割りたい以外の意図を読み取れませんでしたし、気にもしていませんでした。質問の内容から私が勝手にそう思っているだけの話で、実はやりたいことが違うかもしれない、と思ってきました。

というのも、話の感じからfanaさんは、私よりもいろいろな人の意見を聞いてきて、「もう既にプログラムに起こしてしまったよ。そんなことなら初めからもっと詳しく条件を教えてくれればよかったのに、無駄骨だったか――」というようなプログラムとは別の実務経験で鍛えられてきたような厚み（？）が感じられたためです。

[fana]

１つ目のコメント（2021/09/02 10:42 ）の意図は，

単に唐突に「あっているか？」とだけ書かれると判断に困る（判断するためにはコードを完全解読せねばならない）であろうから，
「x座標を基準にして左右方向に関して真っ二つにする　という話であれば，こうすればできるが，そんな話であっているか？」
みたいな形で，何を示しているのかを明確に（コードとは別に）述べる方が望ましいのではないか，というだけの話です．

---

２つ目のコメント（2021/09/02 10:47 ）の方は，実際に私がコードを見てみたらこんな内容だと思われた，という話ですね．
「ROIを左右に真っ二つ」とかいう語があるだけで，質問者がコードを見る際のヒントになり得るかな，と．
（で，ついでに，判断基準や懸念点(?)になりそうな話を示しておいた感じです．）