面向區(qū)域尺度土壤侵蝕建模的地形因子并行計(jì)算方法
【專利摘要】本發(fā)明公開(kāi)了一種面向區(qū)域尺度土壤侵蝕建模的地形因子并行計(jì)算方法，屬于計(jì)算機(jī)和信息【技術(shù)領(lǐng)域】。該方法包括以下過(guò)程：(1)確定數(shù)據(jù)劃分策略，根據(jù)并行節(jié)點(diǎn)數(shù)目，采用行帶緩沖區(qū)劃分方法對(duì)原始地形數(shù)據(jù)進(jìn)行等行劃分，同時(shí)每個(gè)子數(shù)據(jù)塊包含了數(shù)據(jù)區(qū)、緩沖區(qū)以及控制點(diǎn)三個(gè)要素；(2)按照數(shù)據(jù)劃分策略，基于預(yù)處理后的DEM數(shù)據(jù)，進(jìn)行坡度并行計(jì)算、水流方向計(jì)算、匯流累積計(jì)算、河網(wǎng)提取計(jì)算、坡長(zhǎng)計(jì)算以及坡度坡長(zhǎng)因子計(jì)算六步。本發(fā)明在并行算法設(shè)計(jì)過(guò)程中充分考慮了串行算法的計(jì)算特征，針對(duì)局部型算法和全局型算法分別設(shè)計(jì)了不同的并行化策略，有效提高了算法的并行效率，實(shí)現(xiàn)了基于海量地形數(shù)據(jù)的地形因子快速準(zhǔn)確計(jì)算。
【專利說(shuō)明】面向區(qū)域尺度土壤侵蝕建模的地形因子并行計(jì)算方法

【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于計(jì)算機(jī)和信息【技術(shù)領(lǐng)域】，特別涉及一種面向區(qū)域土壤侵蝕建模的地形 因子并行計(jì)算方法。

【背景技術(shù)】
[0002] 土壤侵蝕是全人類共同面對(duì)的重大環(huán)境問(wèn)題，而地形要素對(duì)土壤侵蝕發(fā)生發(fā)展具 有重要的控制作用，也是現(xiàn)有土壤侵蝕模型中不可或缺的因子之一。當(dāng)前學(xué)者普通采用坡 度坡長(zhǎng)因子以表達(dá)地形要素對(duì)土壤侵蝕的影響。坡度坡長(zhǎng)因子包含了坡度因子和坡長(zhǎng)因子 兩部分，其中坡長(zhǎng)的提取是一個(gè)難點(diǎn)。Hickey, Van Remortel，張宏鳴等學(xué)者基于數(shù)字高程 模型產(chǎn)品（^Digital Elevation Models,縮寫(xiě)為DEM)提出并發(fā)展了面向侵蝕過(guò)程的坡長(zhǎng)提 取方法，并已經(jīng)在坡面尺度和流域尺度的土壤侵蝕預(yù)測(cè)中得到廣泛應(yīng)用。
[0003] 近年來(lái)為服務(wù)于區(qū)域水土保持宏觀決策、生態(tài)環(huán)境綜合治理、全球變化等問(wèn)題，在 已有的坡面尺度、流域尺度土壤侵蝕研究的基礎(chǔ)上，區(qū)域尺度土壤侵蝕相關(guān)研究得到越來(lái) 越多的重視，相關(guān)學(xué)者開(kāi)發(fā)了一系列的區(qū)域土壤侵蝕模型。在區(qū)域尺度土壤侵蝕預(yù)測(cè)與模 擬中，由于大區(qū)域高精度數(shù)據(jù)較難獲取，同時(shí)由于計(jì)算機(jī)存儲(chǔ)、計(jì)算等能力的限制，現(xiàn)有研 究通常采用以下兩種方式求取適用于區(qū)域尺度土壤侵蝕定量評(píng)價(jià)的地形指標(biāo)。一是根據(jù)侵 蝕地貌學(xué)理論擬訂替代指標(biāo)，間接表示坡度坡長(zhǎng)因子，如地形起伏度等。二是對(duì)基于中低 分辨率DEM提取的地形因子進(jìn)行尺度變換，以滿足土壤侵蝕建模的需要。然而這兩種方法 都只是權(quán)宜之計(jì)，目前在區(qū)域土壤侵蝕建模中，基于高精度地形數(shù)據(jù)的坡度坡長(zhǎng)因子高效 計(jì)算已成為制約其應(yīng)用與發(fā)展的關(guān)鍵問(wèn)題。
[0004] 隨著地理數(shù)據(jù)采集方式的進(jìn)步，特別是LiDAR、INSAR等技術(shù)的應(yīng)用，大區(qū)域高精 度DEM數(shù)據(jù)越來(lái)越豐富。數(shù)據(jù)來(lái)源的不斷擴(kuò)展更加凸顯了現(xiàn)有計(jì)算方式的落后，運(yùn)用并行 計(jì)算方法實(shí)現(xiàn)對(duì)區(qū)域尺度的坡度坡長(zhǎng)因子提取是解決這一問(wèn)題的有效方法。目前并行計(jì)算 在數(shù)字地形分析領(lǐng)域已經(jīng)得到了廣泛的應(yīng)用，特別是在流域分析方面，基于并行計(jì)算提取 河網(wǎng)、進(jìn)行流域劃分已有大量研究。然而，目前現(xiàn)有針對(duì)坡度坡長(zhǎng)因子的計(jì)算還采用串行方 法，并行計(jì)算在流域分析中的成功運(yùn)用可以借鑒到坡度坡長(zhǎng)因子計(jì)算中。


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

