Najszybsze biblioteki triangulacji Delaunaya dla zestawów punktów 3D


26

Która jest najszybszą biblioteką do przeprowadzania delangijskiej triangulacji zbiorów z milionami punktów 3D? Czy dostępne są również wersje GPU? Z drugiej strony, mając teselację voronoi tego samego zestawu punktów, pomógłby (pod względem wydajności) uzyskać triangulację delaunay?


Czy próbowałeś już implementacji w oprogramowaniu CGAL i TRIANGLE? Oba zawierają algorytmy , które są (teoretycznie) najszybsze dostępne (choć nie równolegle). O(nlog(n))
Paweł

Jonathan Shewchuk ma również wersję strumieniową dla 2D, która może obsłużyć absurdalnie duże zbiory danych, jeśli możesz dodać dodatkowe dane do swojego strumienia. Do czego tego używasz?
Victor Liu

Odpowiedzi:


13

Do obliczania trójwymiarowych triangulacji Delaunaya (tak naprawdę tetraedralizacji), TetGen jest powszechnie stosowaną biblioteką.

Dla Twojej wygody poniżej przedstawiamy, jak długo trzeba obliczać terehedralizację wielu losowych punktów z kostki jednostkowej. Za 100 000 punktów stary Pentium M. zajmuje 4,5 sekundy

Grafika matematyczna

(Dokonano tego za pomocą interfejsu TetGen Mathematiki. Nie wiem, ile to narzuca.)

Jeśli chodzi o twoje drugie pytanie: jeśli już masz teselację Voronoi, to uzyskanie triangulacji Delaunaya jest stosunkowo prostą transformacją .


10

gStar4D to szybki i niezawodny algorytm 3D Delaunay dla GPU. Jest implementowany przy użyciu CUDA i działa na procesorach graficznych NVIDIA.

Algorytm ten, podobnie jak GPU-DT , konstruuje najpierw cyfrowy schemat Voronoi 3D. Jednak w 3D nie można go sprowadzić do triangulacji z powodu problemów topologicznych i geometrycznych. Zamiast tego gStar4D wykorzystuje informacje o sąsiedztwie z tego diagramu do tworzenia gwiazd podniesionych do 4D i wykonuje na nich efektywne rozrzucanie gwiazd na GPU. Po wydobyciu z tego kadłuba uzyskuje się triangulację 3D Delaunay.

Najszybszą implementacją 3D Delaunay jest gDel3D , który jest hybrydowym algorytmem GPU-CPU.

Wykonuje równoległe wstawianie i przerzucanie na GPU. Wynik jest zbliżony do Delaunaya. Następnie naprawia ten wynik przy użyciu konserwatywnej metody rozrzucania gwiazd na procesorze.

Obie te metody są solidne, więc mogą obsługiwać dowolny rodzaj zdegenerowanych danych wejściowych. Mogą obsłużyć miliony punktów, jeśli pamięć GPU jest wystarczająco duża, aby pomieścić pośrednie struktury danych.

Ujawnienie: Jestem autorem tych algorytmów i implementacji :)


Witamy w SciComp.Se, Ashwin! W celu pełnego ujawnienia informacji należy dodać, że jesteś autorem tego oprogramowania (patrz: meta.scicomp.stackexchange.com/a/342/1804 ).
Christian Clason

3

Poleciłbym wypróbowanie CGAL http://www.cgal.org/Manual/latest/doc_html/cgal_manual/Triangulation_3/Chapter_main.html#Section_39.2 , jak zasugerował powyżej Paul. CGAL to solidna i dobrze obsługiwana biblioteka, która istnieje już od dłuższego czasu. Korzystałem z tego z radością w przeszłości, nawet na zestawach punktów z punktami współliniowymi i współpłaszczyznowymi. Nie wiem, czy dziś jest najszybszy, ale z pewnością jest to dobre miejsce, aby zacząć.

Powyższy link zawiera również niektóre wyniki: może zarobić milion punktów w około 10 sekund i 10 milionów w około 1,5 minuty.


Czy mógłbyś też wyjaśnić, dlaczego go polecasz? Czy masz z tym doświadczenie?
Godric Seer,

1

Jeśli masz już schemat voronoi zbioru punktów, to obliczenie triangulacji Delaunaya zajmie ci tylko O ​​(n). Odpowiednio, biorąc pod uwagę punkt voronoi, możesz uzyskać jego trójkąt Delaunay w O (1). Są podwójne, więc staraj się wykorzystywać tę sytuację, kiedy tylko jest to możliwe.


1

Możesz wypróbować oprogramowanie geograficzne, które rozwijam: http://alice.loria.fr/software/geogram/doc/html/index.html

Ma równoległy algorytm, który oblicza ID 14 milionów wierzchołków w mniej niż 19 sekund na Intel Core I7 (dla 1 miliona wierzchołków zajmuje to około 0,8 s)


Witamy w SciComp.SE! W trosce o pełne ujawnienie (i pokazanie, że wiesz, o czym mówisz), powinieneś wspomnieć, że jesteś programistą tego oprogramowania.
Christian Clason,
Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.