File: ClipPoly.c

package info (click to toggle)
acm 5.0-19
  • links: PTS
  • area: main
  • in suites: sarge
  • size: 7,852 kB
  • ctags: 4,792
  • sloc: ansic: 42,427; makefile: 706; cpp: 293; perl: 280; sh: 198
file content (57 lines) | stat: -rw-r--r-- 1,179 bytes parent folder | download | duplicates (9)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57

#include <Alib.h>

#define IsInside(w, pt)	pt.x < w->clip.x1 || pt.x >= w->clip.x2 || \
			pt.y < w->clip.y1 || pt.y >= w->clip.y2

/*
 *  A quick and dirty 2D polygon clipping algorithm.
 */

void
Clip2DPolygon(AWindow * w, Point * pts, int npts, Point * clippedPts, int *nClippedPts)
{

	register int numPts, j, lastj, inside, lastInside;
	int       x1, y1, x2, y2;

	lastj = npts - 1;
	numPts = 0;

	lastInside = IsInside(w, pts[lastj]) ? 1 : 0;

	for (j = 0; j < npts; ++j) {

		if (lastInside == 1) {
			clippedPts[numPts++] = pts[lastj];
		}

		inside = (IsInside(w, pts[j])) ? 1 : 0;

/*
 *  Fast, but not accurate.  This method does not handle polygons in corners
 *  well.  A more accurate method would be to clip by each edge independently.
 */

		if ((inside ^ lastInside) == 1) {
			x1 = pts[lastj].x;
			y1 = pts[lastj].y;
			x2 = pts[j].x;
			y2 = pts[j].y;
			ClipLine(w, &x1, &y1, &x2, &y2);
			if (x2 == pts[lastj].x && y2 == pts[lastj].y) {
				clippedPts[numPts].x = x1;
				clippedPts[numPts++].y = y1;
			}
			else {
				clippedPts[numPts].x = x2;
				clippedPts[numPts++].y = y2;
			}
		}

		lastInside = inside;
		lastj = j;
	}

	*nClippedPts = numPts;
}