一種海陸交錯帶氣象等值線自動繪制方法
【專利摘要】本發(fā)明公開了一種海陸交錯帶氣象等值線自動繪制方法。包括如下步驟：1）從多維海洋氣象數(shù)據(jù)中讀取數(shù)據(jù)，并獲取海陸分界數(shù)據(jù)，進行疊加分析分離出海洋區(qū)域數(shù)據(jù)；2）設(shè)定判斷噪點的閾值，并對每個像元進行判斷，若像元為噪點，對該像元進行校正，得到海洋區(qū)域有效數(shù)值范圍；3）根據(jù)海洋區(qū)域的有效數(shù)值范圍得到全部區(qū)域的一個劃分；4）分別構(gòu)建海洋區(qū)域和其他區(qū)域的間隔表，設(shè)置各個區(qū)域不同疏密的間隔；5）合并海洋區(qū)域和其他區(qū)域的間隔；6）利用MATLAB批量讀取數(shù)據(jù)，并將間隔應(yīng)用于MATLAB等值線繪制指令，實現(xiàn)自動繪制。本發(fā)明繪制的等值線可有效突出海洋區(qū)域要素的變化，并兼顧陸地區(qū)域氣象要素的表現(xiàn)，可應(yīng)用于海陸交錯帶氣象專題圖的制作等領(lǐng)域。
【專利說明】一種海陸交錯帶氣象等值線自動繪制方法

【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及海洋氣象等值線繪制的【技術(shù)領(lǐng)域】。尤其涉及一種海陸交錯帶氣象等值 線自動繪制方法。

【背景技術(shù)】
[0002] 等值線是制圖對象中數(shù)值相等的各點連成的平滑曲線，它在水文、地質(zhì)、海洋、大 氣等諸多領(lǐng)域都有非常重要的用途，比如氣壓場中的等壓線圖、大地測繪中的等高線圖、地 質(zhì)中的礦體等值線圖等等。目前等值線繪制和填充的方法已經(jīng)比較成熟，但在海洋氣象領(lǐng) 域存在問題。海洋氣象數(shù)據(jù)包含陸地和海洋區(qū)域，要素數(shù)值在海洋上變化平緩，而在陸地上 變化顯著。采用傳統(tǒng)的等值線間隔設(shè)置方法會使等值線集中在陸地區(qū)域，在海洋區(qū)域則缺 少可分辨的等值線，而其正是海洋氣象工作者重點關(guān)注的區(qū)域。因此，如何自動地設(shè)置海洋 氣象等值線間隔，從而提高繪制要素的信息量，突出興趣區(qū)域，是一個亟待解決的問題。
[0003] 對于不等間距等值線，目前有兩種間隔設(shè)置的方法：一是繪制人員根據(jù)具體數(shù)據(jù) 手動設(shè)置；二是基于面積統(tǒng)計的間隔設(shè)置。手動設(shè)置法是利用軟件繪制完等值線后，根據(jù)效 果不斷調(diào)整間隔，直到出來比較滿意的結(jié)果。這種方法往往可以得到很好的效果，但是繪制 過程麻煩，耗時耗力。面積統(tǒng)計的方法在一定程度上解決了人工繪制的問題，可以實現(xiàn)業(yè)務(wù) 化批量處理，但是應(yīng)用于海洋氣象這類特殊數(shù)據(jù)時，繪制出的等值線效果無法滿足專業(yè)要 求。比如大氣等壓線的間隔通常需要設(shè)置為特定數(shù)值，如2. 5百帕或者4百帕，而面積統(tǒng)計 的方法只根據(jù)等面積原則劃分，無法解決這類特殊需求。
[0004] 針對以上問題，本發(fā)明提出一種海陸交錯帶氣象等值線自動繪制方法，首先將原 始要素與海陸分界數(shù)據(jù)做疊加分析，提取出海洋區(qū)域，然后進行斑點去噪計算得到海洋區(qū) 域數(shù)值范圍，再結(jié)合映射表選擇合適的間隔值，合并海洋和陸地區(qū)域獲得不等間距間隔，最 后將間隔應(yīng)用于MATLAB等值線繪制指令實現(xiàn)自動繪制?；谠摲椒梢詫崿F(xiàn)海陸交錯帶 氣象等值線的自動繪制，并且能夠突出海洋區(qū)域信息。


【發(fā)明內(nèi)容】

