《图形处理》课件

《图形处理》课程课件从基础算法到前沿应用·计算机图形学完整教学体系Contents课程目录从图形生成原理到前沿渲染技术，系统构建计算机图形学知识体系。01绪论：计算机图形学概述02基本图形生成原理与算法03图形几何变换与投影04曲线曲面理论与建模05图像处理技术基础06前沿应用与真实感渲染CHAPTER01绪论：计算机图形学概述定义、发展历程与应用领域全景CHAPTER01·基础概念计算机图形学的定义与学科定位计算机图形学是研究利用计算机生成、处理和显示图形的学科，核心是将抽象数据转化为可视化图形。它与图像处理互为逆过程，与计算机视觉、CAD/CAM等学科深度交叉，是连接数学理论与视觉应用的桥梁学科。学科定义研究如何通过计算机将数字数据转化为可视图形输出，涉及建模、渲染、交互等完整流程。区别于图像处理与计算机视觉，图形学核心方向是"从数据到图"。从数据到图交叉学科关系与计算几何共享曲线曲面理论基础，矩阵运算是图形变换的数学核心。与CAD/CAM形成上下游：图形学提供显示与交互技术，CAD提供工程建模需求。矩阵运算核心研究内容基础层涵盖图形生成算法、几何变换、裁剪与填充；高级层包括真实感渲染与物理模拟。交互技术贯穿始终，从早期光笔到现代VR手势识别持续推动学科发展。渲染·模拟ComputerGraphics·Timeline计算机图形学发展历程从1950年代CRT简单绘图到实时光线追踪，硬件能力的跃迁始终是图形学发展的核心驱动力。早期CRT显示设备·计算机历史博物馆藏品1950S·萌芽期MIT旋风计算机首次使用CRT输出图形，1963年IvanSutherland开发Sketchpad系统奠定交互式图形学基础1960–70S·奠基期Bresenham直线算法、Cohen-Sutherland裁剪算法等经典算法提出，光栅扫描显示器逐步取代矢量显示器1980–90S·成熟期OpenGL与DirectX图形API标准化，SGI工作站推动三维图形工业化，纹理映射与光照模型理论体系完善21STCENTURY·普适期GPU可编程管线架构革命，实时光线追踪进入消费级产品（NVIDIARTX），云渲染与移动端图形技术快速普及ApplicationDomains计算机图形学核心应用领域计算机图形学已从学术实验室走向产业核心，在游戏娱乐、工业设计、医学影像、科学可视化、虚拟现实和数据呈现六大领域发挥关键作用，其技术渗透率随GPU算力提升和实时渲染能力增强而持续扩大。娱乐与创意产业游戏引擎实现实时光线追踪，电影特效依赖物理渲染与粒子模拟技术从手绘到全CG渲染，图形学算法的进步直接决定视觉表现力的上限Unity·Unreal工业与工程设计CAD/CAM系统将图形学与工程计算结合，支持概念设计到数控加工全流程数字孪生利用三维重建与实时渲染，实现物理世界的虚拟镜像数字孪生医学与科学领域CT/MRI三维重建将断层扫描数据转化为可交互三维解剖模型，辅助手术规划科学计算可视化将流体力学、气象模拟等抽象数据场转化为直观图形CT·MRIGraphicsSystem计算机图形系统架构与核心组件现代图形系统采用CPU+GPU异构架构，通过管线流水线实现高效图形生成；帧缓冲区作为连接计算与显示的核心桥梁，直接影响画面的分辨率与流畅度。01图形管线三阶段：应用阶段（CPU处理场景逻辑）、几何阶段（顶点变换与裁剪）和光栅化阶段（像素着色与输出）02帧缓冲区双缓冲：存储每个像素的颜色与深度值，前后缓冲区交替避免画面撕裂03GPU并行架构：数千个流处理器，专为同时处理百万级像素的矩阵运算和着色器程序优化04图形API标准化：OpenGL/DirectX/Vulkan屏蔽底层硬件差异，实现跨平台运行NVIDIAGPU—大规模并行流处理器架构，专为图形管线与像素着色优化Chapter02基本图形生成原理与算法直线、圆、多边形的扫描转换与区域填充Rasterization·LineDrawing直线扫描转换算法直线扫描转换是将数学上的连续直线映射到离散像素网格的过程。DDA算法直观但依赖浮点运算，Bresenham算法通过整数误差项累加完全消除浮点运算，成为光栅化直线的标准算法。DDA数值微分算法根据直线斜率选择主步进方向，每步在主方向递增1个像素，另一方向按斜率增量累加并取整确定像素位置。实现简单直观，易于理解其几何原理。每步涉及浮点加法与取整运算，在早期硬件上计算开销较大，不适合实时批量绘制场景。精度累积误差可能影响长线绘制质量。核心特征：浮点累加运算Bresenham中点画线算法引入决策变量d，通过判断中点与直线的位置关系确定下一像素，仅使用整数加减和比较操作，完全无需浮点运算。初始决策变量和增量值可预先计算，循环体内仅需一次加法和一次条件判断，效率比DDA提升约3-5倍，成为工业标准。核心特征：整数误差项累加SCANCONVERSION圆的扫描转换：中点画圆法圆的扫描转换利用圆的八对称性，仅需计算1/8圆弧即可推导完整圆周像素位置。中点画圆法通过构造判别式d判断中点相对圆的位置来决定下一像素，全程整数运算，是图形学中"利用几何对称性降低计算复杂度"的经典范例。01八对称性原理圆上一点(x,y)可通过坐标变换映射到其余七个对称点，仅需计算第一象限中从(0,R)到(R/√2,R/√2)的弧段1/8弧段02判别式构造取候选像素中点M(x+1,y-0.5)代入圆方程F(x,y)=x²+y²-R²，根据d=F(M)的正负决定取上方还是下方像素d=F(M)03增量递推优化d的更新仅涉及整数加法（d_new=d+2x+3或d+2x-2y+5），避免了每步重新计算圆方程的乘法开销纯整数运算04算法推广同样的中点判别思想可扩展到椭圆生成，但因椭圆不具有八对称性需分两段处理且增量公式更复杂椭圆扩展REGIONFILLING多边形区域填充算法区域填充是将封闭边界内部像素着色的过程。种子填充适合交互式点选但受递归深度限制，扫描线填充通过活化边表实现高效逐行光栅化，两者各有优劣，按场景选用。种子填充·基本原理从指定种子点出发，按四连通或八连通规则向邻域递归扩展，逐像素着色直到遇到边界颜色停止。4/8-Connected种子填充·局限与改进递归深度与区域面积成正比，大面积区域可能导致栈溢出；改进方案使用扫描线种子填充减少递归层数。StackOverflowRisk扫描线填充·核心机制逐条扫描线处理，利用边的连贯性维护活化边表（AET），仅需对跨越当前扫描线的边做求交和排序。ActiveEdgeTable扫描线填充·优化策略增量更新交点x坐标（x_new=x_old+1/m）避免重复求交，配合奇偶或非零环绕规则判断填充区间。ParityRuleGRAPHICSPRIMITIVES线宽处理与线型控制线宽和线型是图形生成的实用细节，直接影响工程制图的可读性。线宽处理需在垂直于线条方向扩展像素宽度，线型控制则通过模式掩码实现虚线、点划线等样式。线宽处理方法沿线条法线方向复制像素行适用于水平与垂直线，或用圆形、方形笔刷沿路径逐步移动以支持任意角度线条绘制。法线扩展斜线粗细问题由于像素网格离散性，45°斜线的视觉宽度会比水平线细约30%，需根据斜率动态调整像素扩展量来补偿。≈30%线型实现机制定义位模式掩码（如虚线为11110000循环），沿线条路径逐像素匹配掩码位决定是否绘制，支持虚线、点线等样式。位掩码字符生成基础点阵字体通过位图模板映射像素，矢量字体通过贝塞尔曲线轮廓光栅化，TrueType是两者的融合方案。TrueTypeAlgorithmComparison二维图形裁剪算法图形裁剪是在显示前过滤掉窗口外不可见部分的过程，直接影响渲染效率。Cohen-Sutherland算法通过区域编码快速排除可见/不可见线段，Liang-Barsky算法将裁剪转化为一维参数区间求交，两者各有优势。RegionCodingCohen-Sutherland编码裁剪01将裁剪窗口周围区域按上下左右四个方向编码为4位二进制数，端点编码的按位与非零则线段完全在窗外可立即丢弃02按位或为零则线段完全在窗内直接保留；否则选取窗外端点，计算与窗口边界的交点后替换端点并重复判断4-BITREGIONCODEParametricLiang-Barsky参数化裁剪01将线段表示为参数方程P(t)=P₀+t(P₁−P₀)，t∈[0,1]，把四条边界约束转化为关于t的不等式组，求解t的有效区间02只需最多四次除法和比较即可确定可见段，比Cohen-Sutherland减少了求交次数，对大量线段的批量裁剪效率提升显著O(4)PARAMETRICCLIPCHAPTER03图形几何变换与投影从二维变换到三维投影的数学基础与矩阵实现GEOMETRICTRANSFORMS二维基本几何变换二维基本变换包含平移、旋转和缩放三种操作，均可表示为矩阵乘法形式。引入齐次坐标后，原本需要加法表示的平移也能统一为矩阵乘法，使复合变换可以通过矩阵连乘实现，奠定了现代图形管线的数学基础。平移变换将点(x,y)移动到(x+tx,y+ty)，在齐次坐标下表示为3×3矩阵，tx和ty位于矩阵第三列。3×3矩阵缩放变换各坐标分量乘以比例因子(sx,sy)，相等时为等比缩放，不等时为拉伸变形，缩放中心默认为原点。sx,sy旋转变换绕原点旋转θ角，新坐标由三角函数线性组合得出，旋转矩阵行列式恒为1，保持面积不变。det=1复合变换通过矩阵连乘实现：先平移后旋转等价于组合矩阵M=R·T，变换顺序不同结果不同（矩阵乘法不交换）。M=R·THomogeneousCoordinates齐次坐标与变换统一表示齐次坐标通过增加一个维度将n维空间点表示为n+1维向量，使平移变换从加法运算升级为矩阵乘法，实现了平移、旋转、缩放三类变换的统一矩阵表示。这一数学工具是现代图形管线所有变换操作的理论基础，被OpenGL和DirectX等标准深度集成。01二维点(x,y)的齐次坐标表示为(x,y,1)，变换矩阵为3×3；三维点(x,y,z)的齐次坐标为(x,y,z,1)，变换矩阵为4×4。这种统一的维度扩展使得所有变换操作都能在相同的数学框架下进行。3×3/4×4MatrixDimensionExtension02统一性优势：平移矩阵T、旋转矩阵R、缩放矩阵S均为同维度方阵，复合变换只需矩阵连乘M=S·R·T即可一步完成。这种表示方法消除了传统方法中向量加法与矩阵乘法的混合运算，大幅简化了图形管线中的变换计算流程。M=S·R·TUnifiedTransformation03缩放不变性：(x,y,w)与(x/w,y/w,1)表示同一个二维点，w分量编码深度信息，这一性质在透视投影中至关重要。齐次坐标的这一特性使得透视除法能够自然地融入图形渲染管线，实现了从观察空间到裁剪空间的平滑过渡。w=DepthPerspectiveDivision04逆变换便捷性：平移逆矩阵取反平移量，旋转逆矩阵等于转置（正交矩阵性质），复合变换逆矩阵为各矩阵逆的反序乘积。这种结构化的逆变换计算使得相机变换、法线变换等图形学核心操作能够高效实现。OrthogonalPropertyInverseTransformation3DTRANSFORMS三维几何变换与旋转表示三维变换在二维基础上增加了Z轴维度和绕三轴旋转，4×4齐次变换矩阵成为标准工具。三维旋转的表示是工程难点：欧拉角直观但存在万向节锁问题，四元数避免了锁死且支持平滑插值，是游戏引擎和动画系统中旋转插值的首选方案。014×4齐次变换矩阵三维基本变换使用4×4齐次矩阵：平移、缩放结构与二维类似，旋转分为绕X/Y/Z三轴的独立旋转矩阵Rx、Ry、Rz。TranslationRotationScaling02欧拉角与万向节锁欧拉角用三个角度(α,β,γ)描述旋转序列，直观易懂但在俯仰角±90°时出现万向节锁，丢失一个旋转自由度。Roll-Pitch-YawGimbalLock⚠03四元数与Slerp插值四元数q=w+xi+yj+zk用四个分量表示旋转，避免万向节锁且支持Slerp球面线性插值，是动画关键帧平滑过渡的标准方案。QuaternionSlerp✓04观察变换ViewTransform观察变换将相机从世界坐标对齐到标准观察坐标系，本质是对整个场景施加一个逆相机变换矩阵。LookAtViewMatrixProjectionTransform投影变换：平行投影与透视投影投影变换是将三维场景映射到二维视平面的降维过程。平行投影保持平行性和比例关系适合工程制图，透视投影通过近大远小的视觉效果实现真实感。平行投影投影面物体投影方向投影线互相平行，物体大小不随距离变化透视投影视点投影面近远投影线汇聚于灭点，产生近大远小效果Parallel投影线互相平行，物体大小不随距离变化，保持平行性和比例关系，适合需要精确量取尺寸的三视图和轴测图Orthographic正投影是投影线垂直于投影面的特殊情况，三视图（正视、侧视、俯视）均为正投影，是机械制图的标准表达方式Perspective投影线汇聚于观察点（灭点），产生近大远小效果，平行线在远处交汇，符合人眼和相机的成像规律Matrix透视投影矩阵将视锥体变换为标准立方体[-1,1]³，z值经透视除法后编码为深度缓冲所需的非线性深度值CHAPTER04曲线曲面理论与建模Bezier、B样条与NURBS的数学原理与工程应用COMPUTERGRAPHICS曲线基础：插值与逼近曲线设计分为插值（精确过点）和逼近（整体趋近）两大范式。插值方法如拉格朗日多项式在数据点多时出现龙格现象导致剧烈振荡，逼近方法如Bezier曲线通过控制点定义凸包约束保证曲线光滑可控，更适合交互式设计场景。拉格朗日插值要求曲线精确通过所有给定数据点，多项式阶数随点数增长，高阶时边缘出现龙格现象——剧烈振荡RUNGEPHENOMENONBezier曲线曲线不必精确过点，由Bernstein基函数加权控制点生成，始终位于控制点凸包内部CONVEXHULL三次样条插值通过分段三次多项式保证C²连续性，在过点精度与曲线光滑度之间取得平衡C²CONTINUITY交互式设计优势逼近方法控制点少且调整直观，曲线天然光滑无振荡，适合CAD与动画自由形状设计CAD&ANIMATIONMATHEMATICALFOUNDATIONSBezier曲线的数学原理与性质全局控制特性与高阶数值不稳定是Bezier曲线的主要局限，催生了B样条的改进方案。数学定义n次Bezier曲线由Bernstein基函数加权控制点构成，公式表达为P(t)=ΣBi,n(t)·Pi，其中基函数Bi,n(t)=C(n,i)·ti·(1-t)n-i，参数t∈[0,1]实现从起点到终点的连续映射。Bernstein基函数·参数化表示核心性质Bezier曲线具有三大核心性质：端点插值保证曲线必过首末控制点；凸包性确保曲线始终位于控制点构成的凸包内部；仿射不变性意味着对控制点进行仿射变换等价于对曲线直接变换，这一特性在几何建模中至关重要。端点插值·凸包性·仿射不变deCasteljau算法该算法通过n层递推线性插值计算曲线上任意一点，递推公式为Pi(k)=(1-t)·Pi(k-1)+t·Pi+1(k-1)。相比直接计算Bernstein多项式，该几何构造方法具有更好的数值稳定性，也是曲线细分和求交运算的基础。递推线性插值·数值稳定全局控制局限Bezier曲线存在两个根本缺陷：一是全局性控制，修改任一控制点将影响整条曲线形状，无法实现局部编辑；二是高次曲线易产生龙格现象，出现不必要的波动和振荡。因此工程实践中通常将次数限制在3次以内，复杂形状则采用分段拼接策略。≤3次约束·分段拼接方案B-SPLINECURVESB样条曲线：局部控制与高阶连续性B样条曲线通过引入节点向量和分段基函数，实现了控制点的局部支撑性——修改单个控制点仅影响附近k段曲线。同时保持曲线阶数与段数独立，可用大量控制点构建复杂形状而不升高多项式次数，成为CAD系统中最主流的曲线表示方法。B样条基函数与定义Cox-deBoor递推基函数定义在节点向量U={u₀,u₁,...,uₘ}上，k次B样条基函数Nᵢ,ₖ(u)通过Cox-deBoor递推公式逐层计算，具有局部支撑性。B样条基函数与定义区间非零·精确局部编辑每个基函数Nᵢ,ₖ仅在区间[uᵢ,uᵢ₊ₖ₊₁)上非零，因此移动控制点Pᵢ仅影响该区间对应的曲线段，实现精确的局部编辑。对比Bezier的优势局部控制vs全局变化移动一个控制点仅改变附近k段曲线形状，其余部分完全不变；Bezier修改任何点则全局变化。对比Bezier的优势阶数独立·复杂形状构建曲线阶数k由用户选定而非由控制点数决定，可在保持低次多项式的同时用大量控制点构建复杂形状。GEOMETRICMODELING·几何建模核心NURBS曲线与非均匀有理B样条NURBS（非均匀有理B样条）通过为每个控制点引入权重因子wi，在B样条基础上获得了精确表示圆锥曲线的能力。"非均匀"指节点间距可变，"有理"指基函数加权后做有理分式运算。NURBS是ISO工业标准，被CATIA、Rhino等主流CAD系统采用为核心几何表示方法。01NURBS曲线定义P(u)=Σwi·Ni,k(u)·Pi/Σwi·Ni,k(u)，权重wi增大时曲线被拉向对应控制点Pi，wi=0时该点失效wi02精确表示圆锥曲线通过特定权重值（如圆弧端点权重为1、中点为cos(θ/2)）可精确表示圆、椭圆、抛物线，而普通B样条只能近似cos(θ/2)03非均匀节点向量允许节点重复（重节点处曲线降低连续性），k重节点处曲线插值端点，实现Bezier曲线作为NURBS的特例k-fold04工业标准地位STEP国际标准（ISO10303）将NURBS定义为产品几何数据交换的标准格式，是CAD/CAM互操作性的基础保障ISO10303PARAMETRICSURFACES参数曲面与曲面建模方法参数曲面将曲线的单参数t扩展为双参数(u,v)，通过张量积方法将Bezier、B样条和NURBS曲线自然推广到曲面形式。01·基本原理参数曲面基本原理参数曲面P(u,v)将两个参数方向的控制点网格通过张量积基函数加权求和，每个(u,v)对应三维空间中一个点等参数线：固定u变化v得到v方向曲线族，固定v变化u得到u方向曲线族，两组曲线交织形成曲面的网格可视化02·构造与连续性曲面构造与连续性Coons曲面通过四条边界曲线的混合函数构造，Gordon曲面进一步融合两组交叉曲线族，适合从已有边界曲线构造曲面曲面连续性：位置连续G⁰、切线连续G¹、曲率连续G²，汽车A级曲面要求至少G²连续以保证反射光带平滑无断裂NURBS曲面在汽车A级曲面设计中的反射检测CHAPTER05图像处理技术基础图像增强、滤波、边缘检测与频域分析ImageDigitization图像数字化：采样与量化图像数字化通过采样（空间离散化）和量化（幅度离散化）将连续模拟信号转化为数字矩阵。采样率决定空间分辨率，量化位数决定灰度精度。Sampling·采样ΔxΔy空间采样间隔决定像素网格密度。采样不足时高频细节丢失，产生摩尔纹混叠伪影。Quantization·量化256级8位量化产生256级灰度。位数越高灰度越精细，但数据量成倍增长。空间采样01将连续图像在x和y方向等间隔取点，采样间隔Δx和Δy决定像素数量，采样不足时高频细节丢失产生混叠伪影（摩尔纹）02奈奎斯特采样定理：采样频率≥2倍信号最高频率才能无失真重建，实践中通常取2.5–3倍作为安全余量灰度量化03将采样点的连续灰度值映射为有限离散级数，8位量化为256级灰度，10位为1024级，位数越多灰度越精细但数据量成倍增长04量化不足（如仅用4位/16级灰度）会在灰度缓变区域产生伪轮廓效应，可通过抖动技术（dithering）在低位数量化下改善视觉质量SpatialDomainFiltering空间域图像滤波空间域滤波通过卷积核对图像局部邻域进行加权运算，平滑滤波（均值、高斯）抑制噪声但会模糊边缘，锐化滤波（拉普拉斯、Sobel）增强边缘细节但会放大噪声。滤波器的选择本质上是在去噪与保边之间寻找平衡，这是图像处理的核心权衡。常用空间域滤波器对比滤波器类型典型卷积核主要效果局限性均值滤波3×3全1/9均匀平滑去噪边缘模糊严重高斯滤波加权递减核平滑且保留部分边缘大尺度边缘仍模糊中值滤波取邻域中值有效去除椒盐噪声破坏细线条结构拉普拉斯[0,1,0;1,−4,1;0,1,0]各向同性边缘增强放大噪声需配合平滑Sobel算子Gx/Gy方向梯度方向性边缘检测边缘较粗非单像素i不同滤波器在去噪与保边之间各有取舍，实际应用中常组合使用以达到最优效果FrequencyDomainAnalysis频域分析与傅里叶变换频域分析通过傅里叶变换将图像从空间域转换到频率域，低频对应平滑区域，高频对应边缘细节，FFT算法将复杂度从O(N²)降至O(NlogN)。二维离散傅里叶变换将M×N图像f(x,y)变换为频域F(u,v)，幅度谱反映各频率成分强度，相位谱保留空间结构信息DFT·F(u,v)频域滤波流程对图像做DFT→乘以滤波器传递函数H(u,v)→逆DFT回到空间域，等价于空间域卷积运算卷积定理滤波器过渡与振铃理想低通在截止频率处产生振铃效应（Gibbs现象），巴特沃斯和高斯滤波器通过平滑过渡带消除振铃Gibbs现象同态滤波分离处理将照射分量（低频·光照不均）与反射分量（高频·纹理细节）分离，可独立调节实现图像增强照射·反射边缘检测技术边缘检测通过计算图像梯度或二阶导数定位亮度突变位置。Canny检测器通过高斯平滑、非极大值抑制和双阈值追踪实现了细、准、连的最优边缘检测，至今仍是工业视觉系统中最常用的算法。梯度算子与拉普拉斯Sobel梯度算子用3×3卷积核分别计算水平和垂直梯度Gx、Gy，梯度幅值反映边缘强度，方向θ=arctan(Gy/Gx)。Sobel算子计算高效，是实时图像处理的首选方法。3×3卷积核梯度算子与拉普拉斯拉普拉斯二阶算子二阶导数检测灰度零交叉点定位边缘，对噪声更敏感但可检测更精细的纹理变化。常与高斯滤波组合使用形成LoG算子。零交叉点Canny最优边缘检测器四步检测流程高斯平滑去噪→Sobel计算梯度→非极大值抑制细化为单像素宽边缘→双阈值追踪连接真实边缘。流程严谨，输出质量高。四步流程Canny最优边缘检测器三指标最优折中在信噪比、定位精度和单响应三个指标上达到最优，双阈值机制有效区分真实边缘和噪声。被OpenCV等库作为标准实现。工业首选CHAPTER06前沿应用与真实感渲染光线追踪、全局光照、物理模拟与科学可视化RenderingPipeline真实感渲染与光照模型真实感渲染从Phong经验模型发展到基于物理的渲染(PBR)再到路径追踪，核心目标是更精确地模拟光线与材质的物理交互。01Phong光照模型：将表面光照分解为环境光、漫反射（Lambert余弦定律）和镜面反射（高光指数衰减）三项叠加，计算高效但物理精度有限02PBR（基于物理的渲染）：采用微表面模型描述粗糙表面，用BRDF双向反射分布函数和能量守恒约束确保材质在不同光照条件下表现一致03光线追踪：从相机出发反向追踪光线，在表面交点处递归生成反射、折射、阴影光线，自然处理全局光照效果04实时光追里程碑：NVIDIARTX系列引入RTCore硬件加速光线-三角形求交，结合DLSS深度学习超采样，在60fps下实现消费级实时光线追踪光线追踪渲染的真实感效果·全局光照与物理材质交互VolumeVisualization科学计算可视化技术科学计算可视化将仿真和观测产生的海量数值数据转化为直观图形，帮助科学家理解复杂现象。体绘制直接渲染三维数据场无需提取表面，MarchingCubes算法从标量场中提取等值面生成三角网格。随着数据规模爆炸式增长，可视化技术在