同じような処理を２度実行したい　python opencv

画像はkuro1.pngで１つの画像を２回に分けて処理したいです

１つ目の処理と２つの目の処理は少しずつ違っているのですが
ラベリング処理は両方行います

しかし下記のようなプログラムではエラーが出ます
どう直せば良いでしょうか

import cv2
import numpy as np
from matplotlib import pyplot as plt



（（（この部分は後から使いたいとおもって書いているだけ
# 変換前後の対応点を設定
p_original = np.float32([[247,388], [805,408], [17, 833], [1063, 860]])
p_trans = np.float32([[0,0], [500,0], [0,690], [500,690]])
 
# 変換マトリクスと射影変換
M = cv2.getPerspectiveTransform(p_original, p_trans)
））



src= cv2.imread("kuro1.png")


#色抽出
dst1 = cv2.inRange(src, (3,3,3), (10,8,10))
dst2 = cv2.blur(dst1, ksize=(5,5))
cv2.imshow("b",dst2)

ret,thresh1 = cv2.threshold(dst2,127,255,cv2.THRESH_BINARY)

# ラベリング処理(詳細版)
label = cv2.connectedComponentsWithStats(thresh1)

### 追加した部分1
stats = label[2]
area = stats[:, cv2.CC_STAT_WIDTH] * stats[:, cv2.CC_STAT_HEIGHT]
top2_idx = area.argsort()[-3:-1]  # 2番目、3番めに面積が大きいラベル (1番目は背景なので除く)

# オブジェクト情報を項目別に抽出
data = label[2]
center = label[3]

# ラベリング結果書き出し用に二値画像をカラー変換
color_src = cv2.cvtColor(thresh1, cv2.COLOR_GRAY2BGR)

# オブジェクト情報を利用してラベリング結果を表示
for i in top2_idx:
    # 各オブジェクトの外接矩形を赤枠で表示
    x0 = data[i][0]
    y0 = data[i][1]
    x1 = data[i][0] + data[i][2]
    y1 = data[i][1] + data[i][3]
    cv2.rectangle(color_src, (x0, y0), (x1, y1), (0, 0, 255))




    # 各オブジェクトの重心座標をに黄文字で表示
    cv2.putText(color_src, "X: " + str(int(center[i][0])), (x1 - 10, y1 + 15), cv2.FONT_HERSHEY_PLAIN, 1, (0, 255, 255))
    


  

















#「背景」となる画像の取り込み（グレースケール）
img_src01 = cv2.imread("zero1.png", 0)
 
#「差分」をもった画像の取り込み（グレースケール）
img_src02 = cv2.imread("kuro1.png", 0)
 
#「背景差分」計算用オブジェクトの作成　
bgObj = cv2.bgsegm.createBackgroundSubtractorMOG()
 
#差分となっている「前景領域」に対してマスクをかける
fgmask = bgObj.apply(img_src01)
fgmask = bgObj.apply(img_src02)

# グレースケール化
gray = cv2.cvtColor(fgmask, cv2.COLOR_BGR2GRAY)

# ラベリング処理(詳細版)
label2 = cv2.connectedComponentsWithStats(gray)

### 追加した部分1
stats2 = label2[2]
area2 = stats2[:, cv2.CC_STAT_WIDTH] * stats2[:, cv2.CC_STAT_HEIGHT]
top2_idx2 = area2.argsort()[-3:-1]  # 2番目、3番めに面積が大きいラベル (1番目は背景なので除く)

# オブジェクト情報を項目別に抽出
data2 = label2[2]
center2 = label2[3]

# ラベリング結果書き出し用に二値画像をカラー変換
color_src2 = cv2.cvtColor(gray, cv2.COLOR_GRAY2BGR)

# オブジェクト情報を利用してラベリング結果を表示
for i in top2_idx2:
    # 各オブジェクトの外接矩形を赤枠で表示
    x0 = data2[i][0]
    y0 = data2[i][1]
    x1 = data2[i][0] + data2[i][2]
    y1 = data2[i][1] + data2[i][3]
    cv2.rectangle(color_src2, (x0, y0), (x1, y1), (0, 0, 255))

    # 各オブジェクトの重心座標をに黄文字で表示
    cv2.putText(color_src2, "Y: " + str(int(y1)), (x1 - 10, y1 + 30), cv2.FONT_HERSHEY_PLAIN, 1, (0, 255, 255))



cv2.imshow('haikei',color_src2)



cv2.waitKey(0)
cv2.destroyAllWindows()

エラー文

    gray = cv2.cvtColor(fgmask, cv2.COLOR_BGR2GRAY)
cv2.error: OpenCV(4.1.1) /Users/travis/build/skvark/opencv-python/opencv/modules/imgproc/src/color.simd_helpers.hpp:92: error: (-2:Unspecified error) in function 'cv::impl::(anonymous namespace)::CvtHelper<cv::impl::(anonymous namespace)::Set<3, 4, -1>, cv::impl::(anonymous namespace)::Set<1, -1, -1>, cv::impl::(anonymous namespace)::Set<0, 2, 5>, cv::impl::(anonymous namespace)::SizePolicy::NONE>::CvtHelper(InputArray, OutputArray, int) [VScn = cv::impl::(anonymous namespace)::Set<3, 4, -1>, VDcn = cv::impl::(anonymous namespace)::Set<1, -1, -1>, VDepth = cv::impl::(anonymous namespace)::Set<0, 2, 5>, sizePolicy = cv::impl::(anonymous namespace)::SizePolicy::NONE]'
> Invalid number of channels in input image:
>     'VScn::contains(scn)'
> where
>     'scn' is 1

このようなエラーが出るのですが全くわかりません泣
アドバイスよろしくお願いします