我正在尝试对对数比率矩阵进行归一化,因此,为了做到这一点,我想找到矩阵的最大值。但是我得到了无限,这是不可能的。
我编写的代码:
import imageio as im
import numpy as np
imagepath1 = 'Andasol_09051987.jpg'
imagepath2 = 'Andasol_09122013.jpg'
image1 = im.imread(imagepath1)
image2 = im.imread(imagepath2)
Ds = np.abs(image1 - image2)
Dl = np.abs(np.log(image2+1)-np.log(image1+1))
Dsmax = Ds.max()
Dsmin = Ds.min()
Ds = ((Ds - Dsmax)/(Dsmax - Dsmin))*255
Dlmax = np.amax(Dl)
Dlmin = Dl.min()
Dl = ((Dl - Dlmax)/(Dlmax - Dlmin))*255
对于减法Ds部分,一切正常,但Dl部分不起作用。 Dlmax的值是无限的。
用于计算对数比
Dl = np.abs(np.log(image2+1)-np.log(image1+1))
它有一个警告RuntimeWarning:除以日志中遇到的零 Dl = np.abs(np.log(image2 + 1)-np.log(image1 + 1))
我真的想避免被0除,这就是为什么我向每个像素加1。
两个图像都是灰度图像,因此每个像素的值从[0 255]开始变化