center to center distance between objects: Need advise
2 views (last 30 days)
Abdul Hannan Qureshi on 30 Apr 2022
I am trying to calculate center to center distance between lines, i.e., between intersection points. I am working on this image
I want to achieve this
However, following steps I have assumed like
- detect centerlines in the image.
- detect intersection points.
- calculate horizantal distances between intersecting points.
- calculate vertical distances between intersecting points.
I am a bit confused about functions to be applied here. If anyone can advice which fucntions will be helpful to be to achieve above outcome OR if there is some better methodlogy to achieve above shared outcome. Kindly advise.
Another way is to get two images, one of all the vertical lines and one of all the horizontal lines. Use imopen() with the proper structuring element - either a vertical or horizontal rectangle.
Then threshold and get the centerlines with bwskel()
Then get a list of (x,y) for each line.
Then compare every vertical line with every horizontal line to find the distance that is min. This is where they cross. You can either use pdist2() or use a loop with sqrt(). This will give you the (x,y) intersection point for every line.
After that it's just some looping and bookkeeping to get the distances between adjacent points.
Again, if you can't figure it out, let me know.
More Answers (1)
What I'd do is what I actually do in a very similar situation, though I'm looking for 5 points on the perimeter not 5. Here is the algorithm I think you can use
- binarize the image to find dark rectangles mask = imbinarize();
- call mask = bwareafilt(mask, [smallest, largest]) to get rid of small blobs and the super large one touching the border and get ONLY the valid squares.
- call imfill(mask, 'holes');
- Find centroid with props = regionprops(mask, 'Centroid', 'Image')
- Loop over all blobs and for each blob, props(k).Image, find boundaries with bwtraceboundary. Start on the left middle edge so you have a known starting point.
- Find distance of all boundary points to that blob's centroid and plot them. distances = sqrt((xCentroid-x)^2+(yCentroid2-y)^2)
- Find peaks with findpeaks(distances, .....). Specify the MinPeakDistance as half the approximate width of the square. If there are more than 4 peaks found, take the 4 with the highest prominence value. This gives you the index of the corners in the boundary coordinates. You'll have 1 is the upper left, 2 = upper right, 3 = lower right, and 4 = lower left.
- Compute distances like d12 = sqrt((x2-x1)^2+(y2-y1)^2), and similar for d23, d34, and d41. Stored the distances for each blob.
If you can't code it up from the pseudocode I've given so far, then let me know and I'll finish it.