DCMTK BMP转多帧影像

BOOL CDcmAssit::MultibmpConvertToDCM(const char * filename , const char * outputfilename )
{
    OFCondition status;

    DcmFileFormat fileformat;
    DcmDataset* mydatasete=fileformat.getDataset();
    DicomUtils::AddDicomElements((DcmDataset*&)mydatasete);
    Uint16 rows,cols,samplePerPixel,bitsAlloc,bitsStored,highBit,pixelRpr,planConf,pixAspectH,pixAspectV;
    OFString photoMetrInt;
    Uint32 length;
    E_TransferSyntax ts;
    char* mydata=new char[1024*1024*10];
    memset(mydata,0,sizeof(char)*1024*1024*10);
    char* tmpData=mydata;

    //拆分文件名
    OFList<OFString> inputFiles;
    if(filename != NULL)
    {
        char * bufFileNames = new char[strlen(filename)+1];
        strcpy_s(bufFileNames,strlen(filename)+1,filename);
        char * p1,*p2;
        p1 = p2= bufFileNames;
        while(1)
        {
            p1 = strchr(p2,'|');
            if(p1 !=NULL)
            {
                *p1='\0';
                inputFiles.push_back(p2);
                p2 = ++p1;
            }
            else
            {
                inputFiles.push_back(p2);
                break;
            }
        }
        delete[] bufFileNames;
    }

    OFListIterator(OFString) if_iter = inputFiles.begin();
    OFListIterator(OFString) if_last = inputFiles.end();
     OFBool ignoreName;
     const char * currentfilename = NULL;
    while (if_iter != if_last)
    {
        ignoreName = OFFalse;
        currentfilename = (*if_iter).c_str();
        if (OFStandard::fileExists(currentfilename))
        {
            OFString offilename = (*if_iter).c_str();
            I2DBmpSource* bmpSource=new I2DBmpSource();
            bmpSource->setImageFile(offilename);
            char *pixData = NULL;
            bmpSource->readPixelData(rows,cols,samplePerPixel,photoMetrInt,bitsAlloc,bitsStored,highBit,pixelRpr,planConf,pixAspectH,pixAspectV,pixData,length,ts);
            memcpy(tmpData,pixData,length);
            tmpData+=length;
            delete bmpSource;

        }
        ++if_iter;
    }


    mydatasete->putAndInsertUint16(DCM_SamplesPerPixel,samplePerPixel);
    mydatasete->putAndInsertString(DCM_NumberOfFrames,"4");
    mydatasete->putAndInsertUint16(DCM_Rows,rows);
    mydatasete->putAndInsertUint16(DCM_Columns,cols);
    mydatasete->putAndInsertUint16(DCM_BitsAllocated,bitsAlloc);
    mydatasete->putAndInsertUint16(DCM_BitsStored,bitsStored);
    mydatasete->putAndInsertUint16(DCM_HighBit,highBit);
    mydatasete->putAndInsertUint8Array(DCM_PixelData,(Uint8*)mydata,4*length);
    mydatasete->putAndInsertOFStringArray(DCM_PhotometricInterpretation,photoMetrInt);
    mydatasete->putAndInsertUint16(DCM_PlanarConfiguration,0);
    status=fileformat.saveFile(outputfilename,ts);
    if(status.bad())
    {
        return FALSE; 
    }
    return TRUE;

}

void DicomUtils::AddDicomElements(DcmDataset*& dataset)  
{  
    //构建测试数据  
  
    /*  添加患者信息  */  
    dataset->putAndInsertUint16(DCM_AccessionNumber,0);  
    dataset->putAndInsertString(DCM_PatientName,"zssure",true);  
    dataset->putAndInsertString(DCM_PatientID,"2234");  
    dataset->putAndInsertString(DCM_PatientBirthDate,"20141221");  
    dataset->putAndInsertString(DCM_PatientSex,"M");  
  
    /*  添加Study信息   */  
    dataset->putAndInsertString(DCM_StudyDate,"20141221");  
    dataset->putAndInsertString(DCM_StudyTime,"195411");  
    char uid[100];  
    dcmGenerateUniqueIdentifier(uid,SITE_STUDY_UID_ROOT);  
    dataset->putAndInsertString(DCM_StudyInstanceUID,uid);  
    dataset->putAndInsertString(DCM_StudyID,"1111");  
  
  
    /*  添加Series信息  */  
    dataset->putAndInsertString(DCM_SeriesDate,"20141221");  
    dataset->putAndInsertString(DCM_SeriesTime,"195411");  
    memset(uid,0,sizeof(char)*100);  
    dcmGenerateUniqueIdentifier(uid,SITE_SERIES_UID_ROOT);  
    dataset->putAndInsertString(DCM_SeriesInstanceUID,uid);  
    /*  添加Image信息   */  
    dataset->putAndInsertString(DCM_ImageType,"ORIGINAL\PRIMARY\AXIAL");  
    dataset->putAndInsertString(DCM_ContentDate,"20141221");  
    dataset->putAndInsertString(DCM_ContentTime,"200700");  
    dataset->putAndInsertString(DCM_InstanceNumber,"1");  
    dataset->putAndInsertString(DCM_SamplesPerPixel,"1");  
    dataset->putAndInsertString(DCM_PhotometricInterpretation,"MONOCHROME2");  
    dataset->putAndInsertString(DCM_PixelSpacing,"0.3\0.3");  
    dataset->putAndInsertString(DCM_BitsAllocated,"16");  
    dataset->putAndInsertString(DCM_BitsStored,"16");  
    dataset->putAndInsertString(DCM_HighBit,"15");  
    dataset->putAndInsertString(DCM_WindowCenter,"600");  
    dataset->putAndInsertString(DCM_WindowWidth,"800");  
    dataset->putAndInsertString(DCM_RescaleIntercept,"0");  
    dataset->putAndInsertString(DCM_RescaleSlope,"1");  
  
  
}