使用C语言实现二维,三维绘图算法(2)-解析曲面的显示

最新推荐文章于 2024-09-14 13:42:52 发布
最新推荐文章于 2024-09-14 13:42:52 发布
阅读量621 收藏
点赞数
文章标签： c/c++

---- 引言----

每次使用OpenGL或DirectX写三维程序的时候, 都有一种隔靴搔痒的感觉, 对于内部的三维算法的实现不甚了解. 其实想想, Win32中既然存在画线画点函数, 利用计算机图形学的知识, 我们用可以用纯C调用Win32实现三维绘图, 完全不用借助OpenGL和DirectX, 这有重复造轮子的嫌疑, 但是自己动手实现一遍, 毕竟也是有意义的.

[效果演示]

原始效果(100条浮动曲线)

加密以后的效果(200条浮动曲线)

[浮动水平线法绘图过程]

固定一个y值按步长变换给定一个x值, 从而可计算出平面截线一个点的z坐标值. 将改点投影到xoy平面上, 然后再变换到屏幕上. 如果是曲线端点要填充边界值. 接着检验此点的可见性,并用1表示上方可见, 0表示不可见, -1表示下方可见. 可见性检测就是用当前点的y值与上下浮动水平线数组中相应的元素值进行比较,y值大于上水平线数组中元素值或小于下水平线数组中元素值, 则当前点可见, 否则不可见. 往下再计算同一平面截线的另一点, 和上面点一样, 先投影到坐标平面上, 再变换到屏幕上. 先前的点叫紧前点, 当前的点为当前点. 紧前点和当前点的可见性主要有下面一些可能情形:

[编程实现要点]

曲面函数的定义

float SurfaceFun(float X, float Y)
{
    float w1, w2, w3, FV;
    w1=4*(X-2)*(X-2) + (Y-4)*(Y-4) - 1;
    w2=(X-5)*(X-5)/9 + 4*(Y-2)*(Y-2) - 1;
    w3=(X-5)*(X-5)/9 + 4*(Y-6)*(Y-6) - 1;
    
    if(w1>85) w1=85;
    if(w2>85) w2=85;
    if(w3>85) w3=85;    
    
    FV=w1*w1*exp(-w1) + w2*w2*exp(-w2) + w3*w3*exp(-w3);
    
    return(FV);
    
}


绘制曲面函数

