How to plot a point in 3D data displayed in volume Viewer?

I have a 3D medical image.

I want to mark a certain point in that image and display the 3D image with volumeViewer.

The point has x,y,z coordinates.

How can I plot this point into the 3D image and display both in volume Viewer?

[Dunja Gorup]

I have a same question.

Logical indexing is not working.

E.g.

I=logical(zeros(size(Your_volume_size)));

I(XYZCoordinates_Array)=true; %not working

volshow(I); %displays only a line, not a shape, or point cloud. ptcloud, volumes and lines are not supported to display by volshow with logical masks. The answer below does not do anything, even though it was accepted.

Thank you very much

[Rik]

What value does XYZCoordinates_Array have? And what do you interpret as not working?

[Dunja Gorup]

XYZCoordinates_Array was positive integers array (:,:,:).

I(XYZ)=true; %does not preform as expected for e.g. XYZ=[10,13,9;12,11,8;1,2,3;3,4,5;1,7,5];

I(XYZ(1,:))=true % does not produce the same result as I(10,13,9)=true;

The solution was achieves by sub2index as below.

[Rik]

The reason why this is the case, is that the indexing does not produce a comma separated list. You can produce such a list with a struct array and with a cell array:

XYZ=[10,13,9;12,11,8;1,2,3;3,4,5;1,7,5];
XYZ(1,:) % normal array

ans = 1×3

10 13 9

tmp = num2cell(XYZ);
tmp{1,:} % comma separated list

ans = 10

ans = 13

ans = 9

Answer 1

[Selva Karna]

To view 3d Volume:

clc

clear all;

close all;

your volume=V;

volshow(V)

3Point view:

% Make mask 3d based on your 3d points

len=length(your 3d points);

% Create Mask

3d_mask=zeros(your_3d points(size));

3d_mask(:,:,:)=your_3d_mask_holes;

volshow(3d_mask)

[Dunja Gorup]

How is that accepted answer?

[Rik]

Perhaps the OP read the answer with some extra formatting (and replaced the invalid variable name):

clc
clear all;
close all;

your volume=V;

volshow(V)

3Point view:

% Make mask 3d based on your 3d points
len=length(your 3d points);
% Create Mask
3d_mask=zeros(your_3d points(size));
3d_mask(:,:,:)=your_3d_mask_holes;
volshow(3d_mask)

[Dunja Gorup]

Following solution worked by subindexing:

AL=zeros(50,50,50);
XYZ=[10,13,9;12,11,8;1,2,3;3,4,5;1,7,5];

ind = sub2ind(size(AL), XYZ(:,1), XYZ(:,2), XYZ(:,3));
AL(ind) = true;

volshow(AL);

[Dunja Gorup]

However, this mask is sometimes offset to the original image from which the voxel values were derived. Probably due to some transformation?