Plot a heatmap from a matix

Question

Ahmed Madhun 2019-10-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/487459-plot-a-heatmap-from-a-matix

编辑： Ahmed Madhun 2019-10-26

I implemented the minutia heat map presented in this paper (page 7): https://arxiv.org/pdf/1909.09901.pdf

The result is that i have 6 matrixes and i am looking for a way to plot these as shown in the paper (same page), like this:

Basically, the low values is presented in dark, and it gets lighter when it increase.

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

darova 2019-10-26

在 MATLAB Online 中打开

Make range for color axis the same for each figure

min = % minimum of 6 matrix
max = % maximum of 6 matrix
caxis([min max])

Ahmed Madhun 2019-10-26

编辑：Ahmed Madhun 2019-10-26

在 MATLAB Online 中打开

Sorry didn't get you comment: how to apply that for my code ?

% Function to create the HeatMap matrix of a minutia set in "TestData" file
function HeatMapCreator()
    % Start plotting minutia 
    [X, Y, A] = GetMinData("TestData");
    % Minutia 1 has X1, Y1, and A1 => X and Y is the position and A is the angle between 0-359
    % Count Minutiea
    n = length(X);
    % set Width and Hight
    W = 600;
    H = 750;
    K = 6;
    
    % Create the matixes
    MM = cell(K, 1);
    M = zeros(H,W);
    
    % Define the gussian paramater
    GP = 2*2^2;
    
    % Go throw each pixel
    for k = 1:6
        for i = 1 : W
           for j = 1 : H
              Hijk = 0; 
              for t = 1 : n
                Xt = X(t);
                Yt = Y(t);
                At = A(t);
                %Calculate Cs
                ED = sqrt((i-Xt)^2+(j-Yt)^2);
                Cs = exp(-ED/GP);
                %Calculate Co
                DO = At - (2 * k * pi / K);
                if (DO < -pi) || (DO > pi)
                    DO = 2 * pi - DO;
                end                       
                Co = exp(-DO/GP);
                Hijk = Hijk + (Cs * Co);
              end
              M(i,j) = Hijk;
           end
        end
        MM{k} = M;
    end
    
    % Show the 6 figures
    for k = 1:6
        figure(k)
        heatmap(MM{k},'Colormap',gray,'GridVisible','off');
    end
end

Comment: I set the Gussian paramater to test weather it works or not.

请先登录，再进行评论。

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

Answer 1

darova 2019-10-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/487459-plot-a-heatmap-from-a-matix#answer_398313

在 MATLAB Online 中打开

You are using H for rows and W for columns

M1 = zeros(H,W);
% ...
for i = 1 : W
    for j = 1 : H
        % ...
        M1(i,j) = Hijk;        % looks like mistake

You can use cells to create 6 matrices automatically

MM = cell(6,1);
M = zeros(H,W);
for k = 1:6
    for i 
        for j
            for t
                % do stuff
            end
            % ...
            M(i,j) = %...
        end
    end
    MM{k} = M;
end

After you found max and min values (global) use loop to visualize

for k = 1:6
    figure(k)
    heatmap(MM{k});
    caxis([minx maxx])
    colormap gray
end

11 个评论
显示 9更早的评论隐藏 9更早的评论

Ahmed Madhun 2019-10-26

编辑：Ahmed Madhun 2019-10-26

在 MATLAB Online 中打开

Not really, no value is now displayed on the heatmap figures at all.

I change my code as you advised:

    minV = 0;
    maxV = 0;
    
    for k = 1:6
        for i = 1 : W
           for j = 1 : H
              Hijk = 0; 
              for t = 1 : n
                Xt = X(t);
                Yt = Y(t);
                At = A(t);
                %Calculate Cs
                ED = sqrt((i-Xt)^2+(j-Yt)^2);
                Cs = exp(-ED/GP);
                %Calculate Co
                DO = At - (2 * k * pi / K);
                if (DO < -pi) || (DO > pi)
                    DO = 2 * pi - DO;
                end                       
                Co = exp(-DO/GP);
                Hijk = Hijk + (Cs * Co);
              end
              M(i,j) = Hijk;
           end
        end
        MM{k} = M;
        
        % Get the maximum and minimum cell from all matrixes
        if minV > min(M(:))
            minV = min(M(:));
        end
        
        if maxV < max(M(:))
            maxV = max(M(:));
        end
    end
    
    for k = 1:6
        figure(k)
        heatmap(MM{k})
        caxis([minV maxV])
        colormap gray
    end

the result is:

However, if i use the old code to display the heatmap as:

heatmap(MM{1},'Colormap',gray,'GridVisible','off');

the result is:

But this is equal for all matrixes!, which is not correct as all matrixes are different!

darova 2019-10-26

编辑：darova 2019-10-26

Can you attach a screenshot or something? The file is too large, can't donwload .pdf (i have low internet speed)

Ahmed Madhun 2019-10-26

The paper is only 8 Mb. You can find it if you search on google for: "Learning a Fixed-Length Fingerprint Representation"

请先登录，再进行评论。

Plot a heatmap from a matix

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

采纳的回答

11 个评论
显示 9更早的评论隐藏 9更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

Community Treasure Hunt

Plot a heatmap from a matix

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

采纳的回答

11 个评论 显示 9更早的评论隐藏 9更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

Community Treasure Hunt

WeChat

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

11 个评论
显示 9更早的评论隐藏 9更早的评论