计算机图形学——多边形区域填充算法

本文主要分析并实现了2种常见的多边形填充算法:逐点法、有序边表法

1.逐点法

在逐点法实现中,用到了两个函数。

其中的实现原理十分简单,不过PointInPolygon(int pointNumOfPolygon, CPoint tarPolygon[], CPoint testPoint)函数的实现比较复杂,但这个函数也十分实用,在后面的三角剖分中也有使用。

算法描述:
(1)将给定多边形输入;

(2)求出多边形的最小包含矩形;

(3)逐点扫描最小矩形的每一点,并判断是否位于多边形内部,从最小点到最大点一次判断,如果在该多边形内部,则将该点上色;

(4)判断位于多边形内部的方法是,过每一点水平向右作射线,与多边形边界求交点,如果交点个数为奇数,则说明该点在多边形内部,偶数则说明在多边形外部。

效果如下:
在这里插入图片描述
实现代码:

//逐点法填充
void PointFillpoly(int pNumbers, CPoint *points, CDC *pDC)
{
 	BoxRect_t polyRect;
 
 	polyRect = getPolygonRect(pNumbers, points);
	m_pDC->MoveTo(polyRect.minX, polyRect.minY);

	m_pDC->LineTo(polyRect.minX, polyRect.maxY);
	m_pDC->LineTo(polyRect.maxX, polyRect.maxY);
	m_pDC->LineTo(polyRect.maxX, polyRect.minY);
	m_pDC->LineTo(polyRect.minX, polyRect.minY);

	CPoint testPoint;
	//从最小点到最大点一次判断,如果在该多边形内部,则进行填充
	for (testPoint.x = polyRect.minX; testPoint.x < polyRect.maxX; testPoint.x++)
		for (testPoint.y = polyRect.minY; testPoint.y < polyRect.maxY; testPoint.y++) {
			if (PointInPolygon(m_pNumbers, m_pAccord, testPoint))
				pDC->SetPixel(testPoint.x, testPoint.y, RGB(255, 255, 255));
		}
}
//得到该多边形的最大、最小的y、x值
BoxRect_t getPolygonRect(int pointNumOfPolygon, CPoint tarPolygon[])
{
	BoxRect_t boxRect;

	boxRect.minX = 50000;
	boxRect.minY = 50000;
	boxRect.maxX = -50000;
	boxRect.maxY = -50000;

	for (int i = 0; i < pointNumOfPolygon; i++) {
		if (tarPolygon[i].x < boxRect.minX) boxRect.minX = tarPolygon[i].x;
		if (tarPolygon[i].y < boxRect.minY) boxRect.minY = tarPolygon[i].y;
		if (tarPolygon[i].x > boxRect.maxX) boxRect.maxX = tarPolygon[i].x;
		if (tarPolygon[i].y > boxRect.maxY) boxRect.maxY = tarPolygon[i].y;
	}
	return boxRect;
}
//判断点是否位于区域内
BOOL PointInPolygon(int pointNumOfPolygon, CPoint tarPolygon[], CPoint testPoint)
{
	if (pointNumOfPolygon < 3)
		return false;

	bool  inSide = FALSE;
	float lineSlope, interSectX;
	int   i = 0, j = pointNumOfPolygon - 1;

	for (i = 0; i < pointNumOfPolygon; i++) {
		if ((tarPolygon[i].y < testPoint.y && tarPolygon[j].y >= testPoint.y ||
			tarPolygon[j].y < testPoint.y && tarPolygon[i].y >= testPoint.y) &&
			(tarPolygon[i].x <= testPoint.x || tarPolygon[j].x <= testPoint.x)) {
			if (tarPolygon[j].x != tarPolygon[i].x) {
				lineSlope = (float)(tarPolygon[j].y - tarPolygon[i].y) / (tarPolygon[j].x - tarPolygon[i].x);
				interSectX = (int)(tarPolygon[i].x + (testPoint.y - tarPolygon[i].y) / lineSlope);
			}
			else
				interSectX = tarPolygon[i].x;
			if (interSectX < testPoint.x)
				inSide = !inSide;
		}
		j = i;
	}

	return inSide;
}

