Python 円の重心計算

Question

前提・実現したいこと

複数の円に対して、重心計算を行いたいです。

発生している問題・エラーメッセージ

円が１つの場合は正常に作動するのですが、円が複数存在する場合、それぞれの円の重心が計算されるのではなく、重心が１つしか計算されません（その座標も正しいものではありません）
構文エラーはありませんが、円の個数と計算される重心の数が合いません。

該当のソースコード

#輪郭検出
contours,hierarchy =cv2.findContours(dst,cv2.RETR_EXTERNAL,cv2.CHAIN_APPROX_SIMPLE)
#面積計算
areas = np.array(list(map(cv2.contourArea,contours)))
# 重心の計算
max_idx = np.argmax(areas)
max_area = areas[max_idx]
result = cv2.moments(contours[max_idx])
x = int(result["m10"]/result["m00"])
y = int(result["m01"]/result["m00"])
coordination=([x],[y])
#CSVファイル作成
with open('Cordinations.csv','a',newline="") as f:
    writer=csv.writer(f)
    writer.writerow([x,y])
    f.close() 

補足情報（FW/ツールのバージョンなど）

使用言語:Python
作業環境:Visual Studio
OpenCV(4.4)

Comments

[meg_]

> 以下のようなエラーが発生しました。
エラーメッセージが掲載されていません。追記してください。

Answer 1

import cv2
import matplotlib.pyplot as plt
import imutils
%matplotlib inline

img=cv2.imread('3.png')
gray=cv2.cvtColor(img,cv2.COLOR_BGR2GRAY)
_,bw = cv2.threshold(gray,120,255,cv2.THRESH_BINARY)
#輪郭を検出。
contours,hierarchy = cv2.findContours(bw, cv2.RETR_EXTERNAL, cv2.CHAIN_APPROX_SIMPLE)
#輪郭を書き込み。
contimg=cv2.drawContours(img,contours,-1,(0,255,0),3)
count = len(contours)
plt.imshow(contimg)


cnts = cv2.findContours(bw.copy(), cv2.RETR_EXTERNAL,cv2.CHAIN_APPROX_SIMPLE)
cnts = imutils.grab_contours(cnts)
# loop over the contours
cxAll=[]
cyAll=[]
for c in cnts:
    # compute the center of the contour
    M = cv2.moments(c)
    cX = int(M["m10"] / M["m00"])
    cY = int(M["m01"] / M["m00"])
    cxAll.append(cX)
    cyAll.append(cY)
    # draw the contour and center of the shape on the image
    cv2.drawContours(img, [c], -1, (0, 255, 0), 2)
    cv2.circle(img, (cX, cY), 7, (255, 0, 0), -1)
    cv2.putText(img, "center", (cX - 20, cY - 20),
        cv2.FONT_HERSHEY_SIMPLEX, 0.5, (255, 0, 0), 2)

# show the image
cv2.imshow("Image", img)
cv2.waitKey(0)
print(cxAll)
print(cyAll)

[309, 159, 347, 102, 210]
[306, 263, 214, 180, 91]

こんにちは、こちらのページを参考に作成しました。
cv2.findContoursで、各ブロブを検出したのちにcv2.momentsで重心を計算→これをForループ内で繰り返し、、、としていけばうまく重心が計算できると思います。上の参考画像では重心が正しく計算できていることをご確認ください。お役に立てば幸いです。よろしくお願いいたします。