画像のrgbそれぞれの分散を求めたい

Question

前提・実現したいこと

画像に対して、全画素のｒの平均・分散、ｇの平均・分散、ｂの平均・分散を算出したいです。平均は求まったのですが、分散ではエラーが発生してしまいます。どうすれば分散を求められますか。

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

C:\Users\Anaconda3\lib\statistics.py:535: RuntimeWarning: overflow encountered in ubyte_scalars
  T, total, count = _sum((x-c)**2 for x in data)
C:\Users\Anaconda3\lib\statistics.py:538: RuntimeWarning: overflow encountered in ubyte_scalars
  U, total2, count2 = _sum((x-c) for x in data)

該当のソースコード

import cv2
import numpy as np
import os
import statistics
import math
from statistics import variance

img = cv2.imread("C:\photo\Lena.jpg", cv2.IMREAD_COLOR)

# 対象範囲を切り出し
boxFromX = 0 #対象範囲開始位置 X座標
boxFromY = 0 #対象範囲開始位置 Y座標
boxToX = 400 #対象範囲終了位置 X座標
boxToY = 225 #対象範囲終了位置 Y座標
# y:y+h, x:x+w　の順で設定
imgBox = img[boxFromY: boxToY, boxFromX: boxToX]

# RGB平均値・分散を出力
# flattenで一次元化しmeanで平均を取得 
m_b = imgBox.T[0].flatten().mean()
m_g = imgBox.T[1].flatten().mean()
m_r = imgBox.T[2].flatten().mean()

# flattenで一次元化しvarianceで分散を取得 
v_b = imgBox.T[0].flatten()
v_g = imgBox.T[1].flatten()
v_r = imgBox.T[2].flatten()

v_b=variance(v_b)
v_g=variance(v_g)
v_r=variance(v_r)

# RGB平均値・分散を取得
print("m_B: %.2f" % (m_b))
print("m_G: %.2f" % (m_g))
print("m_R: %.2f" % (m_r))

print("v_B: %.2f" % (v_b))
print("v_G: %.2f" % (v_g))
print("v_R: %.2f" % (v_r))

試したこと

v_b = imgBox.T[0].flatten().variance()

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

画像Lenaのサイズは400*225です。

Answer 1

8ビット整数で演算しオーバーフローしてしまっているのでimg = cv2.imread(～の後にimg = img.astype(float)にて浮動小数点化して以降を演算してやればよいかと思います。

Comments

[eg88]

うまくいきました。ありがとうございます！