source: issm/trunk-jpl/externalpackages/pcatool/install/dispSVDnb.m@ 21589

Last change on this file since 21589 was 21589, checked in by Eric.Larour, 8 years ago

CHG: new pcatool for EOF analysis.

File size: 1.7 KB
Line 
1% dispSVDnb(typg,CHP1,CHP2,N,expvar) Display SVDs
2%
3% => Display SVDs in white and black
4% typg is plot type you want (see plotm.m).
5% CHP contains SVDs and is of the form: CHP(iSVD,X,Y)
6% N is number of SVD to display
7% expvar contained the explained variance of SVDs
8%
9% Rq: it''s using the plotmNB function
10%
11%================================================================
12
13% Guillaume MAZE - LPO/LMD - December 2004
14% gmaze@univ-brest.fr
15
16function [] = dispSVDnb(typg,CHP1,CHP2,N,expvar)
17
18global ffile1 ffile2
19% --------------------------------------------------
20% AFFICHAGE
21% --------------------------------------------------
22if (typg~=1)&(typg~=7)&(typg~=3)&(typg~=6)
23 CBAROR='horiz';
24 iw=N;jw=2;
25else
26 CBAROR='vert';
27 iw=N;jw=2;
28end
29
30
31f1=figur;clf; hold on
32for iN=1:N
33 subplot(iw,jw,2*iN-1)
34C = squeeze(real(CHP1(iN,:,:)));
35%[cs,h]=plotmNB(C,typg); whos cs h
36[cs,h]=plotmNB(C,typg);
37%clabel(cs,h,'labelspacing',100);
38titre1=strcat('SVD',num2str(iN),':',ffile1);
39titre2=strcat(' (',num2str(expvar(iN)),'%)');
40title(strcat(titre1,titre2));
41
42 subplot(iw,jw,2*iN)
43C = squeeze(real(CHP2(iN,:,:)));
44[cs,h]=plotmNB(C,typg);
45%clabel(cs,h,'labelspacing',100);
46titre1=strcat('SVD',num2str(iN),':',ffile2);
47titre2=strcat(' (',num2str(expvar(iN)),'%)');
48title(strcat(titre1,titre2));
49
50end %for iN
51
52
53
54
55%titre1=strcat(num2str(lat(il1)),'/',num2str(lat(il2)));
56titre2=strcat(' Total explained variance:',num2str(sum(expvar)),'%');
57%titre3=strcat('(Analysed files: ',ffile1,' vs ',ffile2,')');
58%titre=strvcat(titre2,titre3);
59suptitle(titre2);
60
61set(f1,'Position',[378 39 313 647]); % Laptop screen
62%set(f1,'Position',[369 55 316 899]); % Desktop screen
63%set(f1,'Name',strcat(ffile,'<>',titre1,'<> REAL PART'));
Note: See TracBrowser for help on using the repository browser.