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., 2016

def K_factor(fname,llx,urx,lly,ury,lat_h,lon_h): 

    '''
    This function ca