图形学各种算法代码

最新推荐文章于 2022-05-22 10:45:39 发布
最新推荐文章于 2022-05-22 10:45:39 发布
阅读量5.4k 收藏 18
点赞数 2
分类专栏： 图形学
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/damotiansheng/article/details/43345793
版权
目的：实现画直线算法，画圆，画椭圆算法，二维直线裁剪算法，种子扫描填充算法


1. DDA算法，算法思想参考百度百科：
http://baike.baidu.com/link?url=bKEk1Bmu6YO7gdPuQpXokbS2eDV8-X5UQaObYeg1NWyHvLS82d5RBZWfdr-tkHrVKfQWgdbciVf_BLNGh6jsFK


2. 各种算法实现：

 2.1 画直线算法：

//DDA算法
void drawAlgorithmDialog::drawLineDDA( int x0,int y0,int x1,int y1, QPainter* ppainter )
{
    //QPainter painter(this);

    int dx, dy, n, k;
    float xinc, yinc, x, y;
    dx = x1 - x0;
    dy = y1 - y0;

    if( qAbs(dx) > qAbs(dy) )
    {
        n = abs(dx);
    }
    else
    {
        n = abs(dy);
    }

    xinc = (float)dx / n;
    yinc = (float)dy / n;
    x = (float)x0;
    y = (float)y0;

    for( k = 1; k <= n; k++ )
    {
        ppainter->drawPoint( int(x + 0.5), int(y+0.5) );
        x+=xinc;
        y+=yinc;
    }

    return;
}
//中点画线
void drawAlgorithmDialog::drawLineMidPoint(int x0, int y0, int x1, int y1, QPainter *ppainter)
{
    int a, b, d1, d2, d, x, y;
        a = y0-y1;
        b = x1 - x0;
        d = 2 * a + b;
        d1 = 2 * a;
        d2 = 2 * ( a + b );
        x=x0;
        y=y0;
        ppainter->drawPoint( x, y );
        while( x < x1 )
        {
          if( d < 0 )
          {
              x++, y++, d+=d2;
          }
          else
          {
              x++, d += d1;
          }
          ppainter->drawPoint( x, y );
         //drawpixel (x, y, color);
        }
    return;
}
// BresenHam画线
void drawAlgorithmDialog::drawLineBresenHam(int x1, int y1, int x2, int y2, QPainter *ppainter )
{
    int dx = x2 - x1;//△x
    int dy = y2 - y1;//△y
    int p = (2*dy) - dx ; //P = 2△y - △x
    int dobDy = 2* dy ; // 2 △y
    int dobD = 2*(dy - dx) ; // 2(△y - △x)
    int PointX,PointY;
    //设置两个临时用来显示位置的变量
    if( x1 > x2)
    { //判断线段的方向
        PointX = x2;//起始坐标X
        PointY = y2;起始坐标Y
        x2 = x1;
    }
    else
    {
        PointX = x1;//起始坐标X
        PointY = y1;//起始坐标Y
    }
    //达因第一个起始点
    //cout<<"Point: X:"<<PointX<<" Y:"<<PointY<<endl;
    ppainter->drawPoint( PointX, PointY );
    while( PointX < x2)
    {
        PointX++;
        if(p < 0)
        {
            p += dobDy;
        }
        else
        {
            PointY++;
            p += dobD;
        }
            //cout<<"Point: X:"<<PointX<<" Y:"<<PointY<<endl;
        ppainter->drawPoint( PointX, PointY );
    }
    return;
}

2.2 画圆，椭圆： 

// 中点画圆
void drawAlgorithmDialog::MidPointCircle( int x1, int y1, int R, QPainter* ppainter )
{
    int x = 0;
    int y = R;
    double d = 1.25 - R;
    ppainter->drawPoint( x1,y1 );  //绘制圆心
    while( x < y )\
    {
        ppainter->drawPoint( x+x1,y+y1 );
        ppainter->drawPoint( y+x1,x+y1 );
        ppainter->drawPoint( -x+x1,y+y1 );
        ppainter->drawPoint( -y+x1,x+y1 );

        ppainter->drawPoint( -x+x1,-y+y1 );
        ppainter->drawPoint( -y+x1,-x+y1 );
        ppainter->drawPoint( x+x1,-y+y1 );
        ppainter->drawPoint( y+x1,-x+y1 );


        if( d < 0 )
        {
            d += 2 * x + 3;
            x++;
        }
        else
        {
            d += 2 * ( x - y ) + 5;
            x++;
            y--;
        }

        //ppainter->drawPoint( x, y );
    }

    return;
}