2.有序边表法

首先给出有序边表法中定义使用的变量。

int m_Begin, m_End, m_edgeNumbers, m_Scan;       //求交边集指针,有效边数,当前扫描位置
float m_yMax[N], m_yMin[N], m_Xa[N], m_Dx[N];	//每条边y最大值,y最小值,每条边与扫描线x的交点,每条边斜率倒数

基本思想:用水平扫描线从上到下(或从下到上)扫描由多条首尾相连的线段构成的多边形,每根扫描线与多边形的某些边产生一系列交点。将这些交点按照x坐标排序,将排序后的点两两成对,作为线段的两个端点,以所填的颜色画水平直线。

  • 在实现中首先需要对多边形的每一条边进行处理操作,因为是采用水平线扫描填充,所以对于水平线不做处理。对其余边根据端点的y值大小进行排序,保存yMax,yMin,Xa(相对应的x坐标)以及Dx(斜率的倒数)。
  • 在进行扫描的过程中,很重要的一部分操作便是求交边集指针的移动。初始状态位0,在扫描线开始向下移动后,调用Include()函数,检查是否有边进入扫描线交集(即判断所有y最大值大于扫描线当前y值的边线),此时将m_End++,即尾指针向后移动。在Include()函数中,也会调整起始点位置,将Dx调整为位移量。
  • 之后调用UpdateXvalue()函数,判断是否有边退出求交边集。
    • 如果没有边退出,则移动x,并根据x值大小进行排序。
    • 有边退出,更新数组,删除该边,m_Begin++,即头指针向后移动。
  • 调用pFillScan(pDC)函数,进行填充。
  • m_Scan–,回到第二步进行循环,直到m_Begin==m_End。

算法描述:

  • 建立边表的方法:

(1)与x轴平行的边不计入;

(2)多边形的顶点分为两大类:一类是局部极值点,另外一类是非极值点。当扫面线与第一类顶点相交时,应看作两个点;而当扫描线与第二类顶点相交时,应视为一个点,对于极值点则要记录两条边;

(3)扫描线按照y轴从低到高顺次记录;

(4)一条边按照y轴的高低记录;

(5)多条边以x轴递增顺序记录;

  • 算法流程:

1、根据给出的多边形顶点坐标,建立NET表,求出顶点坐标中最大y值ymax和最小y值ymin。

2、初始化AET表指针,使它为空。

3、执行下列步骤直至NET和AET都为空.

(1)如NET中的第y类非空,则将其中的所有边取出并插入AET中;

(2)如果有新边插入AET,则对AET中各边排序;

(3)对AET中的边两两配对,(1和2为一对,3和4为一对,…),

将每对边中x坐标按规则取整,获得有效的填充区段,再填充.

(4)将当前扫描线纵坐标y值递值1;

(5)如果AET表中某记录的ymax=yj,则删除该记录 (因为每条边被看作下闭上开的);

(6)对AET中剩下的每一条边的x递增dx,即x’ = x+ dx .

效果如下:
在这里插入图片描述
实现代码如下:

//有序边表法填充
void Fillpolygon(int pNumbers, CPoint* points, CDC* pDC)
{
	m_edgeNumbers = 0;
	pLoadPolygon(pNumbers, points);   // Polygon Loading, calculates every edge's m_yMax[],m_yMin[],m_Xa[],m_Dx[]

	//求交边集范围,因为数组已经根据y值大小进行边的排序,所以end向后移动即代表有边进入,start向后移动,代表有边退出
	m_Begin = m_End = 0;
	m_Scan = (int)m_yMax[0];          //从顶向下扫描
	Include();                        //检查是否有边进入扫描线
	UpdateXvalue();                   //检查是否有边退出扫描线
	while (m_Begin != m_End) {
		pFillScan(pDC);
		m_Scan--;
		Include();
		UpdateXvalue();
	}
}
void pLoadPolygon(int pNumbers, CPoint* points)
{
	float x1, y1, x2, y2;

	x1 = points[0].x;    y1 = points[0].y + 0.5;
	for (int i = 1; i < pNumbers; i++) {
		x2 = points[i].x;    y2 = points[i].y + 0.5;
		if (abs(int(y2 - y1)) >= 0)	//水平线不做处理
		{
			pInsertLine(x1, y1, x2, y2);
			x1 = x2;      y1 = y2;
		}
		else
			x2 = x1;
	}
}
void pInsertLine(float x1, float y1, float x2, float y2)
{
	int i;
	float Ymax, Ymin;

	Ymax = (y2 > y1) ? y2 : y1;
	Ymin = (y2 < y1) ? y2 : y1;
	i = m_edgeNumbers;
	//根据y值的大小,进行排序插入,大的在前面
	while (i != 0 && m_yMax[i - 1] < Ymax) {
		m_yMax[i] = m_yMax[i - 1];
		m_yMin[i] = m_yMin[i - 1];
		m_Xa[i] = m_Xa[i - 1];
		m_Dx[i] = m_Dx[i - 1];
		i--;
	}
	m_yMax[i] = Ymax;
	m_yMin[i] = Ymin;
	if (y2 > y1) m_Xa[i] = x2;	//根据y大小确定Xa的值,y大的会先于扫描线相交
	else         m_Xa[i] = x1;

	m_Dx[i] = (x2 - x1) / (y2 - y1);	//斜率的倒数
	m_edgeNumbers++;
}
void Include()
{
	//end向后移动,找出所有边最高点y值大于当前扫描线的边,看是否有新的边进入交集
	while (m_End < m_edgeNumbers && m_yMax[m_End] >= m_Scan) {
		//有边进入,调整起始点位置,然后将Dx调整为位移量
		m_Xa[m_End] = m_Xa[m_End] + (-0.5) * m_Dx[m_End];
		m_Dx[m_End] = -m_Dx[m_End];
		m_End++;
	}
}
void UpdateXvalue()
{
	int i, start = m_Begin;

	for (i = start; i < m_End; i++) {
		if (m_Scan > m_yMin[i]) {
			//当前边没有退出,则移动x,然后在进行排序
			m_Xa[i] += m_Dx[i];
			pXsort(m_Begin, i);
		}
		else {
			//有边退出,更新数组,然后begin++
			for (int j = i; j > m_Begin; j--) {
				m_yMin[j] = m_yMin[j - 1];
				m_Xa[j] = m_Xa[j - 1];
				m_Dx[j] = m_Dx[j - 1];
			}
			m_Begin++;
		}
	}
}
void pXsort(int Begin, int i)
{
	float temp;

	while (i > Begin&& m_Xa[i] < m_Xa[i - 1]) {
		temp = m_Xa[i];   m_Xa[i] = m_Xa[i - 1];   m_Xa[i - 1] = temp;
		temp = m_Dx[i];   m_Dx[i] = m_Dx[i - 1];   m_Dx[i - 1] = temp;
		temp = m_yMin[i]; m_yMin[i] = m_yMin[i - 1]; m_yMin[i - 1] = temp;
		i--;
	}
}
void pFillScan(CDC* pDC)
{
	int x, y;
	FillPolyDoc* pDoc = GetDocument();
	for (int i = m_Begin; i < m_End; i += 2) {
		if (pDoc->m_displayMode == 1) {
			pDC->MoveTo(m_Xa[i], m_Scan);
			pDC->LineTo(m_Xa[i + 1], m_Scan);
		}
		else if (pDoc->m_displayMode == 4) {	//图案填充
			y = m_Scan;
			for (int x = m_Xa[i]; x < m_Xa[i + 1]; x++)
				if (m_patternData[y % 12][x % 12])
					pDC->SetPixel(x, y, RGB(255, 0, 0));
		}
	}
}
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值