Move all counts of the image in random direction

Question

Alex Kurek 2016-2-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/268615-move-all-counts-of-the-image-in-random-direction

评论： Alex Kurek 2016-2-19

I have the following problem. I have 8-bit image. It represents photon counts, so each pixel is capable to register up to 255 photons. I need to introduce some random error - each photons trajectory is disturbed, so that finaly instead of e.g. (150, 100) it is regiesterd at (150+random1, 100+random2). Random1 and random2 both have Poisson distirbution and are small, ~2 pixels.

Notice, that is is not enough to displace each pixel in random direction, since each count from each pixel has to be displaced separately.

What is the easiest and fastest way to code it?

Best regards, Alex

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

请先登录，再进行评论。

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

Answer 1

Image Analyst 2016-2-17

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/268615-move-all-counts-of-the-image-in-random-direction#answer_210176

在 MATLAB Online 中打开

You need to scan your original image, and if the value is not zero (meaning that there is a photon in the original CCD well or whatever the image represents) then you need to calculate the new location (row,col) by adding a deltax and deltay taken from a Poisson distribution. Then add the photon to the output location and (optionally) remove it from the input location.

outputImage = zeros(rows, columns);
for col = 1 : columns
    for row = 1 : rows
        if originalImage(row, col) > 1
            newRow = round(row + random('Poisson', pd));
            newCol = round(col + random('Poisson', pd));
            % Increment the number of pixels at this output image location.
            outputImage(newRow, newCol) = outputImage(newRow, newCol) + 1;
            % Decrement the number of photons at this output image location (only if you want to)
            originalImage(row, col) = originalImage(row, col) - 1;
        end
    end
end

See documentation for random() in the Statistics and Machine Learning Toolbox for info on how to get pd, the probability distribution object.

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

Image Analyst 2016-2-18

You could perhaps get all the Poisson distances in 256 arrays in advance instead of getting them one number at a time.

Alex Kurek 2016-2-19

在 MATLAB Online 中打开

Yes, it worked. Like this it is >300x faster:

countNo = sum(sum(imageTest));
randContainer = round(normrnd(0, pd, [1 2*countNo]));
i = 1;
tic
for col = 1 : columns
    for row = 1 : rows
        while imageTestTemp(row, col) > 0
            newRow = row + randContainer(i);
            newCol = col + randContainer(i+1);
            i = i+2;
            if newRow>0 && newRow<columns+1    &&    newCol>0 && newCol<rows+1
                outputImage(newRow, newCol) = outputImage(newRow, newCol) + 1;
            end
            imageTestTemp(row, col) = imageTestTemp(row, col) - 1;
        end
    end
end
toc