How to transfer all the coordinate binary into my image dicom

Dear all,

I was do segmentation in my dicom image. then I got the binary image.

then how to transfer all the coordinate binary into my image dicom? So that only binary 1 appear in my dicom.

I've tried below but gor error.

My dicom image as attached.

clc
clear all
close all
% read dicom image
[spect map]=dicomread('I-13125610N1.dcm');
info = dicominfo('I-13125610N1.dcm');
%gp=info.SliceThickness;
spect=(squeeze(spect));%smooth3
aa=size(spect);aa=aa(3);
imshow3D(spect)
% do segmentation using adaptthresh
T = adaptthresh(spect, 0.000000000000001, "NeighborhoodSize", 27);
% change to binary images
BW = imbinarize(spect,T);
figure, imshow3D(BW)

% transfer all the coordinate binary into my image dicom
p = spect.*BW;

[Image Analyst]

Not sure what you want. Do you want the BW image's white areas to appear in your original image as pure white (or some other uniform color like you'd get with imoverlay)?