hotTs2: Automatic computation of hot pixel Radiometric Surface...

In gowusu/sebkc: The Surface Energy Balance and Crop Coefficient Estimation with R

Automatic computation of hot pixel Radiometric Surface Temperature (Ts)

Description

This function computes temperature of a hot pixel based on NDVI, Surface Temperature, DEM and albedo. It plots histograms indicating the selected values

Usage

hotTs2(Ts, NDVI, albedo = NULL, sunangle = NULL, DEM = NULL,
  extent = "interactive", upper = 0.95, lower = 0.8, NDVI.probs = 0.1,
  plot = TRUE, layout = "portrait", draw = "poly", folder = NULL,
  welev = NULL)

## Default S3 method:
hotTs2(Ts = NULL, NDVI, albedo = NULL, sunangle = NULL,
  DEM = NULL, cluster = 10, extent = "interactive", upper = 0.95,
  lower = 0.8, NDVI.probs = 0.2, plot = TRUE, layout = "portrait",
  draw = "poly", folder = NULL, welev = NULL)

Arguments

Ts	A RasterLayer data that indicates radiometric surface temperature values from remote sensing image, preferably in Kelvin (K). You can also point to raster file on your computer.
NDVI	A RasterLayer data that indicates NDVI (Normalized Difference Vegetation Index) values. You can also point to a file on your computer.
albedo	A RasterLayer data that has albedo values. You can also provide file path or location on your computer
sunangle	sun angle (degree) above the horizon.
DEM	A digital elevation model[m]
extent	A subset of the raster file than include c(xmin, xmax; serastercond row: ymin, ymax). You can set it to "interactive" where you can digitise polygon. You can set to "auto" where it will be automatically generated. See layout parameter below for more details. see extent
upper	The quantile probability at which Ts should be included in the selection project. It ranges from 0-1. The default is 0.8
lower	The quantile probability at which NDVI should be included in the slection project.The default is 0.2
NDVI.probs	The probability of selecting quantile NDVI that is above the threshold. The value of 0.6 means only NDVI values that is above 60 the analyses.
plot	logical:TRUE or FALSE. To indicate if histograms or triangle should be plotted
layout	character. It takes "portrait" or "landscape". If extent is set to "auto" the middle part of the image will be used so that NAs can be avoided. You indicate whether the automatic selection should indicate west-east (landscape) or north-south direction (portrait)
draw	interactive. The method of digitisation if extent is set "auto". It takes "poly" or "rect"
folder	An original directory of the images that contains all landsat 5,7 or 8 bands and metadata. At the moment only Landsat folder is supported
welev	Weather station elevation in meters
cluster	numeric, a value indicating the number of clusters.

Details

The function first divides the selected area into a number of clusters. The mean and standard deviations are computed. The cluster means with the lowest standard deviations and qualifying for upper (Ts) and lower (NDVI) limits (see above) are selected. The provided upper and lower probabilities are also used to select the final Ts.

Value

• Ts: The selected or input Ts

• Tshot: The Temperature of hot pixel

• TX: The average temperature of the hot pixels

• x: The x cordinate of the hot pixel

• y: The y cordinate of the hot pixel

• xyhot: combined x and y coordinates

• candidates: The similar candidates' pixels that can be used

Author(s)

George Owusu

References

Allen, R. G., Burnett, B., Kramber, W., Huntington, J., Kjaersgaard, J., Kilic, A., Kelly, C., & Trezza, R. 2013. Automated Calibration of the METRIC-Landsat Evapotranspiration Process. JAWRA Journal of the American Water Resources Association, 49(3): 563-576.

See Also

coldTs, hotTs, coldT2

Examples

## Not run: 
 folder=system.file("extdata","stack",package="sebkc")
  modhot=hotTs2(folder=folder,welev=170,extent="auto",cluster=3)
  
   #use object of  landsat578
folder=system.file("extdata","stack",package="sebkc")
data=landsat578(data=folder, welev=170)
modfhot=hotTs2(folder=data,welev=170,extent="auto",cluster=3)

#Or Ts = data
modhotTs2=hotTs2(data,welev=170,extent="auto",cluster=3)
albedo=raster(system.file("extdata","albedo.grd",package="sebkc"))
Ts=raster(system.file("extdata","Ts.grd",package="sebkc"))
NDVI=raster(system.file("extdata","NDVI.grd",package="sebkc"))
#simulate  with default parameters without albedo
modhot=hotTs2(Ts=Ts,NDVI=NDVI,extent="auto",cluster=3)
#simulate with albedo but set extent to "auto" and cluster to 7
modhot=hotTs2(Ts=Ts,NDVI=NDVI,albedo=albedo,cluster=7,
extent="auto")

## End(Not run)

gowusu/sebkc documentation built on July 28, 2023, 11:44 a.m.