[0005] 本發(fā)明所要解決的技術(shù)問(wèn)題在于，針對(duì)現(xiàn)有串行算法無(wú)法高效提取面向區(qū)域尺度 土壤侵蝕建模的地形因子，提出了一種多核集群環(huán)境下坡度坡長(zhǎng)因子并行計(jì)算方法，為區(qū) 域尺度的土壤侵蝕預(yù)測(cè)與模擬提供技術(shù)支撐。
[0006] 針對(duì)上述發(fā)明目的，本發(fā)明所采用的技術(shù)方案如下：
[0007] 面向區(qū)域尺度土壤侵蝕建模的地形因子并行計(jì)算方法，包括以下過(guò)程：
[0008] (1)確定數(shù)據(jù)劃分策略，其劃分方式為：根據(jù)并行節(jié)點(diǎn)數(shù)目，采用行帶緩沖區(qū)劃分 方法對(duì)原始地形數(shù)據(jù)進(jìn)行等行劃分，同時(shí)每個(gè)子數(shù)據(jù)塊包含了三個(gè)要素：a)數(shù)據(jù)區(qū)：本計(jì) 算節(jié)點(diǎn)所分配的數(shù)據(jù)區(qū)；b)緩沖區(qū)：每個(gè)子數(shù)據(jù)塊附加了上下兩個(gè)緩沖區(qū)，對(duì)于第一子數(shù) 據(jù)塊只有下緩沖區(qū)，對(duì)于最后一個(gè)數(shù)據(jù)塊只有上緩沖區(qū)；C)控制點(diǎn)：用于定位每個(gè)數(shù)據(jù)塊 的起始讀寫(xiě)位置；
[0009] (2)按照步驟（1)的數(shù)據(jù)劃分策略，基于預(yù)處理后的DEM數(shù)據(jù)，進(jìn)行坡度并行計(jì) 算；
[0010] (3)按照步驟（1)的數(shù)據(jù)劃分策略，基于預(yù)處理后的DEM數(shù)據(jù)，進(jìn)行坡長(zhǎng)并行計(jì)算； 具體包含了以下幾個(gè)子算法：a)基于預(yù)處理后的DEM數(shù)據(jù)，進(jìn)行水流方向并行計(jì)算；b)基 于水流方向計(jì)算結(jié)果進(jìn)行匯流累積并行計(jì)算；c)基于匯流累積計(jì)算結(jié)果和設(shè)定的河網(wǎng)閾 值，進(jìn)行河網(wǎng)提取并行計(jì)算；d)基于水流方向計(jì)算結(jié)果進(jìn)行坡長(zhǎng)并行計(jì)算，同時(shí)此步驟需 要坡度數(shù)據(jù)用于坡度截?cái)鄸鸥衽袛?，河網(wǎng)水系數(shù)據(jù)用于河網(wǎng)截?cái)鄸鸥衽袛啵?br> [0011] (4)按照步驟（1)的數(shù)據(jù)劃分策略，基于坡度和坡長(zhǎng)計(jì)算結(jié)果，進(jìn)行坡度坡長(zhǎng)因子 并行計(jì)算。
[0012]上述步驟（2)具體的計(jì)算流程如下：
[0013] 步驟21.根據(jù)數(shù)據(jù)劃分策略，對(duì)輸入的DEM數(shù)據(jù)進(jìn)行數(shù)據(jù)劃分，并創(chuàng)建結(jié)果數(shù)據(jù)文 件；
[0014] 步驟22.各子節(jié)點(diǎn)遍歷輸入的DEM數(shù)據(jù)，并執(zhí)行坡度算法，計(jì)算完成后將計(jì)算結(jié)果 寫(xiě)入結(jié)果文件。
[0015] 上述步驟（3)具體的計(jì)算流程如下：
[0016] 步驟31.根據(jù)數(shù)據(jù)劃分策略，對(duì)輸入的DEM數(shù)據(jù)進(jìn)行數(shù)據(jù)劃分，并創(chuàng)建結(jié)果數(shù)據(jù)文 件；各子節(jié)點(diǎn)遍歷輸入的DEM數(shù)據(jù)，并執(zhí)行水流方向算法，計(jì)算完成后將計(jì)算結(jié)果寫(xiě)入結(jié)果 文件；
[0017] 步驟32.根據(jù)數(shù)據(jù)劃分策略，對(duì)水流方向矩陣進(jìn)行數(shù)據(jù)劃分，并創(chuàng)建匯流累積結(jié) 果文件；各子節(jié)點(diǎn)遍歷輸入的水流方向矩陣，將源點(diǎn)柵格賦值為〇同時(shí)入棧，非源點(diǎn)柵格賦 值為-m，m為鄰接匯入柵格數(shù)目，遍歷完成后依次取出棧內(nèi)柵格并沿流向計(jì)算其下游柵格 匯流累積值，當(dāng)棧內(nèi)元素為空時(shí)計(jì)算結(jié)束；
[0018] 步驟33.數(shù)據(jù)塊間更新緩沖區(qū)信息，各子數(shù)據(jù)塊已更新后的緩沖區(qū)柵格為源點(diǎn)柵 格依次更新本數(shù)據(jù)區(qū)匯流累積值；一次更新計(jì)算完成后，遍歷各子數(shù)據(jù)塊邊界行數(shù)據(jù)，若全 部柵格完成計(jì)算則計(jì)算結(jié)束，否則需要重復(fù)步驟33;計(jì)算完成后各子節(jié)點(diǎn)將計(jì)算結(jié)果寫(xiě)入 已創(chuàng)建的匯流累積文件；
[0019]步驟34.根據(jù)數(shù)據(jù)劃分策略，對(duì)輸入的匯流累積數(shù)據(jù)進(jìn)行數(shù)據(jù)劃分，并創(chuàng)建結(jié)果 數(shù)據(jù)文件；各子節(jié)點(diǎn)遍歷輸入數(shù)據(jù)，并執(zhí)行河網(wǎng)提取算法，計(jì)算完成后將計(jì)算結(jié)果寫(xiě)入結(jié)果 文件；
[0020] 步驟35.根據(jù)數(shù)據(jù)劃分策略，分別對(duì)河網(wǎng)矩陣、坡度矩陣、水流方向矩陣進(jìn)行數(shù)據(jù) 劃分，并創(chuàng)建坡長(zhǎng)結(jié)果文件；各子節(jié)點(diǎn)遍歷輸入數(shù)據(jù)，根據(jù)侵蝕學(xué)坡長(zhǎng)算法定義分別執(zhí)行柵 格點(diǎn)標(biāo)識(shí)、非累積坡長(zhǎng)計(jì)算、累積坡長(zhǎng)計(jì)算步驟；
[0021] 步驟36.數(shù)據(jù)塊間更新緩沖區(qū)信息，各子數(shù)據(jù)塊以更新后的緩沖區(qū)柵格為源點(diǎn)柵 格依次更新本數(shù)據(jù)區(qū)的坡長(zhǎng)值；一次更新計(jì)算完成后，遍歷各子數(shù)據(jù)塊邊界行數(shù)據(jù)，若全部 柵格完成計(jì)算則計(jì)算結(jié)束，否則需要重復(fù)步驟36 ;計(jì)算完成后各子節(jié)點(diǎn)將計(jì)算結(jié)果寫(xiě)入已 創(chuàng)建的坡長(zhǎng)文件。
[0022] 上述步驟⑷具體的計(jì)算流程如下：
[0023] 步驟41.根據(jù)數(shù)據(jù)劃分策略，對(duì)輸入的坡度、坡長(zhǎng)數(shù)據(jù)進(jìn)行數(shù)據(jù)劃分，并創(chuàng)建結(jié)果 數(shù)據(jù)文件；
[0024] 步驟42.各子節(jié)點(diǎn)遍歷輸入數(shù)據(jù)，并執(zhí)行坡度坡長(zhǎng)因子算法，計(jì)算完成后將計(jì)算 結(jié)果寫(xiě)入結(jié)果文件。
[0025] 本發(fā)明具有以下兩點(diǎn)優(yōu)勢(shì)：
[0026] (1)將現(xiàn)有的小流域尺度的土壤侵蝕模型地形因子計(jì)算方法擴(kuò)展到了區(qū)域尺度， 運(yùn)用并行技術(shù)克服了現(xiàn)有串行算法面對(duì)大區(qū)域海量高精度數(shù)據(jù)處理能力的瓶頸，為區(qū)域尺 度土壤侵蝕模擬提供有力支持。
[0027] (2)并行算法設(shè)計(jì)過(guò)程中充分考慮了串行算法的計(jì)算特征，針對(duì)局部型算法和全 局型算法分別設(shè)計(jì)了不同的并行化策略，有效提高了算法的并行效率。實(shí)現(xiàn)了基于海量地 形數(shù)據(jù)的地形因子快速準(zhǔn)確計(jì)算。

