GPS 数据（GPGGA）导入 ArcGIS Pro

1. 观察原始数据

提供的GPGGA数据，TXT格式，共15列，以逗号分隔：

根据查询的资料（百度百科-GPGGA）显示，字段含义如下：

$GPGGA 语句包括17个字段：语句标识头，世界时间，纬度，纬度半球，经度，经度半球，定位质量指示，使用卫星数量，水平精确度，海拔高度，高度单位，大地水准面高度，高度单位，差分GPS数据期限，差分参考基站标号，校验和结束标记(用回车符 和换行符 )，分别用14个逗号进行分隔。该数据帧的结构及各字段释义如下：

2. 数据处理

以原始数据的一行为例：
$GNGGA,000233.00,4050.24926,N,12026.77689,E,1,09,0.96,165.1,M,3.2,M,*7E

1. 在整个文档的第一行加入对应的字段名称：

语句标识头,世界时间,纬度,纬度半球,经度,经度半球,定位质量指示,使用卫星数量,水平精确度,海拔高度,高度单位,大地水准面高度,高度单位,差分GPS数据期限,差分参考基站标号

2. 为方便进一步处理，建议将文档另存为csv格式，UTF-8编码
3. 使用 Excel 打开 csv，这一步是因为 excel 比较好对数据做统一的简单计算（比如某一列都除以100之类的）
   

2.1 经纬度处理

C和E分别为经纬度，在现有数据基础上除以100：

2.2 时间戳处理

把B列，UTC时间，也除以100，变成小数 2.33。
hhmmss.sss 小时、分钟、秒 这个格式，000233.00 等同于 2分23秒，出于后期编程的需要，老板要求转换为：hh:mm:ss格式，即 2.33 转化为 00:02:23 这个格式转换。我写了个python代码解决这一问题，代码如下：

"""
Created on Tue Oct  8 16:48:50 2019
读取csv，并将其中第二列转化为时间格式 hh:mm:ss
@author: yao
"""

import csv
import os

odata = 'test.csv' #原始数据
data = [] 

if not os.path.exists(odata):
    print('wrong path!')
    
with open(odata) as csvfile:
    csv_reader = csv.reader(csvfile)
    header = next(csv_reader) #读取表头
    for row in csv_reader:
        mm = row[1].split('.')[0] #分钟
        ss = row[1].split('.')[1] #秒
        if len(mm)<2:
            mm = '0' + mm 
        timestamp = '00:' + mm + ':' + ss #格式化为 hh:mm:ss
        row[1] = timestamp
        data.append(row)
        
data.insert(0,header) #重新加入表头

#写入新csv文档
with open('test2.csv','w') as csvfile:
    writer = csv.writer(csvfile)    
    for row in data:
        writer.writerow(row)
        

写好的新 csv 差不多长这样：

数据处理至此告一段落。

PS. 其实前面的经纬度操作也可以用 python 代码完成，但那时老板还没提出改时间戳，所以全程都用 excel 操作的。
PPS. excel 小能手应该可以全程用 excel 完成，完成的话，务必跟我讲一下怎么做的，感谢！

3. GPS 数据导入 ArcGIS Pro

将 csv 直接拖入 Pro 的左边栏 Contents 部分，右键导入的表格，选择 Display XY data，随后指定 X、Y、Z所在的列名，点击运行。

检查运行结果的属性表，没有问题：

4. ArcGIS Pro 中的要素类导出成 shp

在右边栏找到 Catalog，在gdb库里找到找到需要导出成shp格式的文件，右键点击选择导出为 shape file：

导出的 shp 可以发布为 portal 上的服务，或者传给 whoever 需要的人……