做水文分析的同行應(yīng)該都有過(guò)這種經(jīng)歷手里的DEM明明處理得好好的填洼、流向、流量一環(huán)扣一環(huán)眼看就要出河網(wǎng)了結(jié)果ArcGIS Pro里的工具跑了幾秒鐘彈出一個(gè)紅色報(bào)錯(cuò)——999999或者010067任務(wù)直接翻車。更氣人的是同一個(gè)數(shù)據(jù)、同一條流程換臺(tái)電腦跑就是好的換回自己這臺(tái)就掛。如果你也遇到過(guò)這種“玄學(xué)失敗”先別急著懷疑數(shù)據(jù)十有八九是并行處理因子在搗鬼。這篇文章就把這個(gè)坑徹底講透從原理到操作再到排查流程一次說(shuō)清。我最早被這個(gè)問(wèn)題絆倒是在某次幫一個(gè)項(xiàng)目組處理山區(qū)小流域的河網(wǎng)提取。當(dāng)時(shí)數(shù)據(jù)量不算特別大幾百萬(wàn)像元的DEM按理說(shuō)ArcGIS Pro跑起來(lái)應(yīng)該毫無(wú)壓力。但奇怪的是“填洼”工具跑完沒(méi)問(wèn)題“流向”工具一執(zhí)行就報(bào)錯(cuò)而且錯(cuò)誤代碼每次都不同有時(shí)候是999999有時(shí)候是010067連輸出路徑都彈不出來(lái)。折騰了大半天最后在一位老前輩的提醒下把“并行處理因子”改成0問(wèn)題瞬間消失。后來(lái)我才意識(shí)到這玩意兒在ArcGIS Pro 3.x里尤其是水文分析這組工具里翻車率比想象中高得多。這篇文章適合誰(shuí)看凡是需要在ArcGIS Pro里做水文分析、流域劃分、河網(wǎng)提取的從業(yè)者和學(xué)生只要你用的是Pro 2.5以后的版本都可能撞上這個(gè)坑。文章會(huì)把并行處理因子的原理、三種設(shè)置路徑、以及配套的數(shù)據(jù)預(yù)處理和環(huán)境參數(shù)全部拆開(kāi)講爭(zhēng)取讓你下次再遇到水文操作失敗十分鐘內(nèi)就能定位問(wèn)題。1. 先對(duì)號(hào)入座你的失敗是不是這一種1.1 癥狀表現(xiàn)并行處理因子導(dǎo)致的水文操作失敗癥狀非常有辨識(shí)度。最典型的是工具啟動(dòng)后進(jìn)度條動(dòng)了幾下然后毫無(wú)征兆地變紅報(bào)錯(cuò)信息里經(jīng)常帶“ERROR 999999”或“ERROR 010067”有時(shí)候還夾雜著“無(wú)法創(chuàng)建輸出”“無(wú)法寫(xiě)入柵格數(shù)據(jù)集”“進(jìn)程已停止”之類的描述。有個(gè)容易被忽視的細(xì)節(jié)是這種失敗往往出現(xiàn)在“流向”和“流量”這兩個(gè)工具上而“填洼”反而很少翻車。我見(jiàn)過(guò)不少人在網(wǎng)上提問(wèn)說(shuō)填洼跑完了流向老是失敗甚至有人懷疑是DEM數(shù)據(jù)本身的問(wèn)題重新投影、裁剪、修復(fù)柵格折騰一圈依然無(wú)解。如果你也是這種“前面的工具正常后面的工具連續(xù)失敗”的節(jié)奏多半和并行計(jì)算有關(guān)。另外一個(gè)典型特征是失敗不穩(wěn)定。同一個(gè)流程今天跑掛明天又能跑通或者小范圍數(shù)據(jù)沒(méi)事?lián)Q成全流域的大范圍就掛又或者同樣數(shù)據(jù)在另一臺(tái)電腦上沒(méi)問(wèn)題。這種“隨機(jī)性”“環(huán)境相關(guān)性”的表現(xiàn)說(shuō)明問(wèn)題大概率不在數(shù)據(jù)本身而在于運(yùn)行環(huán)境的計(jì)算調(diào)度機(jī)制。1.2 為什么偏偏是水文工具容易中招ArcGIS Pro里的水文分析工具集包括填洼、流向、流量、柵格河網(wǎng)矢量化等這些工具的計(jì)算邏輯本身不算特別復(fù)雜但它們有一個(gè)共同點(diǎn)柵格運(yùn)算過(guò)程中需要頻繁訪問(wèn)鄰域像元并且存在嚴(yán)格的迭代依賴關(guān)系。拿“流向”來(lái)說(shuō)算法要確定每個(gè)像元的最陡下降方向需要反復(fù)比較中心像元與周圍八個(gè)像元的高程關(guān)系這個(gè)過(guò)程中會(huì)產(chǎn)生大量中間結(jié)果和臨時(shí)隨機(jī)訪問(wèn)。并行計(jì)算因子一旦被設(shè)置成高于0的數(shù)值系統(tǒng)會(huì)把運(yùn)算切分成多個(gè)塊交給不同線程處理。問(wèn)題就出在這里水文工具在處理分塊邊界時(shí)鄰域信息的同步很容易出現(xiàn)處理不當(dāng)分塊銜接處的像元計(jì)算就會(huì)出錯(cuò)嚴(yán)重時(shí)直接讓整個(gè)任務(wù)崩潰。我的理解是水文工具在并行化改造上屬于“先天不良”的那一類。ArcGIS Pro逐步推行的并行運(yùn)算優(yōu)先照顧的是柵格分類、插值、條件函數(shù)這類可以天然分塊的工具而水文分析這種強(qiáng)迭代強(qiáng)依賴的工具并行化收益有限反而更容易踩到線程同步的雷區(qū)。1.3 誰(shuí)最容易被這個(gè)坑絆住從我自己接觸到的案例來(lái)看幾類場(chǎng)景特別容易觸發(fā)問(wèn)題。第一類是處理高分辨率大范圍DEM比如30米分辨率覆蓋整個(gè)地級(jí)市的DEM數(shù)據(jù)像元數(shù)量輕松上億這種情況系統(tǒng)很容易切分大量計(jì)算塊并行沖突的概率直線上升。第二類是機(jī)器配置比較好的用戶核心數(shù)越多越容易觸發(fā)ArcGIS Pro默認(rèn)會(huì)根據(jù)CPU核數(shù)自動(dòng)計(jì)算并行因子一個(gè)16核的機(jī)器可能默認(rèn)就給到8甚至更高這反而增加了故障概率。第三類是使用某些特定柵格存儲(chǔ)格式的用戶。ArcGIS Pro原生推薦的云柵格格式在并行寫(xiě)操作上存在一些兼容性問(wèn)題我之前有幾次就是輸出格式選成了云柵格格式導(dǎo)致報(bào)錯(cuò)后來(lái)改成TIFF就一切正常。如果你屬于以上三者之一這篇文章接下來(lái)的內(nèi)容值得認(rèn)真讀完。2. 并行處理因子到底是啥它和水文操作怎么結(jié)下的梁子2.1 工作機(jī)制簡(jiǎn)單說(shuō)并行處理因子英文界面里叫Parallel Processing Factor是ArcGIS地理處理工具的一個(gè)環(huán)境參數(shù)用來(lái)控制工具允許使用多少個(gè)處理器核心來(lái)執(zhí)行計(jì)算。它位于“環(huán)境設(shè)置”選項(xiàng)卡里常見(jiàn)的可填選項(xiàng)有0、50、100以及任意百分比或整數(shù)。這里的邏輯并不復(fù)雜0表示不進(jìn)行并行計(jì)算或者說(shuō)以最低并行度運(yùn)行50表示使用一半的計(jì)算資源100表示全部核心參與默認(rèn)值通常就是100。如果你填一個(gè)具體的數(shù)字比如4那就表示讓工具只使用4個(gè)核心。有些朋友可能會(huì)以為設(shè)成0是“讓程序只用0個(gè)核心”這其實(shí)是誤解。更準(zhǔn)確地說(shuō)0在ArcGIS Pro的這套并行邏輯里代表關(guān)閉并行計(jì)算調(diào)度讓工具用最保守的順序計(jì)算方式執(zhí)行任務(wù)。這個(gè)設(shè)定在早期版本里效果尤其明顯因?yàn)榕f版本水文工具本身的并行分支不成熟關(guān)掉之后反而走的是最穩(wěn)的代碼路徑。2.2 水文分析為什么容易觸發(fā)并行故障并行計(jì)算的大原則是分而治之把整個(gè)柵格按行或按塊切碎分給多個(gè)線程分別計(jì)算最后再拼起來(lái)。聽(tīng)起來(lái)很美好但水文分析有個(gè)硬傷——它的每一步計(jì)算都強(qiáng)依賴鄰域和上游結(jié)果。拿“流量”工具舉例每個(gè)像元的累積流量等于上游所有匯入它的像元的流量之和。這個(gè)計(jì)算存在明顯的先后依賴上游沒(méi)算完下游就不能開(kāi)工。如果并行調(diào)度把下游的塊先派給某個(gè)線程處理而那個(gè)線程又拿不到上游的最新結(jié)果計(jì)算就會(huì)出錯(cuò)或者需要反復(fù)通信等待線程越多這種協(xié)調(diào)成本越高出錯(cuò)概率也越大。更要命的是Sink和Flat區(qū)域的遞歸處理填洼和流向工具里涉及大量迭代修正。當(dāng)數(shù)據(jù)里有大面積平坦區(qū)域或者錯(cuò)誤洼地時(shí)迭代次數(shù)會(huì)激增并行計(jì)算下的分塊邊界處特別容易發(fā)生數(shù)據(jù)不同步。我懷疑這也是為什么“流向”工具在高并行度設(shè)置下最容易莫名其妙的失敗——它正是水文工具中迭代邏輯最重的那一個(gè)。2.3 設(shè)0為什么被視為“萬(wàn)能解藥”把并行處理因子設(shè)為0等于強(qiáng)制工具走單線程、順序處理的舊路徑。雖然性能會(huì)下降一些但對(duì)于水文分析這種計(jì)算量不算特別夸張的任務(wù)來(lái)說(shuō)犧牲一點(diǎn)時(shí)間換穩(wěn)定完全值得。尤其是當(dāng)你處理的數(shù)據(jù)本身不是TB級(jí)別時(shí)并行計(jì)算帶來(lái)的性能提升其實(shí)感知不強(qiáng)。我自己實(shí)測(cè)過(guò)一個(gè)千萬(wàn)級(jí)像元的DEM跑流向并行因子100的時(shí)候跑掛的概率大約三成設(shè)成0以后雖然耗時(shí)多了約四分之一但連續(xù)跑十個(gè)流程全都穩(wěn)穩(wěn)落地。穩(wěn)定性優(yōu)先的原則在水文分析里永遠(yuǎn)是對(duì)的。還有一個(gè)隱藏好處設(shè)成0之后ArcGIS Pro不會(huì)自動(dòng)創(chuàng)建額外的臨時(shí)并行緩存文件這對(duì)磁盤(pán)空間緊張的環(huán)境也是減壓。我之前在某臺(tái)磁盤(pán)只剩幾個(gè)GB的機(jī)器上就是靠這個(gè)參數(shù)救回來(lái)的。3. 實(shí)操三招把并行處理因子改成03.1 第一招在工具對(duì)話框里臨時(shí)生效最直接的辦法就是每次運(yùn)行水文工具時(shí)在工具對(duì)話框右側(cè)展開(kāi)“環(huán)境”選項(xiàng)卡找到“并行處理”這一項(xiàng)把“并行處理因子”從默認(rèn)值改成0。這個(gè)操作只對(duì)當(dāng)前這次運(yùn)行生效不會(huì)影響其他工具和其他流程非常干凈。需要注意一個(gè)小細(xì)節(jié)在部分ArcGIS Pro版本里“并行處理因子”輸入框的值默認(rèn)是空的顯示為“無(wú)”但你點(diǎn)開(kāi)旁邊的下拉箭頭能看到0、50、100等預(yù)設(shè)值。直接選0就行。如果界面顯示的是百分比輸入框輸入0再回車效果一樣。這個(gè)方法適合偶爾遇到問(wèn)題時(shí)臨時(shí)驗(yàn)證比如你懷疑當(dāng)前這個(gè)報(bào)錯(cuò)是不是并行導(dǎo)致的先用臨時(shí)方式驗(yàn)證一下確認(rèn)有效再去改全局設(shè)置避免盲目改配置影響其他正常工作。3.2 第二招在全局地理處理選項(xiàng)里一勞永逸如果你確定自己的機(jī)器跑水文工具經(jīng)常出問(wèn)題那最好直接在ArcGIS Pro的全局設(shè)置里把“并行處理因子”默認(rèn)值改成0。具體路徑是頂部菜單“工程” - “選項(xiàng)” - “地理處理” - “環(huán)境”選項(xiàng)卡往下找到“并行處理”把因子改為0。這樣設(shè)置之后你后續(xù)打開(kāi)的所有地理處理工具環(huán)境里都會(huì)默認(rèn)帶著這個(gè)0值不需要每次手動(dòng)調(diào)也不用擔(dān)心忘記設(shè)置導(dǎo)致再次翻車。這里我個(gè)人的建議是如果你主要是做水文分析、流域提取這類型的工作干脆全局保持0也別心疼如果平時(shí)還要跑一些大數(shù)據(jù)量插值和分類那再單獨(dú)給這類工具臨時(shí)調(diào)回50或者100也不遲。提示全局修改并行處理因子是“木桶策略”只影響新開(kāi)工具默認(rèn)值已經(jīng)在運(yùn)行的任務(wù)不會(huì)受影響。改完以后要新開(kāi)工具對(duì)話框才能看到變化。3.3 第三招腳本與自動(dòng)化批處理場(chǎng)景如果你平時(shí)喜歡用Python腳本或者M(jìn)odelBuilder做水文處理流水線那設(shè)置并行處理因子就得放在代碼里或者模型參數(shù)里。ArcPy里的寫(xiě)法很標(biāo)準(zhǔn)通過(guò)arcpy.env.parallelProcessingFactor賦值即可import arcpy arcpy.env.workspace rD:\水文項(xiàng)目\流域提取.gdb arcpy.env.parallelProcessingFactor 0 arcpy.env.cellSize 10 # 填洼 filled_dem arcpy.sa.Fill(dem_10m) filled_dem.save(fill_dem) # 流向 flow_dir arcpy.sa.FlowDirection(filled_dem) flow_dir.save(flow_dir)注意arcpy.env.parallelProcessingFactor賦值的是一個(gè)字符串0不是數(shù)字0這個(gè)細(xì)節(jié)容易踩坑。賦值之后后續(xù)調(diào)用任何地理處理工具都會(huì)按照這個(gè)環(huán)境值執(zhí)行包括arcpy.sa.Fill、arcpy.sa.FlowDirection、arcpy.sa.FlowAccumulation這些水文工具。用ModelBuilder的朋友也不用愁在模型環(huán)境屬性里同樣可以設(shè)置并行處理因子。打開(kāi)模型屬性切到“環(huán)境”選項(xiàng)卡勾選并行處理并賦值為0整個(gè)模型里所有工具就都統(tǒng)一了。這個(gè)做法適合做標(biāo)準(zhǔn)化的自動(dòng)處理流程省得每次改。3.4 三種方式的對(duì)比為了讓你看得更清楚我把三種方式和適用場(chǎng)景整理成一張表設(shè)置方式生效范圍適用場(chǎng)景注意事項(xiàng)工具對(duì)話框環(huán)境選項(xiàng)卡僅當(dāng)前一次運(yùn)行臨時(shí)驗(yàn)證、快速定位問(wèn)題需要每次手動(dòng)設(shè)置工程選項(xiàng)里的地理處理環(huán)境當(dāng)前工程內(nèi)新建的所有工具本地機(jī)器上長(zhǎng)時(shí)間做水文分析可能影響其他工具性能arcpy.envy / 模型環(huán)境屬性腳本運(yùn)行期或整個(gè)模型批處理、自動(dòng)化流水線腳本里要賦值字符串0三者的優(yōu)先級(jí)也有講究工具對(duì)話框內(nèi)的環(huán)境設(shè)置高于腳本環(huán)境腳本環(huán)境又高于全局工程選項(xiàng)。如果你全局設(shè)成0但某個(gè)工具對(duì)話框里手動(dòng)填了100那實(shí)際運(yùn)行還是會(huì)按100來(lái)執(zhí)行。所以有時(shí)候你明明全局改了0還報(bào)錯(cuò)要么是工具級(jí)設(shè)置覆蓋了要么是腳本里又賦了新值。4. 只改并行因子還不夠附贈(zèng)一套水文操作保平安流程4.1 數(shù)據(jù)處理前的姿態(tài)投影、裁剪、修邊界水文分析對(duì)數(shù)據(jù)有一個(gè)先決條件DEM必須是投影坐標(biāo)系絕對(duì)不能直接用經(jīng)緯度的地理坐標(biāo)系去做填洼和流向。原因很簡(jiǎn)單流向算法里要計(jì)算水平距離和垂直距離的比值來(lái)決定最陡下降方向經(jīng)緯度坐標(biāo)下單位不一致算出來(lái)的坡度方向就是錯(cuò)的后續(xù)所有結(jié)果都不可信。我見(jiàn)過(guò)不少人跳過(guò)投影這一步直接跑流向工具雖然不報(bào)錯(cuò)但出來(lái)的河網(wǎng)支離破碎、流向錯(cuò)亂還以為是并行問(wèn)題。正確的做法是先把原始DEM重投影到適合目標(biāo)區(qū)域的投影坐標(biāo)系比如常見(jiàn)的UTM或者高斯克呂格投影然后按研究區(qū)范圍裁剪并檢查DEM邊界上有沒(méi)有負(fù)異常值或者邊緣噪點(diǎn)。裁剪這一步尤其重要。很多時(shí)候水文工具失敗不是因?yàn)樗惴ǘ且驗(yàn)镈EM在裁切邊界處出現(xiàn)了孤立像元或者NoData區(qū)域?qū)е锣徲蛴?jì)算時(shí)遇到非法值。所以在填洼之前先把DEM的NoData區(qū)域和邊界檢查一遍最好用邊界清洗工具把邊緣的小塊異常像元處理掉。4.2 填洼和流向的配套環(huán)境參數(shù)除了并行處理因子水文工具組的另一個(gè)關(guān)鍵環(huán)境參數(shù)是“像元大小”。我建議在跑水文流程之前顯式地把環(huán)境里的cellSize設(shè)為DEM原始分辨率避免系統(tǒng)自動(dòng)根據(jù)輸出范圍調(diào)整出一個(gè)莫名其妙的值。還有“處理范圍”這個(gè)參數(shù)如果輸出數(shù)據(jù)和輸入數(shù)據(jù)范圍不一致流向和流量結(jié)果很容易在邊緣出現(xiàn)大面積NoData。最好把處理范圍鎖定為DEM數(shù)據(jù)的實(shí)際范圍或者用“與柵格分析相同”選項(xiàng)確保每次輸出都在同一坐標(biāo)系和邊界下對(duì)齊。如果你處理的區(qū)域里存在大量平地或者洼地填洼的“z限制”參數(shù)也要留意。填洼有兩種模式一種是填充到最低邊界出流點(diǎn)為止一種是只填充指定深度以內(nèi)的洼地。在丘陵、平原混合區(qū)域默認(rèn)的填洼模式可能把大片平原填成“湖面”影響后續(xù)流向的準(zhǔn)確性。這個(gè)參數(shù)雖然不是環(huán)境設(shè)置但在失敗排查中經(jīng)常被忽略——因?yàn)椴⑿袌?bào)錯(cuò)把你帶偏了實(shí)際上有時(shí)候是填洼太狠導(dǎo)致下游計(jì)算溢出。4.3 輸出格式與目錄的選擇ArcGIS Pro 3.x默認(rèn)輸出柵格格式經(jīng)常落到云柵格格式這個(gè)格式在并行讀寫(xiě)時(shí)的穩(wěn)定性讓我踩過(guò)兩次坑一次是流量計(jì)算輸出到云柵格格式直接失敗一次是保存時(shí)卡在“正在寫(xiě)入柵格數(shù)據(jù)集”十幾分鐘沒(méi)反應(yīng)。后來(lái)我的習(xí)慣是水文分析的每一級(jí)輸出都顯式指定文件地理數(shù)據(jù)庫(kù)FileGDB或者干脆用TIF格式。FileGDB里的柵格在順序讀寫(xiě)時(shí)非常穩(wěn)極少出現(xiàn)鎖和臨時(shí)文件沖突TIF則便于后續(xù)在別的軟件里做進(jìn)一步分析。如果你堅(jiān)持要用云柵格格式至少先把并行處理因子設(shè)成0并且確認(rèn)磁盤(pán)空間充足不然臨時(shí)文件爆掉也會(huì)報(bào)錯(cuò)。工作目錄也別選在同步網(wǎng)盤(pán)上。某云端同步目錄的坑我印象很深工具每次運(yùn)行到一半同步客戶端就開(kāi)始把臨時(shí)文件往云端傳鎖文件沖突任務(wù)直接被殺。把工程和中間輸出放在本地純數(shù)據(jù)盤(pán)是最省心的方案。4.4 后臺(tái)地理處理開(kāi)關(guān)64位與32位的區(qū)別ArcGIS Pro默認(rèn)使用64位后臺(tái)地理處理這意味著大多數(shù)地理處理任務(wù)能突破32位內(nèi)存上限。但如果你所在的環(huán)境里后臺(tái)地理處理被關(guān)閉或者意外切換成了32位模式水文工具在處理大范圍DEM時(shí)就會(huì)頻繁內(nèi)存不足、意外崩潰。檢查路徑工程 - 選項(xiàng) - 地理處理確認(rèn)“啟用后臺(tái)處理”是勾選狀態(tài)并且下面的“后臺(tái)處理位數(shù)”是64位。改完設(shè)置后必須重啟ArcGIS Pro才生效。這里還有一個(gè)容易被忽略的內(nèi)存參數(shù)同樣是地理處理選項(xiàng)頁(yè)里的“臨時(shí)柵格數(shù)據(jù)存儲(chǔ)位置”。ArcGIS Pro運(yùn)行水文工具時(shí)會(huì)產(chǎn)生大量臨時(shí)柵格默認(rèn)路徑可能在系統(tǒng)盤(pán)如果C盤(pán)空間不足并行計(jì)算和臨時(shí)數(shù)據(jù)寫(xiě)入會(huì)互相爭(zhēng)搶資源任務(wù)失敗概率大增。我的習(xí)慣是把臨時(shí)柵格路徑指定到大容量數(shù)據(jù)盤(pán)比如D盤(pán)專門建一個(gè)“臨時(shí)柵格緩存”文件夾實(shí)測(cè)能減少不少莫名其妙的報(bào)錯(cuò)。5. 常見(jiàn)報(bào)錯(cuò)速查與實(shí)戰(zhàn)排查經(jīng)驗(yàn)5.1 高頻報(bào)錯(cuò)對(duì)照表報(bào)錯(cuò)內(nèi)容可能原因處理建議ERROR 999999通用錯(cuò)誤原因復(fù)雜先改并行因子為0再檢查磁盤(pán)空間和范圍一致性ERROR 010067無(wú)法創(chuàng)建輸出柵格數(shù)據(jù)集檢查輸出路徑切換輸出格式必要時(shí)重命名輸出名內(nèi)存不足/進(jìn)程終止32位后臺(tái)處理或內(nèi)存配置過(guò)低確認(rèn)64位后臺(tái)處理開(kāi)啟減少同時(shí)運(yùn)行的大任務(wù)工具運(yùn)行后輸出全為NoData處理范圍不一致或投影坐標(biāo)系錯(cuò)誤統(tǒng)一處理范圍重投影DEM為投影坐標(biāo)系在云柵格格式寫(xiě)入時(shí)卡死并行寫(xiě)入與格式兼容性問(wèn)題輸出換為FileGDB或TIF同時(shí)并行因子設(shè)0這些報(bào)錯(cuò)里ERROR 999999是最氣人的因?yàn)樗旧硎且粋€(gè)“萬(wàn)能錯(cuò)誤”本身不告訴你真正的原因。我自己遇到999999時(shí)排查順序是固定的先改并行因子為0不行再看磁盤(pán)空間再不行看輸出路徑是否包含中文或特殊字符最后才懷疑數(shù)據(jù)本身。這個(gè)順序踩了很多次坑才總結(jié)出來(lái)。5.2 一個(gè)真實(shí)的排查過(guò)程有一次幫人在ArcGIS Pro 3.0上處理一個(gè)流域的河網(wǎng)提取到了“流向”這一步連續(xù)三次報(bào)ERROR 999999。第一次我先嘗試把并行因子設(shè)成0結(jié)果還是報(bào)錯(cuò)于是我的排查鏈條繼續(xù)走。第二步我檢查了工作目錄發(fā)現(xiàn)工程放在了某云盤(pán)同步目錄下立刻把整個(gè)工程復(fù)制到本地純數(shù)據(jù)盤(pán)重新跑流向成功了但“流量”又卡住了。接下來(lái)我懷疑是中間臨時(shí)文件的問(wèn)題把菜單選項(xiàng)里的臨時(shí)柵格位置也改到數(shù)據(jù)盤(pán)再把流量輸出格式從云柵格格式換成TIF整個(gè)過(guò)程立刻順暢了。那次經(jīng)歷讓我深刻意識(shí)到并行因子設(shè)0是第一步但光設(shè)0不能解決所有環(huán)境類問(wèn)題整個(gè)運(yùn)行鏈路的穩(wěn)定性要一起檢查。5.3 避坑心得設(shè)置完還要留個(gè)心眼有些朋友把并行因子設(shè)成0后跑通了小范圍數(shù)據(jù)就放心地拿著大面積數(shù)據(jù)直接開(kāi)跑結(jié)果半夜起來(lái)看發(fā)現(xiàn)又掛了。這種情況多半不是并行因子的問(wèn)題而是數(shù)據(jù)量太大導(dǎo)致臨時(shí)文件占用過(guò)多空間。我的建議是大范圍高分辨率的水文分析拆成區(qū)塊處理是常態(tài)。先按子流域或者方格網(wǎng)把大范圍DEM切成幾塊分別跑完水文流程再通過(guò)柵格拼接工具合并。表面上看多了一步操作但每一塊的失敗風(fēng)險(xiǎn)整體降低不少反而比一次性跑完全部更省時(shí)間。另外ArcGIS Pro有時(shí)候改完環(huán)境設(shè)置后新開(kāi)工具對(duì)話框才能生效。如果你在現(xiàn)有工具對(duì)話框里改了環(huán)境值但沒(méi)重新打開(kāi)工具用的還是舊值這個(gè)問(wèn)題非常隱蔽。我的做法是每次調(diào)完并行因子就把當(dāng)前工具對(duì)話框關(guān)掉重新打開(kāi)確?？吹降木褪菍?shí)際用的。還有一個(gè)容易忽視的小細(xì)節(jié)如果你在ArcGIS Pro的Python窗口里直接執(zhí)行工具窗口本身有自己的環(huán)境狀態(tài)有時(shí)候你之前測(cè)試代碼時(shí)設(shè)過(guò)env值但之后手動(dòng)改了全局設(shè)置Python窗口里的舊值還在生效。最簡(jiǎn)單的辦法是把Python窗口關(guān)掉重新打開(kāi)或者顯式在代碼開(kāi)頭重新賦值并行因子。最后再補(bǔ)一個(gè)經(jīng)驗(yàn)ArcGIS Pro的版本差異很大2.x和3.x在并行處理的穩(wěn)定性和選項(xiàng)文案上都有變化。每次升級(jí)之后水文工具第一次跑之前我都習(xí)慣性打開(kāi)環(huán)境看一下并行因子的默認(rèn)值別指望系統(tǒng)自己記住你的偏好。6. 把并行處理因子日?；癁椤澳J(rèn)動(dòng)作”我在實(shí)際工作中已經(jīng)養(yǎng)成了條件反射只要下一步要跑水文分析相關(guān)工具先看并行處理因子是不是0。這不是矯枉過(guò)正而是在無(wú)數(shù)次凌晨排查報(bào)錯(cuò)里換來(lái)的教訓(xùn)?？赡苡腥藭?huì)擔(dān)心并行因子設(shè)0會(huì)讓大范圍柵格處理慢到不可接受。其實(shí)在我的數(shù)據(jù)規(guī)模下——一般也就一億像元以內(nèi)——差異并不明顯。水文分析工具本身并不是原生并行收益很高的工具穩(wěn)定運(yùn)行比節(jié)省幾分鐘重要得多。如果你的數(shù)據(jù)確實(shí)到了需要并行才能跑完的級(jí)別我建議你也別死磕水文工具本身先去想想是不是應(yīng)該對(duì)DEM做分幅、投影重采樣或者降采樣數(shù)據(jù)處理策略比工具參數(shù)更能解決問(wèn)題。這里再分享一個(gè)我后來(lái)經(jīng)常用的小技巧把整個(gè)水文流程的“環(huán)境參數(shù)”存成ArcGIS Pro的一個(gè)模板工程。打開(kāi)ArcGIS Pro后直接從這個(gè)模板工程啟動(dòng)所有環(huán)境設(shè)置都是預(yù)設(shè)好的包括并行處理因子0、臨時(shí)柵格位置、輸出坐標(biāo)系、像元大小連默認(rèn)輸出目錄都幫你固定好。這樣你每次拿到新數(shù)據(jù)只需要加載DEM、按流程跑工具不需要再花時(shí)間檢查環(huán)境設(shè)置。我用這個(gè)辦法之后水文分析的第一次成功率從原來(lái)的六成提到了九成以上剩下的一成基本是數(shù)據(jù)本身的質(zhì)量問(wèn)題。做水文分析這些年越來(lái)越覺(jué)得這行比拼的早就不是“會(huì)不會(huì)用工具”而是“出了問(wèn)題找不找得到原因”。一個(gè)人對(duì)著滿屏報(bào)錯(cuò)毫無(wú)頭緒的時(shí)候能快速定位到是環(huán)境的哪個(gè)參數(shù)在攪局才是真正的核心競(jìng)爭(zhēng)力。如果你現(xiàn)在正被ArcGIS Pro水文操作的某個(gè)旋渦拖著走不妨先把所有“聰明的并行計(jì)算”關(guān)掉讓工具老老實(shí)實(shí)一步一步算結(jié)果反而會(huì)給你驚喜。