逐步解释% matlab mobile传感器数据，从手机导入到本地即可 clear;close all;clc; load('sensorlog_20230601_123433.mat'); lat = Position.latitude; lon = Position.longitude; altitude = Position.altitude; timestamp = Position.Timestamp; spd = Position.speed; nums = length(lat); nBins = 10; binSpacing = (max(spd) - min(spd))/nBins; binRanges = min(spd):binSpacing:max(spd)-binSpacing; % 添加下确界 binRanges(end+1) = inf; % |histc| 确定值落入哪一个bin [~, spdBins] = histc(spd, binRanges); lat = lat'; lon = lon'; spdBins = spdBins'; % 创建一个地理形状矢量，该矢量将线段存储为features s = geoshape(); for k = 1:nBins % 保留与当前bin匹配的经纬度，其余部分保留为NaN，NaN为线段中的中断 latValid = nan(1, length(lat)); latValid(spdBins==k) = lat(spdBins==k); lonValid = nan(1, length(lon)); lonValid(spdBins==k) = lon(spdBins==k); % 保留从当前速度bin转换到另一个速度bin后出现的经纬度使路径连续 transitions = [diff(spdBins) 0]; insertionInd = find(spdBins==k & transitions~=0) + 1; % 预分配空间并插入额外的经纬度 latSeg = zeros(1, length(latValid) + length(insertionInd)); latSeg(insertionInd + (0:length(insertionInd)-1)) = lat(insertionInd); latSeg(~latSeg) = latValid; lonSeg = zeros(1, length(lonValid) + length(insertionInd)); lonSeg(insertionInd + (0:length(insertionInd)-1)) = lon(insertionInd); lonSeg(~lonSeg) = lonValid; % 将经纬度线段添加到地理形状矢量 s(k) = geoshape(latSeg, lonSeg); end wm = webmap('World Imagery'); mwLat = 26.053376; mwLon = 119.187501; name = 'School'; iconDir = fullfile(matlabroot,'toolbox','matlab','icons'); iconFilename = fullfile(iconDir, 'fzu.png'); wmmarker(mwLat, mwLon, 'FeatureName', name, 'Icon', iconFilename); colors = autumn(nBins); wmline(s, 'Color', colors, 'Width', 5); wmzoom(16);

这段 MATLAB 代码的功能是将从手机传感器中获取的位置、速度等信息可视化在地图上，并将速度分成 nBins 个区间，对于每个区间将对应的位置信息连接成线段，最终将所有线段可视化在地图上。

下面是代码的逐行解释：

clear; close all; clc;
load('sensorlog_20230601_123433.mat');

清空工作区，关闭所有窗口，加载存储在 mat 文件中的数据。

lat = Position.latitude;
lon = Position.longitude;
altitude = Position.altitude;
timestamp = Position.Timestamp;
spd = Position.speed;
nums = length(lat);

从存储在 mat 文件中的 Position 变量中获取经度、纬度、高度、时间戳和速度等信息，并计算经度和纬度的数量。

nBins = 10;
binSpacing = (max(spd) - min(spd))/nBins;
binRanges = min(spd):binSpacing:max(spd)-binSpacing;
% 添加下确界
binRanges(end+1) = inf;
% |histc| 确定值落入哪一个bin
[~, spdBins] = histc(spd, binRanges);

将速度分成 nBins 个区间，并确定每个速度值所属的区间。

lat = lat'; lon = lon'; spdBins = spdBins';

转置经度、纬度和速度值，以便后面的操作。

% 创建一个地理形状矢量，该矢量将线段存储为features
s = geoshape();
for k = 1:nBins
    % 保留与当前bin匹配的经纬度，其余部分保留为NaN，NaN为线段中的中断
    latValid = nan(1, length(lat));
    latValid(spdBins==k) = lat(spdBins==k);
    
    lonValid = nan(1, length(lon));
    lonValid(spdBins==k) = lon(spdBins==k);
    
    % 保留从当前速度bin转换到另一个速度bin后出现的经纬度使路径连续
    transitions = [diff(spdBins) 0];
    insertionInd = find(spdBins==k & transitions~=0) + 1;
    
    % 预分配空间并插入额外的经纬度
    latSeg = zeros(1, length(latValid) + length(insertionInd));
    latSeg(insertionInd + (0:length(insertionInd)-1)) = lat(insertionInd);
    latSeg(~latSeg) = latValid;
    
    lonSeg = zeros(1, length(lonValid) + length(insertionInd));
    lonSeg(insertionInd + (0:length(insertionInd)-1)) = lon(insertionInd);
    lonSeg(~lonSeg) = lonValid;
    
    % 将经纬度线段添加到地理形状矢量
    s(k) = geoshape(latSeg, lonSeg);
end

将速度分成的每个区间对应的经纬度信息连接成线段，并将所有线段存储在地理形状矢量 s 中。

wm = webmap('World Imagery');
mwLat = 26.053376; mwLon = 119.187501;
name = 'School';
iconDir = fullfile(matlabroot,'toolbox','matlab','icons');
iconFilename = fullfile(iconDir, 'fzu.png');
wmmarker(mwLat, mwLon, 'FeatureName', name, 'Icon', iconFilename);

创建一个 Web 地图，设置地图底图为 World Imagery，添加一个标记点用于标记地图的中心点。

colors = autumn(nBins);
wmline(s, 'Color', colors, 'Width', 5);
wmzoom(16);

将线段可视化在地图上，并设置线段颜色为彩虹色，线段宽度为 5 像素。最后将地图缩放至级别 16。