Generación de campos aleatorios no correlacionados espacialmente en R

Question

Tengo un MDE (tamaño de celda 10x10 y cells.dim 450x300) y necesito crear 100 cuadrículas de error aleatorias que tengan las mismas dimensiones que el MDE.

Los valores de las cuadrículas de error se generan a partir de una distribución normal que tiene una media de cero y una desviación estándar igual al RMSE del MDE original ~ (0, +-3).

Luego tengo que añadir esas rejillas de error al DEM original y crear 100 realizaciones del MDE original.

Answer 1

Puedes hacer la primera parte de esto en R con bastante facilidad:

 errmats <- replicate(100,
     matrix(rnorm(450*300,mean=0,sd=3),nrow=450),simplify=FALSE)

da una lista de 100 matrices de 450 por 300 con errores aleatorios no correlacionados según lo especificado. (Si se omite el simplify=FALSE En cambio, obtendrá los resultados almacenados como una matriz de 450 por 300 por 100). Como no nos has dicho cómo están almacenados/estructurados tus MDE (modelos digitales de elevación) dentro de R, es un poco difícil decir exactamente cómo añadirlos a tu MDE original.

Un comentario anterior sugiere que tu modelo de error es poco realista: no puedo decirlo, ya que no sé lo que estás tratando de hacer (no hay suficiente contexto). Sin embargo, si decide utilizar estructuras de error más realistas (correlacionadas espacialmente), puede generarlas fácilmente con el programa RandomFields paquete.