// Bresenham画圆
void drawAlgorithmDialog::BresenhamDrawCircle(int x1, int y1, int R, QPainter *ppainter)
{
    int x = 0;
    int y = R;
    int p = 3 - 2 * R;

    ppainter->drawPoint( x1,y1 );  //绘制圆心

    for( ; x <= y; x++ )
    {
        ppainter->drawPoint( x+x1,y+y1 );
        ppainter->drawPoint( y+x1,x+y1 );
        ppainter->drawPoint( -x+x1,y+y1 );
        ppainter->drawPoint( -y+x1,x+y1 );

        ppainter->drawPoint( -x+x1,-y+y1 );
        ppainter->drawPoint( -y+x1,-x+y1 );
        ppainter->drawPoint( x+x1,-y+y1 );
        ppainter->drawPoint( y+x1,-x+y1 );

        if( p >= 0 )
        {
            p += 4 * ( x - y ) + 10;
            y--;
        }
        else
        {
            p += 4 * x + 6;
        }
    }

    return;
}

// 中点画椭圆
void drawAlgorithmDialog::MidPointEllise( int xcenter, int ycenter, int a, int b, QPainter *ppainter) //中点画椭圆法
{
    int x = 0;
    int y = b;
    double d1 = 0;
    double d2 = 0;
    d1 = b * b + a * a * ( -b + 0.25 );
    ppainter->drawPoint( xcenter + x,ycenter + y);
    ppainter->drawPoint( xcenter - x,ycenter + y);
    ppainter->drawPoint( xcenter + x,ycenter - y);
    ppainter->drawPoint( xcenter - x,ycenter - y);

    while( b*b*(x+1) < a*a*(y-0.5) )
    {
        if( d1 < 0 )
        {
            d1 += b*b*(2*x+3);
            x++;
        }
        else
        {
            d1 += ( b*b*(2*x+3) + a*a*(-2*y+2 ) );
            x++;
            y--;
        }

        ppainter->drawPoint( xcenter + x,ycenter + y);
        ppainter->drawPoint( xcenter - x,ycenter + y);
        ppainter->drawPoint( xcenter + x,ycenter - y);
        ppainter->drawPoint( xcenter - x,ycenter - y);
    }

    d2 = b*b*(x+0.5)*(x+0.5) + a*a*(y-1)*(y-1)-a*a*b*b;
    while( y > 0 )
    {
        if( d2 < 0 )
        {
            d2 += b*b*(2*x+2) + a*a*(-2*y+3);
            x++;
            y--;
        }
        else
        {
            d2 += a*a*(-2*y+3);
            y--;
        }

        ppainter->drawPoint( xcenter + x,ycenter + y);
        ppainter->drawPoint( xcenter - x,ycenter + y);
        ppainter->drawPoint( xcenter + x,ycenter - y);
        ppainter->drawPoint( xcenter - x,ycenter - y);
    }
    return;
}

2.3 直线裁剪算法： 

// 二维直线裁剪-区域编码法裁剪算法：
void LineTrimDialog::makecode(double x, double y, int &c)
{
    c = 0;
    if( x < xl )
    {
        c = 1;
    }
    else if( x > xr )
    {
        c = 2;
    }

    if( y < yb )
    {
        c += 4;
    }
    else if( y > yt )
    {
        c += 8;
    }

    return;
}

//xl 窗口左边x坐标
//xr 窗口右边x坐标
//yb 窗口底部y坐标
//yt 窗口顶部y坐标
  //由于坐标原点默认在左上角，故，xr, xl, yt, tb应该为：
 // xl = vertexPoint[0].x();
//  xr = vertexPoint[1].x();
//  yb = vertexPoint[0].y();
//  yt = vertexPoint[1].y();
// vertexPoint[0]为裁剪窗口左上角坐标
// vertexPoint[1]为裁剪窗口右下角坐标

