青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

eryar

PipeCAD - Plant Piping Design Software.
RvmTranslator - Translate AVEVA RVM to OBJ, glTF, etc.
posts - 603, comments - 590, trackbacks - 0, articles - 0

OpenCascade中的Delaunay三角剖分

Posted on 2013-05-26 16:51 eryar 閱讀(28190) 評論(17)  編輯 收藏 引用 所屬分類: 2.OpenCASCADE

Delaunay Triangulation in OpenCascade

eryar@163.com

摘要:本文簡要介紹了Delaunay三角剖分的基礎理論,并使用OpenCascade的三角剖分算法將邊界BRep表示的幾何體進行三角離散化后在OpenSceneGraph中顯示。

關鍵字:Delaunay Triangulation、OpenCascade、OpenSceneGraph

一、 概述

三角剖分是平面剖分中的一個重要課題,在數字圖像處理、計算機三維曲面造型、有限元計算、逆向工程等領域有著廣泛應用。由于三角形是平面域中的單純形,與其他平面圖形相比,其有描述方便、處理簡單等特性,很適合于對復雜區域進行簡化處理。因此,無論在計算幾何、計算機圖形處理、模式識別、曲面逼近,還有有限元網格生成方面有廣泛的應用。

雖然曲線、曲面等有精確的方程來表示,但是在在計算機中,只能用離散的方式來逼近。如曲線可用直線段來逼近,而曲面可用多邊形或三角形來表示。用多邊形網格表示曲面是設計中經常使用的形式,可以根據應用要求選擇網格的密度。利用三角形面片表示的曲面在計算機圖形學中也稱為三角形網格。用三角形網格表示曲面需要解決幾個問題:三角形的產生、描述、遍歷、簡化和壓縮等,這些問題都是計算幾何研究的范疇,相關問題都可以從中找到答案。下圖所示的圓柱和立方體是由OpenCascade生成,使用OpenCascade的算法離散成三角網格后在OpenSceneGraph中顯示的效果。

wps_clip_image-27228

Figure 1.1 Shaded Cylinder and Box

wps_clip_image-22196

Figure 1.2 Mesh generated by OpenCascade

從圖中可以看出,平面的三角形網格效果還不錯,曲面的三角形網格表示只能是近似表示,可以通過提高網格的密度來增加真實性,但相應渲染的數據量就大了。有人說OpenCascade的顯示模塊做得不是很好,上述方法則可以只使用OpenCascade的造型模塊,再結合OpenSceneGraph來對圖形進行顯示。

三維數據交換STL格式文件中保存的都是三角面片的數據,STL文件格式是由美國3D System公司開發,已被工業界認為是目前快速自動成型領域的準標準零件描述文件格式。它對三維實體描述的解釋具有惟一性。幾乎所有的幾何造型系統都提供STL文件數據交換接口。OpenCascade中的數據交換模塊也提供對STL格式的支持,由此可見三角網格在幾何造型系統中的重要性。

Voronoi圖和Delaunay三角剖分的應用領域十分廣泛:幾何建模——用來尋找三維曲面“好的”三角剖分;有限元分析——用來生成“好的”有限元網格;地理信息系統——用來進行空間領域分析;結晶學——用來確定合金的結構;人類學和考古學——用來確定氏族部落、首領權威、居住中心或堡壘等的影響范圍;天文學——用來確定恒星和星系的分布;生物學生態學和林學——用來確定動植物的競爭;動物學——分析動物的領地;統計學和數據分析——用來分析統計聚合;機器人學——用來進行運動軌跡規劃(在存在障礙物的情況下);模式識別——作為尋找物體骨架點的工具;生理學——用來分析毛細作用的領域;氣象學——用來估計區域平均降雨量;市場學——用來建立城市的市場輻射范圍;以及在遙感圖像處理、化學、地理學、地質學、冶金學、數學等學科的應用等。

本文只對OpenCascade中的三角剖分進行簡要介紹,希望對三角剖分在三維幾何造型方面有興趣的朋友可以對其深入研究。水平很有限,文中不當之處歡迎批評指正、指導,聯系郵箱:eryar@163.com

二、 Voronoi圖

Dirichlet于1850年研究了平面點的鄰域問題,Voronoi于1908年將其結果擴展到高維空間。半空間定義Voronoi圖:給定平面上n個點集S,S={p1, p2, …, pn}。定義:

wps_clip_image-1602

