hectornieto / pyTSEB

A python Two Source Energy Balance model for estimation of evapotranspiration with remote sensing data
GNU General Public License v3.0
133 stars 62 forks source link

Can I use the ProcessLocalImage.ipynb to derive ET of urban scenes based on UAV images? #53

Closed HazelSCUT closed 1 year ago

HazelSCUT commented 1 year ago

Dear colleagues,

Recently, I have tried to use ProcessLocalImage.ipynb to process my urban data that were derived from UAV. I used the model of Priestley Taylor. However, the results seems to have many abnormal values (i.e. NaN): The black pixels of the following pictures show NaN, and such areas are almost trees. I don't know whether can I use this method to derive the spatial ET of urban scenes.

sensible heat flux latent heat flux soil heat flux net radiation test area

Thanks.

HazelSCUT commented 1 year ago

By the way, the pictures respectively are: Sensible heat flux, Latent heat flux, Soil heat flux, Net radiation and test area.

hectornieto commented 1 year ago

Usually one should not apply the model to urban environments, especially in thos very patchy conditinos. The model assumes that pixels are not affected by neiboring conditions, which in the case of urban environmnents, buildins are affecting radiaiton and turbulence (eddies). I.e. a rather homogeneity in the scene and medium scale is expected. Your AOI is extremele heterogeneous and buildings around, but alos tiles and walkpaths are affecting the micrometeorology of the green patches. So I would not advice use this model for your purposes. In any case, I do not think the NaN are caused by this issue. So if you want me to dig deeper in what is causing the problem, please share your configuration data and and example (subset) of the inputs you are using

HazelSCUT commented 1 year ago

Hello,

Thanks for your reply. Ok, I understand it.

However, I want to know why NaN values occur. The subset of my dataset is attached. Other configurations are also attached. Please check it.

PS: 1) I did not separately mask each type of vegetation, instead I use the whole image. 2) As the spatial resolution of my data is 7cm, so the fractional vegetation cover was assumed to be 1.0. I don't know whether it is correct. By the way, I have attached the fraction vegetation cover.tif that was calculated based on NDVI. 3) The Green Fraction data (greeneryindex.tif) was calculated based on NDVI. I have also attached NDVIimage.tif that has negative values due to the existence of impermeable surfaces. 4)Other parameters in the configuration were set as default. The configuration file was also attached.

Thanks.

------------------ 原始邮件 ------------------ 发件人: "hectornieto/pyTSEB" @.>; 发送时间: 2023年3月20日(星期一) 晚上7:14 @.>; @.**@.>; 主题: Re: [hectornieto/pyTSEB] Can I use the ProcessLocalImage.ipynb to derive ET of urban scenes based on UAV images? (Issue #53)

Usually one should not apply the model to urban environments, especially in thos very patchy conditinos. The model assumes that pixels are not affected by neiboring conditions, which in the case of urban environmnents, buildins are affecting radiaiton and turbulence (eddies). I.e. a rather homogeneity in the scene and medium scale is expected. Your AOI is extremele heterogeneous and buildings around, but alos tiles and walkpaths are affecting the micrometeorology of the green patches. So I would not advice use this model for your purposes. In any case, I do not think the NaN are caused by this issue. So if you want me to dig deeper in what is causing the problem, please share your configuration data and and example (subset) of the inputs you are using

— Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you authored the thread.Message ID: @.***>

hectornieto commented 1 year ago

Sorry but I cannot find your configuration file nor your attached images.

But the effective resolution that the model should be run is at a spatial scale effectively larger than the main vegetation elements. It does not really make sense to estimate fluxes at 7cm resolution since this process happens at a much coarser scale due to radiation transmission efficiency, tubulent transport and plant hydraulics/phisiology. That is the main reason of the roles of fc and fg together with LAI.

All the best

HazelSCUT commented 1 year ago

Ok, but why have you also used drone images with high spatial resolution for your study? The difference is between homogeneous farmland and highly heterogeneous urban surfaces rather than the spatial resolution. Thanks.

------------------ 原始邮件 ------------------ 发件人: "hectornieto/pyTSEB" @.>; 发送时间: 2023年3月21日(星期二) 凌晨3:08 @.>; @.**@.>; 主题: Re: [hectornieto/pyTSEB] Can I use the ProcessLocalImage.ipynb to derive ET of urban scenes based on UAV images? (Issue #53)

Sorry but I cannot find your configuration file nor your attached images.

But the effective resolution that the model should be run is at a spatial scale effectively larger than the main vegetation elements. It does not really make sense to estimate fluxes at 7cm resolution since this process happens at a much coarser scale due to radiation transmission efficiency, tubulent transport and plant hydraulics/phisiology. That is the main reason of the roles of fc and fg together with LAI.

All the best

— Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you authored the thread.Message ID: @.***>

hectornieto commented 1 year ago

We use VHR thermal imagery from UAV to extract canopy/leaf and soil/background temperatures. In addtiion the 3D cloud point might be used to characterize canopy structure and surface roughness. But the model is applied at an effective resolution that accounts for a whole canopy/soil system.

diviningwater commented 1 year ago

We did conduct a study in urban areas with drones and UAV data using the TSEB model (manuscript is ready for submission to a journal). The major question for us was the spatial generation of LAI. LAI is not a variable that can be downscaled (to cm/pixel), and it is limited to the sampling approach (e,g. defining an area of turfgrass like 2 2 or 3 3 m) where to collect destructive LAI samples with repetition. Therefore any LAI sample area we chose is the final LAI spatial map resolution. As a rule of thumb, the smallest ET spatial estimation (pixel size) is dictated by the coarsest input resolution (in this case LAI).