path='RT-STRUCT.dcm'
info=dicominfo(path);
contours=dicomContours(info);
ROIs=contours.ROIs;
gtv_index=ROIs.Number(strcmp(ROIs.Name, 'GTV'));
plotContour(contours,gtv_index);
xlim=xlim;
ylim=ylim;
zlim=zlim;
x_size=xlim(2)-xlim(1);
y_size=ylim(2)-ylim(1);
z_size=zlim(2)-zlim(1);
points_per_mm=2;
referenceInfo = imref3d([round(x_size*points_per_mm),round(y_size...
*points_per_mm),round(z_size*points_per_mm)],xlim,ylim,zlim);
rtMask = createMask(contours, gtv_index, referenceInfo);
figure;
volshow(rtMask);
I have attached the result of plotContour (as a .fig) and volshow (as both .fig and png as the .fig had some issues loading) to the post, and as you can see the two are quite different.
I offer previews below, here is the contour plot
which as you can see has no flat or straight edges.
However, the volshow rendering has one completely flat side, which the contour plot does not.
Why is this happening and is there any way to change the creation of the mask in volshow to be more accurate? It seems like MATLAB is either doing some autofilling outside the contours or that my xlim, ylim, zlim arguments are causing the mask to be abruply cut off, however I can't see why it would be the latter as my xlim, ylim, zlim comes directly from the boundaries of the contourplot
EDIT: I got some clarity, if instead of
referenceInfo = imref3d([round(x_size*points_per_mm),round(y_size...
*points_per_mm),round(z_size*points_per_mm)],xlim,ylim,zlim);
I just write
referenceInfo = imref3d([128,128,50],xlim,ylim,zlim);
then the resulting figure is much more reasonable
But I still don't quite understand why. I was led to believe that the first argument of imref3d only corresponds to the size of the output matrix, which presumably should not deform the actual shape that gets rendered when createMask is called, so why are the shapes different if all I have changed is the image size but kept world limits the same?
Thank you in advance for any responses
