參考鏈接:https://zhuanlan.zhihu.com/p/477677143
1. 數(shù)據(jù)簡介
土壤質(zhì)地類型數(shù)據(jù)來源于世界土壤數(shù)據(jù)庫(HWSD途样,Harmonized World Soil Database)陕习。該數(shù)據(jù)是空間分辨率為公里的網(wǎng)格數(shù)據(jù)畜隶,提供了各個(gè)格網(wǎng)點(diǎn)的土壤類型颜矿、土壤相位颂翼、土壤理化性狀等信息。
數(shù)據(jù)庫下載地址:https://www.fao.org/soils-portal/soil-survey/soil-maps-and-databases/harmonized-world-soil-database-v12/en/
2. 數(shù)據(jù)下載
FAO SOILS PORTAL下載網(wǎng)址內(nèi)容如上圖妻枕,其中全球土壤數(shù)據(jù)查看器只能查看全球土壤信息焦辅,并不能提供下載服務(wù)。我們只需要下載后面兩個(gè)鏈接即可婚惫。
下載解壓后文件內(nèi)容如下氛赐,包括一個(gè)BIL格式的全球柵格數(shù)據(jù),和HWSD的數(shù)據(jù)庫文件先舷。
3. 數(shù)據(jù)處理
下載的HWSD.mdb文件可以用Microsoft Access 或 MDB viwer plus軟件打開查看并導(dǎo)出xls文件艰管,也可以直接在arcmap中打開。
查看官方手冊蒋川,可知MU_GLOBAL是提供數(shù)據(jù)庫文件和柵格文件連接的ID牲芋。
SU_SYM90和CODE是在FAO90分類系統(tǒng)下的土壤類型代碼。
官方手冊說明2這里給大家附上數(shù)據(jù)說明:數(shù)據(jù)來源于聯(lián)合國糧農(nóng)組織(FAO)和維也納國際應(yīng)用系統(tǒng)研究所(IIASA)所構(gòu)建的世界土壤數(shù)據(jù)庫捺球。中國境內(nèi)數(shù)據(jù)源為第二次全國土地調(diào)查南京土壤所所提供的1:100萬土壤數(shù)據(jù)缸浦。該數(shù)據(jù)可為建模者提供模型輸入?yún)?shù),農(nóng)業(yè)角度可用來研究生態(tài)農(nóng)業(yè)分區(qū)氮兵,糧食安全和氣候變化等裂逐。數(shù)據(jù)格式:grid柵格格式,投影為WGS84泣栈。采用的土壤分類系統(tǒng)主要為FAO-90卜高。
土壤屬性表主要字段包括:
SU_SYM90(FAO90土壤分類系統(tǒng)中土壤名稱);
SU_SYM85(FAO85分類);
T_TEXTURE(頂層土壤質(zhì)地);
DRAINAGE(19.5);
REF_DEPTH(土壤參考深度);
AWC_CLASS(19.5);
AWC_CLASS(土壤有效水含量);
PHASE1: Real (土壤相位);
PHASE2: String (土壤相位)秩霍;
ROOTS: String (到土壤底部存在障礙的深度分類)篙悯;
SWR: String (土壤含水量特征);
ADD_PROP: Real (土壤單元中與農(nóng)業(yè)用途有關(guān)的特定土壤類型)铃绒;
T_GRAVEL: Real (碎石體積百分比)鸽照;
T_SAND: Real (沙含量);
T_SILT: Real (淤泥含量)颠悬;
T_CLAY: Real (粘土含量)矮燎;
T_USDA_TEX: Real (USDA土壤質(zhì)地分類);
T_REF_BULK: Real (土壤容重)赔癌;
T_OC: Real (有機(jī)碳含量)诞外;
T_PH_H2O: Real (酸堿度)
T_CEC_CLAY: Real (粘性層土壤的陽離子交換能力);
T_CEC_SOIL: Real (土壤的陽離子交換能力)
T_BS: Real (基本飽和度)灾票;
T_TEB: Real (交換性鹽基)峡谊;
T_CACO3: Real (碳酸鹽或石灰含量)
T_CASO4: Real (硫酸鹽含量);
T_ESP: Real (可交換鈉鹽);
T_ECE: Real (電導(dǎo)率)既们。
其中以T_開頭屬性字段表示上層土壤屬性(0-30cm)濒析,以S_開頭屬性字段表示下層土壤屬性(30-100cm)。
具體屬性值代表何意義請(qǐng)參考文件夾下說明文檔.pdf及數(shù)據(jù)庫.mdb啥纸。
使用arcmap打開HWSD.bil文件
由于原始柵格沒有空間參考信息号杏,所以需要定義一下投影,選擇WGS_1984地理坐標(biāo)系(可根據(jù)需要修改斯棒,我下到的好像有)盾致。
為柵格創(chuàng)建屬性表
打開arctoolbox→data management tools →raster → build raster attribute table
連接屬性表
右鍵圖層 →joins and relates →join
點(diǎn)擊第二行那個(gè)文件夾,點(diǎn)擊HWSD.mdb文件荣暮,找到HWSD_DATA雙擊
連接字段選擇MU_GLOBAL庭惜,點(diǎn)擊OK就連接好了
單因子?xùn)鸥駡D層導(dǎo)出
使用3D analyst tools工具包
3D analyst tools → raster reclass → lookup
之后選擇相應(yīng)的字段提取相應(yīng)的因子圖層
這樣就導(dǎo)出了柵格圖層,在導(dǎo)出完成之后穗酥,在arcmap中打開該圖層右鍵