[R-es] Interpolación IDW

[Javier Marcuzzi]

Estimado David

Miro su correo, sin poder reproducir como dice nuestro compañero de lista, es complicado responder, pero por lo que es su correo, hay dos puntos, el primero tiene que ver con la posibilidad que usted esté intentando trabajar con.dos cosas que piensa que son iguales pero son distintas, por eso el mensaje de error, cuándo define algo lo hace de tal forma que se crean dos cosas distintas, algo que a mí me ayuda es usar str( lo que quiero ver), el segundo punto es que es una tesis, su director de tesis debería resolver el problema.

Javier Marcuzzi

> El 21 feb. 2023, a las 15:16, David Camilo Gomez Medina <dcgomezme using unal.edu.co> escribió:
> 
> Emilio, muchas gracias. Tuve que rehacer el código para entender mejor
> todo, pero me sale el siguiente error:
> 
> 
> *Error in as.data.frame.default(x[[i]], optional = TRUE) :   cannot coerce
> class ‘structure("RasterLayer", package = "raster")’ to a data.frame*
> 
> El código es el siguiente:
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> *for (i in 1:nrows_mt) {    data_it_mt <-
> as.numeric(as.vector(data_y_mt[i,]))    it_data_mt <- idw(data_it_mt,
> sta_mt, grid, p = 1)    grid <- data.frame(grid, valores = it_data_mt[ ,
> 1])    coordinates(grid) <- ~x+y    gridded(grid) <- T    grid <-
> raster(grid, "valores")    projection(grid) <- crs("+init=epsg:9377")
> path <- 'C:/Users/David Gomez/Desktop/TESIS/DATOS/RESULTS/'    numeration
> <- as.character(i)    name <- "MEAN_TEMP"    exten <- ".tif"
> writeRaster(grid, paste(path, name, numeration, exten), overwrite = T)}*
> 
> No entiendo muy bien a qué parte hace énfasis el error.
> 
> Quedo atento, gracias.
> 
> Saludos.
> 
> On Mon, 20 Feb 2023 at 02:11, Emilio L. Cano <emilopezcano using gmail.com <mailto:emilopezcano using gmail.com>> wrote:
> 
>> Hola,
>> Es difícil responder con un ejemplo no reproducible y sin el mensaje de
>> error, pero te diría que intentes primero comprobar que puedes ajustar un
>> único modelo fuera del bucle.
>> 
>> Dicho esto, parece que estás pasando como primer argumento de la función
>> un “string”, el nombre de la variable, y debe ser una fórmula. Tal vez
>> as.formula(current_var, “~1”) te funcione.
>> 
>> En cuanto a los NA, puedes usar en cada iteración current_df <-
>> tidyr::drop_na(current_var) y usar current_df en vez de sta_mt
>> 
>> Saludos,
>> Emilio
>> 
>> 
>> 
>> El 17 feb 2023, a las 19:56, David Camilo Gomez Medina <
>> dcgomezme using unal.edu.co> escribió:
>> 
>> Hola, buen día a todos.
>> 
>> He estado mirando cómo interpolar unos datos anuales con el método IDW,
>> pero no consigo aún realizarlo.
>> 
>> idw_models <- list()
>> 
>> for (i in 1:ncol(data_y_mt)) {
>>  #Extract the current variable
>>  current_var <- names(data_y_mt)[i]
>> 
>>  #Create an IDW model for the current variable
>>  idw_model <- idw(current_var, sta_mt, grid, p = 1)
>> 
>>  idw_models[[i]] <- idw_model
>> }
>> 
>> Ese es el código que llevo, pero me siento muy perdido. La variable*
>> data_y_mt* es el dataframe donde tengo los valores, *sta_mt* son las
>> coordenadas de los puntos, *grid* es la malla coordenada.
>> 
>> Me gustaría también eliminar los datos con valores NA. Adjunto imagen del
>> dataframe.
>> 
>> Quedo muy atento, gracias.
>> 
>> <image.png>
>> 
>> *Aviso legal:* El contenido de este mensaje y los archivos adjuntos son
>> confidenciales y de uso exclusivo de la Universidad Nacional de Colombia.
>> Se encuentran dirigidos sólo para el uso del destinatario al cual van
>> enviados. La reproducción, lectura y/o copia se encuentran prohibidas a
>> cualquier persona diferente a este y puede ser ilegal. Si usted lo ha
>> recibido por error, infórmenos y elimínelo de su correo. Los Datos
>> Personales serán tratados conforme a la Ley 1581 de 2012 y a nuestra
>> Política de Datos Personales que podrá consultar en la página web
>> www.unal.edu.co <http://www.unal.edu.co/>. Las opiniones, informaciones, conclusiones y cualquier
>> otro tipo de dato contenido en este correo electrónico, no relacionados con
>> la actividad de la Universidad Nacional de Colombia, se entenderá como
>> personales y de ninguna manera son avaladas por la Universidad.
>> _______________________________________________
>> R-help-es mailing list
>> R-help-es using r-project.org <mailto:R-help-es using r-project.org>
>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>> 
>> 
>> 
> 
> -- 
> *Aviso legal:* El contenido de este mensaje y los archivos adjuntos son 
> confidenciales y de uso exclusivo de la Universidad Nacional de Colombia. 
> Se encuentran dirigidos sólo para el uso del destinatario al cual van 
> enviados. La reproducción, lectura y/o copia se encuentran prohibidas a 
> cualquier persona diferente a este y puede ser ilegal. Si usted lo ha 
> recibido por error, infórmenos y elimínelo de su correo. Los Datos 
> Personales serán tratados conforme a la Ley 1581 de 2012 y a nuestra 
> Política de Datos Personales que podrá consultar en la página web 
> www.unal.edu.co <http://www.unal.edu.co/> <http://www.unal.edu.co/>.* *Las opiniones, informaciones, 
> conclusiones y cualquier otro tipo de dato contenido en este correo 
> electrónico, no relacionados con la actividad de la Universidad Nacional de 
> Colombia, se entenderá como personales y de ninguna manera son avaladas por 
> la Universidad.
> 
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> R-help-es mailing list
> R-help-es using r-project.org <mailto:R-help-es using r-project.org>
> https://stat.ethz.ch/mailman/listinfo/r-help-es


	[[alternative HTML version deleted]]