我已经使用评估程序库存储了两个栅格图像(.tif格式),并绘制了这些栅格图像之间的散点图。当我使用“ abline”功能创建回归线时,该线为我提供了斜率和截距。但是我想在x轴上为最小的点生成另外两条倾斜线,为最大的点生成另一条倾斜线,但是我不能。
帮我这样做。
这是我使用R创建的散点图。
我想用两个倾斜线和方程式创建这种类型的图。
lst<- "F:/LC08_L1TP_148043_20180423_20180502_01_T1/Output/ArcGIS/LST.tif"
lst1 <- raster(lst)
plot(lst1)
summary(lst1)
ndvi<- "F:/LC08_L1TP_148043_20180423_20180502_01_T1/Output/ArcGIS/NDVI.tif"
ndvi1 <- raster(ndvi)
plot(ndvi1)
summary(ndvi1)
plot(ndvi1, lst1)
s <- stack(ndvi1, lst1)
names(s) <- c('ndvi1', 'lst1')
#xyplot(ndvi1~lst1, data = s, alpha = 1)
s1 <- data.frame(na.omit(values(s)))
#s2 <- lm(lst1~ndvi1, data = s1) #may give error
s1_1 <- sampleRegular(s, 1000000)
s1_2 <- data.frame(na.omit(s1_1))
s1_2 <- lm(lst1~ndvi1, data = s1_2)
xyplot(lst1~ndvi1, data = s, alpha = 1)
abline(s1_2, col = "red")
#equation of trendline
paste('y =', coef(s1_2)[[2]], '* x', '+', coef(s1_2)[[1]])