intlut() alternative for LUT consisting of doubles

Question

Stefan 2012-8-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/46062-intlut-alternative-for-lut-consisting-of-doubles

Is there a fast alterative to intlut() accepting an array of doubles as LUT?

Problem is, I want to convert grayscale values (uint8) of an image using a LUT consisting of 256 doubles.

Example:

lut = (1:256)+0.1; % upshift by 0.1
A = uint8(randi(256,50));
% current approach
lutMat = repmat(reshape(lut, 1, 1, numel(lut)), ...
        size(A));
idx = sub2ind( ...
  size(lutMat), ...
  repmat((1:size(A, 1))', 1, size(A, 2)), ...
  repmat((1:size(A, 2)) , size(A, 1), 1), ...
  double(A));
B = lutMat(idx);

The above code works fine, but takes a lot of time and needs too much memory

3 个评论
显示 1更早的评论隐藏 1更早的评论

Stefan 2012-8-20

在 MATLAB Online 中打开

I forgot to add the line

 lut = repmat(reshape(lut, 1, 1, numel(lut)), ...
          size(A));

just after

% current approach

I changed the code in the original question, accordingly.

Image Analyst 2012-8-20

You're still having the problem? Didn't my code do what you want?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Jan 2012-8-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/46062-intlut-alternative-for-lut-consisting-of-doubles#answer_56579

在 MATLAB Online 中打开

I'm not sure if I understand the question. Does this help:

B = lut(A);

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Image Analyst 2012-8-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/46062-intlut-alternative-for-lut-consisting-of-doubles#answer_56585

编辑：Image Analyst 2012-8-18

在 MATLAB Online 中打开

If you're just adding a constant value to the lut-transformed values, like 0.1, you can just do this:

lut = uint8(1:256); % upshift by 0.1
A = uint8(randi(256,50));
B = double(intlut(A, lut)) + 0.1;

If it's not some constant added value, then you can do it like this:

lut = 1000 * rand(256, 1); % A Random reassignment between 0 and 1000.
A = uint8(randi(256,50));
B2 = zeros(size(A), class(lut));
for inputValue = uint8(0) : uint8(255)
   indexes = (A == inputValue);
   B2(indexes) = lut(inputValue + 1);
end