Извлеките данные об изоповерхности из данных об объеме
Изоповерхность является 3-D поверхностным представлением точек с равными значениями в 3-D распределении данных. isosurface функция вычисляет и чертит поверхность соединительными точками постоянного значения в объеме пробела.
Создайте 3-D декартовы координаты при помощи meshgrid и создайте данные об объеме V. Затем создайте изоповерхность в объеме соединительными точками в значении.0001. Постройте изоповерхность.
[x,y,z] = meshgrid([-3:0.25:3]); V = x.*exp(-x.^2 -y.^2 -z.^2); isosurface(x,y,z,V,1e-4);
Создайте изоповерхность, где изозначение 0.0001. Возвратите местоположения поверхностей и вершин для получившейся поверхности в структуре.
[x,y,z] = meshgrid([-3:0.25:3]); V = x.*exp(-x.^2 -y.^2 -z.^2); s = isosurface(x,y,z,V,1e-4)
s = struct with fields:
vertices: [1693x3 double]
faces: [3348x3 double]
Создайте изоповерхность, где изозначение.0001. Возвратите местоположения поверхностей и вершин для получившейся поверхности как отдельные массивы.
[x,y,z] = meshgrid([-3:0.25:3]); V = x.*exp(-x.^2 -y.^2 -z.^2); [faces,verts] = isosurface(x,y,z,V,1e-4)
faces = 3348×3
1 2 3
1 4 2
2 4 5
4 6 5
5 6 7
6 8 7
9 10 11
9 12 10
10 12 13
12 14 13
⋮
verts = 1693×3
0.2500 -0.2500 -2.7918
0.2181 -0.2500 -2.7500
0.2500 -0.4366 -2.7500
0.2500 0 -2.8092
0.2049 0 -2.7500
0.2500 0.2500 -2.7918
0.2181 0.2500 -2.7500
0.2500 0.4366 -2.7500
0.5000 -0.7500 -2.7935
0.4165 -0.7500 -2.7500
⋮
Создайте изоповерхность, где изозначение 0.0001. Вычислите местоположения поверхностей и вершин как структура. Затем постройте поверхность и данные о вершине путем передачи данных о структуре patch функция. Настройте представление графика, измените поверхностные цвета и задайте пользовательскую подсветку.
[x,y,z] = meshgrid([-3:0.25:3]); V = x.*exp(-x.^2 -y.^2 -z.^2); s = isosurface(x,y,z,V,1e-4);
p = patch(s); isonormals(x,y,z,V,p) view(3); set(p,'FaceColor',[0.5 1 0.5]); set(p,'EdgeColor','none'); camlight; lighting gouraud;
Визуализируйте данные о потоке, но цветовой код поверхность, чтобы указать на величину вдоль оси X. Во-первых, вычислите цветные данные путем определения вектора, содержащего скалярное значение для каждой вершины в изоповерхности. Цвета сопоставлены с текущей палитрой.
[x,y,z] = meshgrid([-3:0.25:3]); V = x.*exp(-x.^2 -y.^2 -z.^2); [faces,verts,colors] = isosurface(x,y,z,V,1e-4,x);
Затем постройте изоповерхность с вычисленными цветными данными путем передачи данных функции закрашенной фигуры. Установите угол камеры и соотношение сторон и измените палитру.
patch('Vertices',verts,'Faces',faces,'FaceVertexCData',colors,... 'FaceColor','interp','EdgeColor','interp') view(3) colormap copper
Если вы хотите задать пользовательское поле зрения камеры, подсветку или различную карту цветов поверхности, используйте patch функция, чтобы отобразить на графике данные, возвращенные в isosurface. Можно передать структуру, созданную isosurface непосредственно к patch команда. Например:
s = isosurface(X,Y,Z,V,isovalue); patch(s)
В качестве альтернативы можно передать отдельные поверхности и массивы вершин, чтобы исправить путем определения имен свойства. Например:
[faces,verts] = isosurface(X,Y,Z,V,isovalue) ; patch('Faces',faces,'Vertices',verts);
isonormals | shrinkfaces | smooth3 | contour | meshgrid | contourslice