GEE python: RUSLE土壤侵蚀模型的代码
目录
简介
代码1
代码2
代码3
代码4
代码5
代码6
简介
该模型完全由 Python 代码编写,需要借助 Python API for GEE 对其进行转换,以便在 GEE 上运行。 模型由 3 个脚本组成: Input_File.py:从 GEE 中获取输入数据并导出为 tif 文件,然后读取这些文件并将输出转换为数组。 完成所需的处理(重投影、选择测试区域等)。 Adj_RUSLE.py: 这是主要模型,根据文献中的不同方法计算 RUSLE 因子 Output_File.py: 此外,侵蚀模型还需要导入 2 个函数才能运行:erosivity.so 和 erosivity.f90 scaling_func.py
代码1
# -*- coding: utf-8 -*-
import numpy as np
from netCDF4 import Dataset
from sys import argv#code Adjusted RUSLE model as presented in Naipal et al., 2015., 2016def K_factor(fname,llx,urx,lly,ury,lat_h,lon_h): '''This function ca