Warning, /verification/fizhi-cs-aqualev20/scripts/get_new_TUV.m is written in an unsupported language. File is not indexed.
view on githubraw file Latest commit 7ae8fb32 on 2006-04-03 20:55:29 UTC
7ae8fb32b5 Andr*0001 % load COS & SIN of rotation angle:
0002 fid=fopen('proj_cs32_2uEvN.bin','r','b');
0003 nx = 192;
0004 ny = 32;
0005 uvEN=fread(fid,nx*ny*2,'real*8'); fclose(fid);
0006 uvEN=reshape(uvEN,[prod(size(uvEN))/2 2]);
0007 AngleCS=uvEN(:,1);
0008 AngleSN=uvEN(:,2);
0009 CS = reshape(AngleCS,[nx ny]);
0010 SN = reshape(AngleSN,[nx ny]);
0011 for fi = 1:6
0012 is = (fi - 1)*32;
0013 nTUu(:,:,fi) = CS(is+[1:32],1:32);
0014 nTUv(:,:,fi) = -SN(is+[1:32],1:32);
0015 nTVu(:,:,fi) = SN(is+[1:32],1:32);
0016 nTVv(:,:,fi) = CS(is+[1:32],1:32);
0017 end
0018
0019 load ../set_001/TUV.mat
0020
0021 oTUu = TUu;
0022 oTUv = TUv;
0023 oTVu = TVu;
0024 oTVv = TVv;
0025
0026 clear TUu TUv TVu TVv
0027
0028 figure(1)
0029 subplot(1,1,1)
0030 for fi = 1:6
0031
0032 iif = (fi-1)*4;
0033 subplot(6,4,iif+1), surf(oTUu(:,:,fi)), view(2), shading flat
0034 subplot(6,4,iif+2), surf(oTUv(:,:,fi)), view(2), shading flat
0035 subplot(6,4,iif+3), surf(oTVu(:,:,fi)), view(2), shading flat
0036 subplot(6,4,iif+4), surf(oTVv(:,:,fi)), view(2), shading flat
0037
0038 end
0039
0040 figure(2)
0041 subplot(1,1,1)
0042 for fi = 1:6
0043
0044 iif = (fi-1)*4;
0045 subplot(6,4,iif+1), surf(nTUu(:,:,fi)), view(2), shading flat
0046 subplot(6,4,iif+2), surf(nTUv(:,:,fi)), view(2), shading flat
0047 subplot(6,4,iif+3), surf(nTVu(:,:,fi)), view(2), shading flat
0048 subplot(6,4,iif+4), surf(nTVv(:,:,fi)), view(2), shading flat
0049
0050 end
0051
0052 TUu = nTUu;
0053 TUv = nTUv;
0054 TVu = nTVu;
0055 TVv = nTVv;
0056
0057 save TUV_from_proj_cs32_2uEvN TUu TUv TVu TVv
0058
0059
0060