// 区域编码法裁剪 
void LineTrimDialog::Cohen_Sutherland( double x0, double y0, double x2, double y2, QPainter *pPainter, int i  )
{
   // pPainter->drawLine( 100, 100, 300, 300 );
   // QMessageBox::warning(this,tr("Warning"),"cohen",QMessageBox::Yes);
   // update();
    int c = 0;
    int c1 = 0;
    int c2 = 0;
    double x = 0;
    double y = 0;
    makecode( x0, y0, c1 );
    makecode( x2, y2, c2 );
    while( c1 != 0 || c2 != 0 )
    {
        if( c1&c2 == 1 )
        {
            return;
        }
        c = c1;
        if( c == 0 )
        {
            c = c2;
        }
        if( c&1 == 1 )
        {
            y = y0 + ( y2 - y0 ) * ( xl - x0 ) / ( x2 - x0 );
            x = xl;
        }
        else if( c&2 )
        {
            y = y0 + ( y2 - y0 ) * ( xr - x0 ) / ( x2 - x0 );
            x = xr;
        }
        else if( c&4 )
        {
            x = x0 + ( x2 - x0 ) * ( yb - y0 ) / ( y2 - y0 );
            y = yb;
        }
        else if( c&8 )
        {
            x = x0 + ( x2 - x0 ) * ( yt - y0 ) / ( y2 - y0 );
            y = yt;
        }
        if( c == c1 )
        {
            x0 = x;
            y0 = y;
            makecode( x, y, c1 );
        }
        else
        {
            x2 = x;
            y2 = y;
            makecode( x, y, c2 );
        }
    }
    /*
    if( x0 - (int)(x0) > 0.5 ) //防止double转int丢失较大精度
    {
        x0++;
    }
    if( y0 - (int)(y0) > 0.5 )
    {
        y0++;
    }
    if( x2 - (int)(x2) > 0.5 )
    {
        x2++;
    }
    if( y2 - (int)(y2) > 0.5 )
    {
        y2++;
    }*/
   // QMessageBox::warning(this,tr("Warning"),"t344st",QMessageBox::Yes);
   // LinePointArray[i] = QPoint( x0, y0 ); //会把x0,y0,x2,y2转换为int，会丢失精度
  //  LinePointArray[i+1] = QPoint( x2, y2 );
    //在这里直接绘制，可以避免由double转换为int时丢失的精度
    pPainter->drawLine( x0, y0, x2, y2 );  //向上偏移一下，方便看到裁剪后的直线
    update();
    return;
}

// 二维直线裁剪-梁友栋-Barsky算法：
bool LineTrimDialog::cansee(double q, double d, double &t0, double &t1)
{
    double r;
    if( q < 0 )
    {
        r = d / q;
        if( r > t1 )
        {
            return false;
        }
        else if( r > t0 )
        {
            t0 = r;
        }
    }
    else if( q > 0 )
    {
        r = d / q;
        if( r < t0 )
        {
            return false;
        }
        else if( r < t1 )
        {
            t1 = r;
        }
        else if( d < 0 )
        {
            return false;
        }
    }
    return true;
}
void LineTrimDialog::L_Barsky(double x0, double y0, double x2, double y2, QPainter *pPainter, int i )
{
    //QMessageBox::warning(this,tr("Warning"),"barsky",QMessageBox::Yes);
    double t0, t1, deltax, deltay;
    t0 = 0.0;
    t1 = 1.0;
    deltax = x2 - x0;

    if( !cansee( -deltax, x0 - xl, t0, t1 ) )
    {
        return;
    }

    if( !cansee( deltax, xr - x0, t0, t1 ) )
    {
        return;
    }

    deltay = y2 - y0;

    if( !cansee( -deltay, y0 - yb, t0, t1 ) )
    {
        return;
    }

    if( !cansee( deltay, yt - y0, t0, t1 ) )
    {
        return;
    }

    x2 = x0 + t1 * deltax;
    y2 = y0 + t1 * deltay;
    x0 = x0 + t0 * deltax;
    y0 = y0 + t0 * deltay;

   // LinePointArray1.push_back( DPoint( x0, y0 ) );
   // LinePointArray1.push_back( DPoint( x2, y2 ) );

    pPainter->drawLine( x0, y0, x2, y2 );
    update();

   /* if( x0 - (int)(x0) > 0.5 ) //防止double转int丢失较大精度
    {
        x0++;
    }

    if( y0 - (int)(y0) > 0.5 )
    {
        y0++;
    }

    if( x2 - (int)(x2) > 0.5 )
    {
        x2++;
    }

    if( y2 - (int)(y2) > 0.5 )
    {
        y2++;
    }*/

    //pPainter->drawLine( x0, y0, x2, y2 );

    //LinePointArray[i] = QPoint( x0, y0 );
    //LinePointArray[i+1] = QPoint( x2, y2 );

    return;
}

2.4 种子扫描填充算法

