|
No offense taken. I was serious - if I did some stupid crossposting thing, I really didn't mean to and I'd love to learn what I did so I could avoid it in the future.
|
|
|
|
|
The cross-posting reference was about the original post (not you) - Cpallini's reply linked to the an identical post on another forum that had some useful discussion.
What i'd like to know is why folk want to solve this particular problem?
Peter
"Until the invention of the computer, the machine gun was the device that enabled humans to make the most mistakes in the smallest amount of time."
|
|
|
|
|
Ah - okay. Sorry to be so dense. I don't often interact with these forums (fora?) so assumed I'd made some sort of boneheaded blunder.
Don't know about a motive. Could just be "recreational math". I think the original poster was under the impression that there was a single unique solution which suggests to me that he/she really was trying to solve a "real world" problem, but who knows?
I do know that my estimation of the simplicity of finding intersections of quadrics hasn't been good. There's a similar problem in the Fortune algorithm for Voronoi diagrams where you have to find the intersection of a couple of parabolas and I went in with the assumption that it would be a trivial problem. I was way off. Eventually, I just dumped everything into Mathematica and let it solve the problem. Every description I've ever seen of this algorithm passes this problem off as "find the interesection of the parabolas" as if it was some of throw away detail. You'd think somebody would stop and note that this is a damnedly difficult problem and here is the (rather ugly)solution. Oh well - if only the world were run the way I think it ought to be...
Just FYI, the problem is essentially that you're given two focii of parabolas with a common horizontal line as the directrix for both parabolas. Find the (two) intersections between the parabolas. Seemed simple. Wasn't. I'm still not sure of any terribly "intuitive" way of solving this. No obvious transform to do on this one since the directrix is already horizontal.
|
|
|
|
|
Now here's my chance to overestimate the simplicity of a problem, surely the equation of a parabola with horizontal directrix is y = P(x) where P(x) is a quadratic, so for two such parabolas you need to solve
P1(x) = P2(x)
which is a 'simple' quadratic. I can see that the coefficients may be messy, but straightforward (?).
Peter
"Until the invention of the computer, the machine gun was the device that enabled humans to make the most mistakes in the smallest amount of time."
|
|
|
|
|
Yes, that's the kind of thinking I had for both of these questions - equating two quadratics shouldn't be that tough. It's the fact that you've got two quadratic equations in two unknowns in both cases that complicate matters. In the Fortune case, the equations are:
(ys-y)^2 = (x1-x)^2 + (y1-y)^2 = (x2-x)^2 + (y2-y)^2
where ys is the y coordinate of the sweepline and (x1,y1), (x2,y2) are the two points. Using one equation, you can get y as a quadratic function of x. Substituting that into the other equation gives a quartic equation in x and you have similar problems to the ones you deal with in the original problem of this thread. I'm glad I'm not the only one who can be misled by these problems! What's more, like I said, I don't see any nice intuitive transforms or any other easy way, geometric or otherwise, of coming up with the solution for this one. Certainly doesn't mean there isn't one - in fact I suspect that there is such a thing - I just haven't seen it yet. When I got the solution from Mathematica, I was able to break it into lots of redundant pieces to improve performance which suggests to me that there ought to be a way. Perhaps if I meditated on the solution long enough, I could see it, but I haven't spent the time.
Minus some exceptional cases, here's the code I ended up with (directly from my code):
double a1 = 1 / (2 * (pt1.Y - ys));
double a2 = 1 / (2 * (pt2.Y - ys));
double da = a1 - a2;
double s1 = 4 * a1 * pt1.X - 4 * a2 * pt2.X;
double dx = pt1.X - pt2.X;
double s2 = 2 * Math.Sqrt(2 * (2 * a1 * a2 * dx * dx - da * (pt1.Y - pt2.Y)));
double m = 0.25 / da;
double xs1 = m * (s1 + s2);
double xs2 = m * (s1 - s2); It's been a long time since I did this. I think the rounding was done to avoid some numerical instability problems.
xs1 and xs2 are the x coordinates of the two solutions. For my purposes, I didn't really care about the y coordinates.
Darrell
|
|
|
|
|
Oops - okay - I see what you mean. Doh! I think you're right. I'm sorry to be so dim witted! Since you end up with the y^2 terms dropping out, you can solve for y, set the resulting two quadratics in x equal to each other for one quadratic in x and solve. Man, I better just stop before I stick my foot in my mouth again!
Darrell
|
|
|
|
|
Peter
"Until the invention of the computer, the machine gun was the device that enabled humans to make the most mistakes in the smallest amount of time."
|
|
|
|
|
Peter,
cp9876 wrote: as are others (except for those claiming 4 solutions)
Unless the triangle has the same length legs (ac=bc), there will be 4 solutions. Nothing in the stated problem indicated any orientation. For example, line ab can be taken as a vertical line starting with a at the bottom and b at the top, or, flip that upside down, then point c can be on either side of that vertical line, in either orientation of ab. Essentially, you need to draw two circles of radius ac, one at one end of the diameter ab and one at the other end of the diameter ab. There will be 4 points of intersection of these two circles with the circle drawn from the midpoint of ab with radius of ab/2. There is also no constraint on whether ac is the longest or shortest leg.
Dave.
|
|
|
|
|
Chesnokov Yuriy wrote: There is ABC triangle randomly rotated on XY plane. A,B coords are known, ACB angle is 90'. AC length is known. How to find C coords????
Member 4194593 wrote: Nothing in the stated problem indicated any orientation.
We seem to be solving different problems, I assumed that as the poster says that the coordinates of A and B are known that he knows which are which, in particular which coordinates are the end of the side AC. i.e. we are given (xa, ya), (xb, yb) and the length |AC|.
If he really doesn't know which is which, he can easily get the other solutions by swapping the A and B coords.
Peter
"Until the invention of the computer, the machine gun was the device that enabled humans to make the most mistakes in the smallest amount of time."
|
|
|
|
|
Peter,
You are correct as usual. You cannot draw the circle of length AC at B, it has to be drawn at A, thus only two solutions, one on either side of line AB.
One solution could be to start at the origin and call the point A, drop down distance AC on the y axis and call the point C, move to the right (the mirrored solution would be obtained by moving to the left) a distance sqrt(AB^2-ac^2) and call that point B. Now rotate ABC about A until the line AB matches the original slope of the given line AB, then translate ABC to position A at the origin to the original given position of A.
Dave.
|
|
|
|
|
Dave,
Neat constructional idea, but the algebra looks nasty to work out the sin and cosine of that angle for the rotation,
Peter.
Peter
"Until the invention of the computer, the machine gun was the device that enabled humans to make the most mistakes in the smallest amount of time."
|
|
|
|
|
I would like to check if a segment and a polygon have some intersections. Could you plz indicate me some algorithms to solve it. I knew Ray Casting algorithm and I wonder if there are other algorithms ?
|
|
|
|
|
|
The above algorithm won't detect the case where the segment is entirely contained in the polygon, without intersecting any edge. Let me know if you need to detect this case too.
|
|
|
|
|
Plz show me if that is the most optimal algorithm ? I really want to know the most optimal algorithm to solve this problem. So, would you mind helping me ?
|
|
|
|
|
If it isn't optimal it's pretty close. The calculation to determine if two segments intersect isn't time consuming. The only case it doesn't handle, as I've mentioned, is when the segment is entirely contained within the polygon.
modified on Friday, October 17, 2008 9:55 AM
|
|
|
|
|
OK. Thanks so much for replying
|
|
|
|
|
Alan,
Does it handle the case where the intersection is at the node between two line segments of the polygon?
Dave.
|
|
|
|
|
It does if your algorithm to detect line-segment intersection will detect a segment intersecting the end point of another segment.
|
|
|
|
|
Alan,
My question was about the suggested algorithm you gave above and whether it returned an intersection for a crossing at a node.
Dave.
|
|
|
|
|
Yes it does; if you follow the link there's a Java applet that lets you draw two line segments and tells you if they intersect. You can draw two segments that intersect at an end, and it does detect the intersection.
|
|
|
|
|
Alan,
Thank you for the answer.
I did start to try the link. First it wanted to install a plugin in FF, OK, I'll go with that. That failed, and the next thing it wanted to do was install a Sun DLL. I decided to quit at that point because the next thing it would want to do would be enable Java Script which I have disabled.
I'm happy that that solution will work if I ever need it.
Dave.
|
|
|
|
|
Depends on your situation. You should probably do a bounding box check first since it's simple, quick and will eliminate a lot of lines.
There's an entire book essentially on this subject:
http://www.amazon.com/Real-Time-Collision-Detection-Interactive-Technology/dp/1558607323/ref=sr_1_1?ie=UTF8&s=books&qid=1227945892&sr=1-1[^]
Some of the things you need to think about:
Are your polygons convex or can they be concave?
Are you testing fixed polygons against a series of lines or does the polygon change with each test?
Are your lines more likely to miss than not?
BSP trees are worth looking into if you really want to squeeze every last drop of performance. They're a bit complex, though and mostly good for lots of lines testing against a single polygon.
|
|
|
|
|
Iam student have aproject about retreiving image
,so Iwant agood algorithm about analysis images to find the signature of it and uses this step to retreive...
plz plz plz help me in a shortest possible time ..
thanks alot..
|
|
|
|
|
luna2010 wrote: plz plz plz help me in a shortest possible time ..
Done.
BTW: Google & CodeProject's article search engine are your best friends (reading the forum guidelines will be a plus).
If the Lord God Almighty had consulted me before embarking upon the Creation, I would have recommended something simpler.
-- Alfonso the Wise, 13th Century King of Castile.
This is going on my arrogant assumptions. You may have a superb reason why I'm completely wrong.
-- Iain Clarke
[My articles]
|
|
|
|
General 
News 
Suggestion 
Question 
Bug 
Answer 
Joke 
Praise 
Rant 
Admin 
Submit an article or tip
