Plotting the iso-surface in a specific range of values

Hi,

I have a 3D matrix and I want to show the isosurface for a range of values, such as the values greater than 0.2, not for a single value.

Could anybody help me to revise my code?

data = smooth3(porosity_median,'box',5);
p1 = patch(isosurface(data, 0.1 ), 'FaceColor','red','EdgeColor','none');
view(3); axis vis3d tight
camlight; lighting phong

Answer 1

[José M. Requena Plens]

You need ISOCAPS with 'enclose' option.

Using yor code.

To plot values greater than 0.1:

data = smooth3(porosity_median,'box',5);
p1 = patch(isosurface(data, 0.1 ), 'FaceColor','red','EdgeColor','none'); 
p2 = patch(isocaps(data,0.1,'enclose','above'),'FaceColor','red','EdgeColor','none');
view(3); axis vis3d tight
camlight; lighting phong

To plot values lower than 0.1:

data = smooth3(porosity_median,'box',5);
p1 = patch(isosurface(data, 0.1 ), 'FaceColor','red','EdgeColor','none'); 
p2 = patch(isocaps(data,0.1,'enclose','below'),'FaceColor','red','EdgeColor','none');
view(3); axis vis3d tight
camlight; lighting phong

And I recommend that you use the ISONORMALS to better graphics.

isonormals(data,p1)

[hamed]

Thanks Jose.

I get this error when I use isonormals instead of isosurface:

index exceeds matrix dimensions.

Error in isonormals (line 75)

n(:,1)=interp3(x, y, z, nx, verts(:,1), verts(:,2), verts(:,3));

[José M. Requena Plens]

You should only apply it to the isosurface patch 'p1'.

Example:

data = smooth3(randn(50,50,50),'box',5);
p1 = patch(isosurface(data, 0.1 ), 'FaceColor','red','EdgeColor','none');
isonormals(data,p1)
p2 = patch(isocaps(data,0.1,'enclose','above'),'FaceColor','red','EdgeColor','none');
view(3); axis vis3d tight
camlight; lighting phong

Do you still having error?

[hamed]

I dont get error, but the output does not make sense for my data. I set to get the isosurfaces for above 0.3 and I know that the maximum value in my data is 0.27! but I get this figure, for all values I set:

using this command:

data = smooth3(porosity_median,'box',5);
p1 = patch(isosurface(data, 0.3 ), 'FaceColor','red','EdgeColor','none'); 
isonormals(data,p1)
p2 = patch(isocaps(data,0.1,'enclose','above'),'FaceColor','red','EdgeColor','none');
view(3); axis vis3d tight
camlight; lighting phong

[José M. Requena Plens]

The value you want to represent must be defined in both patches.

...
p1 = patch(isosurface(data, 0.3 ), 'FaceColor','red','EdgeColor','none'); 
...
p2 = patch(isocaps(data,0.3,'enclose','above'),'FaceColor','red','EdgeColor','none');
...

[hamed]

Thanks Jose.

Answer 2

[darova]

Use for loop to create several value