// 种子扫描填充算法 (x,y)为种子点
void CZhztchView::boundaryfill4(int x, int y, int boundarycolor, int newcolor)
{
	/* 
	int color;
	CClientDC dc(this);      //获取客户区设备描述表
	color=dc.GetPixel(x,y);
	
	if(color!=newcolor&&color!=boundarycolor)
	{
		dc.SetPixel(x,y,newcolor);
		boundaryfill4(x,y+1,boundarycolor,newcolor);
		boundaryfill4(x,y-1,boundarycolor,newcolor);
		boundaryfill4(x-1,y,boundarycolor,newcolor);
		boundaryfill4(x+1,y,boundarycolor,newcolor);
	} */

	int color;
	//color=dc.GetPixel(x,y);

	CClientDC dc(this);      //获取客户区设备描述表
	using namespace std;
	int x0 = 0;
	int xl = 0;
	int xr = 0;
	int y0 = 0;
	int xid = 0;
	bool flag = false;
	int count = 0;
	stack<Point> s;
	Point p;
	s.push( Point( x, y ) );

	while( !s.empty() ) //判断栈是否为空
	{
		p = s.top();
		s.pop();
		x = p.x;
		y = p.y;
		dc.SetPixel(x,y,newcolor);
		//painter->drawPoint( x, y );
		// SetPixel( hdc, x, y, 16711935 );
		// repaint();
		x0 = x + 1;

		while( dc.GetPixel( x0, y ) != boundarycolor && dc.GetPixel( x0, y ) != newcolor )// 填充种子右方像素
		{
			//painter->drawPoint( x0, y );
			dc.SetPixel(x,y,newcolor);
			//SetPixel( hdc, x, y, 16711935 );
			// repaint();
			x0++;
			//repaint();

		}

		xr = x0 - 1; //最右边像素位置
		x0 = x - 1;

		while( dc.GetPixel( x0, y ) != boundarycolor && dc.GetPixel( x0, y ) != newcolor ) //填充种子左边像素
		{
			//painter->drawPoint( x0, y );
			dc.SetPixel(x,y,newcolor);
			//SetPixel( hdc, x, y, 16711935 );
			// repaint();
			x0--;
		}

		xl = x0 + 1; //最左边像素位置
		y0 = y;

		for( int i = 1; i >= -1; i -= 2 )
		{
			x0 = xr;
			y = y0 + i;
			while( x0 >= xl )
			{
				flag = false;

				while( ( dc.GetPixel( x0, y ) != boundarycolor ) && ( dc.GetPixel( x0, y ) != newcolor ) && ( x0 >= xl ) )
				{
					if( !flag )
					{
						flag = true;
						xid = x0;
					}
					x0--;
				}

				if(flag)
				{
					s.push( Point( xid, y ) );
					flag = false;
				}

				while( ( dc.GetPixel( x0, y ) == boundarycolor ) || ( dc.GetPixel( x0, y ) == newcolor ) && ( x0 >= xl ) )
				{
					x0--;
				}
			}
		}

	}
	

}


