How to plot many images in a loop

Question

Tsitsos Dimitris 2021-3-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/775852-how-to-plot-many-images-in-a-loop

评论： Rik 2021-3-18

Hello,

I am working with a Matlab code where I want to run a loop. In each loop, I will do some processing and then show an image based on the processed data. When I go to the next loop, I want the previously stored image to be closed and generate a new image based on the latest data.

A sum up from my code is this :

for ww= 0.4:0.2:2.6
    **********************image processing code*********************************
                  t1=1-ww*(dark_channel1/A); % part that connect ww with the image processing code
    **********************image processing code*********************************
out = pic;
figure;
imshow (out)
imwrite(out,sprintf('%d.jpg',ww));
End;

What I expect is many different images depend on the "ww", but the output that I take looks like in each iteration each image sticks to the previous one.

I think that somehow I can’t separate the output (image), and after each iteration I don't reset my loop. As a result I display a new image but stack on the previous. It is like in each iteration the new image sticks to the previous.

Any suggestions??

Thank you in advance

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

请先登录，再进行评论。

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

Answer 1

ANKUR KUMAR 2021-3-17

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/775852-how-to-plot-many-images-in-a-loop#answer_650602

在 MATLAB Online 中打开

for i=1:5
    imwrite(randi(255,250,205)/255,sprintf('Image_%d.png',i))
end

10 个评论
显示 8更早的评论隐藏 8更早的评论

Tsitsos Dimitris 2021-3-17

编辑：Rik 2021-3-18

在 MATLAB Online 中打开

sure!

