personas - equipos gps portatiles

Llenar los valores de GPS faltantes en R

Question

Estoy mirando datos de un dispositivo que pasa gran parte del tiempo bajo el agua. Cuando sale a la superficie, obtiene una corrección de GPS (lat y lon) y luego se hunde (perdiendo su corrección de GPS) y continúa recopilando datos nuevamente hasta su próxima aparición. Esto da como resultado una gran cantidad de valores de NA en los datos para los valores latitudinales y no coinciden con las lecturas de datos subacuáticos.

Me gustaría crear un gráfico de la cortina de los datos, pero necesitaré algunos valores lat / lon interpolados para hacer una representación aproximada de los datos subsuperficiales en un mapa 3D.

¿Cómo puedo completar algunos valores linealmente interpolados para las muchas NA que ocurren entre las correcciones de GPS en R? Están en formato numérico DD.DDDDD.

Se puede encontrar un ejemplo de los datos en: http://modata.ceoe.udel.edu/public/gps_example_data.csv

Answer 1

na.approx() en el paquete del zoo hace justo lo que estás buscando.

Con sus datos, haga algo como:

df <- read.csv("http://modata.ceoe.udel.edu/public/gps_example_data.csv", 
               header=T)

library(zoo)

df2 <- 
transform(df, 
    m_gps_lat_dec = na.approx(m_gps_lat_dec, m_present_time, na.rm=FALSE),
    m_gps_lon_dec = na.approx(m_gps_lon_dec, m_present_time, na.rm=FALSE))

Consulte también esta publicación de , que incluye un buen ejemplo trabajado (y algunos comentarios útiles de un autor del paquete del zoo , para arrancar).