Main Content

Cette page a été traduite par traduction automatique. Cliquez ici pour voir la dernière version en anglais.

Visualisez le chemin parcouru dans les cartes de la NASA

Cet exemple utilise :

Ouvrir le script

Cet exemple montre comment visualiser les données de latitude et de longitude stockées dans un flux du canal ThingSpeak ™ sur une carte raster fournie par la NASA. Vous pouvez récupérer une carte raster à partir d'un service de cartographie Web (WMS) et tracer le chemin parcouru par le navire Mary Maersk.

Lire les données

Le canal ThingSpeak 73734 stocke les données de localisation du navire Mary Maersk. Lisez les données des 10 derniers jours à l'aide de la fonction thingSpeakRead.

data = thingSpeakRead(73734,'DateRange',...
    [datetime('now')-caldays(10),datetime('now')],'outputFormat','table');

Traiter les données brutes

Traitez les données brutes de latitude et de longitude avant de les tracer sur la carte. Stockez les informations de latitude et de longitude dans des variables distinctes lat et lon. Utilisez la fonction geoquadline pour trouver les limites du quadrilatère géographique le plus étroit possible qui délimite une ligne reliant les sommets avec des coordonnées géographiques spécifiées par lat et lon. Pour rendre la carte plus lisible, augmentez les limites de latitude et de longitude de 10 unités.

lon = data.Longitude;
lat = data.Latitude;
[latlim,lonlim] = geoquadline(lat,lon);
buf = 10;
[latlim,lonlim] = bufgeoquad(latlim,lonlim,buf,buf);

Récupérer la carte raster de la NASA

Utilisez le service de cartographie Web de la NASA pour produire des cartes de données raster spatialement référencées. La fonction wmsfind recherche dans les entrées du serverurl de la base de données WMS la chaîne de caractère 'nasa'.

nasa = wmsfind('nasa','SearchField','serverurl');

Affinez la recherche de la couche de carte raster requise dans la base de données de la NASA à l'aide de la fonction refine.

layer = refine(nasa,'bluemarbleng','SearchField','layername', ...
    'MatchType','exact');

Appelez la fonction wmsread pour lire la carte raster .

[A,R] = wmsread(layer(1),'Latlim',latlim,'Lonlim',lonlim);

Générer la carte

Utilisez la carte raster renvoyée par la fonction worldmap pour générer un tracé du chemin parcouru par le Mary Maersk au cours des 10 derniers jours.

geoshow(A,R);
geoshow(lat,lon);
xlabel('Longitude');
ylabel('Latitude');
axis tight;

Voir aussi

Fonctions

Sites web externes