Matlab中按元素进行矩阵除法仅返回二进制值

Question

我正在Matlab中读取两个图像作为矩阵，并想要创建第三个图像，其计算公式为：

(im1-im2)./(im1+im2)

但是，我相信Matlab正在舍入按元素划分的值（我当然不希望这样做）。 例如， im1-im2的第一值为32 ， im1+im2的第一值为70 。 因此，结果矩阵的第一个值应该是32/70=0.4571 ，但是它是0 。

如何执行此矩阵运算，将最终除法结果值保留在最终矩阵中？

Answer 1

这取决于您的图像数据及其加载方式。 加载有例如imread的8位图像存储为整数矩阵（ uint8 ）。 位深度较高的图像将使用uint16或更多存储。

因此，您可以通过将矩阵转换为double并先将其转换回适当的位深度来处理图像，也可以直接将其转换为所需的位深度并接受舍入，因为无论如何都需要在结束。

im1 = double( imread('fig1.png') );
im2 = double( imread('fig2.png') ); 

result = uint32( (im1-im2)./(im1+im2) );

Answer 2

im1和im2最有可能是uint8或uint16 ，将它们转换为double并执行此操作。 请注意，将它们隐式转换为*1.0无效。

所以，像

im1 = double(im1);
im2 = double(im2);
thingYouWant = (im1 - ...;

请注意，转换回uint会再次为您提供1（值大于0.5）或0（值小于0.5）。 要将其转换回uint8 ，您需要执行以下操作：

resultInUint8 = uint8(thingYouWant * 255); % Linear scale from 0 to 255

要注意的另一件事是，如果im2大于im1 ，则double(im1 - im2)结果为0。 因此，您可能最终将图像的一半设为0。在进行任何操作之前最好进行转换。