当前位置：首页 > news >正文

Learning vtkjs之ImageMarchingCubes

news 来源：原创 2025/5/1 12:30:07

体积等值面处理

介绍

vtkImageMarchingCubes - 对体积进行等值面处理

给定一个指定的等值，使用Marching Cubes算法生成一个等值面。

效果

在这里插入图片描述
新建了一个球，对比一下原始的（透明的）和ISO的效果

核心代码

参数部分

const updateParam = (type, obj) => {const { sample, sphere, sphereSource, mCubes, renderWindow } =context.current;let value;switch (type) {case "volumeResolution":value = obj.volumeResolution;sample.setSampleDimensions(value, value, value);renderWindow.render();setState({ ...state, ...obj });break;case "isoValue":value = obj.isoValue;mCubes.setContourValue(value);console.log(value);renderWindow.render();setState({ ...state, ...obj });break;case "sphereRadius":value = obj.sphereRadius;sphere.setRadius(value);  renderWindow.render();setState({ ...state, ...obj });break;case "computeNormals":value = obj.computeNormals;mCubes.setComputeNormals(!!value);renderWindow.render();setState({ ...state, ...obj });break;case "mergePoints":value = obj.mergePoints;mCubes.setMergePoints(!!value);renderWindow.render();setState({ ...state, ...obj });break;default:break;}};

主体流程

const fullScreenRenderer = vtkFullScreenRenderWindow.newInstance({background: [0, 0, 0],rootContainer: vtkContainerRef.current,});const renderer = fullScreenRenderer.getRenderer();const renderWindow = fullScreenRenderer.getRenderWindow();const actor = vtkActor.newInstance();renderer.addActor(actor);const mapper = vtkMapper.newInstance();actor.setMapper(mapper);// Build pipelineconst radius = 1;const sampleDimension = 50;const sphere = vtkSphere.newInstance({center: [0.0, 0.0, 0.0],radius: radius,});// value of implicit function // x = x1 - center// x*x + y*y + z*z - R*R = valueconst sample = vtkSampleFunction.newInstance({implicitFunction: sphere,sampleDimensions: [sampleDimension, sampleDimension, sampleDimension],modelBounds: [-radius, radius, -radius, radius, -radius, radius],});// Create the isosurface by contourValueconst mCubes = vtkImageMarchingCubes.newInstance({ contourValue: 0.0 });// Connect the pipeline propermCubes.setInputConnection(sample.getOutputPort());mapper.setInputConnection(mCubes.getOutputPort());context.current = {sample,mCubes,sphere,renderer,renderWindow,};// 同步sphereconst sphereActor = createSphere();renderer.addActor(sphereActor);renderer.resetCamera();renderWindow.render();