1 |
% load COS & SIN of rotation angle: |
2 |
fid=fopen('proj_cs32_2uEvN.bin','r','b'); |
3 |
nx = 192; |
4 |
ny = 32; |
5 |
uvEN=fread(fid,nx*ny*2,'real*8'); fclose(fid); |
6 |
uvEN=reshape(uvEN,[prod(size(uvEN))/2 2]); |
7 |
AngleCS=uvEN(:,1); |
8 |
AngleSN=uvEN(:,2); |
9 |
CS = reshape(AngleCS,[nx ny]); |
10 |
SN = reshape(AngleSN,[nx ny]); |
11 |
for fi = 1:6 |
12 |
is = (fi - 1)*32; |
13 |
nTUu(:,:,fi) = CS(is+[1:32],1:32); |
14 |
nTUv(:,:,fi) = -SN(is+[1:32],1:32); |
15 |
nTVu(:,:,fi) = SN(is+[1:32],1:32); |
16 |
nTVv(:,:,fi) = CS(is+[1:32],1:32); |
17 |
end |
18 |
|
19 |
load ../set_001/TUV.mat |
20 |
|
21 |
oTUu = TUu; |
22 |
oTUv = TUv; |
23 |
oTVu = TVu; |
24 |
oTVv = TVv; |
25 |
|
26 |
clear TUu TUv TVu TVv |
27 |
|
28 |
figure(1) |
29 |
subplot(1,1,1) |
30 |
for fi = 1:6 |
31 |
|
32 |
iif = (fi-1)*4; |
33 |
subplot(6,4,iif+1), surf(oTUu(:,:,fi)), view(2), shading flat |
34 |
subplot(6,4,iif+2), surf(oTUv(:,:,fi)), view(2), shading flat |
35 |
subplot(6,4,iif+3), surf(oTVu(:,:,fi)), view(2), shading flat |
36 |
subplot(6,4,iif+4), surf(oTVv(:,:,fi)), view(2), shading flat |
37 |
|
38 |
end |
39 |
|
40 |
figure(2) |
41 |
subplot(1,1,1) |
42 |
for fi = 1:6 |
43 |
|
44 |
iif = (fi-1)*4; |
45 |
subplot(6,4,iif+1), surf(nTUu(:,:,fi)), view(2), shading flat |
46 |
subplot(6,4,iif+2), surf(nTUv(:,:,fi)), view(2), shading flat |
47 |
subplot(6,4,iif+3), surf(nTVu(:,:,fi)), view(2), shading flat |
48 |
subplot(6,4,iif+4), surf(nTVv(:,:,fi)), view(2), shading flat |
49 |
|
50 |
end |
51 |
|
52 |
TUu = nTUu; |
53 |
TUv = nTUv; |
54 |
TVu = nTVu; |
55 |
TVv = nTVv; |
56 |
|
57 |
save TUV_from_proj_cs32_2uEvN TUu TUv TVu TVv |
58 |
|
59 |
|
60 |
|