【專利附圖】

【附圖說(shuō)明】
[0028] 圖1是本發(fā)明中數(shù)據(jù)劃分示意圖；
[0029] 圖2是全局型算法數(shù)據(jù)依賴示意圖，圖2(a)是水流方向矩陣，圖2(b)是計(jì)算源點(diǎn) 柵格示意圖，圖2(c)是數(shù)據(jù)分塊導(dǎo)致的計(jì)算未完成示意圖；
[0030] 圖3是本發(fā)明中緩沖區(qū)更新計(jì)算策略示意圖，圖3(a)是進(jìn)程間更新緩沖區(qū)示意 圖，圖3(b)是緩沖區(qū)更新后各進(jìn)程分別計(jì)算示意圖；
[0031] 圖4是本發(fā)明中算法并行化策略示意圖；
[0032] 圖5是本發(fā)明中并行化流程示意圖；
[0033] 圖6是實(shí)施例中輸入DEM數(shù)據(jù)；
[0034] 圖7是實(shí)施例中輸出的坡度坡長(zhǎng)因子結(jié)果。

【具體實(shí)施方式】
[0035] 下面結(jié)合附圖和具體實(shí)施例，對(duì)本發(fā)明作進(jìn)一步詳細(xì)說(shuō)明。
[0036] 一、數(shù)據(jù)劃分策略
[0037] 本發(fā)明采用的是一種條帶劃分方式，其基本思想是根據(jù)并行算法的進(jìn)程數(shù)，將原 始數(shù)據(jù)劃分為等數(shù)據(jù)行的子數(shù)據(jù)塊，每個(gè)子數(shù)據(jù)塊主要包含三個(gè)元素：（1)數(shù)據(jù)區(qū)：本計(jì)算 節(jié)點(diǎn)所能分配的數(shù)據(jù)區(qū)域；（2)控制點(diǎn)：用于定位每個(gè)數(shù)據(jù)塊的起始讀寫(xiě)位置；（3)緩沖區(qū)： 每個(gè)子數(shù)據(jù)塊附加了上下兩個(gè)緩沖區(qū)（對(duì)于第一子數(shù)據(jù)塊只有下緩沖區(qū)，最后一個(gè)數(shù)據(jù)塊 只有上緩沖區(qū)）。緩沖區(qū)一方面可解決跨數(shù)據(jù)訪問(wèn)數(shù)據(jù)的問(wèn)題，同時(shí)也可用于進(jìn)程間計(jì)算結(jié) 果的傳遞。（如圖1所示）
[0038] 二、算法并行策略
[0039] 本發(fā)明涉及的六個(gè)子算法，根據(jù)其計(jì)算特征可分為局部型算法和全局型算法兩 類。對(duì)于局部型算法，每一個(gè)柵格的計(jì)算都是相對(duì)獨(dú)立的，即計(jì)算過(guò)程中柵格之間的計(jì)算結(jié) 果不具有依賴性。對(duì)于局部型算法并行化，數(shù)據(jù)劃分后各子數(shù)據(jù)塊之間不存在數(shù)據(jù)依賴，因 此不需要數(shù)據(jù)塊間的通信。而對(duì)于全局型算法，數(shù)據(jù)劃分后各子數(shù)據(jù)塊之間需要信息傳遞 以消除依賴關(guān)系。
[0040] 局部型算法可進(jìn)一步劃分為：（1)點(diǎn)運(yùn)算，如河網(wǎng)提取和因子計(jì)算兩類算法，每個(gè) 柵格計(jì)算是完全獨(dú)立的；（2)窗口運(yùn)算，如坡度和水流方向兩類算法，每個(gè)柵格計(jì)算時(shí)需要 開(kāi)辟一個(gè)分析窗口，在本發(fā)明中采用的是3*3窗口。具體實(shí)施時(shí)，局部型算法的并行化策略 如下所示：
[0041] (1)程序初始化，獲取進(jìn)程數(shù)及算法相關(guān)參數(shù)。
[0042] (2)根據(jù)進(jìn)程數(shù)拆分DEM數(shù)據(jù)。
[0043] (3)讀取數(shù)據(jù)，各進(jìn)程根據(jù)定位信息讀取本節(jié)點(diǎn)數(shù)據(jù)及緩沖區(qū)數(shù)據(jù)。
[0044] (4)算法計(jì)算，根據(jù)具體算法對(duì)各進(jìn)程的數(shù)據(jù)進(jìn)行計(jì)算。
[0045] (5)寫(xiě)出數(shù)據(jù)，根據(jù)定位點(diǎn)信息，各進(jìn)程分別將計(jì)算結(jié)果寫(xiě)入結(jié)果文件。
[0046] (6)釋放內(nèi)存，并行計(jì)算任務(wù)完成。
[0047] 局部型算法的并行化策略的偽代碼如表1所示：
[0048]表1局部型算法的并行化策略的偽代碼

