画素の総数が知りたい（マルチスケール）

Question

以下のコードで２つのマルチスケール画像を比較して、一致しない画素値の総数を出力します。
２つの画像はグレースケール画像なので、cv2.imreadmultiの引数を0にしているのですが、
これであっていますか？ANY_COLORの時と値が変わりません。

import cv2
import numpy as np

# 画像を比較し、一致していない箇所を表示する

# 元画像
# 単一画像
#img_original = cv2.imread("a/569_OUT.jpg")

#マルチスケール画像
r, img_original = cv2.imreadmulti('MVI354_32/Target/00065_00096.tif', [], 0)

# 比較対象画像(単一画像)
#img_comparison1 = cv2.imread("a/570_OUT.jpg")
#マルチスケール画像
r, img_comparison1 = cv2.imreadmulti('test3_3D/32/out02.tiff', [], 0)

#比較画像の画素値が一致しない総数を出力
print(sum(np.sum(img_original[i] != img_comparison1[i]) for i in range(len(img_comparison1))))


'''
# 画素がどのくらい一致しているかを確認する(一致している画素数を出力)
print(np.count_nonzero(img_original == img_comparison1))
print(np.count_nonzero(img_original == img_comparison2))
'''

Answer 1

２つの画像はグレースケール画像なので、cv2.imreadmultiの引数を0にしているのですが、

これであっていますか？ANY_COLORの時と値が変わりません。

合っていると思います。

理由は以下の通りです。

公式のFlags used for image file reading and writingによると、

定数	値	説明
cv::IMREAD_GRAYSCALE	0	If set, always convert image to the single channel grayscale image (codec internal conversion).
cv::IMREAD_ANYCOLOR	4	If set, the image is read in any possible color format.

とあります。

つまり、
グレースケール画像-->IMREAD_GRAYSCALE(0)-->グレースケール画像
グレースケール画像-->IMREAD_ANYCOLOR(4)--->グレースケール画像（最適な条件で読んでも元のグレースケールのまま）

マルチスケール画像-->IMREAD_GRAYSCALE(0)-->グレースケール画像
マルチスケール画像-->IMREAD_ANYCOLOR(4)--->マルチケール画像（最適な条件で読み込むと元のマルチスケール）

---

もし意図と違う回答でしたら、実際の画像をアップロードした上でコメントください。

Comments

[tiroha]

ありがとうございます。
画像は訳あって、公開することができません。すいません。
RGBではなくちゃんとグレースケールで呼び出しができているということでしょうか？
実際の出力が132286です。
8フレーム重ねた画像で、サイズは256×256。白黒画像。
フレーム数で割った値が
132286/8=16535
となります。

[退会済みユーザー]

> RGBではなくちゃんとグレースケールで呼び出しができているということでしょうか？

元がグレースケールであれば、「マルチスケールに変換(例えばimg = cv2.cvtColor(img,cv2.COLOR_GRAY2BGR))する処理をしない限り」読み込んだ画像はグレースケールのはずです。

どうしても心配であれば、以下のようにすると、(height,width,channel)を確認できます。
print(img_original.shape)
print(img_comparison1 .shape)

グレースケールだとchannel情報がないので(height,width)の形で返ってくるか、特殊なフォーマットであれば(height,width,1)のように返ってくるはずです。

参考：
> 画像は訳あって、公開することができません。
代わりに縦横チャネル数、数値の型は実物の画像に沿ったものにして、ダミー画像をnp.ones((height,width,channel),np.uint8)*128みたいな感じで作れます。

[tiroha]

ありがとうございます。
そうなんですね。
教えていただいた、２つのprint分を記述して実行すると、以下のエラーが出ます。
１枚のみの画像（マルチスケール画像ではない）ものだとうまく出力されるのですが、マルチスケール画像だとうまくいかないみたいです。
Traceback (most recent call last):
File "test_image_multi.py", line 20, in <module>
print(img_original.shape)
AttributeError: 'list' object has no attribute 'shape'

[退会済みユーザー]

マルチスケール画像はたぶん[img,img]のように画像がリストに格納されるっぽいですね。
やり方は汚いですが、とりあえず以下のようにやると中身はみえるはずです。
[print(img_single.shape) for img_single in multiscale_img]

[tiroha]

ありがとうございます。
実行すると、以下のように出てきます。
File "test_image_multi.py", line 20
print(img_single.shape) for img_single in multiscale_img
^
SyntaxError: invalid syntax
画像は２枚ともマルチスケール画像なので、img_single.shapeが出てくるのは違和感があるのですが、どういうことでしょうか？

[退会済みユーザー]

[]で囲ってください

[tiroha]

ありがとうございます。
[]で括ると、以下のエラー文が出てきます。
Traceback (most recent call last):
File "test_image_multi.py", line 20, in <module>
[print(img_single.shape) for img_single in multiscale_img]
NameError: name 'multiscale_img' is not defined

[tiroha]

[print(img_original.shape) for img_original in img_comparison1]
であっていますかね？
実行すると、
(256,256)
というのが１１行出力されます。

[退会済みユーザー]

multiscale_imgは、img_originalとimg_comparison1に該当します。

どちらのshapeを確認したいときはmultiscale_imgをimg_originalで回すのと、img_comparison1で回すのをして確認してください。

[tiroha]

ありがとうございます。そうなんですね。
img_comparison1の時は（256,256）が16行。
img_originalの時は（256,256）が16行。
出力されました。
これってチャネル数が省略されているので、ちゃんとグレースケールで呼び出せているということであっていますよね？

[退会済みユーザー]

そうなります、グレースケールだと(256,256)がxyとして吐き出されます。ちゃんとグレースケールになっているとみてよいです。

[tiroha]

ありがとうございます。