void DrawSurface()
{
    int Xe, Ye, Ln, Pt, XPre, YPre, XCur, YCur, Xi, Yi;
    int *pi, LimY, VisCur, VisPre;
    float X, Y, Z;
    
    LimY=GetWindowHeight();
    SetLineColor(BLUE); 

    for(Ln=0; Ln<=LNo; ++Ln)
    {
        Y=Y2-Ln*IncY;
        X=X1;
        Z=SurfaceFun(X,Y);
        CalcuProject(X, Y, Z);
        XPre = 0.5 + (XProj-F1)*EchX + C1;
        YPre = 0.5 + (YProj-F3)*EchY + C3;
        
        FillEdge(XPre, YPre, Xd, Yd);
        
        VisPre = VisibilityTest(XPre, YPre);
        
        for(Pt=0; Pt<=PNo; ++Pt)
        {
            X=X1+Pt*IncX;
            Z=SurfaceFun(X,Y);
            CalcuProject(X, Y, Z);
            XCur = 0.5 + (XProj-F1)*EchX + C1;
            YCur = 0.5 + (YProj-F3)*EchY + C3;
            VisCur = VisibilityTest(XCur, YCur);
            
            if( (HMax[XCur]==0) || (HMin[XCur]==LimY) ) VisCur = VisPre;
            
            if(VisCur == VisPre)
            {
                if( (VisCur==1) || (VisCur==-1) )
                {
                    if(0<=XCur)
                        PlotLine(XPre, LimY-60-YPre, Xi, LimY-60-YCur);
                    else if(0<=YCur)
                        PlotLine(XPre, LimY-60-YPre, XPre, LimY-60-YCur);
                    else
                        PlotLine(Xi, LimY-60-YPre, XPre, LimY-60-YPre);
                    
                    HorizonInc(XPre, YPre, XCur, YCur);
                }
            }
            else    // VisCur!=VisPre
            {
                if(VisCur==0)
                {
                    if(VisPre == 1)
                    {
                        pi = Inter(XPre, YPre, XCur, YCur, HMax);
                        Xi = *pi;
                        Yi = *(pi+1);
                    }
                    else
                    {
                        pi = Inter(XPre, YPre, XCur, YCur, HMin);
                        Xi = *pi;
                        Yi = *(pi+1);
                    }
                    
                    if(0<=Xi)
                        PlotLine(XPre, LimY-60-YPre, Xi, LimY-60-Yi);
                    else if(0<=Yi)
                        PlotLine(XPre, LimY-60-Yi, XPre, LimY-60-Yi);
                    else
                        PlotLine(XPre, LimY-60-YPre, XPre, LimY-60-YPre);
                    
                    HorizonInc(XPre, YPre, Xi, Yi);
                }
                else    
                {
                    if(VisCur == 1)
                    {
                        if(VisPre == 0)  
                        {
                            pi = Inter(XPre, YPre, XCur, YCur, HMax);
                            Xi = *pi;
                            Yi = *(pi+1);  
                            
                            if(0<=Xi)
                                PlotLine(Xi, LimY-60-Yi, XCur, LimY-60-YCur);       
                            else if(0<=Yi)
                                PlotLine(XCur, LimY-60-YCur, XCur, LimY-60-YCur);
                            else
                                PlotLine(XCur, LimY-60-YCur, XCur, LimY-60-YCur);
                            
                            HorizonInc(Xi, Yi, XCur, YCur);   
                        }
                        else   
                        {
                            pi = Inter(XPre, YPre, XCur, YCur, HMin);
                            Xi = *pi;
                            Yi = *(pi+1);  
                            
                            if(0<=Xi)
                                PlotLine(XPre, LimY-60-YPre, Xi, LimY-60-Yi);       
                            else if(0<=Yi)
                                PlotLine(XPre, LimY-60-YPre, XPre, LimY-60-Yi);
                            else
                                PlotLine(XPre, LimY-60-YPre, XPre, LimY-60-YPre);
                            
                            HorizonInc(XPre, YPre, Xi, Yi);    
                            
                            pi = Inter(XPre, YPre, XCur, YCur, HMax);
                            Xi = *pi;
                            Yi = *(pi+1);  
                            
                            if(0<=Xi)
                                PlotLine(Xi, LimY-60-YCur, XCur, LimY-60-YCur);       
                            else if(0<=Yi)
                                PlotLine(XCur, LimY-60-Yi, XCur, LimY-60-YCur);
                            else
                                PlotLine(XCur, LimY-60-YCur, XCur, LimY-60-YCur);
                            
                            HorizonInc(Xi, Yi, XCur, YCur);
                        }
                    }
                    else    // VisCur!=0, VisCur!=1
                    {
                        if(VisPre == 0)  
                        {
                            pi = Inter(XPre, YPre, XCur, YCur, HMin);
                            Xi = *pi;
                            Yi = *(pi+1);  
                            
                            if(0<=Xi)
                                PlotLine(Xi, LimY-60-YCur, XCur, LimY-60-YCur); 
                            else if(0<=Yi)
                                PlotLine(XCur, LimY-60-Yi, XCur, LimY-60-YCur);
                            else
                                PlotLine(XCur, LimY-60-YCur, XCur, LimY-60-YCur);
                            
                            HorizonInc(Xi, Yi, XCur, YCur);   
                        }
                        else    // VisCur!=0, VisCur!=1, VisPre!=0
                        {
                            pi = Inter(XPre, YPre, XCur, YCur, HMax);
                            Xi = *pi;
                            Yi = *(pi+1);  
                            
                            if(0<=Xi)
                                PlotLine(XPre, LimY-60-YPre, Xi, LimY-60-Yi);      
                            else if(0<=Yi)
                                PlotLine(XPre, LimY-60-YPre, XPre, LimY-60-Yi);
                            else
                                PlotLine(XPre, LimY-60-YPre, XPre, LimY-60-YPre);
                            
                            HorizonInc(XPre, YPre, Xi, Yi);    
                            
                            pi = Inter(XPre, YPre, XCur, YCur, HMin);
                            Xi = *pi;
                            Yi = *(pi+1);  
                            
                            if(0<=Xi)
                                PlotLine(Xi, LimY-60-Yi, XCur, LimY-60-YCur);      
                            else
                                PlotLine(XCur, LimY-60-YCur, XCur, LimY-60-YCur);
                            
                            HorizonInc(Xi, Yi, XCur, YCur);
                        }                        
                    }
                }
            }
            
            VisPre = VisCur;
            XPre   = XCur;
            YPre   = YCur;
        }
        
        FillEdge(XCur, YCur, Xg, Yg);
    }
}