【權(quán)利要求】
1. 面向區(qū)域尺度土壤侵蝕建模的地形因子并行計(jì)算方法，其特征在于，包括以下過(guò) 程： (1) 確定數(shù)據(jù)劃分策略，其劃分方式為：根據(jù)并行節(jié)點(diǎn)數(shù)目，采用行帶緩沖區(qū)劃分方法 對(duì)原始地形數(shù)據(jù)進(jìn)行等行劃分，同時(shí)每個(gè)子數(shù)據(jù)塊包含了三個(gè)要素：a)數(shù)據(jù)區(qū)：本計(jì)算節(jié) 點(diǎn)所分配的數(shù)據(jù)區(qū)；b)緩沖區(qū)：每個(gè)子數(shù)據(jù)塊附加了上下兩個(gè)緩沖區(qū)，對(duì)于第一子數(shù)據(jù)塊 只有下緩沖區(qū)，對(duì)于最后一個(gè)數(shù)據(jù)塊只有上緩沖區(qū)；c)控制點(diǎn)：用于定位每個(gè)數(shù)據(jù)塊的起 始讀寫(xiě)位置； (2) 按照步驟（1)的數(shù)據(jù)劃分策略，基于預(yù)處理后的DEM數(shù)據(jù)，進(jìn)行坡度并行計(jì)算； (3) 按照步驟（1)的數(shù)據(jù)劃分策略，基于預(yù)處理后的DEM數(shù)據(jù)，進(jìn)行坡長(zhǎng)并行計(jì)算；其 中包含了水流方向計(jì)算、匯流累積計(jì)算、河網(wǎng)提取計(jì)算以及坡長(zhǎng)計(jì)算四個(gè)子計(jì)算流程； (4) 按照步驟（1)的數(shù)據(jù)劃分策略，基于坡度和坡長(zhǎng)計(jì)算結(jié)果，進(jìn)行坡度坡長(zhǎng)因子并行 計(jì)算。
2. 根據(jù)權(quán)利要求1所述面向區(qū)域尺度土壤侵蝕建模的地形因子并行計(jì)算方法，其特征 在于，步驟（2)具體的計(jì)算流程如下： 步驟21.根據(jù)數(shù)據(jù)劃分策略，對(duì)輸入的DEM數(shù)據(jù)進(jìn)行數(shù)據(jù)劃分，并創(chuàng)建結(jié)果數(shù)據(jù)文件； 步驟22.各子節(jié)點(diǎn)遍歷輸入的DEM數(shù)據(jù)，并執(zhí)行坡度算法，計(jì)算完成后將計(jì)算結(jié)果寫(xiě)入 結(jié)果文件。
3. 根據(jù)權(quán)利要求1所述面向區(qū)域尺度土壤侵蝕建模的地形因子并行計(jì)算方法，其特征 在于，步驟（3)具體的計(jì)算流程如下： 步驟31.根據(jù)數(shù)據(jù)劃分策略，對(duì)輸入的DEM數(shù)據(jù)進(jìn)行數(shù)據(jù)劃分，并創(chuàng)建結(jié)果數(shù)據(jù)文件； 各子節(jié)點(diǎn)遍歷輸入的DEM數(shù)據(jù)，并執(zhí)行水流方向算法，計(jì)算完成后將計(jì)算結(jié)果寫(xiě)入結(jié)果文 件； 步驟32.根據(jù)數(shù)據(jù)劃分策略，對(duì)水流方向矩陣進(jìn)行數(shù)據(jù)劃分，并創(chuàng)建匯流累積結(jié)果文 件；各子節(jié)點(diǎn)遍歷輸入的水流方向矩陣，將源點(diǎn)柵格賦值為〇同時(shí)入棧，非源點(diǎn)柵格賦值 為_(kāi)m，m為鄰接匯入柵格數(shù)目，遍歷完成后依次取出棧內(nèi)柵格并沿流向計(jì)算其下游柵格匯 流累積值，當(dāng)棧內(nèi)元素為空時(shí)計(jì)算結(jié)束； 步驟33.數(shù)據(jù)塊間更新緩沖區(qū)信息，各子數(shù)據(jù)塊已更新后的緩沖區(qū)柵格為源點(diǎn)柵格依 次更新本數(shù)據(jù)區(qū)匯流累積值；一次更新計(jì)算完成后，遍歷各子數(shù)據(jù)塊邊界行數(shù)據(jù)，若全部柵 格完成計(jì)算則計(jì)算結(jié)束，否則需要重復(fù)步驟33 ;計(jì)算完成后各子節(jié)點(diǎn)將計(jì)算結(jié)果寫(xiě)入已創(chuàng) 建的匯流累積文件； 步驟34.根據(jù)數(shù)據(jù)劃分策略，對(duì)輸入的匯流累積數(shù)據(jù)進(jìn)行數(shù)據(jù)劃分，并創(chuàng)建結(jié)果數(shù)據(jù) 文件；各子節(jié)點(diǎn)遍歷輸入數(shù)據(jù)，并執(zhí)行河網(wǎng)提取算法，計(jì)算完成后將計(jì)算結(jié)果寫(xiě)入結(jié)果文 件； 步驟35.根據(jù)數(shù)據(jù)劃分策略，分別對(duì)河網(wǎng)矩陣、坡度矩陣、水流方向矩陣進(jìn)行數(shù)據(jù)劃 分，并創(chuàng)建坡長(zhǎng)結(jié)果文件；各子節(jié)點(diǎn)遍歷輸入數(shù)據(jù)，根據(jù)侵蝕學(xué)坡長(zhǎng)算法定義分別執(zhí)行柵格 點(diǎn)標(biāo)識(shí)、非累積坡長(zhǎng)計(jì)算、累積坡長(zhǎng)計(jì)算步驟； 步驟36.數(shù)據(jù)塊間更新緩沖區(qū)信息，各子數(shù)據(jù)塊以更新后的緩沖區(qū)柵格為源點(diǎn)柵格依 次更新本數(shù)據(jù)區(qū)的坡長(zhǎng)值；一次更新計(jì)算完成后，遍歷各子數(shù)據(jù)塊邊界行數(shù)據(jù)，若全部柵格 完成計(jì)算則計(jì)算結(jié)束，否則需要重復(fù)步驟36 ;計(jì)算完成后各子節(jié)點(diǎn)將計(jì)算結(jié)果寫(xiě)入已創(chuàng)建 的坡長(zhǎng)文件。
4.根據(jù)權(quán)利要求1所述面向區(qū)域尺度土壤侵蝕建模的地形因子并行計(jì)算方法，其特征 在于，步驟（4)具體的計(jì)算流程如下： 步驟41.根據(jù)數(shù)據(jù)劃分策略，對(duì)輸入的坡度、坡長(zhǎng)數(shù)據(jù)進(jìn)行數(shù)據(jù)劃分，并創(chuàng)建結(jié)果數(shù)據(jù) 文件； 步驟42.各子節(jié)點(diǎn)遍歷輸入數(shù)據(jù)，并執(zhí)行坡度坡長(zhǎng)因子算法，計(jì)算完成后將計(jì)算結(jié)果 寫(xiě)入結(jié)果文件。
【文檔編號(hào)】G06F19/00GK104392147SQ201410754205
【公開(kāi)日】2015年3月4日 申請(qǐng)日期:2014年12月10日 優(yōu)先權(quán)日:2014年12月10日
【發(fā)明者】湯國(guó)安, 劉凱, 江嶺, 蔣如喬, 黃驍力 申請(qǐng)人:南京師范大學(xué)