[0005] 本發(fā)明的目的是為克服現(xiàn)有技術(shù)存在的問題，提供一種海陸交錯帶氣象等值線自 動繪制方法。
[0006] -種海陸交錯帶氣象等值線自動繪制方法包括如下步驟：
[0007] 1)從多層多變量的海洋氣象原始數(shù)據(jù)文件中讀取出待繪制等值線的一層數(shù)據(jù)，再 獲取與該數(shù)據(jù)格式一致的海陸分界數(shù)據(jù)，利用海陸分界數(shù)據(jù)對原始數(shù)據(jù)做疊加分析，分離 出海洋區(qū)域數(shù)據(jù)；
[0008] 2)設(shè)定判斷噪點的閾值，遍歷數(shù)據(jù)中的每個像元，若當(dāng)前像元的數(shù)值與周圍相連 的八個像元的方差減去整體的方差大于給定的閾值，則認(rèn)為該像元為噪點，對該像元進行 校正，以此得到海洋區(qū)域有效數(shù)值范圍；
[0009] 3)根據(jù)海洋區(qū)域的有效數(shù)值范圍得到全部區(qū)域的一個劃分；
[0010] 4)分別構(gòu)建海洋區(qū)域和其他區(qū)域的間隔表，根據(jù)間隔表設(shè)置各個區(qū)域不同疏密的 間隔；
[0011] 5)合并海洋區(qū)域和其他區(qū)域的間隔；
[0012] 6)利用MATLAB批量讀取數(shù)據(jù)，并將步驟5)合并后的間隔應(yīng)用于MATLAB等值線繪 制指令，實現(xiàn)自動繪制。
[0013] 所述的步驟1)包括：從一個多層多變量的海洋氣象原始數(shù)據(jù)文件中讀取出待繪 制等值線的一層數(shù)據(jù)，該層數(shù)據(jù)是由η行m列的二維規(guī)則格網(wǎng)構(gòu)成，再獲取與該數(shù)據(jù)格式一 致的η行m列的海陸分界數(shù)據(jù)，海陸分界數(shù)據(jù)包含與其相對應(yīng)的二值數(shù)據(jù)，每個像元都引用 相同的地理位置，海陸分界二值數(shù)據(jù)中0代表其他，1代表海洋，利用海陸分界數(shù)據(jù)對原始 數(shù)據(jù)做疊加分析，分離出海洋區(qū)域數(shù)據(jù)。
[0014] 所述的步驟2)包括：設(shè)定一個判斷噪點的閾值，遍歷數(shù)據(jù)中的每個像元，若當(dāng)前 像元的數(shù)值f(i，j)與周圍相連的八個像元的方差σ 2減去整體的方差&大于給定的閾值， 則認(rèn)為該像元為噪點，該方法可用以下公式表示：
[0015] ^Ii8=i[gi-f(U)]2-^>￡ (1)
[0016] 式中：gi為f (i，j)點周圍鄰點的數(shù)值；
[0017] 當(dāng)確認(rèn)某像元為噪點時，用以下公式計算的結(jié)果替換原值，
[0018] = (1 -P)f(i,j) + Ρ[;Σ?=ι&] (2)
[0019] 式中：f' (i，j)為像元的新亮度；Ρ為該像元是噪點的概率，
[0020] 經(jīng)以上校正后的海洋區(qū)域數(shù)據(jù)可以得到海洋區(qū)域有效數(shù)值范圍。
[0021] 所述的步驟3)包括：得到整體數(shù)值范圍，再去除海洋區(qū)域的計算結(jié)果，以此得到 其他區(qū)域的數(shù)值范圍；海洋區(qū)域數(shù)值范圍Q〇= [min0，max0]可以將其他區(qū)域數(shù)值范圍分為 兩個部分，上界QU = (minLpmaxLj和下界QL2 = [minL2，maxL2)，它們是全部區(qū)域的范圍 QA = [minA，maxA]的一個劃分；
[0022] QA = QL! U Q0 U QL2 (3)
[0023] 式中QLp Q0與QL2兩兩交集都為空；
[0024] 若其他數(shù)值全部大于海洋區(qū)域，則QQ為空集，若其他數(shù)值全部小于海洋區(qū)域，則 QL2為空集。
[0025] 所述的步驟4)包括：對于求出的海洋區(qū)域和其他區(qū)域的數(shù)值范圍，構(gòu)建間隔表設(shè) 置各個區(qū)域不同疏密的間隔，間隔表區(qū)域字段分為興趣區(qū)域和其他區(qū)域，不同區(qū)域的間隔 設(shè)置原則不同，興趣區(qū)域間隔較小，對應(yīng)的等值線也較密，其他區(qū)域間隔相對較大，繪制的 等值線也較稀疏；同一類區(qū)域因為數(shù)值范圍的不同設(shè)置不同間隔，所述的數(shù)值范圍關(guān)注的 是范圍差所屬的區(qū)間，海洋區(qū)域的數(shù)值范圍差為maxO-minO。
[0026] 所述的步驟5)包括：通過步驟4)構(gòu)建的間隔表的映射，得到以下間隔設(shè)置，海洋 區(qū)域數(shù)值范圍間隔01〇=[111；[110，;[1^61^310，1]^1〇]，其它區(qū)域數(shù)值范圍上界( >)11^1=(111；[111^， interval!^，maxLj 和下界 QIL2 = [minL2, intervalL2，maxL2)，合并海洋區(qū)域和其他區(qū)域 的間隔，得到整體區(qū)域的間隔QIA，
[0027] QIA = QILi U QIO U QIL2 (4)
[0028] 式中QIL^QIO和QIL2兩兩交集都為空。
[0029] 所述的步驟6)包括：利用MATLAB批量讀取數(shù)據(jù)，調(diào)用MATLAB等值線繪制指令 contour或contourf，將合并后的間隔QIA應(yīng)用于該指令，實現(xiàn)自動批量繪制。
[0030] 本發(fā)明與現(xiàn)有技術(shù)相比具有的有益效果是：本發(fā)明通過設(shè)置不等間距間隔，確保 海洋區(qū)域有更密集的等值線，能重點展現(xiàn)出海洋區(qū)域的數(shù)值變化，使得在海洋區(qū)域中原本 較為稀疏的等值線分布得更加科學(xué)、豐富；而在次重要的陸地區(qū)域，原本較為密集的等值線 相對趨于稀疏，以降低該區(qū)域等值線的混亂程度；海洋區(qū)域和陸地區(qū)域等值線所體現(xiàn)的整 體信息更加完整，能提供更有意義的參考價值。

【專利附圖】

【附圖說明】：
[0031] 圖1是海陸交錯帶氣象等值線自動繪制方法；
[0032] 圖2是陸地表面氣壓數(shù)值分布；
[0033] 圖3是海洋表面氣壓數(shù)值分布；
[0034] 圖4是等值線繪制成果。

【具體實施方式】：
[0035] 如圖1所示，海陸交錯帶氣象等值線自動繪制方法包括如下步驟：
[0036] 1)從多層多變量的海洋氣象原始數(shù)據(jù)文件中讀取出待繪制等值線的一層數(shù)據(jù)，再 獲取與該數(shù)據(jù)格式一致的海陸分界數(shù)據(jù)，利用海陸分界數(shù)據(jù)對原始數(shù)據(jù)做疊加分析，分離 出海洋區(qū)域數(shù)據(jù)；
[0037] 2)設(shè)定判斷噪點的閾值，遍歷數(shù)據(jù)中的每個像元，若當(dāng)前像元的數(shù)值與周圍相連 的八個像元的方差減去整體的方差大于給定的閾值，則認(rèn)為該像元為噪點，對該像元進行 校正，以此得到海洋區(qū)域有效數(shù)值范圍；
[0038] 3)根據(jù)海洋區(qū)域的有效數(shù)值范圍得到全部區(qū)域的一個劃分；
[0039] 4)分別構(gòu)建海洋區(qū)域和其他區(qū)域的間隔表，根據(jù)間隔表設(shè)置各個區(qū)域不同疏密的 間隔；
[0040] 5)合并海洋區(qū)域和其他區(qū)域的間隔；
[0041] 6)利用MATLAB批量讀取數(shù)據(jù)，并將步驟5)合并后的間隔應(yīng)用于MATLAB等值線繪 制指令，實現(xiàn)自動繪制。
[0042] 所述的步驟1)包括：從一個多層多變量的海洋氣象原始數(shù)據(jù)文件中讀取出待繪 制等值線的一層數(shù)據(jù)，該層數(shù)據(jù)是由η行m列的二維規(guī)則格網(wǎng)構(gòu)成，再獲取與該數(shù)據(jù)格式一 致的η行m列的海陸分界數(shù)據(jù)，海陸分界數(shù)據(jù)包含與其相對應(yīng)的二值數(shù)據(jù)，每個像元都引用 相同的地理位置，海陸分界二值數(shù)據(jù)中0代表其他，1代表海洋，利用海陸分界數(shù)據(jù)對原始 數(shù)據(jù)做疊加分析，分離出海洋區(qū)域數(shù)據(jù)。
[0043] 所述的步驟2)包括：設(shè)定一個判斷噪點的閾值，遍歷數(shù)據(jù)中的每個像元，若當(dāng)前 像元的數(shù)值f(i，j)與周圍相連的八個像元的方差σ 2減去整體的方差^大于給定的閾值， 則認(rèn)為該像元為噪點，該方法可用以下公式表示：
[0044] 金ZiUa - f(i, j)]2 - 7 > ε ( 1 )
[0045] 式中：gi為f (i，j)點周圍鄰點的數(shù)值；
[0046] 當(dāng)確認(rèn)某像元為噪點時，用以下公式計算的結(jié)果替換原值，
[0047] f'oj) = (1 - p)f〇,j) + p glilia] ( 2 )
[0048] 式中：f' (i, j)為像元的新亮度；P為該像元是噪點的概率，
[0049] 經(jīng)以上校正后的海洋區(qū)域數(shù)據(jù)可以得到海洋區(qū)域有效數(shù)值范圍。
[0050] 所述的步驟3)包括：得到整體數(shù)值范圍，再去除海洋區(qū)域的計算結(jié)果，以此得到 其他區(qū)域的數(shù)值范圍；海洋區(qū)域數(shù)值范圍Q〇= [min0，max0]可以將其他區(qū)域數(shù)值范圍分為 兩個部分，上界QU = (minLpmaxLj和下界QL2 = [minL2，maxL2)，它們是全部區(qū)域的范圍 QA = [minA，maxA]的一個劃分；
[0051] QA = QL! U Q0 U QL2 (3)
[0052] 式中QLp Q0與QL2兩兩交集都為空；
[0053] 若其他數(shù)值全部大于海洋區(qū)域，則QQ為空集，若其他數(shù)值全部小于海洋區(qū)域，則 QL2為空集。
[0054] 所述的步驟4)包括：對于求出的海洋區(qū)域和其他區(qū)域的數(shù)值范圍，構(gòu)建間隔表設(shè) 置各個區(qū)域不同疏密的間隔，間隔表區(qū)域字段分為興趣區(qū)域和其他區(qū)域，不同區(qū)域的間隔 設(shè)置原則不同，興趣區(qū)域間隔較小，對應(yīng)的等值線也較密，其他區(qū)域間隔相對較大，繪制的 等值線也較稀疏；同一類區(qū)域因為數(shù)值范圍的不同設(shè)置不同間隔，所述的數(shù)值范圍關(guān)注的 是范圍差所屬的區(qū)間，海洋區(qū)域的數(shù)值范圍差為maxO-minO。
[0055] 所述的步驟5)包括：通過步驟4)構(gòu)建的間隔表的映射，得到以下間隔設(shè)置，海洋 區(qū)域數(shù)值范圍間隔01〇=[111；[110，;[1^61^310，1]^1〇]，其它區(qū)域數(shù)值范圍上界( >)11^1=(111；[111^， interval!^，maxLj 和下界 QIL2 = [minL2, intervalL2，maxL2)，合并海洋區(qū)域和其他區(qū)域 的間隔，得到整體區(qū)域的間隔QIA，
[0056] QIA = QILi U QIO U QIL2 (4)
[0057] 式中QIL^QIO和QIL2兩兩交集都為空。
[0058] 所述的步驟6)包括：利用MATLAB批量讀取數(shù)據(jù)，調(diào)用MATLAB等值線繪制指令 contour或contourf，將合并后的間隔QIA應(yīng)用于該指令，實現(xiàn)自動批量繪制。
[0059] 實施例：
[0060] (1)利用海陸分界數(shù)據(jù)對等值線原始數(shù)據(jù)做疊加分析，分離出海洋部分的數(shù)據(jù) [0061] 以2013年5月23日12點0時次的印度洋區(qū)域表明氣壓數(shù)值預(yù)報產(chǎn)品為例。經(jīng)度 范圍30。E?125° E，緯度范圍20。S?30。N，分辨率為476*251。通過海陸分界數(shù)據(jù) 的疊加分析，分離出海洋和陸地部分?jǐn)?shù)據(jù)。圖2為陸地區(qū)域表面氣壓數(shù)值預(yù)報產(chǎn)品數(shù)據(jù)的 統(tǒng)計圖，對該數(shù)據(jù)的分析可知，陸地區(qū)域的氣壓值范圍為600-1100百帕，其分布成鋸齒形， 從850-1020區(qū)間連續(xù)有幾個峰值，向兩邊逐步遞減。圖3為海洋區(qū)域表面氣壓數(shù)值預(yù)報產(chǎn) 品數(shù)據(jù)的統(tǒng)計圖，海洋區(qū)域的氣壓值集中在1000-1020區(qū)間，其余區(qū)間只有少量的噪點。
[0062] (2)采用八鄰點法去除海洋區(qū)域數(shù)據(jù)的噪點，獲取有效的數(shù)值范圍
[0063] 設(shè)定判斷噪點的閾值為1，遍歷數(shù)據(jù)中的每個像元，若當(dāng)前像元的數(shù)值f(i，j)與 周圍相連的八個像元的方差σ 2減去整體的方差&大于給定的閾值，則認(rèn)為該像元為噪 點。經(jīng)校正后的海洋區(qū)域數(shù)據(jù)可以得到海洋區(qū)域有效數(shù)值范圍。
[0064] (3)根據(jù)海洋區(qū)域的有效數(shù)值范圍得到全部區(qū)域的一個劃分
[0065] 計算得到整體數(shù)值范圍為600-1100百帕，海洋區(qū)域范圍為1000-1020百帕，以此 得到全部區(qū)域的一個劃分為[600,1000] U (1000,1020) U [1020,1100]。
[0066] (4)利用專家經(jīng)驗分別構(gòu)建海洋區(qū)域和其他區(qū)域的間隔表，設(shè)置間隔
[0067] 根據(jù)印度洋區(qū)域表面氣壓分布的特定，設(shè)計了映射間隔表，如表2所示。
[0068] 表 1
[0069]

