function Fusiona(imCT,trozo,SUV,HU_MIN,HU_MAX,gp)
% tam1=size(imCT,1);
% tam2=size(imCT,2);
% aux=zeros(tam1,tam2);
% aux=imCT;
color = zeros(512,512,3);
color(:,:,3) = (imCT-min(min(imCT)))./(max(max(imCT))-min(min(imCT)));
color(:,:,2) = (trozo-min(min(trozo)))./(max(max(trozo))-min(min(trozo)));
if(gp==0)
color(:,:,1) = (trozo > SUV) & (imCT > HU_MIN) & (imCT < HU_MAX);
end
% imtool(color)
% %%
% color = zeros([size(imPETscaled,1) size(imPETscaled,2) 3]);
% posinic = round((CT.ImagePositionPatient(1:2) - PETscaled.ImagePositionPatient(1:2))./CT.PixelSpacing);
% color(posinic(1):posinic(1)+size(imCT,1)-1,posinic(2):posinic(2)+size(imCT,2)-1 ,3) = (imCT-min(min(imCT)))./(max(max(imCT)-min(min(imCT))));
% color(:,:,2)= (imPETscaled-min(min(imPETscaled)))./(max(max(imPETscaled)-min(min(imPETscaled))));
% imtool(color);
%
%
% h = figure (100);
% h.ToolBar= 'none';
% h.Position = [1 1 1024 1024];
% subplot('Position', [0.5 0.5 0.5 0.5]); image(im2uint8(color(:,:,3))); axis('image'); axis('off');
% subplot('Position', [0.0 0.5 0.5 0.5]); image(im2uint8(color(:,:,2))); axis('image'); axis('off');
% subplot('Position', [0.5 0.0 0.5 0.5]); imshow(color,[]);
alpha = 0.33;
% subplot('Position', [0.0 0.0 0.5 0.5]); image(im2uint8(color(:,:,3)));
% hold on; subplot('Position', [0.0 0.0 0.5 0.5]); hPT = image(uint16(im2uint8(color(:,:,2)))+256);hold off;
% set(hPT, 'alphadata',alpha*ones(size(color(:,:,3))))
% axis('image'); axis('off');
% subplot('Position', [0.5 0.0 0.5 0.5]); image(im2uint8(color(:,:,3)));
% hold on; subplot('Position', [0.5 0.0 0.5 0.5]); hPT = image(uint16(im2uint8(color(:,:,2)))+256);hold off;
% alpha = 0.25;
% set(hPT, 'alphadata',alpha*ones(size(color(:,:,3))))
%
% hold on; subplot('Position', [0.5 0.0 0.5 0.5]); hBAT = image(uint16(im2uint8(color(:,:,1)))+512);hold off;
% alpha = 0.5;
%
% set(hBAT, 'alphadata',alpha*ones(size(color(:,:,3))))
%
% axis('image'); axis('off');
% commonCMAP = [gray(256); hot(256); zeros(255,3) ; [0 0 1]];
PET_BAT = im2uint8(color(:,:,2));
PET_BAT(PET_BAT == 255) = 254;
PET_BAT(color(:,:,1) == 1) = 255;
%h = figure(100);
%subplot('Position', [0.5 0.0 0.5 0.5]);
image(im2uint8(color(:,:,3)));
hold on;
%subplot('Position', [0.5 0.0 0.5 0.5]);
hPT = image(uint16(PET_BAT)+256);
hold off;
axis('image'); axis('off');
set(hPT,'alphadata',alpha*ones(size(color(:,:,3))));
commonCMAP = [gray(256); hot(256)];
commonCMAP(512,:) = [0 0 1];
colormap(commonCMAP);