Eine Delaunay-Triangulation filtern

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Tanja86

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 10.06.2009, 10:43 Titel: Eine Delaunay-Triangulation filtern

Hi,

mein Problem ist folgendes: ich hab versucht eine Punktewolke mithilfe von Delaunay-Triangulation zu visualisieren.
Kann mir jemand helfen die ganzen Strukturen von denjenigen Dreiecken zu bereinigen,
deren Normalen orthogonal zu x-y-Ebene sind(sind an den Stellen fn(:,3)==1 oder -1).

Code:

dt = DelaunayTri(X,Y,Z);
[tri Xb] = freeBoundary(dt);
tr = TriRep(tri, Xb);
P = incenters(tr);
fn = faceNormals(tr);
trisurf(tri,Xb(:,1),Xb(:,2),Xb(:,3), ...
'FaceColor', 'cyan', 'faceAlpha', 0.8);
axis equal;
hold on;

quiver3(P(:,1),P(:,2),P(:,3), ...
fn(:,1),fn(:,2),fn(:,3),0.5, 'color','r');
hold off;

Funktion ohne Link?

Bin für jede Hilfe dankbar.
MfG Tanja

Tria.fig

Beschreibung:

Download

Dateiname:

Tria.fig

Dateigröße:

236.18 KB

Heruntergeladen:

441 mal

Maddy

Ehrenmitglied


	Beiträge: 494

	Anmeldedatum: 02.10.08

	Wohnort: Greifswald

	Version: ---

Verfasst am: 10.06.2009, 11:58 Titel:

Ich kenn mich mit dem Verfahren nicht aus, aber wenn du weißt dass die Flächennormale die entsprechende Werte hat, dann suchst du dir die Zeilen raus in denen die Bedingung nicht gilt und plottest eben nur diese.

Code:

zeilen=find(fn(:,3)~=1 & fn(:,3)~=-1);

Funktion ohne Link?

_________________

>> why
The computer did it.

Tanja86

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 10.06.2009, 12:17 Titel:

So ähnlich hab ich das auch versucht, aber das lasst fn zu einem
6400*1-vektor schrumpfen, es sollte aber weiterhin eine 6700*3-Matrix bleiben,
nur halt ohne die Zeilen die an der dritten Stelle ungleich Null sind.

Was mich auch noch irritiert, sind alle Werte in der fn-Matrix zwischen -1 und 1 und nach dem

Code:

zeilen=find(fn(:,3)~=1 & fn(:,3)~=-1) ;

Funktion ohne Link?

besteht der Vektor zeilen aus Werten zwischen 20 und 6700.

Danke für die schnelle Antwort.
MfG Tanja

Tanja86

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 10.06.2009, 12:27 Titel:

Sry hab meinen Denk-Fehler erkannt.

Der Vektor "zeilen" besteht dann aus Indizies der Zeilen mit den Werten ungleich -1 und 1.

Danke nochmal für die schnelle Antwort.
MfG Tanja

Maddy

Ehrenmitglied


	Beiträge: 494

	Anmeldedatum: 02.10.08

	Wohnort: Greifswald

	Version: ---

Verfasst am: 10.06.2009, 12:27 Titel:

Du sollst ja fn nicht schrumpfen lassen, sondern dir die Zeilen ausgeben lassen, die deiner Bedingung entsprechen (also eine flächennormale ohne 1 bzw. -1 in der 3. Komponente)

Wenn du nun

Code:

fn(zeilen,:)

Funktion ohne Link?

nutzt, hast du nur die Flächennormalen die nicht 1 bzw. -1 sind. Das selbe kannst du dann auch für alle anderen Matrizen anwenden. Der Vektor "zeilen" gibt dir also die Indizes vor, welche deinen Wünschen entsprechen.

edit: gleichzeitig Wink

_________________

>> why
The computer did it.


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen

Impressum | Nutzungsbedingungen | Datenschutz | FAQ |

RSS

Hosted by:

Copyright © 2007 - 2024 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.