專(zhuān)利名稱(chēng):基于動(dòng)態(tài)插幀技術(shù)的曲面過(guò)程時(shí)空模擬方法
技術(shù)領(lǐng)域:
本發(fā)明涉及一種曲面時(shí)空模擬方法,特別是一種基于動(dòng)態(tài)插幀技術(shù)的曲面過(guò)程時(shí)空模擬方法,典型應(yīng)用有地面沉降動(dòng)態(tài)模擬、地下水變化過(guò)程等動(dòng)態(tài)模擬。
背景技術(shù):
我們知道,世界上90%的信息是通過(guò)視覺(jué)系統(tǒng)感知的,為了使信息表達(dá)更加形象直觀,易于理解,圖形成為傳遞信息的最主要媒體之一,計(jì)算機(jī)圖形學(xué)因此得到了很大的發(fā)展。尤其隨著現(xiàn)代科學(xué)技術(shù)的發(fā)展,三維圖形可視化技術(shù)逐步發(fā)展起來(lái),并成功應(yīng)用在各行各業(yè)。三維圖形圖像可以形象直觀地表達(dá)客觀世界,三維動(dòng)態(tài)模擬技術(shù)可以實(shí)現(xiàn)對(duì)難以描述的一些動(dòng)態(tài)過(guò)程進(jìn)行直觀形象的模擬。例如,在對(duì)某一個(gè)區(qū)域歷年的地面沉降情況進(jìn)行模擬時(shí),以往一般的做法是將歷年的沉降值以二維圖表的形式表達(dá)出來(lái),這種表達(dá)形式只能給人一種數(shù)值上的變化感知, 沒(méi)有地面的形態(tài)表達(dá),更沒(méi)有沉降動(dòng)態(tài)的模擬過(guò)程。目前利用三維技術(shù),用曲面模型表達(dá)地面的高程情況,用多個(gè)曲面模型模擬歷年的地面沉降的多個(gè)狀態(tài),加上時(shí)間維度,就可以模擬出地面沉降的過(guò)程,這種動(dòng)態(tài)的表達(dá)方式不僅可以直觀地感到數(shù)值上的變化趨勢(shì),還較逼真地模擬了地面沉降動(dòng)態(tài)變化過(guò)程。然而,目前在三維動(dòng)態(tài)過(guò)程模擬技術(shù)中,常常由于模型變化的無(wú)規(guī)則性,以及模型樣本觀測(cè)值的稀缺等原因,要想做出較準(zhǔn)確的、實(shí)時(shí)的時(shí)空動(dòng)態(tài)過(guò)程模擬比較困難。我們知道,在3D動(dòng)畫(huà)中,一般實(shí)現(xiàn)的方法是通過(guò)關(guān)鍵幀,一幀一個(gè)畫(huà)面來(lái)實(shí)現(xiàn)一個(gè)過(guò)程的模擬。因?yàn)槿说娜庋鄣姆磻?yīng)時(shí)間是1Λ4秒,所以當(dāng)畫(huà)面數(shù)一秒鐘達(dá)到M幀以上時(shí),人眼是分不出來(lái)具體的每個(gè)畫(huà)面的,看到的是連續(xù)的流暢的變化過(guò)程。所以在三維動(dòng)態(tài)過(guò)程模擬中,也需要將關(guān)鍵幀設(shè)置在M幀以上,才能達(dá)到流暢的顯示效果。但是原始數(shù)據(jù)量一般沒(méi)有那么多, 另外模型較多對(duì)內(nèi)存要求也更高。
發(fā)明內(nèi)容
本發(fā)明目的在于解決上述現(xiàn)有技術(shù)的不足,尤其是模型樣本觀測(cè)值的稀缺的問(wèn)題而提出一種基于動(dòng)態(tài)插幀技術(shù)的曲面過(guò)程時(shí)空模擬方法,可以實(shí)現(xiàn)流暢平滑的時(shí)空模擬過(guò)程。本方法突破常規(guī)對(duì)觀測(cè)信息以數(shù)據(jù)或二維圖表的形式表達(dá)的方式,而采用曲面時(shí)空變化的畫(huà)面表達(dá)觀測(cè)信息。本方法是由原始的數(shù)據(jù)構(gòu)建模擬的初始曲面,再通過(guò)插值計(jì)算,構(gòu)建出每個(gè)插幀時(shí)刻的曲面空間模型并顯示,從而實(shí)現(xiàn)曲面時(shí)空動(dòng)態(tài)過(guò)程的模擬。為了實(shí)現(xiàn)本發(fā)明目的所采用的技術(shù)方案是提出一種基于動(dòng)態(tài)插幀技術(shù)的曲面過(guò)程時(shí)空模擬方法,通過(guò)計(jì)算機(jī)實(shí)施,包括以下步驟
步驟一,開(kāi)始模擬,構(gòu)建模擬初始曲面根據(jù)初始的已有數(shù)據(jù)利用普通三維曲面構(gòu)建的方法先構(gòu)建模擬的初始曲面;
步驟二,初始化曲面模擬的參數(shù)初始化參數(shù)包括整個(gè)模擬演示的時(shí)間和每秒鐘插入的幀數(shù)步驟三,動(dòng)態(tài)插幀構(gòu)面,包括以下幾個(gè)步驟(1)、計(jì)算時(shí)間值;(2)、插值計(jì)算所有點(diǎn)的 Z值;(3)、將所有點(diǎn)的Z值更新;(4)、重新構(gòu)面;
所述計(jì)算時(shí)間值是計(jì)算當(dāng)前插幀時(shí)刻距離模擬初始時(shí)刻的時(shí)間值;
所述的插值計(jì)算所有點(diǎn)的Z值是根據(jù)原始觀測(cè)信息計(jì)算所有點(diǎn)的Z值;
所述的將所有點(diǎn)的Z值更新是將步驟(2)計(jì)算的所有點(diǎn)的Z值更新到構(gòu)建曲面的程序
中;
所述的重新構(gòu)面是在所有點(diǎn)的Z值更新后,由更新后的Z值重新構(gòu)建曲面; 步驟四,將插幀構(gòu)建的曲面進(jìn)行三維顯示,曲面三維顯示之后繼續(xù)轉(zhuǎn)向步驟三,即繼續(xù)重復(fù)進(jìn)行插幀構(gòu)建曲面,直至完成最后一個(gè)曲面模型的顯示,即模擬過(guò)程結(jié)束。本發(fā)明在曲面模擬過(guò)程中,在模擬之前獲取的原始數(shù)據(jù)可以是地面離散觀測(cè)點(diǎn)以及每一年觀測(cè)得到的地面高程數(shù)據(jù),也可以是每一年地面的曲面模型。如果獲取的數(shù)據(jù)是原始的采樣信息,則由離散的觀測(cè)點(diǎn)和觀測(cè)最初期的地面高程信息構(gòu)建一個(gè)初始的曲面模型,可以利用普通三維曲面構(gòu)建的方法構(gòu)建即可。如果獲取的數(shù)據(jù)已經(jīng)是曲面模型,則可以直接利用這個(gè)原始的曲面為初始曲面模型。本發(fā)明在曲面模擬過(guò)程中,是在刷幀的過(guò)程中完成曲面的構(gòu)建與更新,這樣做可以避免占用計(jì)算機(jī)內(nèi)存,有效地解決了對(duì)計(jì)算機(jī)內(nèi)存需求大的要求。本發(fā)明所述的動(dòng)態(tài)插幀構(gòu)面的過(guò)程中,插值計(jì)算所有點(diǎn)Z值的方法步驟如下
①、原始采樣的離散點(diǎn)不變,根據(jù)時(shí)間間隔ΔT內(nèi)空間上的變化值Δ S,計(jì)算這一時(shí)間間隔內(nèi)的平均變化速度V= Δ S/ Δ T ;
②、計(jì)算出所有時(shí)間間隔內(nèi)的平均變化速度Vi=Δ Si/ Δ Ti,將所有的各段時(shí)間間隔內(nèi)的平均速度作為該段時(shí)間間隔內(nèi)的初始速度,即Ttl至T1時(shí)間間隔內(nèi)的空間變化值為S1-I, 平均速度Vch1= (S1 - S0)/ (T1 - Ttl),則將Ttl時(shí)刻的變化速度Vtl設(shè)為平均速度Vch1, V0 = Vch1= (S1- S0)/ (T1- Tq),其他時(shí)刻的速度同理處理;
③、由得到每個(gè)時(shí)刻的變化速度算出每段時(shí)間間隔內(nèi)的空間變化的加速度^=AVi / Δ Ti,得到每個(gè)插幀時(shí)刻的變化速度,設(shè)Tk為一插幀時(shí)刻,如果Tk處于Ttl至T1時(shí)間段, 則Tk時(shí)刻的速度為Vk = V0+ ai* (Tk- T0),計(jì)算出K點(diǎn)空間上的變化值,Z值Zk= Z0 + 0. 5* (VVk) (Tk- T0);如果 Tk 處于 T1 至 T2 則時(shí)間段,則 Vk= V1+ a2* (Tk- T1), Z 值 Zk= Z1+ 0.5* (V^Vk) (Tk- T1);處于其他時(shí)間段依次類(lèi)推,即得到新的Z值。式中“*”、“ /”分別表示乘號(hào)、除號(hào)。本發(fā)明所述的動(dòng)態(tài)插幀構(gòu)面的過(guò)程中,采用非均勻插值算法,即時(shí)間變化周期Δ T 內(nèi)的空間變化的速度是動(dòng)態(tài)變化的,可以實(shí)現(xiàn)插幀的平滑過(guò)渡。本發(fā)明的基于動(dòng)態(tài)插幀技術(shù)的曲面過(guò)程時(shí)空模擬方法,與傳統(tǒng)的模擬表達(dá)方法相比具有如下顯著的優(yōu)點(diǎn)
(1)采用本發(fā)明動(dòng)態(tài)插幀技術(shù)的曲面過(guò)程時(shí)空模擬方法,突破以往常規(guī)模型分析方式, 提供一種直觀生動(dòng)的表現(xiàn)形式。( 2 )本發(fā)明實(shí)現(xiàn)了時(shí)間與空間上的同步模擬,模擬效果準(zhǔn)確、實(shí)時(shí),表達(dá)力極強(qiáng),可以清晰地看出某一時(shí)刻相應(yīng)的情況及過(guò)程變化趨勢(shì)。(3)在原始數(shù)據(jù)相對(duì)有限的情況下,本發(fā)明采用動(dòng)態(tài)插幀技術(shù)構(gòu)建出符合需要的插幀曲面,插值算法簡(jiǎn)單,較好的解決了關(guān)鍵幀之間的平滑過(guò)渡問(wèn)題。
(4)本發(fā)明提出了動(dòng)態(tài)插幀構(gòu)面,曲面過(guò)程模擬時(shí),在刷幀的過(guò)程中完成曲面模型的構(gòu)建與更新,解決了對(duì)計(jì)算機(jī)內(nèi)存要求高的難題。(5)本發(fā)明實(shí)現(xiàn)的動(dòng)態(tài)過(guò)程模擬,畫(huà)面流暢、自然。
圖1為本發(fā)明時(shí)空動(dòng)態(tài)模擬過(guò)程示意圖。圖2為本發(fā)明曲面時(shí)空動(dòng)態(tài)模擬流程圖。圖3為本發(fā)明插值計(jì)算Z值的過(guò)程示意圖。圖4為兩種不同插值算法的比較示意圖。
具體實(shí)施例方式下面結(jié)合附圖和具體實(shí)施例對(duì)本發(fā)明的方法作進(jìn)一步的詳述。實(shí)施例1 本發(fā)明提出一種基于動(dòng)態(tài)插幀技術(shù)的曲面過(guò)程時(shí)空模擬方法,所述時(shí)空曲面動(dòng)態(tài)過(guò)程模擬是指在時(shí)間上和空間上能保持同步,空間上曲面的變化過(guò)程能以真實(shí)的三維立體直觀的方式體現(xiàn)。將時(shí)間空間模擬的過(guò)程用圖表示,就是隨著時(shí)間軸,在每一時(shí)間點(diǎn)上有對(duì)應(yīng)的空間模型信息,時(shí)間上連貫起來(lái)就是一個(gè)變化過(guò)程的動(dòng)態(tài)模擬,參見(jiàn)圖1?,F(xiàn)在以模擬某一地區(qū)2000-2009年十年的地面沉降變化情況為例,假設(shè)獲取的原始數(shù)據(jù)是地面離散觀測(cè)點(diǎn)以及每一年觀測(cè)得到的地面高程數(shù)據(jù),采用本發(fā)明的基于動(dòng)態(tài)插幀技術(shù)的曲面過(guò)程時(shí)空模擬方法進(jìn)行模擬,通過(guò)計(jì)算機(jī)實(shí)施模擬地面沉降變化,包括以下步驟
步驟一,開(kāi)始模擬,構(gòu)建模擬初始的曲面。由已知的離散的觀測(cè)點(diǎn)和觀測(cè)最初期的信息,即2000年的地面高程信息構(gòu)建一個(gè)初始的曲面模型,利用普通三維曲面構(gòu)建的方法構(gòu)建即可。這里構(gòu)建模擬初始曲面的目的是,為整個(gè)曲面模擬過(guò)程提供一個(gè)曲面模型構(gòu)建的程序,在下步的模擬過(guò)程中,只需要計(jì)算所有點(diǎn)的Z值來(lái)更新曲面。因?yàn)樗星婺P?,都有空間坐標(biāo),整個(gè)地面沉降變化過(guò)程只是地面高程值的變化,即空間上Z值的變化。步驟二,初始化曲面模擬參數(shù)。初始化參數(shù)包括整個(gè)模擬演示過(guò)程的時(shí)間和每秒鐘插幀數(shù)。假設(shè)初始化2000年至2009年的地面沉降變化整個(gè)模擬演示過(guò)程的時(shí)間是10秒鐘,插幀數(shù)是30,則每一年的變化過(guò)程對(duì)應(yīng)一秒鐘的模擬演示時(shí)間,每秒鐘內(nèi)需要插幀30 次。步驟三,插幀構(gòu)建曲面。插幀構(gòu)建曲面是整個(gè)曲面過(guò)程時(shí)空模擬的核心過(guò)程,具體分以下幾個(gè)步驟(1)、計(jì)算時(shí)間值;(2)、插值計(jì)算所有點(diǎn)的Z值;(3)、將所有點(diǎn)的Z值更新;(4)、重新構(gòu)面;
所述的計(jì)算模擬的時(shí)間值是指計(jì)算當(dāng)前插幀時(shí)刻的時(shí)間距離模擬開(kāi)始時(shí)刻的時(shí)間值, 即當(dāng)前插幀時(shí)刻模擬進(jìn)行了多久,由步驟二得到總的模擬時(shí)間是10秒鐘,則這個(gè)時(shí)間值在 0至10秒之間,包括0秒和10秒;
所述的插值計(jì)算所有點(diǎn)的Z值,這一步是整個(gè)模擬的關(guān)鍵步驟,是體現(xiàn)本發(fā)明的特點(diǎn), 插值算法影響模擬的效果;Z值的計(jì)算在后結(jié)合圖3和圖4專(zhuān)門(mén)詳述;
所述的將所有點(diǎn)的Z值更新,是將上一步插值計(jì)算得到的所有點(diǎn)的Z值更新到構(gòu)建曲面的程序中;所述的重新構(gòu)面,由曲面構(gòu)建程序重新構(gòu)建曲面,這里曲面構(gòu)建的程序即步驟一中的普通三維曲面構(gòu)建的程序。步驟四,將構(gòu)建的曲面進(jìn)行三維顯示,構(gòu)建一個(gè)曲面就刷幀顯示,并接著轉(zhuǎn)向步驟三,繼續(xù)計(jì)算時(shí)間值,進(jìn)行插值計(jì)算并重新構(gòu)建曲面,直至完成最后一個(gè)曲面模型的顯示, 即模擬過(guò)程結(jié)束。最后一個(gè)曲面模型就是2009年的曲面模型。下面結(jié)合圖3和圖4說(shuō)明本發(fā)明的插值計(jì)算Z值方法,其步驟為
①、根據(jù)獲取的歷年地面高程信息,即每年觀測(cè)的地面的沉降量,可以直接計(jì)算出每年沉降的平均速度;如圖3所示,用V T1、T2, T3> T4, T5時(shí)刻分別代表2000年、2001年、2002 年、2003年、2004年、2005年,則對(duì)應(yīng)變化的平均速度為V。、V1, V2, V3> V4, V5 ;
②、對(duì)平均速度做這樣的一個(gè)處理,即將第①步計(jì)算的所有的各段時(shí)間間隔內(nèi)的平均速度作為該段時(shí)間間隔內(nèi)的初始速度,例如,Ttl至T1時(shí)間間隔內(nèi)的平均速度Vtl作為T(mén)tl時(shí)刻的變化速度,T1至T2時(shí)間間隔內(nèi)的平均速度V1作為T(mén)1時(shí)刻的變化速度,其他時(shí)刻的速度同理類(lèi)推,如圖3所示;
③、每個(gè)時(shí)刻的變化速度得到以后,可以算出每段時(shí)間間隔內(nèi)的空間變化的加
速度
屮=Δ Vi/Δ Ti,如 Ei1= (V1-V0) / (T1-T0), a2= (V2-V1)/ (T2-T1),從而可以計(jì)算得到任一插幀時(shí)刻的變化速度,如設(shè)Tk為一插幀時(shí)刻,圖3中Tk時(shí)刻的速度為Vk=V1+ B2^(Tk-T1)0 再經(jīng)計(jì)算得出空間上的變化值,即得到新的Z值。假設(shè)Vtl、\、\、W5時(shí)刻的Z值分別為‘ Z1, Z2, Z3, Z4, Z5,這些時(shí)刻的Z值是已知的原始數(shù)據(jù),即每年觀測(cè)得到的地面高程值。 則Vk時(shí)刻的Z值Zk= Z1+ 0.5* (VVk) (Tk- T1),其它插幀時(shí)刻的Z值計(jì)算同理。本發(fā)明所述的動(dòng)態(tài)插幀,一方面是指插幀構(gòu)建曲面是動(dòng)態(tài)連續(xù)的過(guò)程,另一方面是指插幀構(gòu)建的曲面是動(dòng)態(tài)變化的,并不是均勻的變化,即插值算法采用非均勻插值算法。插值算法有兩種,平均插值和非均勻插值。平均插值即將空間變化速度看成均勻變化的;非均勻插值就是每段時(shí)間內(nèi)空間變化的加速度是均勻的,空間變化速度是動(dòng)態(tài)變化的。圖4將兩種插值算法反應(yīng)到空間變化量上,可以清晰地對(duì)比跳躍與平滑過(guò)渡情況。如果采用平均插值,如圖4中(a),即每一年內(nèi)的沉降變化是均勻的,但在每?jī)蓚€(gè)時(shí)間段的過(guò)度時(shí)刻,沉降變化量上可能相差太大,存在跳躍問(wèn)題。如果采用非均勻插值,如圖4中(b), 即每一年內(nèi)的沉降變化的加速度是均勻的,則每一年內(nèi)的沉降變化速度就是動(dòng)態(tài)變化的, 這樣就可以實(shí)現(xiàn)插幀的平滑過(guò)渡,避免跳躍情況。
權(quán)利要求
1.一種基于動(dòng)態(tài)插幀技術(shù)的曲面過(guò)程時(shí)空模擬方法,通過(guò)計(jì)算機(jī)實(shí)施,其特征在于 包括以下步驟步驟一,開(kāi)始模擬,構(gòu)建模擬初始曲面根據(jù)初始的已有數(shù)據(jù)利用普通三維曲面構(gòu)建的方法先構(gòu)建模擬的初始曲面;步驟二,初始化曲面模擬的參數(shù)初始化參數(shù)包括整個(gè)模擬演示的時(shí)間和每秒鐘插入的幀數(shù)步驟三,動(dòng)態(tài)插幀構(gòu)面,包括以下幾個(gè)步驟(1)、計(jì)算時(shí)間值;(2)、插值計(jì)算所有點(diǎn)的 Z值;(3)、將所有點(diǎn)的Z值更新;(4)、重新構(gòu)面;所述計(jì)算時(shí)間值是計(jì)算當(dāng)前插幀時(shí)刻距離模擬初始時(shí)刻的時(shí)間值;所述的插值計(jì)算所有點(diǎn)的Z值是根據(jù)原始觀測(cè)信息計(jì)算所有點(diǎn)的Z值;所述的將所有點(diǎn)的Z值更新是將步驟(2)計(jì)算的所有點(diǎn)的Z值更新到構(gòu)建曲面的程序中;所述的重新構(gòu)面是在所有點(diǎn)的Z值更新后,由更新后的Z值重新構(gòu)建曲面;步驟四,將插幀構(gòu)建的曲面進(jìn)行三維顯示,曲面三維顯示之后繼續(xù)轉(zhuǎn)向步驟三,即繼續(xù)重復(fù)進(jìn)行插幀構(gòu)建曲面,直至完成最后一個(gè)曲面模型的顯示,即模擬過(guò)程結(jié)束。
2.根據(jù)權(quán)利要求1所述的基于動(dòng)態(tài)插幀技術(shù)的曲面過(guò)程時(shí)空模擬方法,其特征在于 所述的動(dòng)態(tài)插幀構(gòu)面的過(guò)程中,插值計(jì)算所有點(diǎn)Z值的方法步驟如下①、原始采樣的離散點(diǎn)不變,根據(jù)時(shí)間間隔ΔT內(nèi)空間上的變化值Δ S,計(jì)算這一時(shí)間間隔內(nèi)的平均變化速度V= Δ S/ Δ T ;②、計(jì)算出所有時(shí)間間隔內(nèi)的平均變化速度Vi= Δ Si / Δ Ti,將所有的各段時(shí)間間隔內(nèi)的平均速度作為該段時(shí)間間隔內(nèi)的初始速度,即Ttl至T1時(shí)間間隔內(nèi)的空間變化值為S1 -S。,平均速度^廠伐-S0)/ (T1- T0),則將Ttl時(shí)刻的變化速度Vtl設(shè)為平均速度Vch1, V0 =V= (S1- S0)/ (T1- Tq),其他時(shí)刻的速度同理處理;③、由得到每個(gè)時(shí)刻的變化速度算出每段時(shí)間間隔內(nèi)的空間變化的加速度^=AVi / ATi,得到每個(gè)插幀時(shí)刻的變化速度,設(shè)Tk為一插幀時(shí)刻,如果Tk處于Ttl至T1時(shí)間段,則 Tk時(shí)刻的速度為Vk = V0 + B1 * (Tk- T0),計(jì)算出K點(diǎn)空間上的變化值,Z值Z0 + 0.5* (V Vk) (Tk- Ttl);如果 Tk 處于 T1 至 T2 則時(shí)間段,則 Vk= V1+ a2* (Tk - T1), Z 值 Zk= Z1+ 0.5 * (V Vk) (Tk- T1);處于其他時(shí)間段依次類(lèi)推,即得到新的Z值。
3.根據(jù)權(quán)利要求2所述的基于動(dòng)態(tài)插幀技術(shù)的曲面過(guò)程時(shí)空模擬方法,其特征在于 所述的動(dòng)態(tài)插幀構(gòu)面的過(guò)程中,采用非均勻插值算法,即時(shí)間變化周期Δ T內(nèi)的空間變化的速度是動(dòng)態(tài)變化的,可以實(shí)現(xiàn)插幀的平滑過(guò)渡。
全文摘要
本發(fā)明涉及一種基于動(dòng)態(tài)插幀技術(shù)的曲面過(guò)程時(shí)空模擬方法,其步驟為開(kāi)始模擬;構(gòu)建模擬初始曲面;初始化曲面模擬的參數(shù);動(dòng)態(tài)插幀構(gòu)面的過(guò)程;最后將插幀構(gòu)建的曲面進(jìn)行三維顯示。所述的動(dòng)態(tài)插幀構(gòu)面的過(guò)程的步驟中包括計(jì)算時(shí)間值,插值計(jì)算所有點(diǎn)的Z值,將所有點(diǎn)的Z值更新,重新構(gòu)面。在時(shí)間上通過(guò)不斷的計(jì)算時(shí)間值,構(gòu)建出每個(gè)插幀時(shí)刻的曲面空間模型并顯示,直至完成最后一個(gè)曲面模型的顯示。本方法的動(dòng)態(tài)插幀技術(shù)解決了當(dāng)樣本觀測(cè)值相對(duì)有限時(shí)模擬過(guò)程跳躍不流暢的問(wèn)題。本發(fā)明突破了以往對(duì)觀測(cè)信息的表達(dá)方式,采用曲面時(shí)空變化的畫(huà)面表達(dá)觀測(cè)信息,逼真直觀地模擬其變化過(guò)程,特別適用于地下水模擬和地面沉降等曲面變化過(guò)程的模擬。
文檔編號(hào)G06T15/00GK102360514SQ201110319819
公開(kāi)日2012年2月22日 申請(qǐng)日期2011年10月20日 優(yōu)先權(quán)日2011年10月20日
發(fā)明者劉修國(guó), 劉培, 張柱, 易云蕾, 花衛(wèi)華, 鄭坤 申請(qǐng)人:中國(guó)地質(zhì)大學(xué)(武漢)