宣纸上墨水扩散模拟.docVIP

宣纸上墨水扩散模拟 　　摘要： 　　墨水扩散是一个复杂的物理现象，针对墨水在宣纸上的扩散模拟问题，提出了一个基于变系数扩散方程的模拟方法，其扩散系数由宣纸结构和随时间而减少的剩余墨水量决定。模拟分为两个阶段：宣纸结构的模拟和墨水扩散动态过程的模拟。为了模拟宣纸结构，使用一种由权重不同且方向随机的直线段均匀分布而成的加权纤维结构。墨水扩散的动态过程由变系数扩散方程来阐述，为了高效地生成扩散图像，使用CrankNicolson数值方法求解墨水扩散方程，并且预计算纤维结构和动态更新扩散图像。与以往类似的模拟方法相比，该方法能够生成更加自然的扩散边界，并有效地解决了边界过于平滑的问题。实验结果表明该方法能够真实地模拟不同宣纸上的墨水扩散效果。 　　关键词： 　　宣纸结构；墨水扩散；扩散模拟；加权纤维结构；变系数扩散方程 　　0引言 　　中国古代文房???宝中的墨和纸是中国传统书法艺术的主要材料和载体，其中宣纸因具有色泽耐久、吸水能力强、纹理清晰、润墨性好等特点被公认为书画创作的上等纸张，许多传世书画艺术珍品都是在宣纸上创作的。独具魅力的宣纸书法艺术引起许多计算机图形领域的研究人员和学者的兴趣，他们致力于研究书法艺术的计算机模拟和生成方法，其中墨水扩散模拟是重要的基础工作，涉及宣纸结构和墨水扩散特性两方面的内容。 　　纸张是由纸浆和纤维丝构成，要完整地模拟纸张，涉及数量众多的纤维丝和纸浆，空间组合结构复杂，运算量很大，因此现有的方法主要是通过构建抽象纤维结构来模拟纸张及其物理特性。1991年，Guo等[1]用显微镜观察纸结构发现纸张中存在许多随机分布的纤维，提出了模拟纸的二维纤维结构模型：用均匀分布的网格将纸平面划分成多个小区域，每个小区域内随机生成大量相互交错的短直线段以模拟纤维丝。此后，Lee[2]改进了Guo等的工作，用正弦曲线代替直线段作为纤维丝的形状；这类模型实现简单，但是扩散过程中粒子的承载单位是方格矩阵，会造成扩散区域过于规则和扩散边界的走样。2005年余斌等[3]提出了呈现正态分布的纤维点和纤维丝组成的纸张模型，将生成的纤维丝???为传输介质，用Wang等[4]提出的“伪布朗运动”模拟墨水中的碳粒子在纤维集合中的运动；为了模拟宣纸结构，该方法需要生成大量的纤维点以及与之关联的纤维走向和长度，其模拟过程需要大量的运算和存储空间。 　　目前，墨水扩散的模拟工作基本上将其视为碳粒子和水粒子分别做自由运动的结果。Small[5]最早提出应用细胞自动机理论仿真水彩画；Curtis等[6]对Small的模型进行了扩充，引入了更复杂的浅水流动仿真模型模拟水彩画的扩散效果；1995年，Kunii等[7]通过仔细分析墨水扩散过程，发现墨水扩散边缘的颜色灰度与碳粒子在水中的浓度存在明显的关系，因而采用线性微分方程来描述碳粒子和水粒子扩散过程中的浓度变化规律。后来他们改进了线性微分方程[8]，以提高墨水扩散的真实感。通过求解偏微分方程得到的边缘一般比较光滑，因此这类方法模拟的扩散边缘形状光顺不够自然。2003年，Way等[9]通过将纸张划分为多层结构，将墨水视为水粒子和碳粒子的集合，将扩散现象视为纸张的吸水性导致水流运动而碰撞碳粒子的过程，此方法能模拟多种墨水扩散的效果，但是由于其运动过程模拟过于复杂，导致其运算效率较低。2006年，Chu等[10]提出了基于改进的LaticceBolzman算法[11]和三层流体渗透模型来模拟水和碳粒子在宣纸上的吸收、渗透和扩散过程，并在文献[12]提出的毛笔模型的基础上，实现了一个电子创作系统，取得了较好的实验结果，但其过程复杂，运算量大。 　　本文首次把宣纸中的纤维视为阻碍墨水扩散的物质，并赋以权重刻画纤维对墨水扩散阻碍作用的大小，提出了模拟宣纸的加权纤维结构，在此基础上，把纤维结构的权重和宣纸上剩余墨水量相结合作为扩散系数，给出了描述墨水扩散的变系数扩散方程，该方程可用稳定的CrankNicolson方法求解，能够在线性时间内从上一帧图像获得随后的扩散图像。本文的方法综合考虑了宣纸结构和墨水属性，能够真实地模拟墨水扩散的颗粒状边缘，克服了类似模拟方法容易产生边界过于平滑的问题。 　　1加权纤维结构 　　宣纸起源于中国唐朝时期，是一种主要的高级书画用纸。因其具有保存方便、经久不脆和不易褪色的特点，而享有“纸寿千年”之誉，一笔落下，墨韵清晰，层次分明。宣纸独特的渗透和润滑性，使其成为中国书法艺术的主要载体。尽管宣纸上的书法作品风格迥异，但字的边缘均有颗粒状现象，如图1所示，而这种现象在比较光洁的普通纸张上不明显或没有。因此要真实地模拟这种扩散现象，必须考虑宣纸的结构及其物理特性。 　　由于原材料和加工方法不同，宣纸品种丰富，但通过使用高倍率的显微镜观察后发现不同品种的宣纸结构基本