% function mov = loadFileYuv(fileName, width, height, idxFrame)
function [mov,imgRgb] = loadFileYuv(fileName, width, height, idxFrame)
% load RGB movie [0, 255] from YUV 4:2:0 file
fileId = fopen(fileName, 'r');
subSampleMat = [1, 1; 1, 1];
nrFrame = length(idxFrame);
for f = 1 : 1 : nrFrame
% search fileId position
sizeFrame = 1.5 * width * height;
fseek(fileId, (idxFrame(f) - 1) * sizeFrame, 'bof');
% read Y component
buf = fread(fileId, width * height, 'uchar');
imgYuv(:, :, 1) = reshape(buf, width, height).'; % reshape RESHAPE(X,M,N) returns the M-by-N matrix
%whose elements are taken columnwise from X.
%An error results if X does not have M*N elements
% read U component
buf = fread(fileId, width / 2 * height / 2, 'uchar');
imgYuv(:, :, 2) = kron(reshape(buf, width / 2, height / 2).', subSampleMat); % reshape and upsample
% read V component
buf = fread(fileId, width / 2 * height / 2, 'uchar');
imgYuv(:, :, 3) = kron(reshape(buf, width / 2, height / 2).', subSampleMat); % reshape and upsample
% normalize YUV values
% imgYuv = imgYuv / 255;
% convert YUV to RGB
imgRgb = reshape(convertYuvToRgb(reshape(imgYuv, height * width, 3)), height, width, 3);
% imgRgb = ycbcr2rgb(imgYuv);
%imwrite(imgRgb,'ActualBackground.bmp','bmp');
mov(f) = im2frame(imgRgb);
% mov(f).cdata = uint8(imgRgb);
% mov(f).colormap = [];
% imwrite(imgRgb,'ActualBackground.bmp','bmp');
%figure, imshow(imgRgb);
%name = 'ActualBackground.bmp';
%Image = imread(name, 'bmp');
%figure, imshow(Image);
end
fclose(fileId);
答案 0 :(得分:0)
不确定我是否对YUV文件有一些基本的误解,但是如果你像我下面那样编辑函数,你就会在一个名为imgYUV
的变量中为每一帧提供YUV组件。请注意,在此过程中可能会耗尽内存,您可能不想一次性加载电影的所有帧。
function imgYUV = loadFileYuv(fileName, width, height, idxFrame)
% load YUV data from YUV 4:2:0 file
fileId = fopen(fileName, 'r');
subSampleMat = [1, 1; 1, 1];
nrFrame = length(idxFrame);
%# preassign imgYUV. In case we can't keep everything in RAM,
%# it is better that the function crashes here, rather than after
%# having wasted time slowly filling up the memory.
%# Since the images are likely to be of class 'uint8',
%# you can save on a lot of memory by initializing
%# imgYUV as zeros(width/2,height/2,3,nrFrame,'uint8');
imgYUV = zeros(width/2 height/2, 3, nrFrame);
for f = 1 : 1 : nrFrame
%# search fileId position
sizeFrame = 1.5 * width * height;
fseek(fileId, (idxFrame(f) - 1) * sizeFrame, 'bof');
%# read Y component
buf = fread(fileId, width * height, 'uchar');
imgYuv(:, :, 1, f) = reshape(buf, width, height).';
%# read U component
buf = fread(fileId, width / 2 * height / 2, 'uchar');
imgYuv(:, :, 2, f) = kron(reshape(buf, width / 2, height / 2).', subSampleMat); % reshape and upsample
% read V component
buf = fread(fileId, width / 2 * height / 2, 'uchar');
imgYuv(:, :, 3, f) = kron(reshape(buf, width / 2, height / 2).', subSampleMat); % reshape and upsample
end
fclose(fileId);
答案 1 :(得分:0)
我只想改变功能定义:
function [mov,imgYuv] = loadFileYuv(fileName, width, height, idxFrame)
无需返回imgRgb,因为它与mov(numel(idxFrame)).cdata
相同。
指定两个输出变量,您将获得第二个变量中最后一帧的YUV分量。
我会更改函数以检查缓冲区是否为空以避免错误,如果您要求的帧数多于文件所包含的帧数。
...
status = fseek(fileId, (idxFrame(f) - 1) * sizeFrame, 'bof');
if status == -1
error('Cannot read frame %d',idxFrame(f));
end
% read Y component
buf = fread(fileId, width * height, 'uchar');
if isempty(buf)
error('Cannot read frame %d',idxFrame(f));
end
imgYuv(:, :, 1) = reshape(buf, width, height).'; % reshape
% read U component
buf = fread(fileId, width / 2 * height / 2, 'uchar');
if isempty(buf)
error('Cannot read frame %d',idxFrame(f));
end
imgYuv(:, :, 2) = kron(reshape(buf, width / 2, height / 2).', subSampleMat); % reshape and upsample
% read V component
buf = fread(fileId, width / 2 * height / 2, 'uchar');
if isempty(buf)
error('Cannot read frame %d',idxFrame(f));
end
...