PiPj連線的垂直平分面將空間分為兩半,V(Pi)表示比其他點更接近Pi的點的軌跡是n-1個半平面的交,它是一個不多于n-1條邊的凸多邊形域,稱為關聯于Pi的Voronoi多邊形或關聯于Pi的Voronoi域。如下圖所示為關聯于P1的Voronoi多邊形,它是一個四邊形,而n=6.

wps_clip_image-4421

Figure 2.1 n=6時的一種V(p1)

對于點集S中的每個點都可以做一個Voronoi多邊形,這樣的n個Voronoi多邊形組成的圖稱為Voronoi圖,記為Vor(S)。如下圖所示:

wps_clip_image-31980

Figure 2.2 Voronoi diagram for 10 randomly points (Generated by MATLAB)

圖中的頂點和邊分別稱為Voronoi頂點和Voronoi邊。顯然,|S|=n時,Vor(S)劃分平面成n個多邊形域,每個多邊形域V(Pi)包含S中的一個點而且只包含S中的一個點,Vor(S)的邊是S中某點對的垂直平分線上的一條線段或半直線,從而為該點對所在的兩個多邊形域所共有。Vor(S)中有的多邊形域是無界的。

wps_clip_image-18135

Figure 2.3 Ten shops in a flat city and their Voronoi cells