版权声明：本文为博主原创文章，未经博主允许不得转载。

Chinamming
关注
二维插值
xiaoyuxianshenging的博客
08-30 1万+
二者均是常用的二维差值方法，两者的区别是，interp2的插值数据必须是矩形域，即已知数据点（x,y)组成规则的矩阵，或称之为栅格，可使用meshgid生成。而griddata函数的已知数据点（X，Y）不要求规则排列，特别是对试验中随机没有规律采取的数据进行插值具有很好的效果。griddata(X,Y,XI,YI,'v4') v4是一种插值算法，没有具体的名字，原文称为“MATLAB 4 gr
使用C语言实现二维,三维绘图算法(1)-透视投影
Chinamming的专栏
11-15 1074
---- 引言---- 每次使用OpenGL或DirectX写三维程序的时候, 都有一种隔靴搔痒的感觉, 对于内部的三维算法的实现不甚了解. 其实想想, Win32中既然存在画线画点函数, 利用计算机图形学的知识, 我们用可以用纯C调用Win32实现三维绘图, 完全不用借助OpenGL和DirectX, 这有重复造轮子的嫌疑, 但是自己动手实现一遍, 毕竟也是有意义的. [效果演示] 线框效...
基于C#实现的三维显示
03-18
本程序可以实现三维的显示，放大，缩小，翻转等操作，并且还可自行添加贴图，只需将其中的图片路径修改一下。
OpenGL数据可视化（5）绘制三维曲面
最新发布
qq_44741020的博客
09-14 1203
介绍了如何使用OpenGL绘制三维曲面，实现类似MATLAB中surface函数的效果，涉及GL_TRIANGLE_STRIP三角条带、GL_PRIMITIVE_RESTART图元重启等。
3D图象算法
waterflier's BLOG
12-20 3852
 转载自http://www.gameres.com3D简介　　我们首先从坐标系统开始。你也许知道在2D里我们经常使用Ren?笛卡儿坐标系统在平面上来识别点。我们使用二维(X,Y):X表示水平轴坐标,Y表示纵轴坐标。在3维坐标系，我们增加了Z，一般用它来表示深度。所以为表示三维坐标系的一个点，我们用三个参数(X,Y,Z)。这里有不同的笛卡儿三维系统可以使用。但是它们都是左手螺旋或右手螺旋的
