Objectif

Corrigez les données GPS à l'aide de Map Matching API de Mapbox J'ai essayé de créer un outil.

Parmi les API Mapbox, nous utilisons l'API Map Matching, qui est classée comme service de navigation.
L'API Map Matching est une API qui capture les données GPS (latitude, longitude) contenant des erreurs sur le réseau routier et les corrige à une position probable.

Environnement d'exécution

Cet article utilise Python3, l'API Map Matching de Mapbox.
Les jetons d'accès sont requis pour utiliser l'API Map Matching de Mapbox.
Les jetons d'accès seront émis automatiquement lorsque vous enregistrerez un compte sur Site Mapbox.
Vous pouvez utiliser ces jetons d'accès pour appeler diverses API. --Mapbox est pay-as-you-go, veuillez donc consulter la page de tarification de Mapbox pour connaître les tarifs détaillés. De plus, étant donné que le niveau gratuit est défini pour jusqu'à 50000 demandes par mois, vous ne serez pas facturé sauf si vous faites un grand nombre de demandes.

Conditions préalables

Des données d'utilisation

――Pour les données d'utilisation, les données GPS acquises par la montre Garmin portée au poignet pendant le tournoi marathon sont utilisées.

Veuillez préparer les données GPS le cas échéant.

Méthode d'exécution

Paramètres d'entrée

--Entrez les données GPS (latitude, longitude) dans le fichier CSV d'entrée.

Courir

Entrez le fichier CSV d'entrée, le fichier CSV de sortie et api_token (jetons d'accès).
Ne modifiez pas l '"URL de l'API de correspondance de carte".
Exécutez le script Python.

`MapboxMapMatchingAPI.py`


import os
import csv
import json
import datetime
import time
import configparser
import requests

#Spécification du fichier CSV d'entrée
input_csvfile = "Entrée Veuillez saisir le fichier CSV."

#Spécification du fichier CSV de sortie
output_csvfile = "Entrez le fichier CSV de sortie."

#URL de l'API de correspondance de la carte * Ne pas modifier.
api_url = "https://api.mapbox.com/matching/v5/mapbox/driving/"

# Mapbox Access token
api_token = "Veuillez saisir le jeton d'accès."

#Liste de travail
data = []
wklst = []

#Définissez le fichier CSV d'entrée sur la liste de travail
with open(input_csvfile) as f:
    #Ignorer la ligne d'en-tête
    h = next(csv.reader(f))
    for v in csv.reader(f):
        data.append(v)

for i in range(len(data) -1):
    wklst.append([data[i][0],data[i][1],data[i+1][0],data[i+1][1]])

# print(wklst)

#Vérifiez l'existence du fichier CSV de sortie(Jugement de la présence ou de l'absence d'un en-tête)
if(os.path.exists(output_csvfile)):
    header_flg = 0
else:
    header_flg = 1

#Vérifiez l'existence du fichier CSV de sortie(Jugement de la présence ou de l'absence d'un en-tête)
if(os.path.exists(output_csvfile)):
    header_flg = 0
else:
    header_flg = 1

#Demander à mapper l'API de correspondance à partir de la liste de travail et écrire le résultat dans le fichier CSV de sortie
with open(output_csvfile, 'a', encoding='shift_jis') as output_csvfile_fp:

    fieldnames = ['bef_lat','bef_lng','aft_lat','aft_lng']
    csvfile_writer = csv.DictWriter(output_csvfile_fp, fieldnames=fieldnames,lineterminator='\n')
    
    #Sortie de l'en-tête de sortie lors de la création d'un nouveau fichier CSV
    if header_flg == 1:
        csvfile_writer.writeheader()
        header_flg == 0

    print(len(wklst))
    for i in range(len(wklst)):
        # print(wklst[i])
        #Génération de paramètres d'API Map Matching
        origins = wklst[i][1] + ','+ wklst[i][0]
        destinations = wklst[i][3] + ','+ wklst[i][2]
        params = {
            'radiuses':'50;50',
            'access_token': api_token
            }
        #Demande de mappage de l'API de correspondance
        raw_response = requests.get(api_url + origins + ';' + destinations, params)
        parsed_response = json.loads(raw_response.text)
        parsed_response_json = json.dumps(parsed_response, indent=4)
        # print(parsed_response_json)
        breakpoint

        #Extraire des éléments de JSON
        if parsed_response['code'] == 'Ok':
            lng1 = parsed_response['tracepoints'][0]['location'][0]
            lat1 = parsed_response['tracepoints'][0]['location'][1]
            lng2 = parsed_response['tracepoints'][1]['location'][0]
            lat2 = parsed_response['tracepoints'][1]['location'][1]
        else:
            lng1 = None
            lat1 = None
            lng2 = None
            lat2 = None

        #Écrire dans le fichier CSV de sortie
        csvfile_writer.writerow({
            'bef_lat': wklst[i][0],
            'bef_lng': wklst[i][1],
            'aft_lat': lat1,
            'aft_lng': lng1
            })
        
        print(wklst[i][0],wklst[i][1],lat1,lng1)

print(u'Sortie de traitement')

Résultat d'exécution

Les données GPS (latitude, longitude) avant correction se trouvent dans les première et deuxième colonnes.
Les données GPS corrigées (latitude, longitude) seront dans les 3ème et 4ème colonnes.

――Le résultat de l'exécution est visualisé par QGIS comme suit.

Les points bleus sont les données GPS avant correction (latitude, longitude), et les points rouges sont les données GPS après correction (latitude, longitude). ――Le point bleu qui était hors de la route avant la correction est accroché sur la route après la correction, et il est corrigé à une certaine position (point rouge).

Pour la carte de fond, la tuile Institut géographique (carte de couleur claire) du National Land Research Institute est utilisée.

Points à noter

** * Veuillez noter que vous serez facturé si vous utilisez l'API au-delà du niveau gratuit (> _ <) **

[PYTHON] J'ai créé un outil pour corriger les données GPS avec l'API Map Matching de Mapbox (Mapbox Map Matching API)