Я сохранил два растровых изображения (в формате .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]])