(http://en.wikipedia.org/wiki/Voronoi_diagram)

wps_clip_image-27044

Figure 2.4 Voronoi tessellation in a cylinder (Voro++ library: http://math.lbl.gov/voro++/)

Voronoi圖有如下性質:

l n個點的點集S的Voronoi圖至多有2n-5個頂點和3n-6條邊;

l 每個Voronoi點恰好是三條Voronoi邊的交點;

l 設v是Vor(S)的頂點,則圓C(v)內不含S的其他點;

l 點集S中點Pi的每一個最近鄰近點確定V(Pi)的一條邊;

l Voronoi圖的直線對偶圖是S的一個三角剖分;

l 如果Pi,Pj屬于S,并且通過Pi,Pj有一個不包含S中其他點的圓,那么線段PiPj是點集S三角剖分的一條邊,反之亦成立。

三、 Delaunay三角剖分 

1. 二維實數域上的三角剖分

假設V是二維實數域上的有限點集,邊e是由點集中的點作端點構成的封閉線段,E為e的集合,那么該點集V的一個三角剖分T=(V,E)是一個平面圖:

l 除了端點,平面圖中的邊不包含點集中的任何點;

l 沒有相交邊;

l 平面圖中所有的面都是三角面,且所有三角面的合集是點集V的凸包。

2. Delaunay邊

假設E中的一條邊(兩端點a,b),e滿足下列條件,則稱為Delaunay邊:存在一個圓經過a,b兩點,圓內不包含點集V中的任何的點。這一特性又稱為空圓特性。

3. Delaunay三角剖分

如果點集V的一個三角剖分T中只包含Delaunay邊,那么該三角剖分稱為Delaunay剖分。

最近點意義下的Voronoi圖的對偶圖實際上是點集的一種三角剖分,該三角剖分就是Delaunay剖分(表示為DT(S)),其中每個三角形的外接圓不包含點集中的其他任何點。因此,在構造點集的Voronoi圖之后,再作其對偶圖,即對每條Voronoi邊作通過點集中某兩點的垂直平分線,即得到Delaunay三角剖分。

wps_clip_image-18943

Figure 3.1 Delaunay Triangulation (Generated by MATLAB)

再看幾個圖片,加深對Delaunay三角剖分的理解:

wps_clip_image-12125

Figure 3.2 Delaunay Edge 

wps_clip_image-2535

Figure 3.3 Illustrate Delaunay Edge

wps_clip_image-16525

Figure 3.4 Delaunay Edge

4. Delaunay三角剖分的特性

l 1978年Sibson證明了在二維的情況下,在點集的所有三角剖分中,Delaunay三角剖分使得生成的三角形的最小角達到最大(max-min angle)。因為這一特性,對于給定點集的Delaunay三角剖分總是盡可能避免“瘦長”三角形,自動向等邊三角形逼近;

l 局部優化與整體優化(locally optimal and globally optimal);

l Delaunay空洞(cavity)與局部重連(local reconnection);

5. 經典的Delaunay三角剖分算法

目前常用的算法分為幾種,有掃描線法(Sweepline)、隨機增量法(Incremental)、分治法(Divide and Conquer)等。

經典的Delaunay三角剖分算法主要有兩類:Bowyer/Watson算法和局部變換法。

l Bowyer/Watson算法又稱為Delaunay空洞算法或加點法,以Bowyer和Watson算法為代表。從一個三角形開始,每次加一個點,保證每一步得到的當前三角形是局部優化的。以英國Bath大學數學分校Bowyer,Green,Sibson為代表的計算Dirichlet圖的方法屬于加點法,是較早成名的算法之一;以澳大利亞悉尼大學地學系Watson為代表的空外接球法也屬于加點法。加點法算法簡明,是目前應用最多的算法,該方法利用了Delaunay空洞性質。Bowyer/Watson算法的優點是與空間的維數無關,并且算法在實現上比局部變換算法簡單。該算法在新點加入到Delaunay網格時,部分外接球包含新點的三角形單元不再符合Delaunay屬性,則這些三角形單元被刪除,形成Delaunay空洞,然后算法將新點與組成空洞的每一個頂點相連生成一個新邊,根據空球屬性可以證明這些新邊都是局部Delaunay的,因此新生成的三角網格仍是Delaunay的。

wps_clip_image-22333

Figure 3.5 Illustration of 2D Bowyer/Watson algorithm for Delaunay Triangulation

l 局部變換法又稱為換邊、換面法。當利用局部變換法實現增量式點集的Delaunay三角剖分時,首先定位新加入點所在的三角形,然后在網格中加入三個新的連接該三角形頂點與新頂點的邊,若該新點位于某條邊上,則該邊被刪除,四條連接該新點的邊被加入。最后,在通過換邊方法對該新點的局部區域內的邊進行檢測和變換,重新維護網格的Delaunay性質。局部變換法的另一個優點是其可以對已存在的三角網格進行優化,使其變換成為Delaunay三角網格,該方法的缺點則是當算法擴展到高維空間時變得較為復雜。

四、 Delaunay三角剖分在OpenCascade的應用

OpenCascade中網格剖分的包主要有BRepMesh、MeshAlgo、MeshVS,其中,類MeshAlgo_Delaunay使用算法Watson來進行Delaunay三角剖分。從類StlTransfer中的注釋The triangulation is computed with the Delaunay algorithm implemented in package BRepMesh.可以看出包BRepMesh就是Delaunay三角剖分的具體實現。使用方法如下:

BRepMesh::Mesh (aShape, Deflection);

這個函數主要是用來對拓撲形狀進行三角剖分。以下通過將一個圓柱三角剖分為例說明如何將一個拓撲形狀進行三角剖分并將結果進行可視化。

 

/*

*    Copyright (c) 2013 eryar All Rights Reserved. 



*        File    : Main.cpp 

*        Author  : eryar@163.com 

*        Date    : 2013-05-26 

*        Version : 0.1 



*    Description : Use BRepMesh_Delaun class to learn  

*                  Delaunay's triangulation algorithm. 



*/
 

// Open Cascade library. 

#include 
<gp_Pnt.hxx> 

#include 
<gp_Pln.hxx> 

#include 
<BRep_Tool.hxx> 

#include 
<TopoDS.hxx> 

#include 
<TopoDS_Edge.hxx> 

#include 
<TopoDS_Wire.hxx> 

#include 
<TopoDS_Face.hxx> 

#include 
<BRepBuilderAPI_MakeEdge.hxx> 

#include 
<BRepBuilderAPI_MakeWire.hxx> 

#include 
<BRepBuilderAPI_MakeFace.hxx> 

#include 
<BRepPrimAPI_MakeBox.hxx> 

#include 
<BRepPrimAPI_MakeCone.hxx> 

#include 
<BRepPrimAPI_MakeCylinder.hxx> 

#include 
<BRepPrimApI_MakeSphere.hxx> 

#include 
<BRepMesh.hxx> 

#include 
<TopExp_Explorer.hxx> 

#include 
<Poly_Triangulation.hxx> 

#include 
<TShort_Array1OfShortReal.hxx> 

#pragma comment(lib, 
"TKernel.lib"

#pragma comment(lib, 
"TKMath.lib"

#pragma comment(lib, 
"TKBRep.lib"

#pragma comment(lib, 
"TKPrim.lib"

#pragma comment(lib, 
"TKMesh.lib"

#pragma comment(lib, 
"TKTopAlgo.lib"

// OpenSceneGraph library. 

#include 
<osgDB/ReadFile> 

#include 
<osgViewer/Viewer> 

#include 
<osgViewer/ViewerEventHandlers> 

#include 
<osgGA/StateSetManipulator> 

#pragma comment(lib, 
"osgd.lib"

#pragma comment(lib, 
"osgDbd.lib"

#pragma comment(lib, 
"osgGAd.lib"

#pragma comment(lib, 
"osgViewerd.lib"

osg::Node
* BuildShapeMesh(const TopoDS_Shape& aShape) 



    osg::ref_ptr
<osg::Group> root = new osg::Group(); 

    osg::ref_ptr
<osg::Geode> geode = new osg::Geode(); 

    osg::ref_ptr
<osg::Geometry> triGeom = new osg::Geometry(); 

    osg::ref_ptr
<osg::Vec3Array> vertices = new osg::Vec3Array(); 

    osg::ref_ptr
<osg::Vec3Array> normals = new osg::Vec3Array(); 

    BRepMesh::Mesh(aShape, 
1); 

    TopExp_Explorer faceExplorer; 

for (faceExplorer.Init(aShape, TopAbs_FACE); faceExplorer.More(); faceExplorer.Next()) 



        TopLoc_Location loc; 

        TopoDS_Face aFace 
= TopoDS::Face(faceExplorer.Current()); 

        Handle_Poly_Triangulation triFace 
= BRep_Tool::Triangulation(aFace, loc); 

        Standard_Integer nTriangles 
= triFace->NbTriangles(); 

        gp_Pnt vertex1; 

        gp_Pnt vertex2; 

        gp_Pnt vertex3; 

        Standard_Integer nVertexIndex1 
= 0

        Standard_Integer nVertexIndex2 
= 0

        Standard_Integer nVertexIndex3 
= 0

        TColgp_Array1OfPnt nodes(
1, triFace->NbNodes()); 

        Poly_Array1OfTriangle triangles(
1, triFace->NbTriangles()); 

        nodes 
= triFace->Nodes(); 

        triangles 
= triFace->Triangles(); 

for (Standard_Integer i = 1; i <= nTriangles; i++



            Poly_Triangle aTriangle 
= triangles.Value(i); 

            aTriangle.Get(nVertexIndex1, nVertexIndex2, nVertexIndex3); 

            vertex1 
= nodes.Value(nVertexIndex1); 

            vertex2 
= nodes.Value(nVertexIndex2); 

            vertex3 
= nodes.Value(nVertexIndex3); 

            gp_XYZ vector12(vertex2.XYZ() 
- vertex1.XYZ()); 

            gp_XYZ vector13(vertex3.XYZ() 
- vertex1.XYZ()); 

            gp_XYZ normal 
= vector12.Crossed(vector13); 

            Standard_Real rModulus 
= normal.Modulus(); 

if (rModulus > gp::Resolution()) 



                normal.Normalize(); 

}
 

else 



                normal.SetCoord(
0., 0., 0.); 

}
 

            vertices
->push_back(osg::Vec3(vertex1.X(), vertex1.Y(), vertex1.Z())); 

            vertices
->push_back(osg::Vec3(vertex2.X(), vertex2.Y(), vertex2.Z())); 

            vertices
->push_back(osg::Vec3(vertex3.X(), vertex3.Y(), vertex3.Z())); 

            normals
->push_back(osg::Vec3(normal.X(), normal.Y(), normal.Z())); 

}
 

}
 

    triGeom
