c# – 如何创建vtk七角棱镜(7个面有4个顶点,2个面有7个顶点)
内容导读
互联网集市收集整理的这篇技术教程文章主要介绍了c# – 如何创建vtk七角棱镜(7个面有4个顶点,2个面有7个顶点),小编现在分享给大家,供广大互联网技能从业者学习和参考。文章包含8193字,纯文字阅读大概需要12分钟。
内容图文
![c# – 如何创建vtk七角棱镜(7个面有4个顶点,2个面有7个顶点)](/upload/InfoBanner/zyjiaocheng/796/b115767aa4404aaea4fb0254b97ada6e.jpg)
在vtk中,您可以使用InsertNextCell(int type,vtkIdType npts,vtkIdType * ptIds)将六边形棱镜单元添加到vtkUnstructuredGrid.这很好,但我想在一个vtkUnstructuredGrid中包含一个七角形和八角形棱镜.单元格类型未在vtk中定义…
只是为了清除一些疑问,我正在使用c#但c或python示例就足以将其转换为c#.
回到问题,我首先复制并略微修改代码in this example以获得体素(8分和celltype = 11)
vtkPoints points = vtkPoints.New();
points.InsertNextPoint(-1.0, -1.0, -1.0);
points.InsertNextPoint(1.0, -1.0, -1.0);
points.InsertNextPoint(1.0, 1.0, -1.0);
points.InsertNextPoint(-1.0, 1.0, -1.0);
points.InsertNextPoint(-1.0, -1.0, 1.0);
points.InsertNextPoint(1.0, -1.0, 1.0);
points.InsertNextPoint(1.0, 1.0, 1.0);
points.InsertNextPoint(-1.0, 1.0, 1.0);
vtkUnstructuredGrid unstructuredGrid1 = vtkUnstructuredGrid.New();
unstructuredGrid1.SetPoints(points);
int[] ptIds = new int[] { 0, 1, 2, 3, 4, 5, 6, 7 };
IntPtr ptIdsPointer = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(int)) * 8);
Marshal.Copy(ptIds, 0, ptIdsPointer, 8);
unstructuredGrid1.InsertNextCell(11, 8, ptIdsPointer);
Marshal.FreeHGlobal(ptIdsPointer);
// Create a mapper and actor
vtkDataSetMapper mapper = vtkDataSetMapper.New();
mapper.SetInputConnection(unstructuredGrid1.GetProducerPort());
vtkActor actor = vtkActor.New();
actor.SetMapper(mapper);
actor.GetProperty().SetColor(1, 0, 0);
actor.GetProperty().SetPointSize(3);
actor.GetProperty().SetOpacity(1);
actor.GetProperty().SetLineWidth(2);
actor.GetProperty().EdgeVisibilityOn();
Renderer3D.AddActor(actor);
这显示了正确的体素:
但是因为我想要一个七面体,我试图使用VTK_POLYHEDRON调整这个想法,更改告诉类型pf单元添加到非结构化网格的线
unstructuredGrid1.InsertNextCell(11, 8, ptIdsPointer);
to(42代表VTK_POLYHEDRON,11代表VTK_VOXEL)和points = 14
unstructuredGrid1.InsertNextCell(42, 14, ptIdsPointer);
像this example.我尝试过类似c的东西,但我没有得到任何结果
int numberOfVertices = 14;
int numberOfFaces = 9; //how to tell heptagonal prism has 2 faces with 7 vertices and 7 with 4 vertices?
int numberOfFaceVertices = 4;
vtkIdType heptagonPointsIds[14] =
{0, 1, 2, 3, 4, 5, 6
7, 8, 9, 10, 11, 12,13};
//points that only change in z
vtkSmartPointer<vtkPoints> points = vtkSmartPointer<vtkPoints>::New();
points->InsertNextPoint(454.6, 678.534, 1.1009301757813); //id 0
points->InsertNextPoint(454.6, 678.545, 1.1500146484375);//id 1
points->InsertNextPoint(454.6, 678.556, 1.2255187988281);//id 2
points->InsertNextPoint(455.1, 678.556, 1.6478076171875);//id 3
points->InsertNextPoint(455.1, 678.534, 1.7926538085938);//id 4
points->InsertNextPoint(455.1, 678.545, 1.7166479492188);//id 5
points->InsertNextPoint(454.85, 678.534, 1.0092297363281);//id 6
points->InsertNextPoint(454.6, 678.534, 5.1009301757813);//id 7
points->InsertNextPoint(454.6, 678.545, 5.1500146484375);//id 8
points->InsertNextPoint(454.6, 678.556, 5.2255187988281);//id 9
points->InsertNextPoint(455.1, 678.556, 5.6478076171875);//id 10
points->InsertNextPoint(455.1, 678.534, 5.7926538085938);//id 11
points->InsertNextPoint(455.1, 678.545, 5.7166479492188);//id 12
points->InsertNextPoint(454.85, 678.534, 5.0092297363281);//id13
vtkIdType heptagonPrismFace[9][4] = {
{0, 1, 7, 8},
{1, 2, 8, 9},
{2, 3, 9, 10},
{3, 4, 10, 11},
{4, 5, 11, 12},
{5, 6, 12, 13},
{6, 1, 13, 8},
//last two faces how to tel they should be 0 1 2 3 4 5 6 and 7 8 9 10 11 12 13 ?
{0, 1, 2, 3},
{7, 8, 9, 10},
};
vtkSmartPointer<vtkCellArray> heptagonFaces=
vtkSmartPointer<vtkCellArray>::New();
for (int i = 0; i < numberOfFaces; i++)
{
heptagonFaces->InsertNextCell(numberOfFaceVertices, heptagonPrismFace[i]);
}
vtkSmartPointer<vtkUnstructuredGrid> uGrid = vtkSmartPointer<vtkUnstructuredGrid>::New();
uGrid->InsertNextCell(VTK_POLYHEDRON,
numberOfVertices, heptagonPointsIds,
numberOfFaces, heptagonFaces->GetPointer());
uGrid->SetPoints(points);
vtkSmartPointer<vtkDataSetMapper> mapper = vtkSmartPointer<vtkDataSetMapper>::New();
mapper->SetInputConnection(uGrid->GetProducerPort());
//add actor to render and stuff...
我怎么会得到像图像这样的东西:
解决方法:
由于没有七边形和八边形棱镜的定义,您必须手动定义所有面,然后将每个面添加到vtkUnstructuredGrid,例如检查此八角形棱镜定义(这应该让您开始):
vtkPoints points = vtkPoints.New();
//up
points.InsertNextPoint(0.0, 0.0, 0.0 );//0
points.InsertNextPoint(0.3, -0.5, 0.0);//1
points.InsertNextPoint(0.6, -.8, 0.0);//2
points.InsertNextPoint(0.7, -.9, 0.0);//3
points.InsertNextPoint(0.8, -0.8, 0.0);//4
points.InsertNextPoint(1.0, 0.0, 0.0);//5
points.InsertNextPoint(1.0, 1.0, 0.0);//6
points.InsertNextPoint(0.0, 1.0, 0.0);//7
//down
points.InsertNextPoint(0.0, 0.0, 1.0);//8
points.InsertNextPoint(0.3, -0.5, 1.0);//9
points.InsertNextPoint(0.6, -.8, 1.0);//10
points.InsertNextPoint(0.7, -0.9, 1.0);//11
points.InsertNextPoint(0.8, -0.8, 1.0);//12
points.InsertNextPoint(1.0, 0.0, 1.0);//13
points.InsertNextPoint(1.0, 1.0, 1.0);//14
points.InsertNextPoint(0.0, 1.0, 1.0);//15
vtkUnstructuredGrid ugrid = vtkUnstructuredGrid.New();
int[] ptIds; IntPtr ptIdsPointer;
//8 quad faces
ptIds = new int[] { 0, 1, 9, 8 }; ptIdsPointer = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(int)) * ptIds.Length); Marshal.Copy(ptIds, 0, ptIdsPointer, ptIds.Length); ugrid.InsertNextCell(9, ptIds.Length, ptIdsPointer); Marshal.FreeHGlobal(ptIdsPointer);
ptIds = new int[] { 1, 2, 10, 9 }; ptIdsPointer = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(int)) * ptIds.Length); Marshal.Copy(ptIds, 0, ptIdsPointer, ptIds.Length); ugrid.InsertNextCell(9, ptIds.Length, ptIdsPointer); Marshal.FreeHGlobal(ptIdsPointer);
ptIds = new int[] { 2, 3, 11, 10 }; ptIdsPointer = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(int)) * ptIds.Length); Marshal.Copy(ptIds, 0, ptIdsPointer, ptIds.Length); ugrid.InsertNextCell(9, ptIds.Length, ptIdsPointer); Marshal.FreeHGlobal(ptIdsPointer);
ptIds = new int[] { 3, 4, 12, 11 }; ptIdsPointer = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(int)) * ptIds.Length); Marshal.Copy(ptIds, 0, ptIdsPointer, ptIds.Length); ugrid.InsertNextCell(9, ptIds.Length, ptIdsPointer); Marshal.FreeHGlobal(ptIdsPointer);
ptIds = new int[] { 4, 5, 13, 12 }; ptIdsPointer = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(int)) * ptIds.Length); Marshal.Copy(ptIds, 0, ptIdsPointer, ptIds.Length); ugrid.InsertNextCell(9, ptIds.Length, ptIdsPointer); Marshal.FreeHGlobal(ptIdsPointer);
ptIds = new int[] { 5, 6, 14, 13 }; ptIdsPointer = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(int)) * ptIds.Length); Marshal.Copy(ptIds, 0, ptIdsPointer, ptIds.Length); ugrid.InsertNextCell(9, ptIds.Length, ptIdsPointer); Marshal.FreeHGlobal(ptIdsPointer);
ptIds = new int[] { 6, 7, 15, 14 }; ptIdsPointer = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(int)) * ptIds.Length); Marshal.Copy(ptIds, 0, ptIdsPointer, ptIds.Length); ugrid.InsertNextCell(9, ptIds.Length, ptIdsPointer); Marshal.FreeHGlobal(ptIdsPointer);
ptIds = new int[] { 7, 0, 8, 15 }; ptIdsPointer = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(int)) * ptIds.Length); Marshal.Copy(ptIds, 0, ptIdsPointer, ptIds.Length); ugrid.InsertNextCell(9, ptIds.Length, ptIdsPointer); Marshal.FreeHGlobal(ptIdsPointer);
//2 octagonal faces
ptIds = new int[] { 0, 1, 2, 3, 4, 5, 6, 7 }; ptIdsPointer = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(int)) * ptIds.Length); Marshal.Copy(ptIds, 0, ptIdsPointer, ptIds.Length); ugrid.InsertNextCell(9, ptIds.Length, ptIdsPointer); Marshal.FreeHGlobal(ptIdsPointer);
ptIds = new int[] { 8, 9, 10, 11, 12, 13, 14, 15 }; ptIdsPointer = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(int)) * ptIds.Length); Marshal.Copy(ptIds, 0, ptIdsPointer, ptIds.Length); ugrid.InsertNextCell(9, ptIds.Length, ptIdsPointer); Marshal.FreeHGlobal(ptIdsPointer);
ugrid.SetPoints(points);
vtkGeometryFilter geometryFilter = vtkGeometryFilter.New();
geometryFilter.SetInput(ugrid);
vtkActor loActor = new vtkActor();
vtkPolyDataMapper pdm = vtkPolyDataMapper.New();
pdm.SetInputConnection(geometryFilter.GetOutputPort());
loActor.SetMapper(pdm);
loActor.GetProperty().SetPointSize(5);
loActor.GetProperty().SetColor(.9, 0, 0);
Renderer3D.AddActor(loActor);
vtkIdFilter ids = new vtkIdFilter();
ids.SetInputConnection(ugrid.GetProducerPort());
ids.PointIdsOn();
ids.CellIdsOn();
ids.FieldDataOn();
vtkLabeledDataMapper ldm = new vtkLabeledDataMapper();
ldm.SetInputConnection(ids.GetOutputPort());
ldm.SetLabelModeToLabelFieldData();
vtkActor2D pointLabels = new vtkActor2D();
这将导致:
内容总结
以上是互联网集市为您收集整理的c# – 如何创建vtk七角棱镜(7个面有4个顶点,2个面有7个顶点)全部内容,希望文章能够帮你解决c# – 如何创建vtk七角棱镜(7个面有4个顶点,2个面有7个顶点)所遇到的程序开发问题。 如果觉得互联网集市技术教程内容还不错,欢迎将互联网集市网站推荐给程序员好友。
内容备注
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 gblab@vip.qq.com 举报,一经查实,本站将立刻删除。
内容手机端
扫描二维码推送至手机访问。