【權(quán)利要求】
1. 一種海陸交錯帶氣象等值線自動繪制方法，其特征在于包括如下步驟： 1) 從多層多變量的海洋氣象原始數(shù)據(jù)文件中讀取出待繪制等值線的一層數(shù)據(jù)，再獲取 與該數(shù)據(jù)格式一致的海陸分界數(shù)據(jù)，利用海陸分界數(shù)據(jù)對原始數(shù)據(jù)做疊加分析，分離出海 洋區(qū)域數(shù)據(jù)； 2) 設(shè)定判斷噪點的閾值，遍歷數(shù)據(jù)中的每個像元，若當(dāng)前像元的數(shù)值與周圍相連的八 個像元的方差減去整體的方差大于給定的閾值，則認(rèn)為該像元為噪點，對該像元進行校正， 以此得到海洋區(qū)域有效數(shù)值范圍； 3) 根據(jù)海洋區(qū)域的有效數(shù)值范圍得到全部區(qū)域的一個劃分； 4) 分別構(gòu)建海洋區(qū)域和其他區(qū)域的間隔表，根據(jù)間隔表設(shè)置各個區(qū)域不同疏密的間 隔； 5) 合并海洋區(qū)域和其他區(qū)域的間隔； 6) 利用MATLAB批量讀取數(shù)據(jù)，并將步驟5)合并后的間隔應(yīng)用于MATLAB等值線繪制指 令，實現(xiàn)自動繪制。
2. 根據(jù)權(quán)利要求1所述的一種海陸交錯帶氣象等值線自動繪制方法，其特征在于所述 的步驟1)包括：從一個多層多變量的海洋氣象原始數(shù)據(jù)文件中讀取出待繪制等值線的一 層數(shù)據(jù)，該層數(shù)據(jù)是由η行m列的二維規(guī)則格網(wǎng)構(gòu)成，再獲取與該數(shù)據(jù)格式一致的η行m列 的海陸分界數(shù)據(jù)，海陸分界數(shù)據(jù)包含與其相對應(yīng)的二值數(shù)據(jù)，每個像元都引用相同的地理 位置，海陸分界二值數(shù)據(jù)中0代表其他，1代表海洋，利用海陸分界數(shù)據(jù)對原始數(shù)據(jù)做疊加 分析，分離出海洋區(qū)域數(shù)據(jù)。
3. 根據(jù)權(quán)利要求1所述的一種海陸交錯帶氣象等值線自動繪制方法，其特征在于所 述的步驟2)包括：設(shè)定一個判斷噪點的閾值，遍歷數(shù)據(jù)中的每個像元，若當(dāng)前像元的數(shù)值 f(i，j)與周圍相連的八個像元的方差σ2減去整體的方差p大于給定的閾值，則認(rèn)為該像 元為噪點，該方法可用以下公式表示： 魯Z]8=i[gi-f(i,D]2-tT2>￡ (1) 式中：gl為f (i，j)點周圍鄰點的數(shù)值； 當(dāng)確認(rèn)某像元為噪點時，用以下公式計算的結(jié)果替換原值， f'(i, j) = (1 - j) + Ρ [5Σ1Ι1 gi] ( 2 ) 式中：f' (i，j)為像元的新亮度；P為該像元是噪點的概率， 經(jīng)以上校正后的海洋區(qū)域數(shù)據(jù)可以得到海洋區(qū)域有效數(shù)值范圍。
4. 根據(jù)權(quán)利要求1所述的一種海陸交錯帶氣象等值線自動繪制方法，其特征在于所述 的步驟3)包括：得到整體數(shù)值范圍，再去除海洋區(qū)域的計算結(jié)果，以此得到其他區(qū)域的數(shù) 值范圍；海洋區(qū)域數(shù)值范圍Q〇 = [minO, maxO]可以將其他區(qū)域數(shù)值范圍分為兩個部分，上 界QQ = (minL^maxLj和下界QL2 = [minL2，maxL2)，它們是全部區(qū)域的范圍QA = [minA, maxA]的一個劃分； QA = QQ U QO U QL2 (3) 式中QU、Q0與QL2兩兩交集都為空； 若其他數(shù)值全部大于海洋區(qū)域，則QU為空集，若其他數(shù)值全部小于海洋區(qū)域，則QL2為 空集。
5. 根據(jù)權(quán)利要求1所述的一種海陸交錯帶氣象等值線自動繪制方法，其特征在于所述 的步驟4)包括：對于求出的海洋區(qū)域和其他區(qū)域的數(shù)值范圍，構(gòu)建間隔表設(shè)置各個區(qū)域不 同疏密的間隔，間隔表區(qū)域字段分為興趣區(qū)域和其他區(qū)域，不同區(qū)域的間隔設(shè)置原則不同， 興趣區(qū)域間隔較小，對應(yīng)的等值線也較密，其他區(qū)域間隔相對較大，繪制的等值線也較稀 疏；同一類區(qū)域因為數(shù)值范圍的不同設(shè)置不同間隔，所述的數(shù)值范圍關(guān)注的是范圍差所屬 的區(qū)間，海洋區(qū)域的數(shù)值范圍差為maxO-minO。
6. 根據(jù)權(quán)利要求1所述的一種海陸交錯帶氣象等值線自動繪制方法，其特征在于所述 的步驟5)包括：通過步驟4)構(gòu)建的間隔表的映射，得到以下間隔設(shè)置，海洋區(qū)域數(shù)值范圍 間隔〇1〇=[111；[11〇，；[1^61^310，1]^1〇]，其它區(qū)域數(shù)值范圍上界( >)11^1=(111；[111^，；[1^61^311^， maxLj和下界QIL2 = [minL2，intervalL2, maxL2)，合并海洋區(qū)域和其他區(qū)域的間隔，得到 整體區(qū)域的間隔QIA， QIA = QILi U QIO U QIL2 (4) 式中QIQ、QI0和QIL2兩兩交集都為空。
7. 根據(jù)權(quán)利要求1所述的一種海陸交錯帶氣象等值線自動繪制方法，其特征在于所 述的步驟6)包括：利用MATLAB批量讀取數(shù)據(jù)，調(diào)用MATLAB等值線繪制指令contour或 contourf，將合并后的間隔QIA應(yīng)用于該指令，實現(xiàn)自動批量繪制。
【文檔編號】G06F17/50GK104102771SQ201410318436
【公開日】2014年10月15日 申請日期:2014年7月4日 優(yōu)先權(quán)日:2014年7月4日
【發(fā)明者】杜震洪, 劉仁義, 張豐, 林賢輝, 王葉晨梓 申請人:浙江大學(xué)