->setVertexArray(vertices.get()); 

    triGeom
->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::TRIANGLES, 0, vertices->size())); 

    triGeom
->setNormalArray(normals); 

    triGeom
->setNormalBinding(osg::Geometry::BIND_PER_PRIMITIVE); 

    geode
->addDrawable(triGeom); 

    root
->addChild(geode); 

return root.release(); 

}
 

int main(int argc, char* argv[]) 



    osgViewer::Viewer myViewer; 

    osg::ref_ptr
<osg::Group> root = new osg::Group(); 

    root
->addChild(BuildShapeMesh(BRepPrimAPI_MakeCylinder(.61))); 

    myViewer.setSceneData(root); 

    myViewer.addEventHandler(
new osgGA::StateSetManipulator(myViewer.getCamera()->getOrCreateStateSet())); 

    myViewer.addEventHandler(
new osgViewer::StatsHandler); 

    myViewer.addEventHandler(
new osgViewer::WindowSizeHandler); 

return myViewer.run(); 

}
 


結果如下圖所示:

wps_clip_image-31478

Figure 4.1 Cylinder mesh generated by BRepMesh::Mesh

BRepMesh::Mesh是經過封裝的,便于對拓撲形狀進行三角剖分。以下通過一個簡單的例子來說明直接使用BRepMesh_Delaun的方法:

 

