comparison toolboxes/MIRtoolbox1.3.2/MIRToolbox/displot.m @ 0:e9a9cd732c1e tip

first hg version after svn

+function res = displot(x,y,xlab,ylab,t,fp,pp,tp,tv,ch,multidata,pm,ap,rp,cl)
+% graphical display of any data (except mirscalar data) computed by MIRToolbox
+
+%opengl('OpenGLWobbleTesselatorBug',1)
+% Rendering complex patch object may cause segmentation violation and
+% return a tesselator error message in the stack trace. This command
+% enables a workaround of that bug that might work sometimes...
+
+if isempty(y)
+    res = 0;
+    return
+end
+manychannels = 0;
+y1 = y{1};
+if length(y) == 1 
+    y = y{1};
+    if length(x) == 1
+        x = x{1};
+    end
+end
+if isempty(y)
+    res = 0;
+    return
+end
+
+if isstruct(cl)
+    y = cl.centr(:,cl.index);
+end
+
+figure
+fp2 = uncell(fp);
+c = size(y,2);
+lx = size(y,1);
+l = size(y1,3);     % Number of channels
+il = (1-0.15)/l;
+
+if iscell(y)
+    if size(y{1},3) > 1 && size(y{1},1) == 1  
+        for j = 1:c
+            y{j} = reshape(y{j},size(y{j},2),size(y{j},3));
+            y{j} = y{j}';
+        end
+        lx = l;
+        l = 1;
+    end
+else
+    if l > 1 && lx == 1
+        y = reshape(y,size(y,2),size(y,3));
+        y = y';
+        x = repmat(ch,[1 size(y,2) size(y,3)]);
+        lx = l;
+        l = 1;
+    end
+end
+if not(iscell(y)) && l > 20 %&& size(y,3) == 1
+    manychannels = 1;
+    if lx == 1
+        y = reshape(y,[c l])';
+        lx = l;
+        l = 1;
+    else
+        y = reshape(y,[lx l])';
+        fp = reshape(x,[1 lx size(x,3)]);
+        fp = fp(:,:,1);
+        fp2 = fp;
+        c = l;
+        l = 1;
+        x = ch';
+    end
+end
+curve = (not(iscell(y)) && not(isequal(fp2,0)) && size(fp2,2)==1) || ...
+        ...%(iscell(y) && size(y{1}) == 1) || ...
+        c == 1 || (strcmp(xlab,'time (s)') && not(manychannels));
+if curve
+    for i = 1:l
+        if l>1
+            subplot(l,1,l-i+1,'align');
+        end
+        if not(iscell(y))
+            y = {y};
+            x = {x};
+        end
+        col = cell(length(y));
+        for h = 1:length(y)
+            col{h} = zeros(size(y{h},2),3,size(y{h},4));
+            for j = 1:size(y{h},2)
+                if ischar(x{h})
+                    xj = x{h};
+                else
+                    xj = x{h}(:,j);
+                end
+                yj = y{h}(:,j,i,:);
+                lk = size(yj,4);
+                for k = 1:lk
+                    if c > 1
+                        col{h}(j,:,k) = rand(1,3);
+                    else
+                        col{h}(j,:,k) = num2col(k);
+                    end
+                    yk = yj(:,:,:,k);
+                    if not(isempty(yk))
+                        if iscell(xj)
+                            lj = length(xj);
+                            abs = 0:lj-1;
+                        elseif ischar(xj)
+                            abs = 1:length(yk);
+                        else
+                            abs = xj;
+                        end
+                        if length(abs) < 20 && lk == 1 && c == 1
+                            bar(abs,yk);
+                        else
+                            plot(abs,yk,'Color',col{h}(j,:,k));
+                        end
+                        if iscell(xj)
+                            tick = 0:ceil(lj/14):lj-1;
+                            set(gca,'xtick',tick);
+                            set(gca,'xticklabel',xj(tick+1));
+                        elseif ischar(xj)
+                            set(gca,'xticklabel',x);
+                        end
+                        hold on
+                    end
+                end
+            end
+            if length(y) > 1
+                if isempty(cl)% || isempty(cl{1}) ... this because of a bug..
+                    %|| isempty(cl{i})
+                    colr = h;
+                %elseif iscell(cl{i})
+                %    colr = cl{i}(h);
+                elseif length(cl) == 1
+                    colr = cl;
+                else
+                    colr = cl(h); %cl{i}(h);
+                end
+                if not(isempty(x{h}))
+                    rectangle('Position',[x{h}(1),...
+                                          min(min(y{h}(:,:,i))),...
+                                          x{h}(end)-x{h}(1),...
+                                          max(max(y{h}(:,:,i)))-...
+                                          min(min(y{h}(:,:,i)))]+1e-16,...
+                        'EdgeColor',num2col(colr),'Curvature',.1,'LineWidth',1)
+                end
+            end
+        end
+        if not(isempty(multidata))
+            legend(multidata{:},'Location','Best')
+        end
+        for h = 1:length(y)
+            for j = 1:size(y{h},2)
+                xj = x{h}(:,j);
+                yj = y{h}(:,j,i,:);
+                lk = size(yj,4);
+                for k = 1:lk
+                    yk = yj(:,:,:,k);
+                    if not(isempty(yk))
+                        if size(y{h},2) > 1
+                            rectangle('Position',[xj(1),min(yk),...
+                                xj(end)-xj(1),max(yk)-min(yk)]+1e-16,...
+                                'EdgeColor',col{h}(j,:,k))
+                        end
+                        if not(isempty(pp)) && not(isempty(pp{1}))
+                            [ppj order] = sort(pp{h}{1,j,i});
+                            if not(isempty(pm)) && not(isempty(pm{1}))
+                                pmj = pm{h}{1,j,i}(order);
+                                [R C] = find(pmj==k);
+                                if iscell(xj)
+                                    plot(ppj(R,C)-1,yk(ppj(R,C)),'or') 
+                                else
+                                    plot(xj(ppj(R,C)),yk(ppj(R,C)),'or') 
+                                end
+                            else
+                                if lk > 1
+                                    [R C] = find(ppj(:,:,2)==k);
+                                    plot(xj(ppj(R,C)),yk(ppj(R,C)),'or')
+                                else
+                                    plot(xj(ppj),yk(ppj),'or')
+                                end
+                            end
+                            if not(isempty(ap)) && not(isempty(ap{1}))
+                                apj = ap{h}{1,j,i};
+                                plot(xj(apj),yk(apj),'dr') 
+                                for g = 1:length(apj)
+                                    line([xj(ppj(g)),xj(apj(g))],...
+                                         [yk(ppj(g)),yk(apj(g))],...
+                                         'Color','r')
+                                end
+                            end
+                            if not(isempty(rp)) && not(isempty(rp{1}))
+                                rpj = rp{h}{1,j,i};
+                                plot(xj(rpj),yk(rpj),'dr') 
+                                for g = 1:length(rpj)
+                                    line([xj(ppj(g)),xj(rpj(g))],...
+                                         [yk(ppj(g)),yk(rpj(g))],...
+                                         'Color','r')
+                                end
+                            end
+                        end
+                    end
+                end
+            end
+        end
+        if i == l
+            title(t)
+        end
+        if i == 1
+            xlabel(xlab)
+        end
+        if l > 1
+            %if iscell(x)
+            %    num = x{i}(1);
+            %else
+                num = ch(i);
+            %end
+            pos = get(gca,'Position');
+            axes('Position',[pos(1)-.05 pos(2)+pos(4)/2 .01 .01],'Visible','off');
+            text(0,0,num2str(num),'FontSize',12,'Color','r')
+        end
+    end
+else
+    % 2-dimensional image
+    displayseg = 0;
+    if iscell(x) && ischar(x{1})
+        if size(y,4) > 1 || size(y,1) > size(x,1)
+            ticky = 0.5:23.5;
+            tickylab = {'CM','C#M','DM','D#M','EM','FM','F#M',...
+            'GM','G#M','AM','A#M','BM','Cm','C#m','Dm','D#m','Em',...
+            'Fm','F#m','Gm','G#m','Am','A#m','Bm'};
+        else
+            ticky = (1:size(x,1))';
+            tickylab = x;
+            displayseg = 1;
+        end
+        x = (1:size(y,1))';
+    elseif iscell(x) && iscell(x{1}) && ischar(x{1}{1})
+        if size(y{1},4) > 1
+            ticky = 0.5:23.5;
+            tickylab = {'CM','C#M','DM','D#M','EM','FM','F#M',...
+            'GM','G#M','AM','A#M','BM','Cm','C#m','Dm','D#m','Em',...
+            'Fm','F#m','Gm','G#m','Am','A#m','Bm'};
+        else
+            ticky = (1:size(x{1},1))';
+            tickylab = x{1};
+            displayseg = 1;
+        end
+        for i = 1:length(x)
+            x{i} = (1:size(x{i},1))';
+        end
+    else
+        ticky = [];
+    end
+    if iscell(y)
+        displayseg = 1;
+        for i = 1:l
+            if l>1
+                subplot(l,1,l-i+1,'align');
+%                subplot('Position',[0.1 (i-1)*il+0.1 0.89 il-0.02])
+            end
+            hold on
+            %surfplot(segt,x{1},repmat(x{1}/x{1}(end)*.1,[1,length(segt)]));
+            if length(x)==1 && length(y)>1
+                for k = 2:length(y)
+                    x{k} = x{1};
+                end
+            end
+            for j = 1:length(x)
+                if length(x{j}) > 1
+                    if size(x{j},1) == 1 && size(x{j},3) > 1
+                        x{j} = reshape(x{j},size(x{j},2),size(x{j},3))';
+                        mel = 1;
+                    else
+                        mel = 0;
+                    end
+                    xx = zeros(size(x{j},1)*size(y{j},4),1); %,size(x{j},2));
+                    yy = zeros(size(y{j},1)*size(y{j},4),size(y{j},2));
+                    for k = 1:size(y{j},4)
+                        xx((k-1)*size(x{j},1)+1:k*size(x{j},1),1) = x{j}(:,1);
+                        yy((k-1)*size(y{j},1)+1:k*size(y{j},1),:) = y{j}(:,:,i,k);
+                    end
+                    if size(y{j},4) == 1 && not(mel)
+                        xxx = [1.5*xx(1,1)-0.5*xx(2,1);...
+                                    (xx(1:end-1)+xx(2:end))/2;...
+                                    1.5*xx(end,1)-0.5*xx(end-1,1)];
+                    else
+                        xxx = (0:size(yy,1))';
+                    end
+                    if size(fp{j},2) > 1
+                        segt = [fp{j}(1,1),...
+                                mean([fp{j}(2,1:end-1);fp{j}(1,2:end)]),...
+                                fp{j}(2,end)];
+                    else
+                        segt = fp{j}';
+                    end
+                    surfplot(segt,xxx,yy(:,:));
+                    if not(isempty(ticky))
+                        set(gca,'ytick',ticky);
+                        set(gca,'yticklabel',tickylab);
+                    end
+                    set(gca,'YDir','normal')
+                    if (exist('pp') == 1) && not(isempty(pp))
+                        if not(isempty(pp{j}))
+                            for k = 1:length(pp{j})
+                                ppj = pp{j}{k};
+                                if size(ppj,3) == 2
+                                    ppj(:,:,1) = ppj(:,:,1) + (ppj(:,:,2)-1)*size(x,1);
+                                    ppj(:,:,2) = [];
+                                    plot(mean(fp{j}(:,k)),ppj-.5,'+w')
+                                elseif 0 %(exist('pm') == 1) && not(isempty(pm))
+                                    pmj = pm{k}{1,1,i};
+                                    ppj(:,:) = ppj(:,:) + (pmj(:,:)-1)*size(x,1);
+                                    plot(mean(fp{j}(:,k)),ppj-.5,'+w') % fp shows segmentation points
+                                elseif not(isempty(ppj))
+                                    plot(mean(fp{j}(:,k)),xx(ppj),'+k')
+                                end
+                            end
+                        end
+                    end
+                end
+            end
+            if i == l
+                title(t)
+            end
+            if i == 1
+                if manychannels
+                    xlabel(xlab);
+                elseif displayseg
+                    xlabel('time axis (in s.)');
+                else
+                    xlabel('temporal location of beginning of frame (in s.)');
+                end
+            end
+            if l > 1
+                num = ch(i);
+                pos = get(gca,'Position');
+                axes('Position',[pos(1)-.05 pos(2)+pos(4)/2 .01 .01],'Visible','off');
+                text(0,0,num2str(num),'FontSize',12,'Color','r')
+            end
+        end
+        y = y{1};
+    else
+        mel = 0;
+        if iscell(x)
+            x = x{1};
+        end
+        if size(x,1) == 1 && size(y,1) > 1
+            if size(x,3) == 1
+                x = x';
+            else
+                mel = 1;
+            end
+        end
+        if not(isempty(pp)) && iscell(pp{1})
+            pp = pp{1};
+            pm = pm{1};
+        end
+        for i = 1:l
+            if l>1
+                subplot('Position',[0.1 (i-1)*il+0.1 0.89 il-0.02])
+            end
+            xx = zeros(size(x,1)*size(y,4),1); %,size(x,2));
+            yy = zeros(size(y,1)*size(y,4),size(y,2));
+            for k = 1:size(y,4)
+                xx((k-1)*size(x,1)+1:k*size(x,1),1) = x(:,1);
+                yy((k-1)*size(y,1)+1:k*size(y,1),:) = y(:,:,i,k);
+            end
+            if iscell(fp)
+                fp = uncell(fp);
+            end
+            lx = length(xx);
+            if lx < 6
+                curve = 1;
+                ttt = (fp(1,:)+fp(2,:))/2;
+                plot(ttt,yy)
+                if i == 1
+                    legdata = cell(1,lx);
+                    if strcmp(xlab(end),'s')
+                        xlab(end) = [];
+                    end
+                    for j = 1:lx
+                        legdata{j} = [xlab,' ',num2str(j)];
+                    end
+                    legend(legdata)
+                end
+            else
+                ttt = [fp(1,:) 2*fp(1,end)-fp(1,end-1)];
+                if size(y,4) == 1 && not(mel)
+                    xxx = [1.5*xx(1)-0.5*xx(2);...
+                           (xx(1:end-1)+xx(2:end))/2;...
+                           1.5*xx(end)-0.5*xx(end-1)];
+                else
+                    xxx = (0:size(yy,1))';
+                end
+                surfplot(ttt,xxx,yy);
+                if not(isempty(ticky))
+                    set(gca,'ytick',ticky);
+                    set(gca,'yticklabel',tickylab);
+                end
+                set(gca,'YDir','normal');
+            end
+            hold on
+            %if iscell(pp)
+            %    pp = uncell(pp);
+            %    if iscell(pm)
+            %        pm = uncell(pm);
+            %    end
+            %end
+            if (exist('tp') == 1) && not(isempty(tp))
+                tp = tp{i}{1};
+                tv = tv{i}{1};
+                for k = 1:size(tp,1)
+                    prej = 0;
+                    for j = 1:size(tp,2)
+                        if tv(k,j)
+                            if prej% && not(isempty(tp(k,j)))
+                                plot([(ttt(prej)+ttt(prej+1))/2,...
+                                      (ttt(j)+ttt(j+1))/2],...
+                                     [(xxx(tp(k,prej))+xxx(tp(k,prej)+1))/2,...
+                                      (xxx(tp(k,j))+xxx(tp(k,j)+1))/2],...
+                                     'k','LineWidth',1)
+                                plot([(ttt(prej)+ttt(prej+1))/2,...
+                                      (ttt(j)+ttt(j+1))/2],...
+                                     [(xxx(tp(k,prej))+xxx(tp(k,prej)+1))/2,...
+                                      (xxx(tp(k,j))+xxx(tp(k,j)+1))/2],...
+                                     'w+','MarkerSize',10)
+                                plot([(ttt(prej)+ttt(prej+1))/2,...
+                                      (ttt(j)+ttt(j+1))/2],...
+                                     [(xxx(tp(k,prej))+xxx(tp(k,prej)+1))/2,...
+                                      (xxx(tp(k,j))+xxx(tp(k,j)+1))/2],...
+                                     'kx','MarkerSize',10)
+                            end
+                            prej = j;
+                        end
+                    end
+                end
+            elseif (exist('pp') == 1) && not(isempty(pp))
+                if size(pp,3) == l
+                    for j = 1:size(pp,2)
+                        ppj = pp{1,j,i};
+                        if (exist('pm') == 1)
+                            pmj = pm{1,j,i};
+                        else
+                            pmj = [];
+                        end
+                        if iscell(ppj)
+                            ppj = uncell(ppj);
+                            if iscell(pmj)
+                                pmj = uncell(pmj);
+                            end
+                        end
+                        if size(ppj,3) == 2
+                            ppj(:,:,1) = ppj(:,:,1) + (ppj(:,:,2)-1)*size(x,1);
+                            ppj(:,:,2) = [];
+                            %plot((fp(1,j)+fp(2,j))/2,ppj-.5,'+w')
+                            plot((ttt(j)+ttt(j+1))/2,(xxx(ppj)+xxx(ppj+1))/2,'+w')
+                        else
+                            if not(isempty(pmj))
+                                ppj(:,:,:) = ppj(:,:,:) + (pmj(:,:,:)-1)*size(x,1);
+                            end
+                            if not(isempty(ppj))
+                                %plot((segt(j)+segt(j+1))/2,xx(ppj),'+k','MarkerSize',10)
+                                plot((ttt(j)+ttt(j+1))/2,(xxx(ppj)+xxx(ppj+1))/2,'+w','MarkerSize',10)
+                            end
+                        end
+                    end
+                else
+                    for j = 1:size(pp,3)
+                        ppj = pp{1,1,j};
+                        if (exist('pm') == 1)
+                            pmj = pm{1,1,j};
+                        else
+                            pmj = [];
+                        end
+                        if not(isempty(ppj))
+                            %plot((segt(j)+segt(j+1))/2,xx(ppj),'+k','MarkerSize',10)
+                            plot((ttt(ppj)+ttt(ppj+1))/2,(xxx(j)+xxx(j+1))/2,'+w','MarkerSize',10)
+                        end
+                    end
+                end
+            end
+            if i == l
+                title(t)
+            end
+            if i == 1
+                xlabel('time axis (in s.)');
+            end
+            if l > 1
+                if iscell(x)
+                    num = x{i}(1);
+                else
+                    num = ch(i);
+                end
+                pos = get(gca,'Position');
+                hfig = axes('Position',[pos(1)-.05 pos(2)+pos(4)/2 .01 .01],'Visible','off');
+                text(0,0,num2str(num),'FontSize',12,'Color','r')
+            end
+        end
+    end
+end
+if l == 1
+    if curve
+        %if (exist('dbv') == 1) && dbv
+        %    ylabel([ylab ' (logarithmic scale)']);
+        %else
+            ylabel(ylab);
+        %end
+    elseif manychannels
+        ylabel('Channels');
+    else
+        ylabel(xlab)
+    end
+end
+res = 1;
+drawnow