getPscale_urban: calculate Pscale for urbanVPRM

In Timothy-W-Hilton/VPRMLandSfcModel: R implementation of VPRM with parameter estimation

calculate Pscale for urbanVPRM

Description

calculates Pscale according to SI eqn 4 in Hardiman et al, 2017

Usage

getPscale_urban(EVI)

Arguments

EVI	1xN numeric vector; enhanced vegetation index

Value

Pscale term in urbanVPRM GEE equation (eqn SI eqn 4 in Hardiman et al, 2017)

Author(s)

Timothy W. Hilton

References

Hardiman, B. S., Wang, J. A., Hutyra, L. R., Gately, C. K., Getson, J. M., & Friedl, M. A. (2017). Accounting for urban biogenic fluxes in regional carbon budgets. Science of The Total Environment, 592, 366–372. https://doi.org/10.1016/j.scitotenv.2017.03.028

Examples

data(Park_Falls)
pfa_dd <- VPRM_driver_data(name_long="Park Falls",
                           name_short = "US-PFa",
                           lat=45.9459,
                           lon=-90.2723,
                           PFT='MF',
                           tower_date=PFa_tower_obs[['date']],
                           NEE_obs=PFa_tower_obs[['FC']],
                           T=PFa_tower_obs[['TA']],
                           PAR=PFa_tower_obs[['PAR']],
                           date_nir = PFa_refl[['date']],
                           rho_nir=PFa_refl[['nir']],
                           date_swir = PFa_refl[['date']],
                           rho_swir = PFa_refl[['swir']],
                           date_EVI = PFa_evi[['date']],
                           EVI=PFa_evi[['evi']],
                           phen=NA)
pfa_df <- as.data.frame(pfa_dd)
pscale <- getPscale_urban(pfa_df[['EVI']])

Timothy-W-Hilton/VPRMLandSfcModel documentation built on July 29, 2023, 8:43 p.m.