GoogleCodeExporter
commented
9 years ago I'm not sure that your sugest is correct, maybe I need a test case for try to
replicate the error.
BTW, I've check the method in JTS and looks really different from NTS impl
(maybe a
chabge that I've not see in my port).
So please try to change your NTS code, maybe retrieving the code from the trunk
or
change manually using this code:
public static bool IsPointInRing(ICoordinate p, ICoordinate[] ring)
{
if (p == null)
throw new ArgumentNullException("p", "coordinate 'p' is null");
if (ring == null)
throw new ArgumentNullException("ring", "ring 'ring' is null");
int crossings = 0;
/*
* For each segment l = (i-1, i), see if it crosses ray from test point
in positive x direction.
*/
for (int i = 1; i < ring.Length; i++)
{
int i1 = i - 1;
ICoordinate p1 = ring[i];
ICoordinate p2 = ring[i1];
if (((p1.Y > p.Y) && (p2.Y <= p.Y)) || ((p2.Y > p.Y) && (p1.Y <= p.Y)))
{
double x1 = p1.X - p.X;
double y1 = p1.Y - p.Y;
double x2 = p2.X - p.X;
double y2 = p2.Y - p.Y;
/*
* segment straddles x axis, so compute intersection.
*/
double xInt = RobustDeterminant.SignOfDet2x2(x1, y1, x2, y2) /
(y2 - y1);
/*
* crosses ray if strictly positive intersection.
*/
if (0.0 < xInt)
crossings++;
}
}
/*
* p is inside if number of crossings is odd.
*/
if ((crossings % 2) == 1)
return true;
else return false;
}
Let me see if this solve the problem, or please post a test case :)
ThanksOriginal comment by diegogu...@gmail.com on 8 Mar 2008 at 11:25
- Changed state: Accepted
- Added labels: Priority-High
- Removed labels: Priority-Medium
Attachments:
Original issue reported on code.google.com by
LivelyTe...@gmail.comon 8 Mar 2008 at 10:48