damotiansheng
关注
  • 2
    点赞
  • 踩
  • 18
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
计算机图形学全代码,计算机图形学作业参考代码
weixin_31787335的博客
07-26 1120
该楼层疑似违规已被系统折叠隐藏此楼查看此楼代码仅供参考，至于怎么修改，你懂的....#include#include#include#includevoid DDALine(int x0, int y0, int x1, int y1, int color);void MidBresenhamLine(int x0, int y0, int x1, int y1, int color);void...
图像处理python代码（自用）
qq_42450300的博客
05-09 2616
1、二值图求差并计算比例 import cv2 import numpy as np import os from PIL import Image # 需要输入图像的路径 img_path = "./gtFine/train/raw/" mask_path = "./mask/train/" imgROI_path = "./diff/train/" if not os.path.exists(imgROI_path): os.makedirs(imgROI_path) sum = 0 co
Python 图形学 基于OpenGL图形的区域填充，种子递归算法
03-23
计算机图形学作业 单个py文件
【DeeplabV3+ MIoU】DeeplabV3+计算评价指标
计算机视觉方面的点点滴滴，欢迎一起讨论
04-18 5311
文章目录1 什么是miou2 获取分割结果+计算miou3 get_miou_png的代码解析4 compute_mIoU, show_results的代码解析感谢链接 1 什么是miou MIoU全称为Mean Intersection over Union，平均交并比。可作为语义分割系统性能的评价指标(像素级分类嘛)。 miou=1k∑i=1kP∩GP∪Gmiou=\frac{1}{k} \sum_{i=1}^{k} \frac{P \cap G}{P \cup G}miou=k1​i=1∑k​P∪GP
【计算机图形学 】多边形种子填充算法 | OpenGL+鼠标交互
Rover95的博客
10-21 4435
文章目录前言步骤1.点的结构体2.确定线的颜色3.设计画点函数（用于填充颜色的点）4.判断是否为边界的函数5. 扫描线算法实现用队列`queue`实现的 4连通区域的填充算法用栈实现的 扫描线形式的填充算法6. 鼠标响应事件完整代码 前言 实现多边形种子填充算法，分为四连通填充和扫描线填充，（分别由队列和栈实现的）并和鼠标进行交互。 具体原理略过，最后会贴上完整代码，可直接运行。 环境： vs2019,OpenGL的库（可以搜索如何用vs使用OpenGL的库，可以使用vs自带的插件或者其他方法，很方便） 要
图的各种算法
05-24 431
找工作时当练习写的，仅供参考。 头文件： #ifndef GRAPH_H_INCLUDED #define GRAPH_H_INCLUDED #define Max 10000 //adjacent matrix class MGraph{ private: int TotVer; vector Vertics; vector > E
作业_计算机图形学代码_
10-02
在这个特定的“作业_计算机图形学代码_”中，我们关注的是一个使用easyx库实现的一笔绘制金刚石图案的项目。EasyX是一个轻量级的C++图形库，它简化了在Windows上进行图形编程的过程，特别适合教学和初学者使用。 ...
图形学算法（伪代码）.rar_C++图形学_伪代码_图形学_图形学算法
09-23
这个文本文件可能包含了各种图形学算法的伪代码实现，如物体的坐标变换（平移、旋转、缩放）、视口变换、投影变换等。也可能涉及光栅化算法，如扫描线算法或Bresenham算法，以及纹理映射和混合算法。每个算法都用...
计算机图形学画图及代码实现——山农
05-30
在这个主题下，"山农"可能是指一种特定的算法或者一个用于学习和实践计算机图形学的实例。从提供的信息来看，这份资源可能是一个学生在学习过程中完成的期中和期末作业，包含了绘制各种图形的算法说明以及相应的代码...
图形学绘图代码-python
最新发布
04-20
用图形学的知识绘制椭圆,python
计算机图形学源代码（实现各种图形绘制和变换）
11-28
功能强大，模拟直线、圆、各种曲线、变换等
计算机图形学 梁友栋-Barsky 算法
05-06
计算机图形学经典算法 ，希望对大家有所帮助
计算机图形学实验代码（包括图形旋转、图形平移、图形缩放）
05-20
计算机图形学实验代码（包括图形旋转、图形平移、图形缩放） 内含C++源代码和测试图片
电子科技大学计算机图形学3Dview
09-05
电子科技大学研究生计算机图形学project1-3D VIEW完整代码，可直接运行
图形学最常用的几个算法
07-24
几个基本的几何图形学问题,可以启发一下解决问题的思路,分享给大家
计算机图形学_中点画圆算法
03-20
中点画圆算法的matlab实现代码
种子填充算法c语言代码实现,OpenGL绘图实例三之种子填充算法
weixin_32536777的博客
05-23 2169
综述博主研究了一下午加一晚上，终于把种子填充算法实现出来并把机器人填充完毕，路途很艰辛，不过也学到了很多，在此和大家一起分享。吐槽与我不是同学的小伙伴，请自动忽略，我是来吐槽教材的。 在此不得不吐槽一下，不得不说教材实在太坑爹了。对于种子填充算法的后半部分，下一个种子点的寻找过程中，从while(x<=xright)开始，我实在无法搞懂它里面的神逻辑，最初我认为它是对的，后来按照它的思路实现...
opengl填充编程以及案例（机器人）种子填充、polygon
qq_45228525的博客
05-16 1002
1.线框机器人： #include <gl/glut.h> #include <math.h> #define PI 3.1415926 //画矩形，传入的是左下角XY坐标和右上角XY坐标 void glRect(int leftX, int leftY, int rightX, int rightY) { //画封闭曲线 glBegin(GL_LINE_LOOP); //左下角 glVertex2d(leftX, leftY); //右
计算机图形学 算法汇总（期末复习）
hanmo22357的博客
05-22 5940
文章目录一、光栅图形扫描转换算法1.扫描转换一个点1.1.自动取整法1.2.四舍五入法2.扫描转换一条直线2.1.数值微分(DDA)算法2.2.中点画线算法2.3.Bresenham算法3.扫描转换一个圆3.1.中点画圆法3.2.Bresenham画圆法4.扫描转换一个椭圆(Bresenham算法） 一、光栅图形扫描转换算法 1.扫描转换一个点 1.1.自动取整法 算法描述：直接取各个点坐标的整数部分。 算法举例：将点P(1.7,0.8)经过自动取整，转化为点(1,0)。 1.2.四舍五入法 算法描述
计算机图形学复习：光线跟踪算法详解
"光线跟踪算法是计算机图形学中一种用于模拟光在场景中传播和反射的技术，常用于创建逼真的图像。此复习课件主要针对期末复习，涵盖计算机图形学的重要概念和技术，如二维和三维图形生成、观察、对象表示以及真实感...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值