Homec4science
Diffusion ShoulderDB_ImportSegmentMeasureAnalyze (Bashardoust_SystCoord) Browse

R8218/CT/dicomCopier_gui.m
05553cd3e2fcBashardoust_SystCoord

Last Change
Raw File

dicomCopier_gui.m
View Options

function varargout = dicomCopier_gui(varargin)
% DICOMCOPIER_GUI MATLAB code for dicomCopier_gui.fig
% DICOMCOPIER_GUI, by itself, creates a new DICOMCOPIER_GUI or raises the existing
% singleton*.
%
% H = DICOMCOPIER_GUI returns the handle to a new DICOMCOPIER_GUI or the handle to
% the existing singleton*.
%
% DICOMCOPIER_GUI('CALLBACK',hObject,eventData,handles,...) calls the local
% function named CALLBACK in DICOMCOPIER_GUI.M with the given input arguments.
%
% DICOMCOPIER_GUI('Property','Value',...) creates a new DICOMCOPIER_GUI or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before dicomCopier_gui_OpeningFcn gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to dicomCopier_gui_OpeningFcn via varargin.
%
% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one
% instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES
% Edit the above text to modify the response to help dicomCopier_gui
% Last Modified by GUIDE v2.5 31-Oct-2013 09:01:18
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...
'gui_Singleton', gui_Singleton, ...
'gui_OpeningFcn', @dicomCopier_gui_OpeningFcn, ...
'gui_OutputFcn', @dicomCopier_gui_OutputFcn, ...
'gui_LayoutFcn', [] , ...
'gui_Callback', []);
if nargin && ischar(varargin{1})
gui_State.gui_Callback = str2func(varargin{1});
end
if nargout
[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT
% --- Executes just before dicomCopier_gui is made visible.
function dicomCopier_gui_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% varargin command line arguments to dicomCopier_gui (see VARARGIN)
% Choose default command line output for dicomCopier_gui
handles.output = hObject;
% Update handles structure
guidata(hObject, handles);
% UIWAIT makes dicomCopier_gui wait for user response (see UIRESUME)
% uiwait(handles.figure1);
movegui('center')
% --- Outputs from this function are returned to the command line.
function varargout = dicomCopier_gui_OutputFcn(hObject, eventdata, handles)
% varargout cell array for returning output args (see VARARGOUT);
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Get default command line output from handles structure
varargout{1} = handles.output;
% --- Executes on button press in loadbutton.
function loadbutton_Callback(hObject, eventdata, handles)
% hObject handle to loadbutton (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% make clean
if isfield(handles, 'myinfo');
handles = rmfield(handles, 'myinfo');
end
if isfield(handles, 'mydicom');
handles = rmfield(handles, 'mydicom');
end
%[handles.dicomlist, handles.dicomdir] = uigetfile('D:\DICOM\*.*','Select the images to load','MultiSelect', 'on');
[handles.dicomlist, handles.dicomdir] = uigetfile('C:\Users\dewarrat\Documents\current projects\Scapula segmentation\raw\Cases 12.2014 n2\*.*','Select the images to load','MultiSelect', 'on');
if ~handles.dicomdir
return
end
handles.N = length(handles.dicomlist);
set(handles.slider_slice, 'Max', handles.N)
set(handles.slider_slice, 'Value', int32(handles.N/2))
set(handles.slider_slice, 'SliderStep',[1/handles.N 10/handles.N])
set(handles.slider_slice, 'Min', 1)
set(handles.info_text, 'ForegroundColor', 'black', 'FontWeight', 'normal')
% Read dicom data
handles.myinfo(1) = dicominfo(sprintf('%s%s', handles.dicomdir, char(handles.dicomlist(1))));
set(handles.IPP_text, 'String', sprintf('%s (%s%s)', handles.myinfo(1).PatientID, handles.myinfo(1).PatientName.FamilyName(1), handles.myinfo(1).PatientName.GivenName(1)))
set(handles.resolution_text, 'String', sprintf('%f mm', handles.myinfo(1).PixelSpacing(1)))
set(handles.CTdate_text, 'String', sprintf('%s', handles.myinfo(1).AcquisitionDate))
if isfield(handles.myinfo(1), 'ConvolutionKernel')
if strfind(handles.myinfo(1).ConvolutionKernel, 'BONE')
set(handles.radiobutton0, 'Value', 0)
set(handles.radiobutton1, 'Value', 1)
set(handles.radiobutton2, 'Value', 0)
elseif strfind(handles.myinfo(1).ConvolutionKernel, 'STANDARD')
set(handles.radiobutton0, 'Value', 0)
set(handles.radiobutton1, 'Value', 0)
set(handles.radiobutton2, 'Value', 1)
else
set(handles.radiobutton0, 'Value', 1)
set(handles.radiobutton1, 'Value', 0)
set(handles.radiobutton2, 'Value', 0)
end
else
set(handles.radiobutton0, 'Value', 1)
set(handles.radiobutton1, 'Value', 0)
set(handles.radiobutton2, 'Value', 0)
end
% pre-allocate ?
% handles.mydicom = zeros(handles.myinfo(1).Height, handles.myinfo(1).Width, handles.N,'int16');
handles.mydicom(:,:,int32(handles.N/2)) = dicomread(sprintf('%s%s', handles.dicomdir, char(handles.dicomlist(int32(handles.N/2)))));
handles.slice = int32(handles.N / 2);
handles.CTmin = -200;
handles.CTmax = 2000;
handles.bigmax = 3000;
handles.bigmin = -2000;
set(handles.slider_contrast, 'Max', handles.bigmax - handles.bigmin - 1)
set(handles.slider_contrast, 'Min', 0)
set(handles.slider_contrast, 'Value', 900)
set(handles.slider_brightness, 'Max', handles.bigmax)
set(handles.slider_brightness, 'Min', handles.bigmin)
set(handles.slider_brightness, 'Value', 1100)
set(handles.slice_num, 'String', handles.slice)
imshow(handles.mydicom(:,:,handles.slice), [handles.CTmin handles.CTmax])
guidata(hObject,handles)
% --- Executes on slider movement.
function slider_slice_Callback(hObject, eventdata, handles)
% hObject handle to slider_slice (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'Value') returns position of slider
% get(hObject,'Min') and get(hObject,'Max') to determine range of slider
handles.slice = int32(get(hObject,'Value'));
set(handles.slice_num, 'String', handles.slice)
handles.mydicom(:,:,handles.slice) = dicomread(sprintf('%s%s', handles.dicomdir, char(handles.dicomlist(handles.slice))));
imshow(handles.mydicom(:,:,handles.slice), [handles.CTmin handles.CTmax])
guidata(hObject,handles)
% --- Executes during object creation, after setting all properties.
function slider_slice_CreateFcn(hObject, eventdata, handles)
% hObject handle to slider_slice (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: slider controls usually have a light gray background.
if isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor',[.9 .9 .9]);
end
function patient_num_Callback(hObject, eventdata, handles)
% hObject handle to patient_num (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of patient_num as text
% str2double(get(hObject,'String')) returns contents of patient_num as a double
% --- Executes during object creation, after setting all properties.
function patient_num_CreateFcn(hObject, eventdata, handles)
% hObject handle to patient_num (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
% --- Executes on slider movement.
function slider_contrast_Callback(hObject, eventdata, handles)
% hObject handle to slider_contrast (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'Value') returns position of slider
% get(hObject,'Min') and get(hObject,'Max') to determine range of slider
if ~isfield(handles, 'mydicom');
return
end
handles.CTmax = (get(handles.slider_brightness,'Max') + get(handles.slider_brightness,'Min') - get(handles.slider_brightness,'Value')) + (get(hObject,'Max') - get(hObject,'Value'))/2;
handles.CTmin = (get(handles.slider_brightness,'Max') + get(handles.slider_brightness,'Min') - get(handles.slider_brightness,'Value')) - (get(hObject,'Max') - get(hObject,'Value'))/2;
if handles.CTmax > handles.bigmax
handles.CTmax = handles.bigmax;
end
if handles.CTmin < handles.bigmin
handles.CTmin = handles.bigmin;
end
imshow(handles.mydicom(:,:,handles.slice), [handles.CTmin handles.CTmax])
guidata(hObject,handles)
% --- Executes during object creation, after setting all properties.
function slider_contrast_CreateFcn(hObject, eventdata, handles)
% hObject handle to slider_contrast (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: slider controls usually have a light gray background.
if isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor',[.9 .9 .9]);
end
% --- Executes on slider movement.
function slider_brightness_Callback(hObject, eventdata, handles)
% hObject handle to slider_brightness (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'Value') returns position of slider
% get(hObject,'Min') and get(hObject,'Max') to determine range of slider
if ~isfield(handles, 'mydicom');
return
end
handles.CTmax = (get(hObject,'Max') + get(hObject,'Min') - get(hObject,'Value')) + (get(handles.slider_contrast,'Max') - get(handles.slider_contrast,'Value'))/2;
handles.CTmin = (get(hObject,'Max') + get(hObject,'Min') - get(hObject,'Value')) - (get(handles.slider_contrast,'Max') - get(handles.slider_contrast,'Value'))/2;
if handles.CTmax > handles.bigmax
handles.CTmax = handles.bigmax;
end
if handles.CTmin < handles.bigmin
handles.CTmin = handles.bigmin;
end
imshow(handles.mydicom(:,:,handles.slice), [handles.CTmin handles.CTmax])
guidata(hObject,handles)
% --- Executes during object creation, after setting all properties.
function slider_brightness_CreateFcn(hObject, eventdata, handles)
% hObject handle to slider_brightness (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: slider controls usually have a light gray background.
if isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor',[.9 .9 .9]);
end
% --- Executes on button press in copy_button.
function copy_button_Callback(hObject, eventdata, handles)
% hObject handle to copy_button (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
sortPictures = 1;
if isempty(get(handles.patient_num,'String'))
set(handles.info_text, 'ForegroundColor', 'red', 'FontWeight', 'bold', 'String', 'Please give a patient number')
return
else
patient = get(handles.patient_num,'String');
end
if isempty(get(handles.output_dir,'String'))
outputdir = uigetdir;
if outputdir
set(handles.output_dir,'String', outputdir)
end
else
outputdir = get(handles.output_dir,'String');
end
if ~outputdir
set(handles.info_text, 'ForegroundColor', 'red', 'FontWeight', 'bold', 'String', 'Please select an output directoy')
return
end
set(handles.info_text, 'ForegroundColor', 'black', 'FontWeight', 'normal')
% Create name string
newname = [patient, '-', handles.myinfo(1).PatientID, '_', handles.myinfo(1).PatientName.FamilyName(1), handles.myinfo(1).PatientName.GivenName(1)];
if get(handles.radiobutton1,'Value')
newname = [newname, '-bone.'];
elseif get(handles.radiobutton2,'Value')
newname = [newname, '-soft.'];
else
newname = [newname, '.'];
end
% Create output directory
outputdir = [outputdir, '\', patient, '-', handles.myinfo(1).PatientID, '\', 'CT-', patient, '-', handles.myinfo(1).PatientID, '-xxx\dicom\'];
mkdir(outputdir)
% Sort pictures in correct order
if sortPictures == 1
for h = 1:handles.N
infoTemp = dicominfo(sprintf('%s%s', handles.dicomdir, char(handles.dicomlist(h))));
SliceLoc(h,1) = h;
SliceLoc(h,2) = infoTemp.SliceLocation;
end
SliceSorted = sortrows(SliceLoc,2);
end
% Read all data
for i = 1:handles.N
str = sprintf('Copying dicom...\n%d / %d\n', i, handles.N);
set(handles.info_text, 'String', str)
drawnow
% Anonymous
if sortPictures == 1
caseID = SliceSorted(i,1);
else
caseID = i;
end
handles.myinfo(i) = dicominfo(sprintf('%s%s', handles.dicomdir, char(handles.dicomlist(caseID))));
if get(handles.radiobutton1,'Value')
handles.myinfo(i).PatientName.GivenName = [handles.myinfo(i).PatientName.FamilyName(1), handles.myinfo(i).PatientName.GivenName(1), '_bone'];
elseif get(handles.radiobutton2,'Value')
handles.myinfo(i).PatientName.GivenName = [handles.myinfo(i).PatientName.FamilyName(1), handles.myinfo(i).PatientName.GivenName(1), '_soft'];
else
handles.myinfo(i).PatientName.GivenName = [handles.myinfo(i).PatientName.FamilyName(1), handles.myinfo(i).PatientName.GivenName(1)];
end
handles.myinfo(i).PatientName.FamilyName = [patient, '_', handles.myinfo(1).PatientID];
% Read and copy
dicomwrite(dicomread(sprintf('%s%s', handles.dicomdir, char(handles.dicomlist(caseID)))),[outputdir, newname, sprintf('%04d',i), '.dcm'], handles.myinfo(i));
end
set(handles.info_text, 'String', 'Copied')
function output_dir_Callback(hObject, eventdata, handles)
% hObject handle to output_dir (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of output_dir as text
% str2double(get(hObject,'String')) returns contents of output_dir as a double
% --- Executes during object creation, after setting all properties.
function output_dir_CreateFcn(hObject, eventdata, handles)
% hObject handle to output_dir (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
% --- Executes on button press in radiobutton1.
function radiobutton1_Callback(hObject, eventdata, handles)
% hObject handle to radiobutton1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hint: get(hObject,'Value') returns toggle state of radiobutton1
if get(hObject,'Value')
set(handles.radiobutton0,'Value',0)
set(handles.radiobutton2,'Value',0)
else
set(hObject,'Value',1)
end
% --- Executes on button press in radiobutton2.
function radiobutton2_Callback(hObject, eventdata, handles)
% hObject handle to radiobutton2 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hint: get(hObject,'Value') returns toggle state of radiobutton2
if get(hObject,'Value')
set(handles.radiobutton0,'Value',0)
set(handles.radiobutton1,'Value',0)
else
set(hObject,'Value',1)
end
% --- Executes on button press in radiobutton0.
function radiobutton0_Callback(hObject, eventdata, handles)
% hObject handle to radiobutton0 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hint: get(hObject,'Value') returns toggle state of radiobutton0
if get(hObject,'Value')
set(handles.radiobutton1,'Value',0)
set(handles.radiobutton2,'Value',0)
else
set(hObject,'Value',1)
end
c4science · Help