jjzjj

c# - 如何创建vtk七角棱柱(7个面4个顶点和2个面7个顶点)

coder 2024-02-20 原文

vtk 中,您可以使用 InsertNextCell (int type, vtkIdType npts, vtkIdType *ptIds) 将六角棱柱单元格添加到 vtkUnstructuredGrid >。这很好,但我想在 vtkUnstructuredGrid 中包含一个 heptagonaloctagonal 棱镜。单元格类型未在 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);

到(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...

我怎样才能得到像图像这样的东西:

最佳答案

使用 vtkLinearExtrusionFilter 更容易。我写了一些 C++ 代码来生成具有可变边的多边形棱镜。我希望这会有所帮助:

    vtkSmartPointer<vtkVectorText> vecText = 
vtkSmartPointer<vtkVectorText>::New();
vecText->SetText("Text!");
vecText->Update();
vtkSmartPointer<vtkPoints> points =
vtkSmartPointer<vtkPoints>::New();
double Radius=1;
int NumberOfSides=7;
for (int i=0;i<NumberOfSides;i++)
    points->InsertNextPoint(Radius*cos((double)i*2*vtkMath::Pi()/NumberOfSides), Radius*sin((double)i*2*vtkMath::Pi()/NumberOfSides), 0.0); 
// Create the polygon
vtkSmartPointer<vtkPolygon> polygon =
vtkSmartPointer<vtkPolygon>::New();
polygon->GetPointIds()->SetNumberOfIds(NumberOfSides); 

for(int i=0;i<NumberOfSides;i++)
    polygon->GetPointIds()->SetId(i, i);

// Add the polygon to a list of polygons
vtkSmartPointer<vtkCellArray> polygons =
vtkSmartPointer<vtkCellArray>::New();
polygons->InsertNextCell(polygon);

// Create a PolyData
vtkSmartPointer<vtkPolyData> polygonPolyData =
vtkSmartPointer<vtkPolyData>::New();
polygonPolyData->SetPoints(points);
polygonPolyData->SetPolys(polygons);

// Apply linear extrusion 
vtkSmartPointer<vtkLinearExtrusionFilter> extrude = 
vtkSmartPointer<vtkLinearExtrusionFilter>::New();
extrude->SetInputData( polygonPolyData);
extrude->SetExtrusionTypeToNormalExtrusion();
extrude->SetVector(0, 0, 1 );
extrude->SetScaleFactor (2);
extrude->Update(); 
vtkSmartPointer<vtkRenderWindow> renderWindow = 
vtkSmartPointer<vtkRenderWindow>::New();
vtkSmartPointer<vtkPolyDataMapper> mapper=vtkSmartPointer<vtkPolyDataMapper>::New();
vtkSmartPointer<vtkActor> actor=vtkSmartPointer<vtkActor>::New();
mapper->SetInputData( extrude->GetOutput());
actor->SetMapper(mapper);
vtkSmartPointer<vtkRenderer> renderer=vtkSmartPointer<vtkRenderer>::New();
renderer->AddActor(actor);
renderWindow->AddRenderer(renderer);

vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor =  vtkSmartPointer<vtkRenderWindowInteractor>::New();
renderWindow->SetInteractor( renderWindowInteractor);
renderWindow->Render(); 

我运行了代码并得到了这个七棱柱的 3d 形状:

关于c# - 如何创建vtk七角棱柱(7个面4个顶点和2个面7个顶点),我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/35932279/

有关c# - 如何创建vtk七角棱柱(7个面4个顶点和2个面7个顶点)的更多相关文章

  1. ruby - 如何使用 Nokogiri 的 xpath 和 at_xpath 方法 - 2

    我正在学习如何使用Nokogiri,根据这段代码我遇到了一些问题:require'rubygems'require'mechanize'post_agent=WWW::Mechanize.newpost_page=post_agent.get('http://www.vbulletin.org/forum/showthread.php?t=230708')puts"\nabsolutepathwithtbodygivesnil"putspost_page.parser.xpath('/html/body/div/div/div/div/div/table/tbody/tr/td/div

  2. ruby - 如何从 ruby​​ 中的字符串运行任意对象方法? - 2

    总的来说,我对ruby​​还比较陌生,我正在为我正在创建的对象编写一些rspec测试用例。许多测试用例都非常基础,我只是想确保正确填充和返回值。我想知道是否有办法使用循环结构来执行此操作。不必为我要测试的每个方法都设置一个assertEquals。例如:describeitem,"TestingtheItem"doit"willhaveanullvaluetostart"doitem=Item.new#HereIcoulddotheitem.name.shouldbe_nil#thenIcoulddoitem.category.shouldbe_nilendend但我想要一些方法来使用

  3. ruby - 如何在 Ruby 中顺序创建 PI - 2

    出于纯粹的兴趣,我很好奇如何按顺序创建PI,而不是在过程结果之后生成数字,而是让数字在过程本身生成时显示。如果是这种情况,那么数字可以自行产生,我可以对以前看到的数字实现垃圾收集,从而创建一个无限系列。结果只是在Pi系列之后每秒生成一个数字。这是我通过互联网筛选的结果:这是流行的计算机友好算法,类机器算法:defarccot(x,unity)xpow=unity/xn=1sign=1sum=0loopdoterm=xpow/nbreakifterm==0sum+=sign*(xpow/n)xpow/=x*xn+=2sign=-signendsumenddefcalc_pi(digits

  4. python - 如何使用 Ruby 或 Python 创建一系列高音调和低音调的蜂鸣声? - 2

    关闭。这个问题是opinion-based.它目前不接受答案。想要改进这个问题?更新问题,以便editingthispost可以用事实和引用来回答它.关闭4年前。Improvethisquestion我想在固定时间创建一系列低音和高音调的哔哔声。例如:在150毫秒时发出高音调的蜂鸣声在151毫秒时发出低音调的蜂鸣声200毫秒时发出低音调的蜂鸣声250毫秒的高音调蜂鸣声有没有办法在Ruby或Python中做到这一点?我真的不在乎输出编码是什么(.wav、.mp3、.ogg等等),但我确实想创建一个输出文件。

  5. ruby-on-rails - 如何验证 update_all 是否实际在 Rails 中更新 - 2

    给定这段代码defcreate@upgrades=User.update_all(["role=?","upgraded"],:id=>params[:upgrade])redirect_toadmin_upgrades_path,:notice=>"Successfullyupgradeduser."end我如何在该操作中实际验证它们是否已保存或未重定向到适当的页面和消息? 最佳答案 在Rails3中,update_all不返回任何有意义的信息,除了已更新的记录数(这可能取决于您的DBMS是否返回该信息)。http://ar.ru

  6. ruby-on-rails - 'compass watch' 是如何工作的/它是如何与 rails 一起使用的 - 2

    我在我的项目目录中完成了compasscreate.和compassinitrails。几个问题:我已将我的.sass文件放在public/stylesheets中。这是放置它们的正确位置吗?当我运行compasswatch时,它不会自动编译这些.sass文件。我必须手动指定文件:compasswatchpublic/stylesheets/myfile.sass等。如何让它自动运行?文件ie.css、print.css和screen.css已放在stylesheets/compiled。如何在编译后不让它们重新出现的情况下删除它们?我自己编译的.sass文件编译成compiled/t

  7. ruby - 如何将脚本文件的末尾读取为数据文件(Perl 或任何其他语言) - 2

    我正在寻找执行以下操作的正确语法(在Perl、Shell或Ruby中):#variabletoaccessthedatalinesappendedasafileEND_OF_SCRIPT_MARKERrawdatastartshereanditcontinues. 最佳答案 Perl用__DATA__做这个:#!/usr/bin/perlusestrict;usewarnings;while(){print;}__DATA__Texttoprintgoeshere 关于ruby-如何将脚

  8. ruby - 如何指定 Rack 处理程序 - 2

    Rackup通过Rack的默认处理程序成功运行任何Rack应用程序。例如:classRackAppdefcall(environment)['200',{'Content-Type'=>'text/html'},["Helloworld"]]endendrunRackApp.new但是当最后一行更改为使用Rack的内置CGI处理程序时,rackup给出“NoMethodErrorat/undefinedmethod`call'fornil:NilClass”:Rack::Handler::CGI.runRackApp.newRack的其他内置处理程序也提出了同样的反对意见。例如Rack

  9. ruby - 使用 Vim Rails,您可以创建一个新的迁移文件并一次性打开它吗? - 2

    使用带有Rails插件的vim,您可以创建一个迁移文件,然后一次性打开该文件吗?textmate也可以这样吗? 最佳答案 你可以使用rails.vim然后做类似的事情::Rgeneratemigratonadd_foo_to_bar插件将打开迁移生成的文件,这正是您想要的。我不能代表textmate。 关于ruby-使用VimRails,您可以创建一个新的迁移文件并一次性打开它吗?,我们在StackOverflow上找到一个类似的问题: https://sta

  10. ruby - 如何每月在 Heroku 运行一次 Scheduler 插件? - 2

    在选择我想要运行操作的频率时,唯一的选项是“每天”、“每小时”和“每10分钟”。谢谢!我想为我的Rails3.1应用程序运行调度程序。 最佳答案 这不是一个优雅的解决方案,但您可以安排它每天运行,并在实际开始工作之前检查日期是否为当月的第一天。 关于ruby-如何每月在Heroku运行一次Scheduler插件?,我们在StackOverflow上找到一个类似的问题: https://stackoverflow.com/questions/8692687/

随机推荐