surveyPoints: A function to generate georeferenced locations for monitoring...

In soilassessment: Assessment Models for Agriculture Soil Conditions and Crop Suitability

A function to generate georeferenced locations for monitoring soil conditions

Description

This function uses stratified random sampling to generate georeferenced locations for monitoring soil conditions

Usage

surveyPoints(soilmap,scorpan,conditionclass,mapproportion)

Arguments

soilmap	input classified map of soil condition
scorpan	number of scorpan factors that generated teh soil condition map. The range is 1-5
conditionclass	reference class in the soil condition map to be monitored. The class code should be in the map
mapproportion	Proportion in percent of the reference class in the soil condition map to be monitored.

Details

The number of scorpan factors can be assumed but need to be with respect to the soil forming factors.The maximum possible number of factors is 5 irrespective of number of layers in each factor while the minimum number is 1.The soil condition class is the class code in the map which is to be targeted

Value

A spatial points dataframe with projection similar to the soil condition map projection

Author(s)

Christian Thine Omuto

See Also

featureRep, imageIndices, pedoTransfer, classCode

Examples

library(sp)
library(raster)
ec=suitabinput["ec"]
ph=suitabinput["ph"]
soc=nutrindicator["soc"]
clay=textureinput["clay"]
texture=suitabinput["texture"]
newmap=ec
newmap$ph=ph$ph
newmap$ECe=ECconversion1(ec$ec*0.1,texture$texture,"FAO","1:5", soc$soc,clay$clay)
newmap$salt=saltSeverity(newmap$ECe,newmap$ph,0.84,"FAO")
newmap$salineclass=classCode(newmap$salt,"saltseverity")
newmap$salineclass1=as.factor(newmap$salineclass)
spplot(newmap["salineclass"], main="Salinity Code")
summary(newmap$salt)
summary(newmap$salineclass)
salt=raster(newmap["salt"])
salt1=newmap["salt"]
n_points=surveyPoints(salt1,4,11,80)
length(n_points$new)
spplot(salt1, scales=list(draw=TRUE),sp.layout=list("sp.points",n_points,pch=8,col="cyan"))
spplot(salt, scales=list(draw=TRUE),sp.layout=list("sp.points",n_points,pch=8,col="cyan"))

soilassessment documentation built on July 9, 2023, 7:40 p.m.