一、DCMTK編譯（-mD/-mDd模式）

0. 材料：

• DCMTK3.6.2原始碼
• CMake
• VS2015

1. DCMTK原始碼下載

2. CMake下載

3. 編譯

3.1定位原始碼

開啟CMake，選擇DCMTK原始碼所在目錄，編譯結果存放目錄，如下如所示：開啟CMake，選擇DCMTK原始碼所在目錄，編譯結果存放目錄，如下如所示：

3.2 確定編譯器

選擇正確的原始碼目錄和編譯結果存放目錄後，點選【Configuire】，選擇所需編譯器，這裡編譯Win32位，所以選擇Visual Studio 14 2015編譯器，如果需要編譯x64，選擇Visual Studio 14 2015 Win64編譯器，然後點選【確定】.

3.3 修改CMake配置

這裡需要修改許多處，不然會編譯出錯，修改處見下圖：

3.5 生成工程

上步完成後，再次點選【Configure】,在點選【Generate】生成工程，點選【Open Project】,用VS2015開啟構建的工程。

3.6 確認工程構建正確

隨便選擇一個解決方案：
1. 確認是使用字符集為：多位元組字串集。

1. 確認是生成動態連結庫形式。
   
   

確定無誤後，執行編譯。

3.7 編譯

選擇ALL BUILD解決方案，右鍵【生成】開始編譯。

3.8 安裝

3.7步無誤後，選擇INSTALL解決方案，右鍵【生成】，開始安裝，則在安裝目錄下將生成編譯的結果。

3.9 Debug和Release版本各生成一份

通過切換VS2015的解決方案配置，Debug版本和Release版本執行3.7和3.8各一次,則能生成Debug和Release版本各生成一份。

3.10 結果

結果如下圖所示

二、VS2015配置DCMTK，OpenCV

0. 材料

• 編譯好的DCMTK動態連結庫
• OpenCV(這裡是2.4.13版本)

1.配置OpenCV

2. 配置DCMTK

跟配置OpenCV一模一樣，此處略過。

三、*.dcm轉OpenCV格式

程式碼及解釋如下：

#include "dcmtk/config/osconfig.h"
 
    /* make sure OS specific configuration is included first */
#include "dcmtk/ofstd/ofstdinc.h"

#include "dcmtk/dcmdata/dctk.h"
#include "dcmtk/dcmdata/cmdlnarg.h"
#include "dcmtk/ofstd/ofconapp.h"
#include "dcmtk/dcmdata/dcuid.h"       /* for dcmtk version name */
#include "dcmtk/dcmjpeg/djdecode.h"    /* for dcmjpeg decoders */
#include "dcmtk/dcmjpeg/dipijpeg.h"    /* for dcmimage JPEG plugin */

#pragma comment(lib, "ws2_32.lib")      // 包含socket，因為.dcm影象包含網路通訊協議
#pragma comment(lib,"iphlpapi.lib")
#pragma comment(lib, "netapi32.lib")

/*
* dcm2Mat：將.dcm格式影象轉換為OpenCV的Mat格式
* fileName：CT影象的名字
* mat: 儲存轉換後的CT影象
*返回值：成功返回true，失敗返回false
*/

bool dcm2Mat(const char* const fileName, cv::Mat & mat)
{

    DJDecoderRegistration::registerCodecs(); // register JPEG codecs  
    DcmFileFormat fileformat;                                        //DcmFileFormat fileformat;

    if (!(fileformat.loadFile(fileName).good()))
        return false;

    DcmDataset *dataset = fileformat.getDataset();
    // decompress data set if compressed    
    dataset->chooseRepresentation(EXS_LittleEndianExplicit, NULL);
    DcmElement* element = NULL;
    OFCondition result = dataset->findAndGetElement(DCM_PixelData, element);
    if (result.bad() || element == NULL)
        return false;
    unsigned short* pixData;
    result = element->getUint16Array(pixData);


    OFintptr_t rows = 0;                    // 影象行數
    OFintptr_t cols = 0;                    // 影象列數
    dataset->findAndGetLongInt(DCM_Rows, rows);
    dataset->findAndGetLongInt(DCM_Columns, cols);

    if (rows == 0 || cols == 0)
        return false;

    unsigned short *pbuf = new unsigned short[rows*cols];

    memcpy(pbuf, pixData, rows*cols * 2);
    Mat xxx(rows, cols, CV_16SC1, pbuf);
    Mat xx;
    cv::normalize(xxx, xx, 0, 255, NORM_MINMAX, CV_8UC1);  // 畫素值歸一化到0~255

    Mat afull(xx.rows, xx.cols, CV_8UC1, 255);  // 全白影象
    xx = afull - xx;                            // 影象顏色反轉，黑變白，白變黑
    mat = xx.clone();                           // 為什麼不直接用xx？ 因為為了防止記憶體洩漏

    cvtColor(mat, mat, CV_GRAY2RGB);            // 灰度轉彩色，視具體工程而定，如果只需灰度圖，可以不轉換

    delete[]pbuf;                               // 防止記憶體洩漏
    DJDecoderRegistration::cleanup(); // deregister JPEG codecs 
    return true;
}

四、*.dcm檔案頭資訊讀取

// 功能：讀取病人姓名
    OFFilename filename= "ct.1";//這裡改為你的CT影象檔案
    OFCondition status = fileformat.loadFile(filename);
    DcmDataset *dataset = fileformat.getDataset();      // 通過dataset 同時獲取檔案頭資訊和影象資料

    OFString PatientName;
    dataset->findAndGetOFString(DCM_PatientName, PatientName);   // DCM_PatientName:病患姓名
    cout <<"姓名:"<<PatientName<<endl;

其他資訊的關鍵字都定義在標頭檔案“dcdeftag.h”中，如下圖所示：

不清楚CT片的檔案頭到底包含什麼資訊？？ 使用Sante DICOM Viewer 8 檢視CT影象，從【view】——【file Header】中檢視。

五、可能會遇到的程式設計錯誤

• 出現錯誤的提示的函式以WSA開頭，如：LNK2001 無法解析的外部符號 [email protected] CC E:\visualstudio\CC\CC\oflog.lib(oflog.obj),如圖：
  

解決方案：顯示連結ws2_32.lib庫

#pragma comment(lib, "ws2_32.lib")

• 出現錯誤：LNK2001 無法解析的外部符號 [email protected] CC E:\visualstudio\CC\CC\ofstd.lib(ofstd.obj)
  

解決方案：顯示連結iphlpapi.lib庫

#pragma comment(lib,"iphlpapi.lib")

• 出現錯誤：LNK2001 無法解析的外部符號 [email protected] CC E:\visualstudio\CC\CC\dcmdata.lib(dcuid.obj)
  

解決方案：顯示連結netapi32.lib庫

#pragma comment(lib, "netapi32.lib")

六、總結

本文內容包含以下幾個方面：
1. 如何編譯DCMTK；
2. 如何進行*.dcm影象資料的讀取，並轉換為OpenCV的Mat格式；
3. *.dcm檔案的檔案頭資訊分析；
4. 期間遇到的常見錯誤列舉。

最後，祝你在醫學影象處理的道路上越走越好！