MNavigation/DicomLoader.cpp at master · alankzh/MNavigation · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
#include "DicomLoader.h"


#define TESTLOG(X) std::cout << #X  << ": " << X << std::endl

bool DicomLoader::has_resource_(false);

itk::GDCMSeriesFileNames::Pointer DicomLoader::file_parser_(itk::GDCMSeriesFileNames::New());

DicomLoader::ReaderType::Pointer DicomLoader::reader_(ReaderType::New());

itk::GDCMImageIO::Pointer DicomLoader::medicalinfo_parser_(itk::GDCMImageIO::New());

itk::ImageToVTKImageFilter<DicomLoader::ImageType>::Pointer DicomLoader::filter_(itk::ImageToVTKImageFilter<ImageType>::New());

DicomLoader::Resource DicomLoader::resource_;

DicomLoader::DicomLoader()
{
}

DicomLoader::~DicomLoader()
{
    ReleaseData();
}

DicomLoader::ErrorCode DicomLoader::ReadDicomData(std::string path,ProgressObserver::Pointer observer) {

    file_parser_->SetInputDirectory(path);
    auto file_names = file_parser_->GetInputFileNames();
    if (file_names.size() == 0) {
        return kNoFilesFound;
    }

    // ProgressObserver::Pointer observer=ProgressObserver::New();

    reader_->AddObserver(itk::ProgressEvent(),observer);

    reader_->SetFileNames(file_names);
    reader_->SetImageIO(medicalinfo_parser_);
    reader_->Update();

    //患者姓名
    medicalinfo_parser_->GetValueFromTag("0010|0010", resource_.info.patient_name);
    //患者id
    medicalinfo_parser_->GetValueFromTag("0010|0020", resource_.info.patient_ID);
    //患者出生日期
    medicalinfo_parser_->GetValueFromTag("0010|0030", resource_.info.patient_DOB);
    //患者性别
    medicalinfo_parser_->GetValueFromTag("0010|0040", resource_.info.patient_sex);
    //患者年龄
    medicalinfo_parser_->GetValueFromTag("0010|1010", resource_.info.patient_age);
    //检查id
    medicalinfo_parser_->GetValueFromTag("0020|0010", resource_.info.study_ID);
    //检查日期
    medicalinfo_parser_->GetValueFromTag("0008|0020", resource_.info.study_date);
    //检查时间
    medicalinfo_parser_->GetValueFromTag("0008|0030", resource_.info.study_time);
    //检查描述
    medicalinfo_parser_->GetValueFromTag("0008|1030", resource_.info.study_description);
    //检查模态
    medicalinfo_parser_->GetValueFromTag("0008|0060", resource_.info.modality);
    //制造商，然而未查询到此tag
    medicalinfo_parser_->GetValueFromTag("0008|0070", resource_.info.manufacturer);
    //序列实例号
    medicalinfo_parser_->GetValueFromTag("0020|000E", resource_.info.series_instance_UID);
    //检查描述和说明
    medicalinfo_parser_->GetValueFromTag("0008|103E", resource_.info.series_description);
    //看名字描述是建立机构名 未查询到此tag
    medicalinfo_parser_->GetValueFromTag("0008|0080", resource_.info.institution_name);
    //图像的左上角在控件坐标系中xyz坐标，单位是毫米
    medicalinfo_parser_->GetValueFromTag("0020|0032", resource_.info.image_position);
    //层厚
    medicalinfo_parser_->GetValueFromTag("0018|0050", resource_.info.slice_thickness);
    //层间距
    medicalinfo_parser_->GetValueFromTag("0018|0088", resource_.info.spacing_between_slices);


    TESTLOG(resource_.info.patient_name);
    TESTLOG(resource_.info.patient_ID);
    TESTLOG(resource_.info.patient_DOB);
    TESTLOG(resource_.info.patient_sex);
    TESTLOG(resource_.info.patient_age);
    TESTLOG(resource_.info.study_ID);
    TESTLOG(resource_.info.study_date);
    TESTLOG(resource_.info.study_time);
    TESTLOG(resource_.info.study_description);
    TESTLOG(resource_.info.modality);
    TESTLOG(resource_.info.manufacturer);
    TESTLOG(resource_.info.series_instance_UID);
    TESTLOG(resource_.info.series_description);
    TESTLOG(resource_.info.institution_name);
    TESTLOG(resource_.info.image_position);
    TESTLOG(resource_.info.slice_thickness);
    TESTLOG(resource_.info.spacing_between_slices);

    filter_->SetInput(reader_->GetOutput());
    filter_->Update();

    has_resource_ = true;
    resource_.data = filter_->GetOutput();


    int image_numbers = resource_.data->GetDimensions()[2];
    if (image_numbers > 20) {
        return kDicom;
    }
    else {
        return kImage;
    }
}

void DicomLoader::ReleaseData() {

    file_parser_ = itk::GDCMSeriesFileNames::New();
    reader_ = ReaderType::New();
    medicalinfo_parser_ = itk::GDCMImageIO::New();
    filter_ = itk::ImageToVTKImageFilter<ImageType>::New();
    has_resource_ = false;
}

vtkImageData* DicomLoader::GetData() {
    if (HasData()) {
        return resource_.data;
    }
    else {
        return nullptr;
    }
}

const DicomLoader::ResourceInfo* DicomLoader::GetDataInfo() {
    if (HasData()) {
        return &resource_.info;
    }
    else {
        return nullptr;
    }
}

bool DicomLoader::HasData() {
    return has_resource_;
}