/*

*    Copyright (c) 2013 eryar All Rights Reserved. 



*        File    : Main.cpp 

*        Author  : eryar@163.com 

*        Date    : 2013-05-26 

*        Version : 0.1 



*    Description : Use BRepMesh_Delaun class to learn  

*                  Delaunay's triangulation algorithm. 



*/ 

#include 
<BRepMesh_Edge.hxx> 

#include 
<BRepMesh_Delaun.hxx> 

#include 
<BRepMesh_Array1OfVertexOfDelaun.hxx> 

#include 
<TColStd_MapIteratorOfMapOfInteger.hxx> 

#pragma comment(lib, 
"TKernel.lib"

#pragma comment(lib, 
"TKMesh.lib"

int main(int argc, char* argv[]) 



    BRepMesh_Array1OfVertexOfDelaun vertices(
14); 

    vertices.SetValue(
1, BRepMesh_Vertex(00, MeshDS_Free)); 

    vertices.SetValue(
2, BRepMesh_Vertex(10, MeshDS_Free)); 

    vertices.SetValue(
3, BRepMesh_Vertex(11, MeshDS_Free)); 

    vertices.SetValue(
4, BRepMesh_Vertex(01, MeshDS_Free)); 

    BRepMesh_Delaun triangulation(vertices); 

//triangulation.AddVertex(BRepMesh_Vertex(0.5, 0.5, MeshDS_OnSurface)); 

    Handle_BRepMesh_DataStructureOfDelaun meshData 
= triangulation.Result(); 

    std::cout
<<"Iterate Mesh Triangles:"<<std::endl; 

    MeshDS_MapOfInteger::Iterator triDom; 

for (triDom.Initialize(meshData->ElemOfDomain()); triDom.More(); triDom.Next()) 



        Standard_Integer triId 
= triDom.Key(); 

const BRepMesh_Triangle& curTri = meshData->GetElement(triId); 

        Standard_Integer vertexIndex1 
= 0

        Standard_Integer vertexIndex2 
= 0

        Standard_Integer vertexIndex3 
= 0

        Standard_Integer edgeIndex1 
= 0

        Standard_Integer edgeIndex2 
= 0

        Standard_Integer edgeIndex3 
= 0

        Standard_Boolean o1 
= Standard_False; 

        Standard_Boolean o2 
= Standard_False; 

        Standard_Boolean o3 
= Standard_False; 

        curTri.Edges(edgeIndex1, edgeIndex2, edgeIndex3, o1, o2, o3); 

const BRepMesh_Edge& edge1 = meshData->GetLink(edgeIndex1); 

const BRepMesh_Edge& edge2 = meshData->GetLink(edgeIndex2); 

const BRepMesh_Edge& edge3 = meshData->GetLink(edgeIndex3); 

        vertexIndex1 
= (o1? edge1.FirstNode(): edge1.LastNode()); 

        vertexIndex2 
= (o1? edge1.LastNode() : edge1.FirstNode()); 

        vertexIndex3 
= (o2? edge2.LastNode() : edge2.FirstNode()); 

const BRepMesh_Vertex& vertex1 = meshData->GetNode(vertexIndex1); 

const BRepMesh_Vertex& vertex2 = meshData->GetNode(vertexIndex2); 

const BRepMesh_Vertex& vertex3 = meshData->GetNode(vertexIndex3); 

const gp_XY& p1 = vertex1.Coord(); 

const gp_XY& p2 = vertex2.Coord(); 

const gp_XY& p3 = vertex3.Coord(); 

        std::cout
<<"--------"<<std::endl; 

        std::cout
<<p1.X()<<" , "<<p1.Y()<<std::endl; 

        std::cout
<<p2.X()<<" , "<<p2.Y()<<std::endl; 

        std::cout
<<p3.X()<<" , "<<p3.Y()<<std::endl; 

        std::cout
<<"========"<<std::endl; 



return 0




上述程序是以一個正方形為例,使用BRepMesh_Delaun三角剖分的結果為兩個三角形,如下所示:

 

Iterate Mesh Triangles: 

-------- 

1 , 1 

0 , 0 

1 , 0 

======== 

-------- 

1 , 1 

0 , 1 

0 , 0 

======== 

 

以上結果都是二維空間上的,三維空間中的使用方法可以參考類:BRepMesh_FastDiscretFace。這個類說明了如何將一個面進行網格劃分。

五、 結論

Delaunay三角剖分理論在三維幾何造型中還是比較重要的,通過對形狀的三角剖分,不僅可以對其進行可視化,還便于對形狀做進一步的處理,如消隱、光照處理等。通過對OpenCascade中三角剖分算法的使用,以進一步了解三角剖分理論應用及其算法實現。

