422转420 并 从YUV转RGB

fid=fopen('F:/QZG/YUV422toYUV420/Convert/vicap_chn0_1024_768_p422_256.yuv','rb');
outfid=fopen('F:/QZG/YUV422toYUV420/Convert/matlab_out.yuv','wb');

fseek(fid,1024*768*2*5,'bof');

for i=1:10
Y=fread(fid,[1024,768],'uint8');
U=fread(fid,[1024/2,768],'uint8');
V=fread(fid,[1024/2,768],'uint8');

UU=U(:,1:2:768);
VV=V(:,1:2:768);
fwrite(outfid,Y,'uint8');
fwrite(outfid,UU,'uint8');
fwrite(outfid,VV,'uint8');
end

fclose(fid);
fclose(outfid);

Y1=Y';
U1=UU';
V1=VV';


Y2=Y1(221:410,41:400);
U2=U1(111:205,21:200);
V2=V1(111:205,21:200);

figure(1);
imshow(uint8(Y2));title('Y分量');
U_old=U(21:200,221:410)';
V_old=V(21:200,221:410)';
figure(2)
subplot(121);
imshow(uint8(U_old));title('原始U分量');
subplot(122);
imshow(uint8(V_old));title('原始V分量');


figure(3);
subplot(121);
imshow(uint8(U2));title('420采样后U分量');
subplot(122);
imshow(uint8(V2));title('420采样后V分量');


for i=1:size(U2,1)*2
    for j=1:size(U2,2)*2
        U3(i,j)=U2(ceil(i/2),ceil(j/2));
        V3(i,j)=V2(ceil(i/2),ceil(j/2));
    end
end

Y3=Y2;
R=1.164.*(Y3-16)+1.596.*(U3-128);
G=1.164.*(Y3-16)-0.813.*(U3-128)-0.392.*(V3-128);
B=1.164.*(Y3-16)+2.017.*(V3-128);

RGB=uint8(ones(size(R,1),size(R,2),3));
RGB(:,:,1)=R;
RGB(:,:,2)=G;
RGB(:,:,3)=B;

figure(4);
imshow(uint8(RGB));title('RGB显示');
%pixval on;