Function to read in netCDF file, calculate LAI and save it back in a netCDF file.

select_data_to_LAI.py

import xarray as xr
import satellitetools.biophys_xarray as bio
from satellitetools import gee
import geopandas as gpd
import os

def select_data_LAI(pathtoinputfile, outdir):
    
    ds_disk = xr.open_dataset(pathtoinputfile)
    #calculating LAI
    area = bio.run_snap_biophys(ds_disk, "LAI")
    
    timeseries = {}
    timeseries_variable = ["lai"]
    
    if not os.path.exists(outdir):
        os.makedirs(outdir, exist_ok=True)

    # creating a timerseries and saving the netCDF file
    area.to_netcdf(os.path.join(outdir, area.name + ".nc"))
    timeseries[area.name] = gee.xr_dataset_to_timeseries(area, timeseries_variable)