Calcul de l'aire des points discrets dans Matlab

Question

J'ai un ensemble de points discrets en 2D, qui ont besoin de calculer la surface interne, ou simplement la longueur et la hauteur, n'a pas d'importance évidemment. Mais le point de cette question est: pas nécessairement nous avons des bords tranchants qui façonnent la surface interne. J'ai essayé les commandes boundary et inpolygon, mais semble ne sont pas les bonnes pour ce problème. Que diriez-vous des bords que nous n'avons pas des bords droits, comme des arcs ou des demi-cercles.

https://drive.google.com/file/d/0B-DuVkE6qbU3ZnJsbWdQNjVmVkk/view?usp=sharing

Les fichiers fournis à l'intérieur de lien:

• Fichier 1 sont constitués de coordonnées XYZ, plus l'intensité (fourni par les données numérisées laser) Valeurs
• dossier 2 se composent de coordonnées XYZ, plus l'intensité , alors que toutes les valeurs Z et d'intensité sont fixées dans le fichier.
• Les figures 1-1, 1-2, 1-3, 1-4 représentent les données en 3D avec des valeurs d'intensité associées. Différentes valeurs d'intensité affichent différents matériaux. Les flèches blanches présentent la profondeur des différentes parties du segment.
• Les figures 2-1, 2-2 représentent les données en 2D avec une seule valeur d'intensité. En fait, dans cette étape, les données 3D ont été projetées en 2D et toutes les valeurs d'intensité ont changé à une seule valeur, en outre, la partie interne de la fenêtre a été supprimée. La zone de la partie interne est de l'intérêt à calculer.

Ça a été un défi pour moi! Quelqu'un pourrait-il conseiller s'il vous plaît?

Answer 1

Peut-être que vous pouvez utiliser imclose, imdilate ou imerode fonctions afin d'obtenir un « objet » unique et vous pouvez calculer la zone en additionnant toutes les valeurs diferent de 0 ou 255, en fonction de votre mise en œuvre. Je suppose que l'image est en noir et blanc, sinon convertissez-la en premier. documentation Matlab:

• imclose: http://www.mathworks.com/help/images/ref/imclose.html

• im2bw: http://www.mathworks.com/help/images/ref/im2bw.html

Answer 2

Si vos points sont les coordonnées des sommets d'un polygone, vous pouvez utiliser polyarea pour retourner la zone de ce polygone.

Ex:

Area = polyarea(pointsX, pointsY) 

Prenez en compte que l'ordre des points la matière, et ils doivent former un polygone valide (segments ne doit pas traverser)