算法-标准绘图(StdDraw)
热门推荐
大白
08-07 1万+
标准会图库很简单：我们可以将它想象为一个抽象的能够在二维画布上画出点和直线的绘图设备。这个设备能够调用StdDraw中的静态方法画出一些基本的几何图形，包括点、直线、圆、多边形等。 使用标准绘图库的话，需要把algs4.jar放到我工程的依赖中，我存到了H盘 package base.第一章; import edu.princeton.cs.algs4.StdDraw; i...
使用C#三维图形控件进行曲线曲面分析
weixin_30918415的博客
04-09 1403
使用AnyCAD.Net三维图图形控件能够计算曲线的切线、法线、曲率、长度等，能够计算曲面的uv切线、法线、面积等。 代码示例一：曲线分析 Platform.LineStyle lineStyle = new Platform.LineStyle(); lineStyle.SetLineWidth(0.5f); ...
c语言编程b样条曲线,B样条曲线曲面和NURBS曲线曲面C语言算法源程序
weixin_35327395的博客
05-18 926
《B样条曲线曲面和NURBS曲线曲面C语言算法源程序》由会员分享，可在线阅读，更多相关《B样条曲线曲面和NURBS曲线曲面C语言算法源程序(27页珍藏版)》请在人人文库网上搜索。1、学习小结：前面学习了Bezier曲线，B样条基函数和B样条曲线的一些基础知识。掌握关键问题是一条B样条曲线间的多段曲线的光滑连接。因为现在是用多段Bezier曲线来描绘一条B样条曲线，所以问题变为两段Bezier曲线间...
C语言实现三维飞机模型：OpenGL光照与Nurbs曲面
本资源是一份关于使用C语言编写的飞机游戏模型文档，重点介绍了如何通过三维立体图来实现飞机的屏幕绘图。文档以OpenGL函数为基础，探讨了光照在游戏中的应用，以及双二次Nurbs曲面在创建飞机模型表面数据中的作用。...
C 代码 定义一些 1D、2D 和 3D 测试功能 用于检测边缘或不连续性， 并使用 gnuplot（） 显示曲线或曲面.rar
05-27
在这个名为"C代码 定义一些 1D、2D 和 3D 测试功能 用于检测边缘或不连续性，并使用gnuplot（）显示曲线或曲面.rar"的压缩包中，包含的是一组C语言编写的程序，它们专注于检测图像或数据中的边缘和不连续性。...
C语言游戏模型（飞机）编程实现.docx
04-12
通过上述内容，我们可以看到如何使用C语言结合OpenGL库来开发一款具有三维立体效果的飞机游戏。从光照到Nurbs曲面的实现，每一个环节都是实现高质量游戏画面不可或缺的部分。对于希望深入了解游戏开发尤其是图形渲染...
三维图形几何变换算法.
10-28
三维图形几何变换是二维图形几何变换的扩展。在三维空间中，用规范化齐次坐标[x y z 1]表示三维点，变换原理是把齐次坐标点(x, y, z, 1)通过变换矩阵变换成新的齐次坐标点(x’ y’, z’, 1)
c#利用direct库显示3d图控件
07-20
可以显示3d图像，支持鼠标移动视角，自带基础图形绘制，不会使用可以联系qq971321447
C语言-三维图形变换
12-24
C语言-三维图形变换
浮动水平线算法源代码
12-29
这是用浮动水平线算法实现的一个球体模型，可以进行各项参数设置
C语言实现边缘检测与连续性测试工具及可视化
资源摘要信息: 本资源提供了用C语言编写的代码示例，用于实现一维(1D)、二维(2D)和三维(3D)的边缘检测或不连续性检测功能，并能够借助gnuplot工具展示检测结果的曲线或曲面。这些功能在图像处理、计算机视觉、信号...
c# 3d图像显示
weixin_30246221的博客
08-23 2510
1.应用directx图形库进行开发； 2.代码： public class TClass : System.Windows.Forms.Form { /// <summary> /// 设备对象，场景中所有图形对象的父对象 /// </summary> private De...
图形算法资源
weixin_30909575的博客
11-19 146
Computer Graphics: http://kesen.realtimerendering.com Computer Vision: http://www.cvpapers.com http://openaccess.thecvf.com GitXiv: http://www.gitxiv.com/ Prof. Olga Sorkine-Hornung: ...
【C#】21. 输出Excel曲面图
huiwuhuiwu的专栏
03-18 3614
VS 与 Excel的结合！
C语言实现二维数组行列排序算法
"吉林大学计算机科学与技术学院的C语言上机实验题目，涉及二维整型数组的排序问题。实验要求编写程序，实现对数组的每一行和每一列分别进行从小到大的排序，并不断重复这一过程，直到数组不再发生变化。实验报告中...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值