百度360必应搜狗淘宝本站头条
当前位置:网站首页 > 热门文章 > 正文

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

bigegpt 2024-08-30 12:52 6 浏览

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

---- 引言----

每次使用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);
    }
}

原文:https://www.cnblogs.com/chinamming/p/3427442.html

相关推荐

Java 泛型大揭秘:类型参数、通配符与最佳实践

引言在编程世界中,代码的可重用性和可维护性是至关重要的。为了实现这些目标,Java5引入了一种名为泛型(Generics)的强大功能。本文将详细介绍Java泛型的概念、优势和局限性,以及如何在...

K8s 的标签与选择器:流畅运维的秘诀

在Kubernetes的世界里,**标签(Label)和选择器(Selector)**并不是最炫酷的技术,但却是贯穿整个集群管理与运维流程的核心机制。正是它们让复杂的资源调度、查询、自动化运维变得...

哈希Hash算法:原理、应用(哈希算法 知乎)

原作者:Linux教程,原文地址:「链接」什么是哈希算法?哈希算法(HashAlgorithm),又称为散列算法或杂凑算法,是一种将任意长度的数据输入转换为固定长度输出值的数学函数。其输出结果通常被...

C#学习:基于LLM的简历评估程序(c# 简历)

前言在pocketflow的例子中看到了一个基于LLM的简历评估程序的例子,感觉还挺好玩的,为了练习一下C#,我最近使用C#重写了一个。准备不同的简历:image-20250528183949844查...

55顺位,砍41+14+3!季后赛也成得分王,难道他也是一名球星?

雷霆队最不可思议的新星:一个55号秀的疯狂逆袭!你是不是也觉得NBA最底层的55号秀,就只能当饮水机管理员?今年的55号秀阿龙·威金斯恐怕要打破你的认知了!常规赛阶段,这位二轮秀就像开了窍的天才,直接...

5分钟读懂C#字典对象(c# 字典获取值)

什么是字典对象在C#中,使用Dictionary类来管理由键值对组成的集合,这类集合被称为字典。字典最大的特点就是能够根据键来快速查找集合中的值,其键的定义不能重复,具有唯一性,相当于数组索引值,字典...

c#窗体传值(c# 跨窗体传递数据)

在WinForm编程中我们经常需要进行俩个窗体间的传值。下面我给出了两种方法,来实现传值一、在输入数据的界面中定义一个属性,供接受数据的窗体使用1、子窗体usingSystem;usingSyst...

C#入门篇章—委托(c#委托的理解)

C#委托1.委托的定义和使用委托的作用:如果要把方法作为函数来进行传递的话,就要用到委托。委托是一个类型,这个类型可以赋值一个方法的引用。C#的委托通过delegate关键字来声明。声明委托的...

C#.NET in、out、ref详解(c#.net framework)

简介在C#中,in、ref和out是用于修改方法参数传递方式的关键字,它们决定了参数是按值传递还是按引用传递,以及参数是否必须在传递前初始化。基本语义对比修饰符传递方式可读写性必须初始化调用...

C#广义表(广义表headtail)

在C#中,广义表(GeneralizedList)是一种特殊的数据结构,它是线性表的推广。广义表可以包含单个元素(称为原子),也可以包含另一个广义表(称为子表)。以下是一个简单的C#广义表示例代...

「C#.NET 拾遗补漏」04:你必须知道的反射

阅读本文大概需要3分钟。通常,反射用于动态获取对象的类型、属性和方法等信息。今天带你玩转反射,来汇总一下反射的各种常见操作,捡漏看看有没有你不知道的。获取类型的成员Type类的GetMembe...

C#启动外部程序的问题(c#怎么启动)

IT&OT的深度融合是智能制造的基石。本公众号将聚焦于PLC编程与上位机开发。除理论知识外,也会结合我们团队在开发过程中遇到的具体问题介绍一些项目经验。在使用C#开发上位机时,有时会需要启动外部的一些...

全网最狠C#面试拷问:这20道题没答出来,别说你懂.NET!

在竞争激烈的C#开发岗位求职过程中,面试是必经的一道关卡。而一场高质量的面试,不仅能筛选出真正掌握C#和.NET技术精髓的人才,也能让求职者对自身技术水平有更清晰的认知。今天,就为大家精心准备了20道...

C#匿名方法(c#匿名方法与匿名类)

C#中的匿名方法是一种没有名称只有主体的方法,它提供了一种传递代码块作为委托参数的技术。以下是关于C#匿名方法的一些重要特点和用法:特点省略参数列表:使用匿名方法可省略参数列表,这意味着匿名方法...

C# Windows窗体(.Net Framework)知识总结

Windows窗体可大致分为Form窗体和MDI窗体,Form窗体没什么好细说的,知识点总结都在思维导图里面了,下文将围绕MDI窗体来讲述。MDI(MultipleDocumentInterfac...