六、 參考資料

1. 周培德. 計算幾何—算法設計與分析. 清華大學出版社, 2011

2. 李海生. Delaunay三角剖分理論及可視化應用研究. 哈爾濱工業大學出版社, 2010

3. 何援軍. 計算機圖形學. 機械工業出版社, 2010

4. 周元峰, 孫峰, 王文平, 汪嘉業, 張彩明. 基于局部修復的移動數據點Delaunay三角化快速更新方法. 計算機輔助設計與圖形學學報, 2011, 12: 2006-1012

5. http://en.wikipedia.org/wiki/Voronoi_diagram

 

Feedback

# re: OpenCascade中的Delaunay三角剖分  回復  更多評論   

2014-03-13 19:18 by OpenCASCADE->3D
稱呼您一聲E老師,一點都不為過。路漫漫其修遠兮,吾將上下而求索!

# re: OpenCascade中的Delaunay三角剖分  回復  更多評論   

2014-03-13 19:52 by eryar
@OpenCASCADE->3D
老師實不敢當。

后半句共勉之……

# re: OpenCascade中的Delaunay三角剖分  回復  更多評論   

2014-09-24 11:08 by 佚名
博主很厲害,怎么對occ研究這么透。如果occ能在國內普及。博主功不可沒!另請教博主一個問題,如何將一組散點組成nurbs曲面,假設這些點是從已有曲面上得到的,不僅僅包含邊界線上的點,還有面內部點。

# re: OpenCascade中的Delaunay三角剖分  回復  更多評論   

2014-09-24 17:49 by eryar
@佚名
過獎了。

我看OCC的例子程序中有個根據點生成面的功能,使用的類有:
1. Plate_Plate;
2. GeomPlate_Surface;
3. GeomPlate_MakeApprox;

相關實現可以參考示例程序。

# re: OpenCascade中的Delaunay三角剖分  回復  更多評論   

2014-09-24 22:10 by 佚名
@eryar
謝謝您的回復,我好好看看這些算例。

# re: OpenCascade中的Delaunay三角剖分  回復  更多評論   

2014-09-24 22:33 by eryar
@佚名
不客氣。。。

# re: OpenCascade中的Delaunay三角剖分  回復  更多評論   

2014-09-25 17:22 by 佚名
博主您好!
按您的指點我看了occ附帶的例子11,它確實可以由離散點組成樣條曲面,美中不足的是曲面的一部分延伸出了邊界線(由離散的外圍點組成),試著通過調整參數將面限制在點的范圍內,然而毫無效果. 是否有更好的方法,如果我先將這些點連成delaunay三角面,能否由這些面組成樣條曲面?謝謝!


# re: OpenCascade中的Delaunay三角剖分  回復  更多評論   

2014-09-25 17:39 by eryar
@佚名

你好!

擬合曲面的算法還沒有細看。

不過好像OCC有收費的模塊,專門用來從點云數據生成幾何數據,可以看看:
Surfaces from Scattered Points:
http://www.opencascade.org/support/products/ssp/

# re: OpenCascade中的Delaunay三角剖分  回復  更多評論   

2014-09-25 19:49 by 佚名
非常感謝

# re: OpenCascade中的Delaunay三角剖分  回復  更多評論   

2014-11-25 10:52 by xmzzy
請教個問題:
在使用BRepMesh::Mesh (aShape, Deflection);時,編譯提示“Mesh”: 不是“BRepMesh”的成員。
這個怎么處理呢?

# re: OpenCascade中的Delaunay三角剖分  回復  更多評論   

2014-11-25 12:43 by eryar
@xmzzy
你可能用的occ6.8.0,這個版本中去掉BRepMesh::Mesh()函數了,官方給的答復是這個函數只是調用了BRepMesh_IncrementalMesh來生成網格,沒有存在的意義,所以予以去除。
詳細可參考官方論壇:
http://www.opencascade.org/org/forum/thread_26649/?forum=3

6.8.0: BRepMesh::Mesh removed
去除原因及新的用法都有提及,原文如下:
Forum supervisor 2014/11/13 17:30
Dear Thorsen,

Method BRepMesh::Mesh did nothing more than called BRepMesh_IncrementalMesh with the given parameter, so it had no significant meaning and as result it was removed.
To perform meshing user should use BRepMesh_IncrementalMesh directly, for instance:

#include<BRepMesh_IncrementalMesh.hxx>
...
BRepMesh_IncrementalMesh(aShape, 0.001);
...

