CN102831291B - 虚拟绘制中浅色扩散区墨浓度计算方法 - Google Patents

虚拟绘制中浅色扩散区墨浓度计算方法 Download PDF

Info

Publication number
CN102831291B
CN102831291B CN201210235777.7A CN201210235777A CN102831291B CN 102831291 B CN102831291 B CN 102831291B CN 201210235777 A CN201210235777 A CN 201210235777A CN 102831291 B CN102831291 B CN 102831291B
Authority
CN
China
Prior art keywords
diffusion region
ink
dark
black
concentration
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201210235777.7A
Other languages
English (en)
Other versions
CN102831291A (zh
Inventor
侯增选
郭超
孙景华
赵宁
刘文慧
李少武
李万杰
张邦磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN201210235777.7A priority Critical patent/CN102831291B/zh
Publication of CN102831291A publication Critical patent/CN102831291A/zh
Application granted granted Critical
Publication of CN102831291B publication Critical patent/CN102831291B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种虚拟绘制中浅色扩散区墨浓度计算方法,属于虚拟绘制技术领域,扩散区分为深色扩散区和浅色扩散区,将浅色扩散区划分为宽度相等的圆环区,引入墨浓度衰变度参数,给出了浅色扩散区墨浓度计算方法,包括计算深色扩散区墨浓度、计算浅色扩散区墨浓度衰变度、计算浅色扩散区墨浓度等步骤,从而得到浅色扩散区某个圆环区的墨浓度。该方法揭示了浅色扩散区墨浓度变化规律,克服了现有虚绘制方法中浅色扩散区墨浓度呈现均匀分布的局限,增强了虚拟绘制中水墨扩散仿真效果。<!--1-->

Description

虚拟绘制中浅色扩散区墨浓度计算方法
技术领域
本发明涉及虚拟绘制技术领域,特别是涉及虚拟绘制中浅色扩散区墨浓度计算方法。
背景技术
在虚拟绘制中为了达到真实的水墨扩散效果,需要在宣纸模型的基础上进行水墨扩散的仿真。水墨是由微小粒子构成的,它们以宣纸作为载体,在宣纸细胞中进行运动。在国外Kunii等在文献”Tosiyasu L.Kunii,G.V.Nosovskij,T.Hayashi.A Diffusion Model for Computer Animation of Diffuse Ink Painting[J].Computer Animation,1995”中提出水墨扩散是由水扩散开始的,水粒子的扩散带动墨粒子在相邻纸细胞之间运动,并阐述了深色扩散区和浅色扩散区以及扩散边界的产生原理。Zhang等在文献”Zhang Qing,Sato Youetsu,Takahashi JunYa,et al.Simple Cellular Automation-Based Simulation of Ink Behavior and ItsApplication to Suibokuga-Like 3D Rendering of Trees[J].The Journal ofVisualization and Computer Animation,1999,19(1)”中提出一种基于宣纸细胞的自动计算模型来仿真水粒子和墨粒子运动状态的方法。国内的研究成果主要体现在在水墨扩散算法研究方面,提出了吸水力平衡算法、伪布朗运动算法、交换临界值算法等。在此基础上Young等在文献”J.T.Lee.Simulating OrientalBlack-Ink Painting[J].IEEE Computer Graphics and Applications,1999,19(3)”中提出用局部平衡算法和三层宣纸模型来仿真水墨扩散,取得较好的仿真效果。
综合国内外的研究现状发现,在虚拟绘制的水墨扩散仿真中存在浅色扩散区墨浓度呈现均匀分布的局限,无法反映真实的水墨扩散情况,降低了水墨扩散仿真效果。
发明内容
为了克服现有虚绘制方法中浅色扩散区墨浓度呈现均匀分布的局限,本发明提出了虚拟绘制中浅色扩散区墨浓度计算方法。该方法在虚拟绘制中浅色扩散区内引入墨浓度衰变度参数,揭示了浅色扩散区墨浓度变化规律,给出了浅色扩散区墨浓度计算公式,增强了虚拟绘制中水墨扩散仿真效果。
本发明所采用的技术方案是:虚拟绘制中浅色扩散区墨浓度计算方法,包括下述步骤:
(1)计算深色扩散区墨浓度
在虚拟绘制中,水墨扩散区分为深色扩散区和浅色扩散区。深色扩散区扩散范围为圆形区域,半径为L1,单位为纸元边长;浅色扩散区扩散范围为深色扩散区周围的圆环区,其宽度为L2,单位为纸元边长。
深色扩散区内墨浓度用θ表示,在深色扩散区内,墨浓度对显示效果的影响较小。假设深色扩散区的墨浓度呈均匀分布。
θ的计算如公式(1)所示:
θ=M深墨/M                                    (1)
M深墨为深色扩散区墨粒子量,M为深色扩散区水墨粒子量。
(2)引入墨浓度衰变度参数
在浅色扩散区内,墨浓度变化对显示效果的影响较大,其规律是随着扩散范围的扩大,扩散区内墨浓度逐渐降低。墨浓度的衰减是一个连续的过程,在计算机仿真中将浅色扩散区等分为n个宽度为d的圆环区,每个圆环区内部的墨浓度呈均匀分布,相邻两个圆环区的墨浓度产生衰减,引入墨浓度衰变度参数λ,第i个圆环区的墨浓度计算如公式(2)所示:
θ=θ×(1-λ)1                                    (2)其
中:n为正整数,1□i□n。
(3)计算浅色扩散区墨浓度衰变度
设第1圆环区内的墨粒子量为m1,第2圆环区内的墨粒子量为m2,以此类推,第n圆环区内的墨粒子量为mn。m1,m2,……,mn满足如下关系式,其中N为单个纸元面积:
mi=[π(L1+d)2-πL1 2]×i×θ×(1-λ)/N
m2=[π(L1+2d)2-π(L1+d)2]×i×θ×(1-λ)2/N
m3=[π(L1+3d)2-π(L1+2d)2]×i×θ×(1-λ)3/N
mn={π(L1+n×d)2-π[L1+(n-1)×d]2}×i×θ×(1-λ)n/N
由此可以推导出mn  的计算公式为:
n个圆环区墨粒子的总和即为浅色扩散区的墨粒子量M浅墨
m1+m2+……+mn=M浅墨                                (4)
由此可以计算出浅色扩散区的墨浓度衰变度λ。
(4)计算浅色扩散区墨浓度
将λ的值代入公式(2),即可求得浅色扩散区第i个圆环区内的墨浓度计算公式:
θ=θ×(1-λ)1
本发明的有益效果是:
1、本发明在虚拟绘制中浅色扩散区内引入墨浓度衰变度参数,揭示了浅色扩散区墨浓度变化规律,给出了浅色扩散区墨浓度计算公式,克服了现有虚拟绘制方法中浅色扩散区墨浓度呈现均匀分布的局限;
2、本发明增强了虚拟绘制中水墨扩散仿真效果。由图2中看出浅色扩散区的墨浓度逐渐衰减。
附图说明
图1是深色扩散区和浅色扩散区的扩散范围示意图。
图中:1深色扩散区;2浅色扩散区;3深色扩散区边界;4浅色扩散区边界。
图2是虚拟绘制中浅色扩散区示意图。
图3是虚拟绘制中浅色扩散区墨浓度计算流程图。
具体实施方式
本具体实施例中以虚拟绘制中一个点的绘制为例说明浅色扩散区墨浓度计算方法。本具体实施例中深色扩散区扩散半径L1=16纸元边长,浅色扩散区的宽度L2=2纸元边长,宣纸模型吸储度i=3。深色扩散区墨粒子量M深墨=2628墨粒子个数,深色扩散区水墨粒子量M=3582水墨粒子个数,浅色扩散区墨粒子量M浅墨=172墨粒子个数。M深墨、M、M浅墨均为相对量,对仿真效果不产生影响。具体过程如下:
1、计算深色扩散区墨浓度,根据公式(1)求得θ=0.73367;
2、计算浅色扩散区墨浓度衰变度,本具体实施例中将浅色扩散区等分为n=20圆环区,每个圆环区的宽度根据公式(4)求得λ=0.1065;
3、计算浅色扩散区墨浓度,将λ的值代入公式(2)求得浅色扩散区第i个圆环区内的墨浓度计算公式θ=0.73367×0.8935i,例如浅色扩散区第5个圆环区内(i=5)的墨浓度θ=0.4178。

Claims (1)

1.一种虚拟绘制中浅色扩散区墨浓度计算方法,其特征包括如下步骤:
(1)计算深色扩散区墨浓度
在虚拟绘制中,水墨扩散区分为深色扩散区和浅色扩散区;深色扩散区扩散范围为圆形区域,半径为L1,单位为纸元边长;浅色扩散区扩散范围为深色扩散区周围的圆环区,其宽度为L2,单位为纸元边长;
深色扩散区内墨浓度用θ表示;假设深色扩散区的墨浓度呈均匀分布;
θ的计算如公式(1)所示:
θ=M深墨/M  (1)
M深墨为深色扩散区墨粒子量,M为深色扩散区水墨粒子量;
(2)引入墨浓度衰变度参数
将浅色扩散区等分为n个宽度为d的圆环区,每个圆环区内部的墨浓度呈均匀分布,相邻两个圆环区的墨浓度产生衰减,引入墨浓度衰变度参数λ,第i个圆环区的墨浓度计算如公式(2)所示:
θ=θ×(1-λ)i  (2)其中:n为正整数,1≤i≤n;
(3)计算浅色扩散区的墨浓度衰变度参数
设第1圆环区内的墨粒子量为m1,第2圆环区内的墨粒子量为m2,以此类推,第n圆环区内的墨粒子量为mn;m1,m2,……,mn满足如下关系式,其中N为单个纸元面积;由此推导出mn的计算公式为:
n个圆环区墨粒子的总和即为浅色扩散区的墨粒子量M浅墨
m1+m2+……+mn=M浅墨  (4)
由此计算出浅色扩散区的墨浓度衰变度参数λ;
(4)计算浅色扩散区墨浓度
将λ的值代入公式(2),即可求得浅色扩散区第i个圆环区内的墨浓度计算公式:θ=θ×(1-λ)i
CN201210235777.7A 2012-07-09 2012-07-09 虚拟绘制中浅色扩散区墨浓度计算方法 Expired - Fee Related CN102831291B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210235777.7A CN102831291B (zh) 2012-07-09 2012-07-09 虚拟绘制中浅色扩散区墨浓度计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210235777.7A CN102831291B (zh) 2012-07-09 2012-07-09 虚拟绘制中浅色扩散区墨浓度计算方法

Publications (2)

Publication Number Publication Date
CN102831291A CN102831291A (zh) 2012-12-19
CN102831291B true CN102831291B (zh) 2015-04-29

Family

ID=47334426

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210235777.7A Expired - Fee Related CN102831291B (zh) 2012-07-09 2012-07-09 虚拟绘制中浅色扩散区墨浓度计算方法

Country Status (1)

Country Link
CN (1) CN102831291B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1773516A (zh) * 2005-07-08 2006-05-17 天津大学 基于毛笔和宣纸的水墨传输模型与仿真算法
CN102254333A (zh) * 2011-07-25 2011-11-23 浙江大学 基于图像的水墨画画风图像生成方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1773516A (zh) * 2005-07-08 2006-05-17 天津大学 基于毛笔和宣纸的水墨传输模型与仿真算法
CN102254333A (zh) * 2011-07-25 2011-11-23 浙江大学 基于图像的水墨画画风图像生成方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Physical Modeling Based Graphical Simulator of Water-Ink Diffusion and Multi-Stroke Superposition;Meijun SUN et.al;《2010 international conference on Multimedia Technology(ICMT)》;20101029;全文 *
一种应用GA的墨扩散效果仿真新方法;唐棣等;《中国图像图形学报》;20090331;第14卷(第3期);全文 *
基于纸的物理建模的水墨画扩散效果仿真;余斌;《系统仿真学报》;20050930;第17卷(第9期);全文 *
基于纹理合成的中国山水画系统仿真;孙美君等;《系统仿真学报》;20041031;第16卷(第10期);全文 *

Also Published As

Publication number Publication date
CN102831291A (zh) 2012-12-19

Similar Documents

Publication Publication Date Title
CN103679264B (zh) 基于人工鱼群算法的人群疏散路径规划方法
CN104143027B (zh) 一种基于sph算法的流体热运动仿真系统
CN102915552B (zh) 一种可控的火焰动画生成方法及其系统
CN102592297B (zh) 一种基于3d网格的实时水波纹理处理方法
CN102789650B (zh) 一种基于粒子系统的海面航迹并行化仿真方法
CN111768502A (zh) 一种基于gpu加速技术的非结构网格二维洪水模拟系统
CN102156785A (zh) 一种动态布料仿真方法
CN102831275B (zh) 一种3d流体的仿真方法及系统
CN102147932A (zh) 一种基于可移动欧拉网格的模型驱动的烟雾模拟方法
CN106934192A (zh) 一种参数优化的浅水方程模型水体建模方法
CN109948109A (zh) 含截面变化的明渠非恒定流无网格粒子模拟方法
CN106407577A (zh) 一种模拟风力机尾流的改进制动面模型建立方法
CN101615300A (zh) 一种屏幕空间微结构表面对象环境光遮挡方法
CN102193837A (zh) 一种基于前线的包围盒碰撞检测方法
CN102831291B (zh) 虚拟绘制中浅色扩散区墨浓度计算方法
CN106530375A (zh) 一种个性化情绪感染的人群动画生成方法
CN102867336B (zh) 一种基于热力学模型的固体燃烧过程模拟方法
CN103412991B (zh) 一种网箱对水流影响的模拟方法
CN102831292B (zh) 一种基于力反馈技术的虚拟水墨扩散范围计算方法
CN103440684B (zh) 一种将球谐光照技术应用于面绘制的方法
CN103886636B (zh) 基于光线投射步进补偿的实时烟雾渲染算法
CN114282403A (zh) 一种耦合生境适宜模型的高效高精度栖息地模拟方法
Xu et al. Physically based simulation of water currents and waves
CN109684649A (zh) 一种基于地形的风速修订方法及系统
CN103632003A (zh) 一种推土铲内动态土壤仿真建模方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150429

Termination date: 20180709

CF01 Termination of patent right due to non-payment of annual fee