Dear 3dslicer community. I try to write script which would calculate the angle between two segments. Here is the image I have [image] I obtained the segments using simple threshold Then I would need to draw (automatically by python script as it is a part of an extension) lines through the segments and calculate the angle between the two lines. As part of the extension I already have the center coordinates of the segments. do you have any ideas ? if it is needed, the whole script which give me the image is here: https://github.com/AurelienCD/Radiotherapy_Quality_Control/blob/master/Twist_MLC_Alignment_TOMOTHERAPY_QC.py Many thanks, AurelienCD

Segment Statistics module solves this very nicely: it can compute oriented bounding box and principal axes directions for each segment, directly from the segmentation. [image] New Segment Statistics: Oriented Bounding Box, Diameter and More Support New statistics have been added to the Segment Statistics module: Centroid Oriented bounding box Feret diameter Perimeter Roundness Flatness Elongation Principal axes/moments These statistics are not enabled by default. To enable them, click on the “Options” button beside the “Labelmap Statistics” checkbox, and select the statistics that you would like calculated from the popup window. [image] You can read the documentation for the Segment Statistics module here. Oriented bounding box: […

Greetings, Get two points (e.g. along axis of symmetry) from each segment. Two points gives you a vector for each segment. Calculate cross or dot product of that vectors, that gives you either sine or cosine value of the angle. The problem is how to calculate two points along axis of symmetry, or something like that? If the task is not on the plane, then you need three points.

@AurelienCD Hello! I would use the script repository: https://www.slicer.org/wiki/Documentation/Nightly/ScriptRepository#Measure_angle_between_two_markup_lines

@AurelienCD I would consider using Principal Component Analysis, specifically the sklearn implementation. Form an n-by-2 matrix of points in image space, either taking the positions of all points in one segment, or all points on the boundary of a segment (experiment to find which works better for you). Call this matrix X. Use sklearn.decomposition.PCA().fit(X).components_ to get the first principal component. This should be a 2-vector that points along the length of the segment. Measure the angles between these orientation vectors.

Thank you James for your correction of my code on github it help a lot. I am a beginner and comments are still necessary :slight_smile:

Thank you Andras, I will use that new module. thank you everyone for your help. all the best, Aurélien

Calculate angle between two segments

Support

lassoan (Andras Lasso) July 7, 2020, 4:27pm 5

Segment Statistics module solves this very nicely: it can compute oriented bounding box and principal axes directions for each segment, directly from the segmentation.

Topic		Replies	Views
Calculate angle between oriented bounding boxes Development segmentation	2	1647	December 28, 2020
Calculate dihedral angles using plug-ins Support	1	155	May 13, 2024
How to measure the angle? Support	22	3200	April 9, 2020
Angle Planes measurement Support	1	793	September 24, 2019
Measuring Error between desired and current trajectories Support surgical-navigation	3	1056	July 10, 2021

Calculate angle between two segments

Related topics