It is exactly the same what BRepMesh::Mesh did.

Best regards
FSR

# re: OpenCascade中的Delaunay三角剖分  回復  更多評論   

2014-11-25 13:29 by xmzzy
@eryar
十分感謝。大神啊。

# re: OpenCascade中的Delaunay三角剖分  回復  更多評論   

2014-12-11 11:26 by zhx
您好,樓主!我想實現對點集的DELAUNAY三角剖分,點集在三維空間的一個平面上(該平面時任意位置的平面),點可以獲取為gp_pnt類型,請問能不能用BRepMesh_Delaun實現對該點集的剖分呢?
另外,請問下BRepMesh_Array1OfVertexOfDelaun類型的賦值怎么由GP_PNT類的值實現呢 非常感謝樓主!

# re: OpenCascade中的Delaunay三角剖分  回復  更多評論   

2014-12-11 14:52 by zhx
@xmzzy
你好啊 我也是在用OCC,之前用的是OCC6.5.0,現在想換位OCC6.8.0,可是安裝完之后發現6.8.0版本的網頁版的OCC文檔不能查找函數,是為什么呢?

# re: OpenCascade中的Delaunay三角剖分  回復  更多評論   

2014-12-11 17:33 by eryar
@zhx
Hi,

你好!

OpenCASCADE的mesh算法只是針對二維的點集(即參數空間中的點),主要是為了替換掉6.3版本以前用的Triangle庫。

如果你是三維的點且可確定都在一個平面上,可以考慮將這些點投影到這個平面,Delaunay三角剖分后再映射回三維空間。

# re: OpenCascade中的Delaunay三角剖分  回復  更多評論   

2014-12-15 22:25 by eryar
@zhx
你好啊 我也是在用OCC,之前用的是OCC6.5.0,現在想換位OCC6.8.0,可是安裝完之后發現6.8.0版本的網頁版的OCC文檔不能查找函數,是為什么呢?
=================================================

我也試了下,occ6.8.0開發文檔中的search功能的確不能用。
看了下,它是用php來調用javascript的,所以應該需要放到支持php的服務器上運行這個文檔就可以搜索了。

我在occ官網上search是可以的,如下為在occ官網上search類gp_Pnt:
http://dev.opencascade.org/doc/refman/html/search.php?query=gp_Pnt

# re: OpenCascade中的Delaunay三角剖分  回復  更多評論   

