# Intersect() Problem when none of the pixels overlap

3 views (last 30 days)
Hg on 28 Jan 2016
Commented: Hg on 1 Feb 2016
When I used intersect() to find the intersection between A (red) and B (blue), I failed to get any because the pixels don't overlap. What can I do to estimate the intersection point?
Ashish Uthama on 28 Jan 2016
Edited: Ashish Uthama on 28 Jan 2016
How would you define (for a computer) an intersection? For example, what if they just 'touch'? (Image 1 and 2 represent A and B)
1 0 0 2
0 1 2 0
0 1 2 0
1 0 0 2
Hg on 29 Jan 2016
In my case, the A (red) is a line so I can thicken or extend it to get an intersection.

Image Analyst on 28 Jan 2016
There is no one pixel where the overlap occurs. If you'll accept any of those 4 pixel locations as an overlap, then perhaps if you dilated, ANDed, then called bwulterode. Something like (untested)
intImage = imdilate(bw1, true(3)) & imdilate(bw2, true(3));
intPoints = bwulterode(intImage);
bwulterode() has some connectivity options, so you might play around with those.
##### 3 CommentsShow 1 older commentHide 1 older comment
Image Analyst on 29 Jan 2016
Edited: Image Analyst on 29 Jan 2016
That's exactly what the first line of code does. It thickens the lines. Then I AND it to find the intersection. It sounds like you think your method is different somehow, but I think they're the same. What I do in my second line of code, to prevent getting an intersection point that is thicker, is to call bwulterode() which shrinks everything down to a single point.
Another thing you could try is to OR the two binary images together and then use bwmorph() to get branchpoints:
intPoints = bwmorph(bw1 | bw2, 'Branchpoints');
Hg on 1 Feb 2016
Got it. It is exactly what I need. Thanks for the insight!

### Categories

Find more on Morphological Operations in Help Center and File Exchange

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!