Importing 3ddose

[Dominic Rafferty]

Hello,

I am trying to import a 3ddose file (generated by DOSXYZNRC) to an existing planC object using CERR. However, when I attempt to view the scan and dose distribution using the GUI there doesn't seem to be any dose there :( Can someone please suggest what I might be doing incorrectly?

Here is my script

%Open .mat/.mat.bz2 file

close all;

clear all;

fprintf('Please select .mat/.mat.bz2 file containing CERR-imported DICOM series and structure set\n\n');

[filename,path] = uigetfile('*'); %gets file

global planC;

planC = openPlanCFromFile(fullfile(path,filename));

%Open file containing dose distribution file (3ddose format output by
%DOSXYZNRC)

fprintf('Please select .3ddose file containing dose distribution\n\n');

[dosefilename,dosepath] = uigetfile('*'); %gets file

[dose,error,xbound,ybound,zbound] = dicomrt_read3ddose(1,fullfile(dosepath,dosefilename));

% Normalise dose distribution

dose = dose/max(max(max(dose)));

indexS=planC{end};

assocScanNum=1;

assocScanUID=planC{indexS.scan}(assocScanNum).scanUID;

dx = mean(diff(xbound));

dy = mean(diff(ybound));

dz = mean(diff(zbound));

zV = -zbound(end:-1:2) + dz;

regParamsS.horizontalGridInterval = dx;

regParamsS.verticalGridInterval   = -dy;

regParamsS.coord1OFFirstPoint   = xbound(1) + dx/2;

regParamsS.coord2OFFirstPoint   = ybound(end) - dy/2;

regParamsS.zValues  = zV;

dose2CERR(dose,[], 'MC_doseNew','test','test','non CT',regParamsS,'no',assocScanUID);

planC{indexS.dose}.doseArray = flipdim(planC{indexS.dose}.doseArray,3);

save_planC(planC,[], 'passed', 'test.bz2'C);

[ap...@mskcc.org]

Hi, I think there might be a mismatch in coordinate specification. Can you provide sample data?

--
You received this message because you are subscribed to the Google Groups "CERR FORUM" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cerr-forum+...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cerr-forum/94b941e7-9e9d-49a9-b829-9299aae2f432n%40googlegroups.com.

*** Only open attachments or links from trusted senders. Report phishing to inf...@mskcc.org ***

 

=====================================================================

Please note that this e-mail and any files transmitted from
Memorial Sloan Kettering Cancer Center may be privileged, confidential,
and protected from disclosure under applicable law. If the reader of
this message is not the intended recipient, or an employee or agent
responsible for delivering this message to the intended recipient,
you are hereby notified that any reading, dissemination, distribution,
copying, or other use of this communication or any of its attachments
is strictly prohibited. If you have received this communication in
error, please notify the sender immediately by replying to this message
and deleting this message, any attachments, and all copies and backups
from your computer.

[JO]

Hi.

Did you ever solve this? If yes, what was the solution?

For similar workflow, I'm using lines below, where seems to be some differences compared to yours (I'm doing flipdim before this block):

global planC

indexS = planC{end};

assocScanNum=1;

assocScanUID=planC{indexS.scan}(assocScanNum).scanUID ;

dx = mean(diff(xbound));

dy = mean(diff(ybound));

dz = mean(diff(zbound));

zV = -zbound(end:-1:2) + dz/2;-------------------------------------------------------------------difference here

regParamsS.horizontalGridInterval = dx;

regParamsS.verticalGridInterval   = -dy;

regParamsS.coord1OFFirstPoint   = xbound(1) + dx/2;

regParamsS.coord2OFFirstPoint   = -ybound(1) - dy/2;-------------------------------------------------------------------difference here

regParamsS.zValues  = zV;

regParamsS.doseUnits = 'GRAYS';

regParamsS.dateWritten = date;

dose2CERR(DxyzABS,[], 'MC_doseNew','test','test','non-CT',regParamsS,'no',assocScanUID);-------------------------------------------------------------------difference here ('non-CT')

[JO]

Hello Dominic.

Did you ever solve this? If yes, what was the solution?

Br,

JO 

[Dominic Rafferty]

Hi Jo,

Apologies for late response I am not receiving alerts on this query. Yes I did manage to solve the issue and I pretty much did the same as you. Thanks for your help

[JO]

Great!

Do you have any comments to those differences I have indicated in my previous post?

BR,

JO