将avi视频转换为多帧的dicom图

网络上有很多将多帧dicom图转换为avi的小工具.
但是反转换的工具却没有找到.
可能是因为需要自己写meta数据.
正好要用到, 所以用matlab实现了这个功能.

注意:
1. 要有和目标输出dicom图相同元数据的参考dicom图.
2. 如果没有参考dicom, 需要自己写相应的元数据.
3. 多帧dicom图的数据是按照4-D存储的, 格式为(width, height, 1, frame_num).
4. matlab读avi数据是uint8, dicom图数据是uint16. 需要一个扩展.
5. 这种生成的dicom图比较简单, 不能用于医学诊断.

matlab代码:

%% 把avi转录成dcm
clear;close all;clc;
metadata=dicominfo('ref.dcm'); % 参考的dicom文件的元数据
refmov=dicomread('ref.dcm');% 参考dicom的图像文件
xyloObj = VideoReader('录像7.avi');

% 帧数如果和ref.dcm不一样, 需要更改元数据的参数
nFrameOut=20; 
if nFrameOut>xyloObj.NumberOfFrames
    nFrameOut=xyloObj.NumberOfFrames;
end

% 多帧的dicom图是按照4-D矩阵存储的.
mov=nan(960,960,1,nFrameOut);
% 一次读一帧
for k = 1 : nFrameOut
    framedata=imresize(rgb2gray(read(xyloObj, k)),[960,960]);
    % uint8扩展为uint16
    mov(:,:,1,k) = uint16(uint16(framedata)*16);
end

% 写成多帧dicom,'CreateMode'设置为'copy'
metadata.NumberOfFrames=nFrameOut;
dicomwrite(uint16(mov),'录像7.dcm',metadata,'CreateMode','copy','MultiframeSingleFile',true);