一、制作土壤數(shù)據(jù)的思想

作用：保存每個網(wǎng)格的基本參數(shù)和各種土壤參數(shù)。包括網(wǎng)格編號（要和植被參數(shù)文件的網(wǎng)格編號一一對應(yīng)），經(jīng)緯度等，VIC就是從這里獲取基本的地理信息的。

VIC模型的土壤參數(shù)文件，一共包含53個參數(shù)，每一個參數(shù)都要獨立制做?。?！

二、土壤質(zhì)地編號的輸出

數(shù)據(jù)準(zhǔn)備：土壤數(shù)據(jù)（寒區(qū)早區(qū)科學(xué)數(shù)據(jù)中心下載的基于世界土壤數(shù)據(jù)庫（HWSD）的中國土壤數(shù)據(jù)集）、氣象數(shù)據(jù)（國家氣象科學(xué)數(shù)據(jù)中心）、研究區(qū)網(wǎng)格數(shù)據(jù)、矢量數(shù)據(jù)

軟件：ArcGIS Pro

步驟：按掩膜提取-柵格轉(zhuǎn)點-點轉(zhuǎn)柵格-打開HWSD_DATA-柵格轉(zhuǎn)點-數(shù)據(jù)匹配-重采樣-柵格轉(zhuǎn)點-數(shù)據(jù)匹配-導(dǎo)出屬性表

1.按掩膜提?。禾崛⊙芯繀^(qū)土壤數(shù)據(jù)

2.柵格轉(zhuǎn)點

3.點轉(zhuǎn)柵格

4.打開HWSD.mdb文件中的HWSD_DATA

此處備注：由于我使用的ArcGIS Pro不支持個人地理數(shù)據(jù)庫，即.mdb格式文件，所以我又下載安裝了ArcGIS10.8及mdb批量轉(zhuǎn)gdb工具，將mdb數(shù)據(jù)庫轉(zhuǎn)為gdb數(shù)據(jù)庫。

5.柵格轉(zhuǎn)點：將研究區(qū)原始土壤數(shù)據(jù)轉(zhuǎn)點處理

6.將轉(zhuǎn)點后柵格中的gird_code值與HWSD_DATA數(shù)據(jù)庫中的MU_GLOBAL值進(jìn)行匹配

7.重采樣：將第3步點轉(zhuǎn)柵格輸出數(shù)據(jù)進(jìn)行重采樣至原始土壤數(shù)據(jù)分辨率

8.柵格轉(zhuǎn)點（第7步結(jié)果）

9.將第8步得到數(shù)據(jù)的pointid值與第6步得到數(shù)據(jù)的pointid值進(jìn)行匹配

10.導(dǎo)出屬性表

三、土壤參數(shù)文件的制備

軟件：ArcGIS Pro、RStudio

步驟：①在RStudio中，輸出上層、下層土壤質(zhì)地編號及其對應(yīng)的土壤參數(shù)、站點年均降水量；②在ArcGIS Pro中，插值掩膜轉(zhuǎn)點處理后，導(dǎo)出年均降水量、網(wǎng)格經(jīng)緯度、網(wǎng)格平均高程、及初始土壤含水量等參數(shù)的屬性表；③在RStudio中，匯總土壤參數(shù)文件需要的所有參數(shù)數(shù)據(jù)。

1.RStudio：打開所需土壤數(shù)據(jù)處理代碼，修改輸入數(shù)據(jù)（ArcGIS導(dǎo)出的數(shù)據(jù)表）、輸出數(shù)據(jù)位置，運行即可（代碼在此不詳細(xì)展示）?

（1）分別輸出上層、下層土壤質(zhì)地編號（代碼1_XJ_soil_T/S）

（2）分別輸出上層、下層質(zhì)地編號對應(yīng)的土壤參數(shù)（代碼2_Res_soil_param_T/S）：飽和導(dǎo)水率（Ksat ）、田間持水量（Wc）、凋萎含水率（Wp ）、水利傳導(dǎo)率（Expt）

（3）按站點求年均降水量（代碼3_year_prec）

2.ArcGIS Pro

（1）年均降水：將年均降水?dāng)?shù)據(jù)進(jìn)行插值處理（分辨率為0.0833333），并掩膜至研究區(qū)范圍，轉(zhuǎn)點處理后，輸出屬性表。

①導(dǎo)入數(shù)據(jù)，顯示XY數(shù)據(jù)

②插值處理（反距離權(quán)重/克里金）

③按掩膜提取

④柵格轉(zhuǎn)點：將研究區(qū)降水柵格數(shù)據(jù)

⑤導(dǎo)出屬性表

（2）網(wǎng)格經(jīng)緯度

①研究區(qū)網(wǎng)格轉(zhuǎn)點處理

②在屬性表中添加字段，右鍵添加字段選擇計算幾何

③導(dǎo)出屬性表

（3）初始土壤含水量、平均土壤氣溫、土壤中的石英含量（init_moist、Avg_T、Quartz）

①導(dǎo)入全球5分土壤數(shù)據(jù)，顯示XY數(shù)據(jù)，并導(dǎo)出點數(shù)據(jù)

②根據(jù)研究區(qū)網(wǎng)格數(shù)據(jù)，選點進(jìn)行插值

③插值處理（反距離權(quán)重/克里金）：Z字段選擇Field19-21、26、31-33，分別輸出柵格

④按掩膜提?。悍謩e對soil_moist1、soil_moist2、soil_moist3、Avg_T、QUARTZ_1、QUARTZ_2、QUARTZ_3做提取處理

⑤柵格轉(zhuǎn)點：分別對soil_moist1_E、soil_moist2_E、soil_moist3_E、Avg_T_E、QUARTZ_1_E、QUARTZ_2_E、QUARTZ_3_E進(jìn)行轉(zhuǎn)點處理

⑥導(dǎo)出屬性表（分別導(dǎo)出soil_moist1_P、soil_moist2_P、soil_moist3_P、Avg_T_P、QUARTZ_1_P、QUARTZ_2_P、QUARTZ_3_P的屬性表）

（4）平均高程（根據(jù)dem填挖數(shù)據(jù)制作）

①重采樣：與模型分辨率一致

②按掩膜提取

③柵格轉(zhuǎn)點

④導(dǎo)出屬性表

3.RStudio（4_Sum_soil_param）

利用前期制備好的數(shù)據(jù)（應(yīng)該有12個數(shù)據(jù)文件），更改“4_Sum_soil_param”程序中對應(yīng)的各數(shù)據(jù)文件地址，運行程序，匯總土壤參數(shù)文件需要的所有參數(shù)數(shù)據(jù)

在運行前，記得逐一檢查數(shù)據(jù)文件內(nèi)容及地址，以免出現(xiàn)如下及其他錯誤

最近在家冒著被批評的風(fēng)險一心投入學(xué)習(xí)中（不怎么干家務(wù)，經(jīng)常坐在電腦前且熬夜），這部分學(xué)習(xí)了兩天，效率不高，預(yù)計還有一半內(nèi)容，盡快學(xué)完，加油?。?！

此處分享一個在解決“mdb數(shù)據(jù)庫轉(zhuǎn)gdb數(shù)據(jù)庫”問題時，遇到的一件趣事~

起因是去TB下單買了一個處理工具，店家一直不發(fā)下載地址，大概過了兩小時（我已經(jīng)在WX找到了免費版），來問我還要不要，于是發(fā)生了如下對話~~~