Displaying 1 result from an estimated 1 matches for "mu_x0".
Did you mean:
mu_x
2023 Nov 08
1
Problem in R code
...sults <- data.frame() #A data frame to store the results
for(i in 1:8)
{
LST_city <- extract(LST, c(lon[i],lat[i]), fun = mean, buffer = 10000,
na.rm = TRUE) #error
}
# Fit a 2D Gaussian surface to the LST data
x <- coordinates(LST)[,1]
y <- coordinates(LST)[,2]
z <- values(LST)
mu_x0 <- mean(x)
mu_y0 <- mean(y)
sigma_x0 <- sd(x)/2
sigma_y0 <- sd(y)/2
amp0 <- max(z)
opt <- optim(c(mu_x0, mu_y0, sigma_x0, sigma_y0, amp0), sse) #error 2
#Calculate the footprint of SUHI effect (FP) by the Gaussian surface
FP <- which(gaussian2d(x, y, opt$par, opt$par, opt$...