一、數(shù)據(jù)來(lái)源介紹

（一）行政邊界數(shù)據(jù)

本文所用到的河北唐山行政邊界數(shù)據(jù)來(lái)源于中國(guó)科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心（https://www.resdc.cn/Default.aspx）。

（二）降水量數(shù)據(jù)

本文所用到的降水量數(shù)據(jù)來(lái)源于國(guó)家地球系統(tǒng)科學(xué)數(shù)據(jù)中心（http://www.geodata.cn/）。

該網(wǎng)站注冊(cè)完成后是沒(méi)有數(shù)據(jù)的下載權(quán)限的，需要在完善個(gè)人信息后才可以下載數(shù)據(jù)，信息越完善每日可下載的數(shù)據(jù)量也就越多，認(rèn)證到哪一步視自己需求決定。

另外該網(wǎng)站下方還有許多數(shù)據(jù)中心和數(shù)據(jù)庫(kù)官方鏈接，對(duì)于正在急于尋求數(shù)據(jù)的伙伴們非常的友好。

*本教程僅供參考，不足之處歡迎斧正。

二、數(shù)據(jù)預(yù)處理

（一）行政邊界數(shù)據(jù)預(yù)處理

河北唐山行政邊界數(shù)據(jù)存放于中國(guó)地市行政邊界數(shù)據(jù)圖層當(dāng)中，我們可以利用ArcGIS提取出我們需要的地市數(shù)據(jù)。記得備份喔，備份以后在用到其他省份的地市行政邊界數(shù)據(jù)時(shí)就不用在從網(wǎng)站中下載了。

我們將中國(guó)地市行政邊界數(shù)據(jù)圖層加載到ArcGIS，因?yàn)槲覀冎恍枰颖碧粕揭粋€(gè)地市的數(shù)據(jù)，我們可以先利用識(shí)別工具識(shí)別其信息，可以了解到其FID為311。

然后右鍵圖層，打開(kāi)圖層的屬性表，利用【按屬性選擇】工具選中河北唐山，接下來(lái)可以將選中的要素導(dǎo)出或者反選將其余要素刪掉，這里建議是導(dǎo)出比較好，選中后右鍵圖層，選擇數(shù)據(jù) -> 導(dǎo)出數(shù)據(jù)即可。

?

（二）降水量數(shù)據(jù)預(yù)處理

網(wǎng)站下載下來(lái)的年降水量數(shù)據(jù)是一個(gè)TIFF文件，其包含20個(gè)波段，每個(gè)波段都包含一個(gè)年份的全國(guó)年降水量，由于體量太大，這里在進(jìn)行裁剪處理時(shí)可以右鍵裁剪選擇批量處理。

工具位置：數(shù)據(jù)管理工具 -> 柵格 -> 柵格處理 -> 裁剪

?

網(wǎng)站下載下來(lái)的逐月降水量數(shù)據(jù)集是一個(gè)NC格式的文件，沒(méi)有辦法直接加載到ArcMAP中，需要用創(chuàng)建NetCDF柵格圖層工具來(lái)打開(kāi)。

工具位置：多維工具 -> 創(chuàng)建NetCDF柵格圖層

但若是像下圖中出現(xiàn)無(wú)法讀取變量、X維度、Y維度這種情況，就需要考慮到可能是因?yàn)閿?shù)據(jù)的擴(kuò)展名雖然是.nc，但實(shí)際上是HDF格式的文件，可以在文件夾中復(fù)制一個(gè)副本，修改其擴(kuò)展名為.hdf，這樣就可以直接在ArcMAP中打開(kāi)了；若是依然無(wú)法打開(kāi)，可以換一臺(tái)電腦再次嘗試使用創(chuàng)建NetCDF柵格圖層工具。

?

添加數(shù)據(jù)時(shí)雙擊【pre_2020 - 副本.hdf】可以看到該數(shù)據(jù)集包含有12個(gè)子數(shù)據(jù)集，由此我們可以推測(cè)這12個(gè)數(shù)據(jù)集所包含的數(shù)據(jù)為12個(gè)月份，每個(gè)月的月降水量。和年降水量數(shù)據(jù)一樣的原因，同樣進(jìn)行批量裁剪處理。

三、降雨侵蝕力因子簡(jiǎn)介

降雨侵蝕力因子其實(shí)是用來(lái)描述降雨對(duì)土壤侵蝕的潛在能力大小的，就是降雨的沖刷對(duì)土壤的侵蝕效應(yīng)。

到目前為止，關(guān)于降雨侵蝕力因子R的計(jì)算方法非常的繁多，不過(guò)這些計(jì)算方法總結(jié)下來(lái)，核心思路其實(shí)大概可以分為兩種：一種是將降雨的動(dòng)能引入模型來(lái)計(jì)算，因?yàn)榻涤甑膭?dòng)能越大，說(shuō)明沖刷的力度越大；另一種是直接引用月降雨量或年降雨量來(lái)計(jì)算，降雨量越大對(duì)土壤的侵蝕能力就越大。

在這里主要給大家介紹如何運(yùn)用月降水量和年降水量來(lái)計(jì)算。

其實(shí)，降雨侵蝕力因子的計(jì)算可以根據(jù)所獲取的降雨資料分別采用不同的計(jì)算公式（由于逐日的降雨量數(shù)據(jù)資料獲取權(quán)限較高，與之相關(guān)的計(jì)算公式這里暫不做介紹，想了解的GiSer可以在知網(wǎng)上多查閱下相關(guān)資料）：

四、按公式計(jì)算

（一）計(jì)算公式（3）

利用柵格計(jì)算器來(lái)進(jìn)行柵格圖層的疊加分析，在柵格計(jì)算器中按公式添加數(shù)據(jù)圖層，在輸入框中輸入如下公式：

工具位置：Spatial Analyst工具 -> 地圖代數(shù) -> 柵格計(jì)算器

1.?0.067?*?Power(("c2001_Clip"+"c2002_Clip"+"c2003_Clip"+"c2004_Clip"+"c2005_Clip"+"c2006_Clip"+"c2007_Clip"+"c2008_Clip"+"c2009_Clip"+"c2010_Clip"+"c2011_Clip"+"c2012_Clip"+"c2013_Clip"+"c2014_Clip"+"c2015_Clip"+"c2016_Clip"+"c2017_Clip"+"c2018_Clip"+"c2019_Clip"+"c2020_Clip")?/?200,1.627)??

?

這里計(jì)算多年平均降雨量時(shí)可以看到是除以了200，20年的數(shù)據(jù)為什么不是除以20呢，原因在于官網(wǎng)的數(shù)據(jù)摘要中介紹到降雨量單位為0.1mm，而公式中的單位為1mm。

（二）計(jì)算公式（4）

利用柵格計(jì)算器來(lái)進(jìn)行柵格圖層的疊加分析，在柵格計(jì)算器中按公式添加數(shù)據(jù)圖層，在輸入框中輸入如下公式：

?

五、成果展示

樹(shù)谷資料庫(kù)資源大全（2月9日更新）