BSL: Bare Soil Line

Description Usage Arguments Details Value Author(s) References Examples

View source: R/BSL.R

Description

Finds Bare Soil Line (BSL) and maximum vegetation point.

Usage

1
BSL(band3, band4, method = "quantile", ulimit = 0.99, llimit = 0.005, maxval = 255)

Arguments

band3

File name or image file (matrix, data frame, or SpatialGridDataFrame) for Landsat band 3 DN (red).

band4

File name or image file (matrix, data frame, or SpatialGridDataFrame) for Landsat band 4 DN (NIR).

method

Either "quantile" or "minimum" – describes way in which soil line is identified.

ulimit

Upper limit for quantile of band ratios (ulimit < 1).

llimit

Lower limit for quantile of band ratios (llimit > 0).

maxval

Maximum value for band data; default of 255 for Landsat 5 and 7.

Details

Finding the BSL requires identifying the lowest NIR values for each level of red. The quantile method takes the lowest set of points, those with a NIR/red ratio less than the llimit-th quantile. The minimum value method takes the lowest NIR value for each level of red. However they are found, these points with low NIR for their red values are used in a major axis regression to find the Bare Soil Line. This function also identifies the full canopy point (maximum vegetation), by using the ulimit to identify the top points, with NIR/red ratio greater than the ulimit-th quantile, and with high NIR values. Red or NIR values of 255 (saturated sensor) are omitted when calculating the BSL.

Value

BSL

Regression coefficients for the Bare Soil Line

top

band 3 and band 4 values for the full canopy point

Author(s)

Sarah Goslee

References

Maas, S. J. & Rajan, N. 2010. Normalizing and converting image DC data using scatter plot matching. Remote Sensing 2:1644-1661.

Examples

1
2
3
4
5
6
	data(nov3)
	data(nov4)
	nov.bsl <- BSL(nov3, nov4)
	plot(as.vector(as.matrix(nov3)), as.vector(as.matrix(nov4)))
	abline(nov.bsl$BSL, col="red")
	points(nov.bsl$top[1], nov.bsl$top[2], col="green", cex=2, pch=16)

Example output

Loading required package: rgdal
Loading required package: sp
rgdal: version: 1.5-18, (SVN revision 1082)
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 3.0.4, released 2020/01/28
Path to GDAL shared files: /usr/share/gdal
GDAL binary built with GEOS: TRUE 
Loaded PROJ runtime: Rel. 6.3.1, February 10th, 2020, [PJ_VERSION: 631]
Path to PROJ shared files: /usr/share/proj
Linking to sp version:1.4-4
To mute warnings of possible GDAL/OSR exportToProj4() degradation,
use options("rgdal_show_exportToProj4_warnings"="none") before loading rgdal.

landsat documentation built on Dec. 19, 2019, 5:06 p.m.

Related to BSL in landsat...