create 3D image from coordinates
显示 更早的评论
I have got a 3D nifti image (img) that is a breast MRI. For the same image, there are NX3 coordinates (x,y,z) in the file coordsV2f. I need to substitute these coordinates to new values that are in the file pred.txt, in order to highlight some pxels and show them marked in the original 3D image. It works and now I have my zeros image with some highlighted pixel. How can I visualize the marked pixels on the original one (img) and scroll down the slices?
many thanks !!!
img = niftiread([fCompletePath '_image']);
load (['19012302a_coordsV2f']);
x = RCS(:,1);
y = RCS(:,2);
z = RCS(:,3);
new_img = zeros(size(img));
k = 1;
for slice = 1 : length(z)
f = z(slice); %3rd dimension of RCS
newX= (x(k)); % get x
newY= (y(k)); % get y
new_img(newX,newY,f) = pred(k); %assign new value on a zeros image
k = k+1;
end
imshow(new_img(:,:,5),'DisplayRange',[])
save new_img %save in workspace
niftiwrite(new_img, [fName '_prediction.nii.gz']); %save as .nii
采纳的回答
Image Analyst
2019-12-24
0 个投票
You can extract each slice of the image in a for loop and call imshow(), then call hold on, and plot() to plot a red spot (or whatever) over the changed pixels. Then at the bottom of the loop call questdlg() to let the user see the image and spots before proceeding with the next image.
2 个评论
France
2019-12-28
Image Analyst
2019-12-28
You certainly would not want to do that. Let's say img is a 3-D RGB image with 3 color channels. Let's say each channel was a megapixel. So you'd be calling imshow() a million times since n would go from 1 to a million. To plot a red dot over the changed pixels stored in x and y, you'd do
imshow(img);
hold on;
plot(x, y, 'r.', 'MarkerSize', 2);
hold off;
questdlg() asks the user a question. It takes a string as a user prompt, and some strings for response options the user can pick. Like
promptMessage = sprintf('Do you want to Continue processing,\nor Quit processing?');
titleBarCaption = 'Continue?';
buttonText = questdlg(promptMessage, titleBarCaption, 'Continue', 'Quit', 'Continue');
if contains(buttonText, 'Quit')
return;
end
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Image Processing Toolbox 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
