How to replace 0s between 1s with 1s "filling in" meshgrid-derived XYZ matrices?

Question

Daniel Bridges 2016-7-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/294274-how-to-replace-0s-between-1s-with-1s-filling-in-meshgrid-derived-xyz-matrices

回答： KSSV 2016-7-12

DICOM ContourData gives points representing organ structures. (Originally a vector for each slice, [x1 y1 z1 x2 y2 z2 ...], I have created three vectors grouping x y z values.)

figure; scatter3(xcontour,ycontour,zcontour,'.')
>> title('PTV1 DICOM Coordinates')

I have another matrix with RTDOSE values. I have created a meshgrid (X,Y,Z) giving the dose's DICOM coordinates and, I think, derived another matrix from it that changes these PTV1 coordinates into 1s with the empty space between them into 0s:

xcontourpointslocation = ismember(X,xcontour);
ycontourpointslocation = ismember(Y,ycontour);
zcontourpointslocation = ismember(Z,zcontour);

So I think now this figure can be thought to represent in "DICOM space" where those 1s are. If I multiply it by the dose matrix, I expect to thereby " extract" the dose at these locations. Yet I obviously want the dose "inside" the structure as well -- i.e., I must convert all those 0s in xcontourpointslocation etc "between the 1s" into 1s to "fill in" the structure.

What MATLAB commands do this elegantly? I have an idea to use the find command with conditioned loops or for loops, but I expect there to be a simpler way, e.g. a command like makesolid(matrix), because this seems straightforward or commonplace.

While hoping you will tell me such a method, I have found some related discussion:

I will continue working on this problem and post my answer if I solve it before hearing from you.

3 个评论
显示 1更早的评论隐藏 1更早的评论

José-Luis 2016-7-8

Would you be happy assuming your organ is convex?

Daniel Bridges 2016-7-12

I want to access all the points corresponding to the interior of that volume and the points on its surface in a corresponding matrix containing dose. My thinking is that if I can use a same-sized meshgrid with all the points on surface & interior being 1, otherwise 0, then I can achieve this through element-wise matrix multiplication.

I am trying to avoid assumptions, since ideally the program should be robust to work with any structure. Moreover, structures may be complicated (e.g. a mixture of concave and convex) as their shape is a function of both anatomy and disease-spread.

请先登录，再进行评论。

请先登录，再回答此问题。