2016-03-12 14:33 by loveabcc
@佚名
你好:不知道你能收到我的回復不,我想實現的是和你是一樣的東西,就是點云轉換成Nurbs,你說的OCC例子11是哪個呢,還有現在你找到OCC的方法了嗎,如果可以分享下。謝謝兄弟
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            亚洲国产精品久久精品怡红院| 亚洲一区二区三区四区五区黄| 欧美在线看片| 亚洲一区二区三区午夜| 国产精品视频第一区| 午夜视频在线观看一区二区| 亚洲尤物在线| 亚洲大胆在线| 亚洲人成人99网站| 国产精品www.| 美腿丝袜亚洲色图| 欧美精品日韩精品| 欧美一区二区三区在线看| 欧美中在线观看| 亚洲理论在线观看| 亚洲免费人成在线视频观看| 精品成人一区二区三区| 亚洲精品三级| 韩国av一区| 亚洲久久视频| 黄色一区二区在线| 亚洲精选在线观看| 一区二区三区在线看| 日韩一二三区视频| 精品91在线| 一区二区三区免费观看| 悠悠资源网亚洲青| 亚洲一区二区三区视频| 亚洲精品国久久99热| 亚洲女与黑人做爰| 99v久久综合狠狠综合久久| 午夜一区在线| 中文在线不卡视频| 蜜桃久久av一区| 欧美一区二区三区免费观看视频 | 一区二区三区四区五区视频| 午夜精品网站| 亚洲天堂男人| 欧美激情视频给我| 麻豆freexxxx性91精品| 国产精品一区久久| 夜色激情一区二区| 亚洲伦伦在线| 美女免费视频一区| 久久婷婷丁香| 国产小视频国产精品| 一区二区三区精品视频在线观看| 亚洲精品社区| 久久天堂国产精品| 久久久久亚洲综合| 国产日韩精品视频一区| 亚洲午夜伦理| 亚洲一区二区三区影院| 欧美精品一级| 亚洲精品在线一区二区| 日韩午夜电影| 欧美激情按摩在线| 亚洲国产导航| 亚洲日本欧美| 欧美黑人国产人伦爽爽爽| 欧美成人一区二区三区| 在线观看欧美日韩| 免费看亚洲片| 91久久精品国产91性色| 亚洲国产精品ⅴa在线观看| 久久这里只有| 亚洲第一区在线| 日韩视频不卡中文| 欧美日韩三级一区二区| 夜夜嗨网站十八久久| 亚洲淫片在线视频| 国产精品一区二区你懂的| 午夜精品福利在线| 久久嫩草精品久久久精品| 精品成人一区二区三区四区| 久久亚洲精品一区二区| 亚洲国产成人在线| 亚洲视频欧美视频| 国产美女扒开尿口久久久| 午夜在线不卡| 美女免费视频一区| 亚洲精品小视频| 欧美日韩综合视频| 午夜精品成人在线| 亚洲成人在线网站| 亚洲一区二区三区精品动漫| 国产精品自拍一区| 久久综合一区| 一区二区三区高清| 久久性色av| 亚洲天堂网在线观看| 国产欧美综合在线| 免费人成精品欧美精品| 一区二区av| 免费不卡在线观看| 亚洲视频综合在线| 精品不卡一区二区三区| 欧美日韩一区二区高清| 欧美一区二区三区免费视| 欧美激情免费观看| 欧美亚洲色图校园春色| 亚洲成人资源| 国产老女人精品毛片久久| 久久综合九色| 亚洲欧美综合精品久久成人| 欧美国产一区二区在线观看| 午夜精品理论片| 亚洲欧洲日产国产网站| 国产欧美一区二区三区国产幕精品| 另类专区欧美制服同性| 亚洲欧美综合网| 一本久久精品一区二区| 欧美激情精品久久久久| 久久国内精品自在自线400部| 最新日韩在线视频| 国内揄拍国内精品少妇国语| 欧美亚洲不卡| 欧美区国产区| 免费看成人av| 蜜桃久久精品乱码一区二区| 亚洲欧美色一区| 一级成人国产| 99国产一区| 亚洲国产精品小视频| 麻豆精品在线视频| 久久免费99精品久久久久久| 午夜精品久久久久久久久久久久| 亚洲精品一区在线| 亚洲激情综合| 亚洲国产欧美在线人成| 狠狠做深爱婷婷久久综合一区| 国产精品人人做人人爽| 欧美色欧美亚洲另类七区| 欧美精品啪啪| 欧美日韩免费观看一区二区三区 | 亚洲一区激情| 夜夜狂射影院欧美极品| 亚洲精品一区二区三| 亚洲国产欧美国产综合一区 | 亚洲欧美国产制服动漫| 亚洲一区欧美二区| 亚洲一区二区三区乱码aⅴ蜜桃女| 99综合精品| 正在播放亚洲| 亚洲综合视频网| 午夜视频一区在线观看| 久久国产精品第一页| 久久精品91久久久久久再现| 久久国产精品色婷婷| 久久久噜噜噜久久中文字幕色伊伊| 欧美中文字幕在线播放| 久久久噜噜噜久久中文字免| 免费成人黄色av| 欧美黄色免费网站| 亚洲三级色网| 亚洲特级片在线| 欧美一级一区| 欧美国产一区二区| 欧美日韩在线影院| 国产日韩av在线播放| 娇妻被交换粗又大又硬视频欧美| 亚洲人成亚洲人成在线观看| 99re热精品| 久久精品国产亚洲精品| 欧美国产一区视频在线观看| 一本一本大道香蕉久在线精品| 亚洲专区欧美专区| 美日韩精品免费| 欧美四级在线观看| 国产亚洲精品久久久久婷婷瑜伽| 极品尤物一区二区三区| 亚洲天堂av在线免费观看| 久久精品国产清高在天天线 | 久久久久久亚洲精品杨幂换脸| 蘑菇福利视频一区播放| 99视频国产精品免费观看| 香蕉国产精品偷在线观看不卡| 老司机免费视频一区二区| 国产精品v日韩精品| 在线激情影院一区| 亚洲一区二区在线| 欧美sm视频| 亚洲一区二区三区精品在线| 久久一日本道色综合久久| 国产精品www.| 日韩亚洲视频在线| 久久免费精品视频| 亚洲天堂第二页| 欧美成人蜜桃| 狠狠久久婷婷| 午夜一区不卡| 9国产精品视频| 免费日韩av| 国内精品视频一区| 亚洲女性裸体视频| 亚洲黄色三级| 美女视频一区免费观看| 国产一区二区精品在线观看| 亚洲主播在线播放| 亚洲人成网站色ww在线|