如何移除缺失经纬度的GPS观测行并完成数据预处理？

以下是整合所有需求的完整代码，用管道操作串联所有步骤：

library(lubridate)
library(amt)
library(dplyr) # 管道操作依赖dplyr，需加载

# 读取数据并完成所有处理步骤
data_clean <- read.csv("gpslocations.csv") %>%
  # 移除LATITUDE或LONGITUDE为NA的行
  filter(!is.na(LATITUDE) & !is.na(LONGITUDE)) %>%
  # 合并DATE和TIME为完整的datetime列
  mutate(datetime = ymd_hms(paste(DATE, TIME))) %>%
  # 转换为amt的track对象（GPS原始数据默认是WGS84坐标系，EPSG:4326）
  make_track(.x = LONGITUDE, .y = LATITUDE, .t = datetime, crs = 4326) %>%
  # 转换到目标坐标系（示例为UTM 10N，EPSG:32610，可根据研究区域替换）
  transform_coords(crs_to = 32610)

关键问题修正与步骤说明： #

• 移除缺失坐标行：
  你之前的代码!is.na('LATITUDE')错误在于，它检查的是字符串"LATITUDE"而非数据框中的列。正确做法是用filter()函数，同时筛选出LATITUDE和LONGITUDE都不为NA的行，确保彻底移除无效数据。

• 日期时间标准化：
  用paste()合并DATE和TIME列，再通过ymd_hms()（lubridate工具）将字符串转换为标准POSIXct日期时间格式，方便后续时间序列分析。

• 坐标系转换：

  1. 先用make_track()将数据转为amt包的track对象，明确指定原始GPS数据的坐标系为WGS84（EPSG:4326）；
  2. 再用transform_coords()转换到目标坐标系，示例中用的是UTM 10N（EPSG:32610），可根据你的研究区域替换为对应EPSG编码。

结果验证： #

运行代码后，data_clean将仅保留前4行有效GPS数据，同时包含标准化的日期时间字段和转换后的坐标系数据。

内容的提问来源于stack exchange，提问作者EleMan