CN102760305B - 采用旋转tin网和距离与夹角双权内插法的三维复杂建模方法 - Google Patents
采用旋转tin网和距离与夹角双权内插法的三维复杂建模方法 Download PDFInfo
- Publication number
- CN102760305B CN102760305B CN201110110179.2A CN201110110179A CN102760305B CN 102760305 B CN102760305 B CN 102760305B CN 201110110179 A CN201110110179 A CN 201110110179A CN 102760305 B CN102760305 B CN 102760305B
- Authority
- CN
- China
- Prior art keywords
- reference mark
- point
- triangle
- adjacent
- basis
- 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
Links
Landscapes
- Image Generation (AREA)
Abstract
采用旋转TIN网和距离与夹角双权内插法的三维复杂建模方法,它涉三维建模方法,以解决现有三维复杂建模方法的数据结构中数据间的分散性和交错性不足、呼应配合程度不足,现有算法中没有对相邻数据间组成的相邻三角形平面间夹角的影响予以直接考虑的缺陷。包括下述步骤:在旋转网的节点处设采样点;如果旋转TIN网中某采样点在以该采样点为中心,在基础正方形边长的0.15倍范围内调整后,可作为建模体特征点,则将该采样点调整到新位置;根据建模体特征增设其它特征点;将增设的特征点与其所在的三角形的三个顶点相连,形成局部次级加密网;采用距离与夹角双权内插法进行特征值内插;制作等值线;进行建模体三维立体透视图的制作。
Description
技术领域
本发明涉及三维建模领域的一种三维复杂建模方法。
背景技术
三维建模领域中的三维复杂建模方法首先要有好的数据结构,在此基础上要有好的算法。现有三维复杂建模方法中的数据结构有规则形和不规则形两类,它们的共同不足之处是没有对三维方向上的控制距离分割进行更好地考虑,数据间的分散性和交错性不足、呼应配合程度不足。现有算法中有内插法和近似法两类,它们均没有对相邻数据间组成的相邻三角形平面间夹角的影响予以直接考虑。由于上述两方面存在的不足之处,使三维复杂建模方法的质量和精度仍然有待提高。
发明内容
本发明的目的是提供一种采用旋转TIN网和距离与夹角双权内插法的三维复杂建模方法,以解决现有三维复杂建模方法的数据结构中数据间的分散性和交错性不足、呼应配合程度不足,现有算法中没有对相邻数据间组成的相邻三角形平面间夹角的影响予以直接考虑的缺陷。它包括下述步骤:一、以分辨率的二倍作为基础正方形的边长;二、在基础正方形的四个顶点各布置一个控制点,在基础正方形相邻的两条边上分散布置n1个控制点,n1为自然数,在基础正方形内分散布置n2个控制点,n2为正整数,且上述控制点中任意两点之间的距离在基础正方形边长的0.20~0.85倍之间,基础正方形顶点、边上和内部的控制点以三角形连接,保证相邻的任意两个控制点之间的连线长度在基础正方形边长的0.20~0.85倍之间,基础正方形中所有三角形中所有的内角都应大于18度,并使连线所分割的相邻两三角形的6个内角中最小者最大为原则,最终形成基本单元1(图1);三、把基本单元1逆时针旋转90×n3度形成一次单元体2(图2),n3为整数,把基本单元1逆时针旋转90×n4度形成二次单元体3(图3),n4为整数,把基本单元1逆时针旋转90×n5度形成三次单元体4(图4),n5为整数,;四、用基本单元1、一次单元体2、二次单元体3和三次单元体4拼接成一个正方形的配套单元(图5),配套单元中相邻两个基本单元公共边上只有顶点处的两个控制点时,公共边两侧、分处于两个基本单元之内的三角形顶点处的控制点相连,保证相邻的任意两个控制点之间的连线长度在基础正方形边长的0.20~0.85倍之间,配套单元中所有三角形中所有的内角都应大于18度,并使连线所分割的相邻两三角形的6个内角中最小者最大、网络中的控制点具有较好的交错性和分散性为原则,可适当修改基础正方形中的三角形连线;五、以配套单元为复制单位,重复平移复制配套单元,形成覆盖建模区的旋转网,相邻两个配套单元中的相邻两个基本单元的公共边上只有顶点处的两个控制点时,公共边两侧的、分处于两个基本单元之内的距离最近的三角形顶点处的控制点相连,保证相邻的任意两个控制点之间的连线长度在基础正方形边长的0.20~0.85倍之间,配套单元中所有三角形中所有的内角都应大于18度,并使连线所分割的相邻两三角形的6个内角中最小者最大、网络中的控制点具有较好的交错性和分散性为原则,形成初始采样网络(图6),所述初始采样网络是去除基础正方形相邻两个顶点处控制点的连线的所有控制点与相邻控制点之间的连线,初始采样网络中相邻两个控制点间连线的长度在基础正方形边长的0.20~0.85倍范围内;六、依次分别选择配套单元中的1个单元作为采样网络的起点可得到初始采样网络的4个亚方案,在4个亚方案中选择一个;七、在初始采样网的控制点处布置控制点;八、如果初始采样网中某控制点在以该初始控制点为中心,在基础正方形边长的0.15倍范围内调整后,可以作为建模体特征点,则将该控制点调整到新位置,形成动态采样网,动态采样网中所有三角形中所有的内角都应大于18度。初始采样网中的初始控制点也可根据其它实际情况调整,调整范围也是正方形边长的0.15倍范围内,调整后形成的动态采样网中所有三角形中所有的内角都应大于18度;九、根据建模体特征布置其它特征点;十、将九中布置的建模体特征点与其所在三角形的三个顶点相连,形成局部次级加密网;十一、在八中动态采样网中的控制点和九中的建模体特征点处采集平面位置和特征值数据;十二、将控制点位置和特征值及采样网中三角形连接数据输入计算机,进行数字化储存和管理;十三、采用距离与夹角双权内插法进行特征值内插,
(一)内插点在三角形平面内时内插点特征值的求解方法:
见示意图32-a(32-b与其同理),A(xa,ya,za),B(xb,yb,zb),C(xc,yc,zc)是三角形的三个顶点,且za>zb>zc,ΔABJ,ΔBCI,ΔACH是与ΔABC临边的三个三角形,M、N、L、K分别是ΔABJ,ΔBCI,ΔACH和ΔABC的分布中心,D(xd,yd,zd)点在AC边上,zd=zb。BD为等值线。G(xg,yg,zg)点是ABC三角形内欲内插特征值的点,zgo是G点在ΔABC平面内内插的特征值,zg是G点根据ΔABC和ΔABJ,ΔBCI,ΔACH四个三角形的数据内插的特征值。与ABC三角形邻边的三个三角形的分布中心(可根据情况选择算术平均中心、加权平均中心和中位中心)分别为M(xm,ym,zm)、N(xn,yn,zn)和L(xl,yl,zl)。E点在AB边上,F点在AC边上,E,G,F点三点在同一条直线上,特征值相同,即zg=ze=zf。S、F、E、Q四点特征值相同,即zs=zf=ze=zq。P点在BC边上,N、P、U三点特征值相同,即zn=zp=zu。θ1(∠GEQ),θ2(∠GFS),θ3(∠UPN)分别为相邻两三角形同名等值线间的夹角(θ1,θ2,θ3均取小于180°者)。自G点到以G点所在三角形为内核向外的第一圈三角形分布中心的距离分别为d1(GM),d2(GL),d3(GN),……dn。
1、在三角形ABC中用线性内插的方法求出G点的特征值zGO。
2、用算术平均法求出以ΔABC为内核向外的第一圈三角形分布中心M、N、L,……n点的平面坐标。用线性插值的方法分别在M、N、L,……n点所在的平面上求出M、N、L,……n点的特征值,ZM,ZN,ZL,……Zn。
3、用在本三角形内线性内插的方法求出ΔABC和以ΔABC为内核向外的第一圈中三角形间同名等值线端点和转折点Q、E、F、S、U、P、N,……点的平面坐标,zQ=zE=zF=zS=zG,zU=zP=zN,……。
4、求ΔABC和以ΔABC为内核向外的第一圈中三角形间同名等值线间夹角
xEQ=xE-xQ,yEQ=yE-yQ,zEQ=zE-zQ,xEG=xE-xG,yEG=yE-yG,zEG=zE-zG
a1·b1=xEQ·xEG+yEQ·yEG+zEQ·zEG
xFS=xF-xS,yFS=yF-yS,zFS=zF-zS,xFG=xF-xG,yFG=yF-yG,zFG=zF-zG
a2·b2=xFS·xFG+yFS·yFG+zFS·zFG
xPN=xP-xN,yPN=yP-yN,zPN=zP-zN,xPU=xP-xU,yPU=yP-yU,zPU=zP-zU
a3·b3=xPN·xPU+yPN·yPU+zPN·zPU
同理求出ΔABC和以ΔABC为内核向外的第一圈中其它三角形间同名等值线间夹角,……直至
5、求以ΔABC为内核向外的第一圈中三角形分布中心距G点的距离
直至
6、求G点特征值
式中ZG是内插点G的特征值,λ与ε均是正整数,ZGi是用线性插值的方法分别求得的以ΔABC为内核向外的第一圈三角形分布中心M、N、L,……n点的特征值ZM,ZN,ZL,……Zn的统称,当θ=90°、ε=2、λ=1时,因此,原有反距离加权法是该方法在θ=90°、ε=2、λ=1时的特例。
(二)内插点在三角形公共边上时的特征值内插
见示意图32--c,AB是相邻ΔABJ和ΔABC的公共边,内插点G(xg,yg,zg)在AB边上,ΔAJE,ΔBJF,ΔBCI,ΔACH,分别是与ΔABJ和ΔABC的外边相邻的四个三角形,W、D、N、L分别是ΔAJE,ΔBJF,ΔBCI,ΔACH的分布中心。zgo是G点在AB边上线性内插的特征值,zg是G点根据ΔABC、ΔABJ和其周边ΔAJE、ΔBJF、ΔBCI、ΔACH四个三角形的数据内插的特征值。θ1(∠TQU),θ2(∠MKL),θ3(∠GPS),θ4(∠GVR)分别为相邻两三角形同名等值线间的夹角(θ1,θ2,θ3,θ4均取小于180°者)。自G点到以G点所在边为对角线的四边形作为内核向外的第一圈三角形分布中心的距离分别为d1(GD),d2(GL),d3(GN),d4(GW),···dn。
1、在AB边上用线性内插的方法求G点特征值ZGO。
2、用算术平均法求以AB边为对角线的四边形为内核向外的第一圈中三角形分布中心的平面坐标,W、D、N、L,……直至n。
3、求G点分别到以AB边为对角线的四边形为内核向外的第一圈中三角形分布中心W、D、N、L,……直至n的距离
直至
在以AB边为内核向外的第一圈中三角形ΔAJE,ΔBJF,ΔBCI,ΔACH,……中分别用线性内插的方法求分布中心D、N、L、W,……的特征值。
4、用JB和GD两直线方程联立的方法求Q点的平面坐标
5、在ΔBJF中用线性内插的方法求Q点的特征值
6、在所属三角形中用线性内插的方法求V、P、R、S、K、M、T、U,……点的平面坐标
7、分别求ΔABJ、ΔABC和以AB边为对角线的四边形为内核向外的第一圈三角形中与其相邻的三角形同名等高线间夹角
xQT=xQ-xT,yQT=yQ-yT,zQT=zQ-zT,xQU=xQ-xU,yQU=yQ-yU,zQU=zQ-zU
a1·b1=xQT·xQU+yQT·yQU+zQT·zQU
xKL=xK-xL,yKL=yK-yL,zKL=zK-zL,xKM=xK-xM,yKM=yK-yM,zKM=zK-zM
a2·b2=xKL·xKM+yKL·yKM+zKL·zKM
xPS=xP-xS,yPS=yP-ys,zPS=zp-zs,xPG=xP-xG,yPG=yP-yG,zPG=zP-zG
a3·b3=xPS·xPG+yPS·yPG+zPS·zPG
xVW=xV-xW,yVW=yV-yW,zVW=zV-zW,xVG=xV-xG,yVG=yV-yG,zVG=zV-zG
a4·b4=xVW·xVG+yVW·yVG+zVW·zVG
8、求G点特征值
式中ZG是内插点G的特征值,λ与ε均是正整数,ZGi是用线性插值的方法分别求得的以ΔABC为内核向外的第一圈三角形分布中心M、N、L,……n点的特征值ZM,ZN,ZL,……Zn的统称,当θ=90°、ε=2、λ=1时,因此,原有反距离加权法是该方法在θ=90°、ε=2、λ=1时的特例。
十四、按设置的等值线间距离制作等值线;十五、进行建模体三维立体透视图的制作。
发明的优点
一、采用旋转TIN网的优点
1、构网方便、快速、算法简单、便于程序控制;
2、较好地考虑了第三维数据的分布;
3、较好地考虑了小区域内三角形间拓扑关系的互相呼应,配合可以对建模体进行更好地控制的问题。
4、数据结构的难易程度介于规则网和TIN之间,是一种在正方形规则网之上有规律的TIN,比较好处理。数据存储规律性较强,正方形格网交点处的数据点按矩阵存储(用行号和列号隐性表示),正方形格网内的三个数据点根据配套网中四个基本单元的类型分为四种类型处理,即在标明所在格网中一个格网点的行号和列号的基础上,再标明基本单元中的类型(1、2、3、4个基本单元类型中的一种),就可以方便、迅速地计算处格网内数据点的坐标;
5、将旋转TIN与特征线配合使用效果更好;
6、沿各方向上的控制距离相近的程度大;
7、相对密集区和相对稀疏区呈围邻关系,因此,可以将相对密集区的趋势用于指导相对稀疏区,使趋势分析的整体准确度提高
8、有利于三维空间整个建模体序特征和序结构的保持。
二、采用距离与夹角双权内插法进行特征内插的优点
1、推出一个相邻三角形同名等值线间夹角余弦值的参数;
2、以线性内插的方法求得非线性内插的结果,且方法与算法简单、结果合理;
3、有利于三维空间整个建模体序特征和序结构的保持;
4、在相邻两个三角形的公共边处为平缓、光滑过渡。
三、采用旋转TIN网和距离与夹角双权内插法的复杂三维建模方法的优点
将旋转配套TIN网和距离与夹角双权内插法的优点结合起来可以更好地提高三维建模体的质量和精度。
附图说明
图1是实施方式一中基本单元1的示意图,图2是实施方式一中一次单元体2的示意图,图3是实施方式一中二次单元体3的示意图,图4是实施方式一中三次单元体4的示意图,图5是实施方式一中配套单元的示意图,图6是实施方式一中以配套单元为复制单位,进行平移复制的采样网示意图,图7是正方形网的示意图。图8是实施方式二中基本单元1的示意图,图9是实施方式二中一次单元体2的示意图,图10是实施方式二中二次单元体3的示意图,图11是实施方式二中三次单元体4的示意图,图12是实施方式二中配套单元的示意图,图13是实施方式二中以配套单元为复制单位,进行平移复制的采样网示意图。图14是实施方式三中基本单元1的示意图,图15是实施方式三中一次单元体2的示意图,图16是实施方式三中二次单元体3的示意图,图17是实施方式三中三次单元体4的示意图,图18是实施方式三中配套单元的示意图,图19是实施方式三中以配套单元为复制单位,进行平移复制的采样网示意图。图20是实施方式五中基本单元1的示意图,图21是实施方式五中一次单元体2的示意图,图22是实施方式五中二次单元体3的示意图,图23是实施方式五中三次单元体4的示意图,图24是实施方式五中配套单元的示意图,图25是实施方式五中以配套单元为复制单位,进行平移复制的采样网示意图。图26是实施方式四中基本单元1的示意图,图27是实施方式四中一次单元体2的示意图,图28是实施方式四中拼接体A1的示意图,图29是实施方式四中拼接体A2的示意图,图30是实施方式四中配套单元的示意图,图31是实施方式四中以配套单元为复制单位,进行平移复制的采样网示意图,图32-a是相邻三角形中过分布中心的等值线与公共边的交点在三角形内的示意图,图32-b是相邻三角形中过分布中心的等值线与公共边的交点在三角形外的示意图,图32-c是内插点在相邻两三角形公共边上的示意图。
具体实施方式
具体实施方式一:下面结合图8至图13具体说明本实施方式。本实施方式包括下述步骤:一、设分辨率是15mm,以15mm的二倍30mm作为基础正方形的边长;二、在基础正方形的四个顶点各布置一个控制点,在基础正方形内分散布置三个控制点,见示意图8,四个顶点的平面坐标分别为(0,0),(30,0),(30,30)和(0,30),a点为(9,9),b点为(12,22.5),c点为(24,12),上述a、b、c三点连接成三角形,上述的内部三个控制点连接成三角形,然后从基础正方形的四个顶点中选取一个顶点,该顶点距离与其相邻的内部三角形两个顶点的距离之和最小,把选取出来的这个正方形的顶点处的控制点和上述的与其相邻的内部三角形两个顶点处的控制点分别连接,过基础正方形的其余三个顶点分别向基础正方形内部三个控制点中就近的的控制点连线,最终形成基本单元1(图8);三、把基本单元1逆时针旋转90度形成一次单元体2(图9),把基本单元1逆时针旋转180度形成二次单元体3(图10),把基本单元1逆时针旋转270度形成三次单元体4(图11);四、把一次单元体2(图9)设置在基本单元1(图8)的下侧,把三次单元体4(图11)设置在基本单元1(图8)的右侧,把二次单元体3(图10)设置于一次单元体2(图9)的右侧和三次单元体4(图11)的下侧,基本单元1、一次单元体2、二次单元体3和三次单元体4拼接成一个正方形的配套单元(图12),配套单元中相邻两个基本单元公共边两侧、分处于两个基本单元之内的三角形顶点处的控制点相连;五、以配套单元为复制单位,重复平移复制配套单元,形成覆盖建模区的旋转网,相邻两个配套单元中的相邻两个基本单元的公共边两侧的、分处于两个基本单元之内的距离最近的三角形顶点处的控制点相连,形成初始采样网络(图6),所述初始采样网络是去除基础正方形相邻两个顶点处控制点的连线的所有控制点与相邻控制点之间的连线,,初始采样网络中相邻两个控制点间连线的长度在基础正方形边长的0.20~0.85倍范围内;六、依次分别选择配套单元中的1个单元作为采样网络的起点可得到初始采样网络的4个亚方案,在4个亚方案中选择一个;七、在初始采样网的控制点处布置控制点;八、如果初始采样网中某控制点在以该初始控制点为中心,在基础正方形边长的0.15倍范围内调整后,可以作为建模体特征点,则将该控制点调整到新位置,形成动态采样网,动态采样网中所有三角形中所有的内角都应大于18度。初始采样网中的初始控制点也可根据其它实际情况调整,调整范围也是正方形边长的0.15倍范围内,调整后形成的动态采样网中所有三角形中所有的内角都应大于18度;九、根据建模体特征布置其它特征点;十、将九中布置的建模体特征点与其所在三角形的三个顶点相连,形成局部次级加密网;十一、在八中动态采样网中的控制点和九中的建模体特征点处采集平面位置和特征值数据;十二、将控制点位置和特征值及采样网中三角形连接数据数据输入计算机,进行数字化储存和管理;十三、采用距离与夹角双权内插法进行特征值内插,
(一)内插点在三角形平面内时内插点特征值的求解方法:
见示意图32-a(32-by与其同理),A(xa,ya,za),B(xb,yb,zb),C(xc,yc,zc)是三角形的三个顶点,且za>zb>zc,ΔABJ,ΔBCI,ΔACH是与ΔABC临边的三个三角形,M、N、L、K分别是ΔABJ,ΔBCI,ΔACH和ΔABC的分布中心,D(xd,yd,zd)点在AC边上,zd=zb。BD为等值线。G(xg,yg,zg)点是ABC三角形内欲内插特征值的点,zgo是G点在ΔABC平面内内插的特征值,zg是G点根据ΔABC和ΔABJ,ΔBCI,ΔACH四个三角形的数据内插的特征值。与ABC三角形邻边的三个三角形的分布中心(可根据情况选择算术平均中心、加权平均中心和中位中心)分别为M(xm,ym,zm)、N(xn,yn,zn)和L(xl,yl,zl)。E点在AB边上,F点在AC边上,E,G,F点三点在同一条直线上,特征值相同,即zg=ze=zf。S、F、E、Q四点特征值相同,即zs=zf=ze=zq。P点在BC边上,N、P、U三点特征值相同,即zn=zp=zu。θ1(∠GEQ),θ2(∠GFS),θ3(∠UPN)分别为相邻两三角形同名等值线间的夹角(θ1,θ2,θ3均取小于180°者)。自G点到以G点所在三角形为内核向外的第一圈三角形分布中心的距离分别为d1(GM),d2(GL),d3(GN),……dn。
1、在三角形ABC中用线性内插的方法求出G点的特征值zGO。
2、用算术平均法求出以ΔABC为内核向外的第一圈三角形分布中心M、N、L,……n点的平面坐标。用线性插值的方法分别在M、N、L,……n点所在的平面上求出M、N、L,……n点的特征值,ZM,ZN,ZL,……Zn。
3、用在本三角形内线性内插的方法求出ΔABC和以ΔABC为内核向外的第一圈中三角形间同名等值线端点和转折点Q、E、F、S、U、P、N,……点的平面坐标,zQ=zE=zF=zS=zG,zU=zP=zN,……。
4、求ΔABC和以ΔABC为内核向外的第一圈中三角形间同名等值线间夹角
xEQ=xE-xQ,yEQ=yE-yQ,zEQ=zE-zQ,xEG=xE-xG,yEG=yE-yG,zEG=zE-zG
a1·b1=xEQ·xEG+yEQ·yFG+zEQ·zEG
xFS=xF-xS,yFS=yF-yS,zFS=zF-zS,xFG=xF-xG,yFG=yF-yG,zFG=zF-zG
a2·b2=xFS·xFG+yES·yFG+zFS·zFG
xPN=xP-xN,yPN=yP-yN,zPN=zP-zN,xPU=xP-xU,yPU=yP-yU,zPU=zP-zU
a3·b3=xPN·xPU+yPN·yPU+zPN·zPU
同理求出ΔABC和以ΔABC为内核向外的第一圈中其它三角形间同名等值线间夹角,……直至
5、求以ΔABC为内核向外的第一圈中三角形分布中心距G点的距离
直至
6、求G点高程
(二)内插点在三角形公共边上时的特征值内插
见示意图32--c,AB是相邻ΔABJ和ΔABC的公共边,内插点G(xg,yg,zg)在AB边上,ΔAJE,ΔBJF,ΔBCI,ΔACH,分别是与ΔABJ和ΔABC的外边相邻的四个三角形,W、D、N、L分别是ΔAJE,ΔBJF,ΔBCI,ΔACH的分布中心。zgo是G点在AB边上线性内插的特征值,zg是G点根据ΔABC、ΔABJ和其周边ΔAJE、ΔBJF、ΔBCI、ΔACH四个三角形的数据内插的特征值。θ1(∠TQU),θ2(∠MKL),θ3(∠GPS),θ4(∠GVR)分别为相邻两三角形同名等值线间的夹角(θ1,θ2,θ3,θ4均取小于180°者)。自G点到以G点所在边为对角线的四边形作为内核向外的第一圈三角形分布中心的距离分别为d1(GD),d2(GL),d3(GN),d4(GW),···dn。
1、在AB边上用线性内插的方法求G点特征值ZGO。
2、用算术平均法求以AB边为对角线的四边形为内核向外的第一圈中三角形分布中心的平面坐标,W、D、N、L,……直至n。
3、求G点分别到以AB边为对角线的四边形为内核向外的第一圈中三角形分布中心的,W、D、N、L,……直至n的距离
直至
在以AB边为内核向外的第一圈中三角形ΔAJE,ΔBJF,ΔBCI,ΔACH,……中分别用线性内插的方法求分布中心D、N、L、W,……的特征值。
4、用JB和GD两直线方程联立的方法求Q点的平面坐标
5、在ΔBJF中用线性内插的方法求Q点的特征值
6、在所属三角形中用线性内插的方法求V、P、R、S、K、M、T、U,……点的平面坐标
7、分别求ΔABJ、ΔABC和以AB边为对角线的四边形为内核向外的第一圈三角形中与其相邻的三角形同名等高线间夹角
xQT=xQ-xT,yQT=yQ-yT,zQT=zQ-zT,xQU=xQ-xU,yQU=yQ-yU,zQU=zQ-zU
a1·b1=xQT·xQU+yQT·yQU+zQT·zQU
xKL=xK-xL,yKL=yK-yL,zKL=zK-zL,xKM=xK-xM,yKM=yK-yM,zKM=zK-zM
a2·b2=xKL·xKM+yKL·yKM+zKL·zKM
xPS=xP-xS,yPS=yP-ys,zPS=zp-zs,xPG=xP-xG,yPG=yP-yG,zPG=zP-zG
a3·b3=xPS·xPG+yPS·yPG+zPS·zPG
xVW=xV-xW,yVW=yV-yW,zVW=zV-zW,xVG=xV-xG,yVG=yV-yG,zVG=zV-zG
a4·b4=xVW·xVG+yVW·yVG+zVW·zVG
8、求G点特征值
取ε=2、λ=1,
十四、按设置的等值线间距离制作等值线;十五、进行建模体三维立体透视图的制作。
具体实施方式二:见示意图1至示意图6,本实施方式与实施方式一的不同之处是:步骤二中,在基础正方形内分散布置的三个控制点分别为,a(6,6),b(15,24),c(24,15);其它步骤与实施方式一都相同。
具体实施方式三:见图14至图19,本实施方式与实施方式一的不同之处是:步骤二中,在基础正方形的四个顶点各布置一个控制点,在基础正方形内分散布置三个控制点,上述的内部三点连接成三角形,然后在基础正方形右边的中点布置一个控制点,与基础正方形的四个顶点上的控制点和基础正方形内分散布置的三个控制点相加一共八个控制点。示意图14中八个控制点分别为(0,0),(30,0),(30,30),(0,30),a(8.31,13.62),b(15,24),c(20.37,6),d(30,15)。将正方形右边中点处的控制点与内部就近的三角形的两个控制点分别连接,将与正方形右边中点处的控制点同在一条直线上的两个正方形顶点与内部就近的三角形的一个控制点分别连接,将正方形的其余两个顶点与内部就近的三角形的两个控制点分别连接,最终形成基本单元1(图14);步骤三中,把基本单元1逆时针旋转180度形成一次单元体2,把基本单元1逆时针旋转360度形成二次单元体3,把基本单元1逆时针旋转540度形成三次单元体4;四、把一次单元体2设置在基本单元1的下侧,把三次单元体4设置在基本单元1的右侧,把二次单元体3设置于一次单元体2的右侧和三次单元体4的下侧,基本单元1、一次单元体2、二次单元体3和三次单元体4拼接成一个正方形的配套单元(图18)。
具体实施方式四:见图26至图31,本实施方式与实施方式一的不同之处是:步骤二中,在基础正方形的四个顶点各布置一个控制点,在基础正方形内分散布置三个控制点,上述的内部三点连接成三角形,然后在基础正方形右边的中点布置一个控制点,与基础正方形的四个顶点上的控制点和基础正方形内分散布置的三个控制点相加一共八个控制点。图26中八个控制点分别为(0,0),(30,0),(30,30),(0,30),a(6,13.5),b(13.5,24),c(16.5,6),d(30,15)。将正方形右边中点处的控制点与内部就近的三角形的两个控制点分别连接,将与正方形右边中点处的控制点同在一条直线上的两个正方形顶点与内部就近的三角形的一个控制点分别连接,将正方形的其余两个顶点与内部就近的三角形的两个控制点分别连接,最终形成基本单元1;步骤三中,以基本单元1(图26)的右边为轴镜像对称形成一次单元体2(图27),把一次单元体2(图27)接于基本单元1(图26)的右侧形成拼接体A1(图28),把拼接体A1逆时针旋转180度形成拼接体A2(图29);步骤四中、把拼接体A2接于拼接体1的右侧,形成配套单元(图30)。
具体实施方式五:见图20至图25,本实施方式与实施方式一的不同之处是:步骤二中,在基础正方形的四个顶点各布置一个控制点,在基础正方形内分散布置三个控制点,上述的内部三点连接成三角形,然后在基础正方形底边的中点和右边的中点各布置一个控制点,与基础正方形的四个顶点上的控制点和基础正方形内分散布置的三个控制点相加一共九个控制点。示意图20中的九个点分别为(0,0),(30,0),(30,30),(0,30),a(5.13,10.89),b(14.16,24.96),c(20.16,15.9),d(30,15),e(15,0)。将基础正方形底边和右边中点处的控制点向与其内部就近的三角形两个顶点处的控制点分别连线,在上述的四个连线中去掉最长的连线,连接与上述最长的连线交叉的基础正方形顶点与内部三角形顶点,然后将基础正方形四个顶点处的控制点分别与内部就近的三角形一个顶点处的控制点连线,最终形成基本单元1,如图20所示;步骤四中,配套单元中相邻两个基本单元公共边上的控制点重合对接;步骤五中,相邻两个配套单元中的相邻两个基本单元的公共边两侧的、分处于两个基本单元之内的距离最近的三角形顶点处的控制点以最小距离相连,形成初始采样网(图25)。
具体实施方式六:本实施方式与实施方式一、实施方式二、实施方式三、实施方式四或实施方式五的不同在于:分别在实施方式一、实施方式二、实施方式三、实施方式四或实施方式五的步骤四和步骤五之间还包括步骤A,把配套单元的长度或宽度方向缩放,使之成为长方形,内部所有控制点的位置按照缩放比例调整。其它步骤仍与原实施方式相同。
Claims (2)
1.一种采用旋转TIN网和距离与夹角双权内插法的三维复杂建模方法,其特征在于它包括下述步骤:一、确定基础正方形的边长;二、在基础正方形的四个顶点各布置一个控制点,在基础正方形相邻的两条边上分散布置n1个控制点,n1为自然数,在基础正方形内分散布置n2个控制点,n2为正整数,且上述控制点中任意两点之间的距离在基础正方形边长的0.20~0.85倍之间,基础正方形顶点、边上和内部的控制点以三角形连接,保证相邻的任意两个控制点之间的连线长度在基础正方形边长的0.20~0.85倍之间,基础正方形中所有三角形中所有的内角都应大于18度,在由同样的4个相邻点连接成两个相邻三角形的两种方案中,选择相邻两三角形的6个内角中最小者与另外一方案中相邻两三角形的6个内角中最小者相比较大的一个方案,最终形成基本单元1;三、把基本单元1逆时针旋转90×n3度形成一次单元体2,n3为整数,把基本单元1逆时针旋转90×n4度形成二次单元体3,n4为整数,把基本单元1逆时针旋转90×n5度形成三次单元体4,n5为整数;四、用基本单元1、一次单元体2、二次单元体3和三次单元体4拼接成一个正方形的配套单元,配套单元中相邻两个基本单元公共边上只有顶点处的两个控制点时,公共边两侧、分处于两个基本单元之内的三角形顶点处的控制点相连,保证相邻的任意两个控制点之间的连线长度在基础正方形边长的0.20~0.85倍之间,配套单元中所有三角形中所有的内角都应大于18度,在由同样的4个相邻点连接成两个相邻三角形的两种方案中,选择相邻两三角形的6个内角中最小者与另外一方案中相邻两三角形的6个内角中最小者相比较大的一个方案、网络中的控制点具有较好的交错性和分散性为原则,可适当修改基础正方形中的三角形连线;五、以配套单元为复制单位,重复平移复制配套单元,形成覆盖建模区的旋转网,相邻两个配套单元中的相邻两个基本单元的公共边上只有顶点处的两个控制点时,公共边两侧的、分处于两个基本单元之内的距离最近的三角形顶点处的控制点相连,保证相邻的任意两个控制点之间的连线长度在基础正方形边长的0.20~0.85倍之间,配套单元中所有三角形中所有的内角都应大于18度,在由同样的4个相邻点连接成两个相邻三角形的两种方案中,选择相邻两三角形的6个内角中最小者与另外一方案中相邻两三角形的6个内角中最小者相比较大的一个方案、网络中的控制点具有较好的交错性和分散性为原则,形成初始采样网络,所述初始采样网络是去除基础正方形相邻两个顶点处控制点的连线的所有控制点与相邻控制点之间的连线,初始采样网络中相邻两个控制点间连线的长度在基础正方形边长的0.20~0.85倍范围内;六、依次分别选择配套单元中的1个单元作为采样网络的起点可得到初始采样网络的4个亚方案,在4个亚方案中选择一个;七、在初始采样网的控制点处布置控制点;八、如果初始采样网中某控制点在以该初始控制点为中心,在基础正方形边长的0.15倍范围内调整后,可以作为建模体特征点,则将该控制点调整到新位置,形成动态采样网,动态采样网中所有三角形中所有的内角都应大于18度,初始采样网中的初始控制点也可根据其它实际情况调整,调整范围也是正方形边长的0.15倍范围内,调整后形成的动态采样网中所有三角形中所有的内角都应大于18度;九、根据建模体特征布置其它特征点;十、将九中布置的建模体特征点与其所在三角形的三个顶点相连,形成局部次级加密网;十一、在八中动态采样网中的控制点和九中的建模体特征点处采集平面位置和特征值数据;十二、将控制点位置和特征值及采样网中三角形连接数据输入计算机,进行数字化储存和管理;十三、采用距离与夹角双权内插法进行特征值内插,
式中θ1,θ2,θ3,θ4…θn分别为相邻两三角形同名等高线间的夹角,均取小于180°者,ZG是内插点G的特征值,λ与ε均是正整数,ZGi是用线性插值的方法分别求得的以ΔABC为内核向外的第一圈三角形分布中心M、N、L,……n点的特征值ZM,ZN,ZL,……Zn的统称,
内插点G在三角形平面内时,d1,d2,d3,d4,……dn分别是自G点到以G点所在三角形为内核向外的第一圈三角形分布中心的距离,ZGO是内插点G在本三角形中线性内插的高程值,
内插点G在三角形公共边上时,d1,d2,d3,d4,……dn分别是自G点到以G点所在边为对角线的四边形作为内核向外的第一圈三角形分布中心的距离,ZGO是内插点G在公共边上线性内插的高程值;
十四、按设置的等值线间距离制作等值线;十五、进行建模体三维立体透视图的制作。
2.依据权利要求1所述的采用旋转TIN网和距离与夹角双权内插法的三维复杂建模方法,其特征在于在步骤四和步骤五之间还包括步骤A,把配套单元的长度或宽度方向缩放,使之成为长方形,内部所有控制点位置按照缩放比例调整。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110110179.2A CN102760305B (zh) | 2011-04-29 | 2011-04-29 | 采用旋转tin网和距离与夹角双权内插法的三维复杂建模方法 |
PCT/CN2012/000581 WO2012146068A1 (zh) | 2011-04-29 | 2012-05-02 | 一种矿山开采三维仿真系统的设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110110179.2A CN102760305B (zh) | 2011-04-29 | 2011-04-29 | 采用旋转tin网和距离与夹角双权内插法的三维复杂建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102760305A CN102760305A (zh) | 2012-10-31 |
CN102760305B true CN102760305B (zh) | 2015-04-01 |
Family
ID=47054758
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110110179.2A Expired - Fee Related CN102760305B (zh) | 2011-04-29 | 2011-04-29 | 采用旋转tin网和距离与夹角双权内插法的三维复杂建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102760305B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104077810B (zh) * | 2014-06-26 | 2017-09-26 | 中南大学 | 基于tin的被错断断层复原建模方法 |
CN112132967B (zh) * | 2020-10-23 | 2023-04-18 | 中建一局集团第二建筑有限公司 | 一种山区洪水模拟方法、存储介质及设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101114276A (zh) * | 2007-08-28 | 2008-01-30 | 中国地质大学(武汉) | 固体矿床三维可视化储量计算系统及计算方法 |
CN101251930A (zh) * | 2007-08-30 | 2008-08-27 | 中国地质大学(武汉) | 三维地质模型定制式无损任意分割方法 |
CN101345001A (zh) * | 2008-08-28 | 2009-01-14 | 杨小平 | 地下采矿三维仿真模型的制作方法 |
-
2011
- 2011-04-29 CN CN201110110179.2A patent/CN102760305B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101114276A (zh) * | 2007-08-28 | 2008-01-30 | 中国地质大学(武汉) | 固体矿床三维可视化储量计算系统及计算方法 |
CN101251930A (zh) * | 2007-08-30 | 2008-08-27 | 中国地质大学(武汉) | 三维地质模型定制式无损任意分割方法 |
CN101345001A (zh) * | 2008-08-28 | 2009-01-14 | 杨小平 | 地下采矿三维仿真模型的制作方法 |
Non-Patent Citations (2)
Title |
---|
《基于TIN的三维地层建模及可视化技术研究》;熊祖强等;《岩土力学》;20070930;第28卷(第9期);1954-1958 * |
《基于三维TIN的精细表面建模方法》;李逢春等;《计算机应用研究》;20060831(第8期);159-161 * |
Also Published As
Publication number | Publication date |
---|---|
CN102760305A (zh) | 2012-10-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Barr | Superquadrics and angle-preserving transformations | |
CN110516388B (zh) | 基于调和映射的曲面离散点云模型环切刀轨生成方法 | |
CN101254732B (zh) | 生成具有互锁几何形状的无定形图形的片材的方法 | |
Mirante et al. | The radial sweep algorithm for constructing triangulated irregular networks | |
CN109360260A (zh) | 一种三角网格三维模型的切挖重构算法 | |
CN109377561A (zh) | 一种基于共形几何的数模表面网格生成方法 | |
CN107665242A (zh) | 一种区域空间多尺度网格编码方法及装置 | |
US20100162637A1 (en) | Supporting Structure for Freeform Surfaces in Buildings | |
CN100383823C (zh) | 三维扫描的点云孔洞填补方法 | |
CN101976468B (zh) | 一种多分辨率动态地形可视化方法及系统 | |
CN102853803B (zh) | 一种文物破损面积的测定方法 | |
CN102853763A (zh) | 一种不规则表面文物体积的测量方法 | |
CN106055794A (zh) | 一种参数化建立三维管线模型的方法 | |
CN101783023B (zh) | AutoCAD图块与基于路径的点状地图符号交换方法 | |
CN108733892A (zh) | 一种关于混凝土中非凸形颗粒的三维数值构造方法 | |
CN108053477A (zh) | 一种管道内变形的数值处理方法 | |
CN106777822A (zh) | 二维高含量颗粒增强复合材料三相细观模型及建立方法 | |
CN107146287A (zh) | 二维投影图像至三维模型的映射方法 | |
CN102760305B (zh) | 采用旋转tin网和距离与夹角双权内插法的三维复杂建模方法 | |
CN100386779C (zh) | 基于通用图形显示卡的被测体正投影与反投影方法 | |
CN102760377B (zh) | 采用旋转tin网和距离与夹角双权内插法的大比例尺地形图测绘方法 | |
Wang et al. | Hexagonal meshes with planar faces | |
Stupariu et al. | Geometric approaches to computing 3D-landscape metrics | |
CN102760297B (zh) | 采用旋转tin网的多边形影像数字化方法 | |
Sadoc et al. | Modelling of the structure of glasses |
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: 20150401 Termination date: 20180429 |
|
CF01 | Termination of patent right due to non-payment of annual fee |