for ww= 0.4:0.2:2.6
    
    for itr=1:5
        under_50=0;
        for i=1:h
            for j=1:w
                dark_I(i,j)=min(img(i,j,:));
            end
        end
        for i=1:h
            for j=1:w
                if(dark_I(i,j)<50)
                    under_50=under_50+1;
                end
            end
        end
        total=size(img,1)*size(img,2)*size(img,3);
        percent=under_50/total;
        if(percent<m)
            minV=255;
            maxV=0;
            countlabel=0;
            img=double(img);
            
            
            dark_I1=dark_I;
            kh=floor(h/degree);
            kw=floor(w/degree);
            lh=mod(h,degree);
            lw=mod(w,degree);
            
            for op1=1:degree:(kh-1)*degree+1
                for op2=1:degree:(kw-1)*degree+1
                    darkI1_1=double(dark_I1(op1:op1+degree-1,op2:op2+degree-1));
                    label=min(min(darkI1_1));
                    for oop1=1:degree
                        for oop2=1:degree
                            darkI1_1(oop1,oop2)=label;
                        end
                    end
                    countlabel=countlabel+label;
                    if label<minV
                        minV=label;
                    end
                    if label>maxV
                        maxV=label;
                    end
                    %splice
                    if(op2==1)
                        I_linkcol=darkI1_1;
                    else
                        I_linkcol=cat(2,I_linkcol,darkI1_1);
                    end
                end
                if(op1==1)
                    I_linkrow=I_linkcol;
                else
                    I_linkrow=cat(1,I_linkrow,I_linkcol);
                end
            end
            I_row=I_linkrow;
            % splice
            if(lh==0&&lw==0)
                I_row=uint8(I_row);
            end
            if(lh==0&&lw~=0)
                rightlink=double(dark_I(1:h,w-lw+1:w));
                I_row=cat(2,I_row,rightlink);
                I_row=uint8(I_row);
            end
            if(lh~=0&&lw==0)
                downlink=double(dark_I(h-lh+1:h,1:w));
                I_row=cat(1,I_row,downlink);
                I_row=uint8(I_row);
            end
            if(lh~=0&&lw~=0)
                rightlink=double(dark_I(1:h-lh,w-lw+1:w));
                downlink=double(dark_I(h-lh+1:h,1:w));
                I_row=cat(2,I_row,rightlink);
                I_row=cat(1,I_row,downlink);
                I_row=uint8(I_row);
            end
            % =======================================
            dark_channel=double(I_row);
            for opp1=1:degree:(kh-1)*degree+1
                for opp2=1:degree:(kw-1)*degree+1
                    dark_channel1=double(dark_channel(opp1:opp1+degree-1,opp2:opp2+degree-1));
                    %                 numa=klabels(opp1,opp2);
                    %                 A=Anum(numa,1);
                    
                    A=150;
                    t1=1-ww*(dark_channel1/A);
                    %splice
                    if(opp2==1)
                        Tlink=t1;
                    else
                        Tlink=cat(2,Tlink,t1);
                    end
                end
                if(opp1==1)
                    t=Tlink;
                else
                    t=cat(1,t,Tlink);
                end
            end
            %optimize t
            for opp11=1:degree:(kh-1)*degree+1
                for opp22=1:degree:(kw-1)*degree+1
                    tt=double(t(opp11:opp11+degree-1,opp22:opp22+degree-1));
                    dark_II=double(dark_I(opp11:opp11+degree-1,opp22:opp22+degree-1));
                    I_rowR=double(I_row(opp11:opp11+degree-1,opp22:opp22+degree-1));
                    T=tt-0.001*(dark_II-I_rowR).*tt;
                    
                    if(opp22==1)
                        T1=T;
                    else
                        T1=cat(2,T1,T);
                    end
                end
                if(opp11==1)
                    ttt=T1;
                else
                    ttt=cat(1,ttt,T1);
                end
            end
            
            LastT=max(ttt,t0);
            [h1,w1]=size(LastT);
            
            new1=img(:,:,1);
            new2=img(:,:,2);
            new3=img(:,:,3);
            I_1(:,:,1)=double(new1(1:h1,1:w1));
            I_1(:,:,2)=double(new2(1:h1,1:w1));
            I_1(:,:,3)=double(new3(1:h1,1:w1));
            
            
            J1(:,:,1) = uint8((I_1(:,:,1) - (1-LastT)*A)./LastT);
            J1(:,:,2) = uint8((I_1(:,:,2) - (1-LastT)*A)./LastT);
            J1(:,:,3) = uint8((I_1(:,:,3) - (1-LastT)*A)./LastT);
            
            if(lh==0&&lw==0)
                J=J1;
            end
            if(lh==0&&lw~=0)
                right=double(img(1:h,w-lw+1:w,:));
                J=cat(2,J1,right);
                J=uint8(J);
            end
            if(lh~=0&&lw==0)
                down=double(img(h-lh+1:h,1:w,:));
                J=cat(1,J1,down);
                J=uint8(J);
            end
            if(lh~=0&&lw~=0)
                right=double(img(1:h-lh,w-lw+1:w,:));
                down=double(img(h-lh+1:h,1:w,:));
                I_row=cat(2,J1,right);
                J=cat(1,I_row,down);
                J=uint8(J);
            end
        else
            J=img;
        end
        pic=J;
        img=pic;
    end
    out = pic;
    figure;
    imshow (out)
    imwrite(out,sprintf('%d.jpg',ww));
end
Unrecognized function or variable 'h'.

I think that THE MAIN problem is that the loop does not reset in each iteration and the photos look like they are stick to the previous one.

For any questions please let me know!

Tsitsos Dimitris 2021-3-18

编辑：Tsitsos Dimitris 2021-3-18

I undestand!

Give me some time and I will post my updated code here!

Thank you again

Rik 2021-3-18

在 MATLAB Online 中打开

Regarding your email:

There must be a way to reproduce this problem. Try removing code line by line to see what happens when.

As an example:

for i=1:h
    for j=1:w
        dark_I(i,j)=min(img(i,j,:));
    end
end

This sets some values to dark_I. Assuming img contains integers 0-255, you can replace this loop with this:

dark_I(1:h,1:w)=uint8(randi(255,h,w));

It is also possible to create a mat file containing dark_I wrapped in a cell with the values for the first 3 iterations. That should also be enough to replicate the result.

I don't do private consulting, there are other people you can hire for that. You need to find a way to reproduce your issue in fewer lines of code than you posted above. Unless there is a bug caused by a very specific condition (making it hard/impossible to reproduce on my end anyway), it is almost always possible to reduce the cause of such behavior down to code that fits on one screen.

请先登录，再进行评论。