Hi,
As per my understanding, you are trying to render two 3D volumes with outer volume (V2) to be transparent using "volshow" function and facing issues related to blurring.
The following code uses 'AlphaData' property of the "volshow" function to render a transparent outer volume.
% dimensions of the volume
dim = 100;
[x, y, z] = ndgrid(1:dim, 1:dim, 1:dim);
% centers for V1 and V2
center = [dim/2, dim/2, dim/2];
% radii for V1 and V2
radiusV1 = 20;
radiusV2 = 40;
% binary masks for V1 and V2
V1 = sqrt((x - center(1)).^2 + (y - center(2)).^2 + (z - center(3)).^2) <= radiusV1;
V2 = sqrt((x - center(1)).^2 + (y - center(2)).^2 + (z - center(3)).^2) <= radiusV2;
% ensure V1 is fully contained within V2 by combining the masks correctly
V2 = V2 & ~V1;
% create an RGB volume
V = zeros([size(V1), 3]);
% assign V1 to the green channel and V2 to the blue channel
V(:,:,:,2) = V1; % green channel for V1
V(:,:,:,3) = V2; % blue channel for V2
% initialize AlphaData with full transparency
Adata = zeros(size(V1));
% make V2 semi-transparent, but only where V1 is not present
Adata(V2 > 0 & V1 == 0) = 0.01;
% ensure V1 remains fully opaque
Adata(V1 > 0) = 1;
% Display the volume with volshow
h = volshow(V);
viewer = h.Parent;
viewer.BackgroundColor="white";
viewer.BackgroundGradient="off";
% set alpha data
h.AlphaData = Adata;
Output :
Note : The above code was run on R2024a.
For viewing edges of the volume with some transparency, the 'RenderingStyle' property (specifically, 'GradientOpacity' rendering style) seems to be more relevant. Refer to the following documentation links for more information:
The modified code using 'RenderingStyle' property is as follows :
% dimensions of the volume
dim = 100;
[x, y, z] = ndgrid(1:dim, 1:dim, 1:dim);
% centers for V1 and V2
center = [dim/2, dim/2, dim/2];
% radii for V1 and V2
radiusV1 = 20;
radiusV2 = 40;
% binary masks for V1 and V2
V1 = sqrt((x - center(1)).^2 + (y - center(2)).^2 + (z - center(3)).^2) <= radiusV1;
V2 = sqrt((x - center(1)).^2 + (y - center(2)).^2 + (z - center(3)).^2) <= radiusV2;
% ensure V1 is fully contained within V2 by combining the masks correctly
V2 = V2 & ~V1;
% create an RGB volume
V = zeros([size(V1), 3]);
% assign V1 to the green channel and V2 to the blue channel
V(:,:,:,2) = V1; % green channel for V1
V(:,:,:,3) = V2; % blue channel for V2
% display the volume with volshow
h = volshow(V);
viewer = h.Parent;
viewer.BackgroundColor="white";
viewer.BackgroundGradient="off";
% set rendering style
h.RenderingStyle = "GradientOpacity";
Output:
Hope this helps!
