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

first hg version after svn

+function display(m)
+% DISSIMATRIX/DISPLAY display of a dissimilarity matrix
+disp(' ');
+d = get(m,'Data');
+nam = get(m,'Name');
+fp = get(m,'FramePos');
+cha = get(m,'Channels');
+t = get(m,'Title');
+tp = get(m,'TrackPos');
+tv = get(m,'TrackVal');
+pp = get(m,'PeakPos');
+g = get(m,'Graph');
+b = get(m,'Branch');
+for i = 1:length(d)
+    if iscell(d{i})
+        d{i} = d{i}{1};
+    end
+    l = size(d{1},3);     % Number of channels
+    il = (1-0.15)/l;
+    figure
+    for k = 1:l         % For each channel
+        if l>1
+            subplot('Position',[0.1 (k-1)*il+0.1 0.89 il-0.02])
+        end
+        fpi = cell2mat(fp{i});
+        if size(fpi,1) == 2
+            fpi = (fpi(1,:,:,:)+fpi(2,:,:,:))/2;
+        end
+        if strcmp(m.view,'l')
+            h = imagesc(fpi,fpi(1:size(d{i},1))-fpi(1),d{i}(:,:,k));
+        else
+            h = imagesc(fpi,fpi,d{i}(:,:,k));
+        end
+        if not(isempty(g))
+            hold on
+            pi = cell(1,size(g{i}{1},2));
+            for j = 1:size(g{i}{1},2)
+                pi{j} = sort(pp{i}{1}{j});
+            end
+            for j = 1:size(g{i}{1},2)
+                for h = 1:length(g{i}{1}{1,j,l})
+                    next = g{i}{1}{1,j,l}{h};
+                    if length(next)>1
+                        for n = 1:size(next,1)
+                            plot([fpi(j) fpi(next(n,1))],...
+                                 [fpi(pi{j}(h)) - fpi(1), ...
+                                  fpi(pi{next(n,1)}(next(n,2))) - fpi(1)],...
+                                 'k','LineWidth',1)
+                            %plot([fpi(j) fpi(next(n,1))],...
+                            %     [fpi(pi{j}(h)) - fpi(1), ...
+                            %      fpi(pi{next(n,1)}(next(n,2))) - fpi(1)],...
+                            %     'w+','MarkerSize',10)
+                            plot([fpi(j) fpi(next(n,1))],...
+                                 [fpi(pi{j}(h)) - fpi(1), ...
+                                  fpi(pi{next(n,1)}(next(n,2))) - fpi(1)],...
+                                 'kx','MarkerSize',10)
+                        end
+                    end
+                end
+            end
+            bi = b{i}{1}{1};
+            for j = 1:size(bi,1)-1
+                plot([fpi(bi(j,1)) fpi(bi(j+1,1))],...
+                     [fpi(pi{bi(j,1)}(bi(j,2))) - fpi(1), ...
+                      fpi(pi{bi(j+1,1)}(bi(j+1,2))) - fpi(1)],...
+                     'w','LineWidth',1.5)
+                plot([fpi(bi(j,1)) fpi(bi(j+1,1))],...
+                     [fpi(pi{bi(j,1)}(bi(j,2))) - fpi(1), ...
+                      fpi(pi{bi(j+1,1)}(bi(j+1,2))) - fpi(1)],...
+                     'w+','MarkerSize',15)
+                plot([fpi(bi(j,1)) fpi(bi(j+1,1))],...
+                     [fpi(pi{bi(j,1)}(bi(j,2))) - fpi(1), ...
+                      fpi(pi{bi(j+1,1)}(bi(j+1,2))) - fpi(1)],...
+                     'kx','MarkerSize',15)
+            end
+        elseif not(isempty(tp)) && not(isempty(tp{i}))
+            hold on
+            for k = 1:size(tp{i}{1}{1},1)
+                prej = 0;
+                for j = 1:size(tp{i}{1}{1},2)
+                    if tv{i}{1}{1}(k,j)
+                        if prej% && not(isempty(tp(k,j)))
+                            plot([fpi(prej) fpi(j)],...
+                                 [fpi(tp{i}{1}{1}(k,prej)) - fpi(1) ...
+                                  fpi(tp{i}{1}{1}(k,j)) - fpi(1)],...
+                                 'k','LineWidth',1)
+                            plot([fpi(prej) fpi(j)],...
+                                 [fpi(tp{i}{1}{1}(k,prej)) - fpi(1) ...
+                                  fpi(tp{i}{1}{1}(k,j)) - fpi(1)],...
+                                 'w+','MarkerSize',10)
+                            plot([fpi(prej) fpi(j)],...
+                                 [fpi(tp{i}{1}{1}(k,prej)) - fpi(1) ...
+                                  fpi(tp{i}{1}{1}(k,j)) - fpi(1)],...
+                                 'kx','MarkerSize',10)
+                        end
+                        prej = j;
+                    end
+                end
+            end
+        elseif not(isempty(pp)) && not(isempty(pp{i}))
+            hold on
+            for k = 1:length(pp{i}{1})
+                for j = 1:length(pp{i}{1}{k})
+                    plot(fpi(k),fpi(pp{i}{1}{k}(j)) - fpi(1), ...
+                                 'w+','MarkerSize',10)
+                end
+            end
+        end
+        set(gca,'YDir','normal')
+        if k == l
+            title(t)
+        end
+        if k == 1
+            xlabel('temporal location of frame centers (in s.)')
+        end
+        if k == ceil(l/2)
+            if strcmp(m.view,'h')
+                ylabel('relative distance between compared frames (in s.)')
+            elseif strcmp(m.view,'l')
+                ylabel('temporal lag (in s.)')
+            else
+                ylabel('temporal location of frame centers (in s.)')
+            end
+        end
+        if l > 1
+            pos = get(gca,'Position');
+            hfig = axes('Position',[pos(1)-.05 pos(2)+pos(4)/2 .01 .01],...
+                        'Visible','off');
+            text(0,0,num2str(cha{i}(k)),'FontSize',12,'Color','r')
+        end
+    end
+    fig = get(0,'CurrentFigure');
+    disp(['The ',t,' related to file ',nam{i},' is displayed in Figure ',num2str(fig),'.']);
+end
+disp(' ');
+drawnow