Rui Gao: [email protected] | [email protected]

• This python script is used for the AggieAir-image processing to support the TSEB model.
• This repository can be used for the user who has high-resolution data, including red-band optical image, near-infrared-band optical image, DSM image, thermal image, and LAI image (map). The results contain the image to run the TSEB-PT model.
• This repository is not written in a python function since the some results from some cells can help researchers to find more information, and it is easier for users to debug when apply this script for new research objectives.

These involved python functions are used to extract information from the corresponding inputs. They could be found in other repositories in my GitHub. Details for each function can be found there.

• CanopyHeight_Vine: [Vine_Info_Extraction]| canopy height at a subfield scale is generated at the end as one image.
• CanopySoilTemperatureSeparation_Vine: [Vine_Info_Extraction]| canopy and soil temperature (and also the correlation coefficient between NDVI and temperature) at a subfield scale is generated at the end as one image.
• TellResolution: [Rui_Python_Functions_Package]| resolution (2 values) is generated.
• TellTheGeoInfo: [Rui_Python_Functions_Package]| (1) dimension, (2) geographic information, and (3) projection information of the input image are generated.
• FolderCreater: [Rui_Python_Functions_Package]| a python function used to generated a folder.
• WriteTiffData: [Rui_Python_Functions_Package]| put required information and the array into this function, image in "Tiff" format is generated.
• TellExtent: [Rui_Python_Functions_Package]| the north, south, east, and west (4 values) are gained based on the input image.

• file_grids: the folder path of the research grids (fishnet), and the resolution I used is 3.6 meter by 3.6 meter.
  
• file_LAI: the folder path of the LAI image, and the resolution is 3.6 meter by 3.6 meter.
  
• file_RGBN: the folder path of the optical image, and the resolution is 0.1 meter by 0.1 meter.
  
• file_DSM: the folder path of the DSM image, and the resolution is 0.1 meter by 0.1 meter.
  
• file_tmp: the folder path of the thermal image, and the resolution is 0.6 meter by 0.6 meter
  
• threshold_veg: this is a threshold from NDVI to define which pixel is recognized as vegetation index. This threshold varies from flight to flight, and user can get a good guess from the cell called "Generate R, G, B, NIR, and DSM images" and then re-define the threshold when you cannot get the exact threshold.
  
• threshold_soil: this is similar to above one, and this NDVI threshold is used to define the soil pixel.
  
• threshold_veg_height: the shortest height of the vegetation (vine-canopy height for my case).
  
• threshold_thermal_upscale_ratio: the ratio between the square of the LAI pixel and the square of the thermal pixel, in string format. It is "6" for my case.
• resolution_grid: the resolution for the fishnet with only one length (e.g., 3.6 meter).
  
• resolution_rgb_dsm: the resolution for the optical and DSM data (e.g., 0.1 meter).
  
• resolution_tmp: the desired resolution of the temperature image (e.g., 0.6 meter).
• resolution_tmp_2t: the desired resolution for the mid-temperature product to get the separated temperature (canopy and soil) for TSEB-2T model.

• Canopy_Height.tif: the map for the canopy height.
• Canopy_W_H.tif: the ratio between canopy width and canopy height.
• Fractional_Cover.tif: the fractional cover of the canopy vegetation.
• Temperature_K.tif: the temperature (one layer) in degree K.

Except the result called "Temperature_K.tif", others can also be used for TSEB-2T model.

• Temperature_2L_K.tif: canopy (layer 1), soil temperature in degree K (layer 2), and correlation coefficient between NDVI and temperature (layer 3) as one image is generated.

Thanks for the suggestions and guidence from Dr. Torres

Please cite the paper below when you are using this script for your paper work.
Evapotranspiration partitioning assessment using a machine-learning-based leaf area index and the two-source energy balance model with sUAV information