地理折线图（`geoplot`）

地理折线图（geoplot）

06_极坐标与特殊　5 张　MATLAB R2023b　100 mm 宽 @ 300 dpi

这个方法适合什么

在地图底图上画路径，坐标用经纬度。适合：输电线路走向、巡检航线、抢修路径、线路沿程的电气量分布。

电气场景

1　500kV输电线路走向

2　无人机巡检航线

3　10kV馈线路径

4　海底电缆登陆段走向

5　抢修车辆行驶路径

代码

脚本 p6_polar_special.m 中对应 geoplot 的一段，SEL 用于选择只跑某几个方法。

% ---- 9-1 500 kV 输电线路走向 ----
try
    la1 = [22.62 22.66 22.70 22.74 22.79 22.84 22.88 22.91 22.95];
    lo1 = [113.18 113.27 113.36 113.46 113.55 113.66 113.76 113.85 113.93];
    la2 = la1 + 0.052; lo2 = lo1 - 0.008;
    f = ml_fig(PP,100,90); gx = axg(f);
    geoplot(gx,la1,lo1,'-o','LineWidth',1.2,'MarkerSize',3.2,'Color',C(1,:));
    geoplot(gx,la2,lo2,'-s','LineWidth',1.2,'MarkerSize',3.0,'Color',C(2,:));
    geolimits(gx,[22.55 23.06],[113.06 114.05]);
    lgp(gx,PP,{'甲乙Ⅰ回 216 km','甲乙Ⅱ回 218 km'},'northeast');
    save_png(f,FIGD,'500kV输电线路走向');
catch ME, fprintf('[FAIL] 9-1 %s\n',ME.message); close all; end

% ---- 9-2 无人机巡检航线 ----
try
    la = [22.62 22.66 22.70 22.74 22.79 22.84 22.88 22.91 22.95];
    lo = [113.18 113.27 113.36 113.46 113.55 113.66 113.76 113.85 113.93];
    laOut = la - 0.020; laBack = la - 0.055;          % 往返两侧各偏出线路 2~6 km
    laRoute = [laOut laBack]; loRoute = [lo fliplr(lo)];
    f = ml_fig(PP,100,90); gx = axg(f);
    geoplot(gx,la,lo,'--','LineWidth',1.0,'Color',C(5,:));
    geoplot(gx,laRoute,loRoute,'-o','LineWidth',1.2,'MarkerSize',3.0,'Color',C(2,:));
    geoscatter(gx,22.62,113.18,50,C(1,:),'filled');
    geolimits(gx,[22.53 23.02],[113.10 114.02]);
    lgp(gx,PP,{'线路中心线','巡检航迹','巡检起点'},'northeast');
    save_png(f,FIGD,'无人机巡检航线');
catch ME, fprintf('[FAIL] 9-2 %s\n',ME.message); close all; end

% ---- 9-3 10 kV 馈线路径 ----
try
    laM = [22.66 22.68 22.70 22.71 22.73 22.74 22.76 22.77];
    loM = [113.24 113.28 113.31 113.35 113.38 113.42 113.45 113.49];
    laB1 = [22.70 22.72 22.73]; loB1 = [113.31 113.30 113.26];
    laB2 = [22.73 22.75 22.76]; loB2 = [113.38 113.40 113.43];
    f = ml_fig(PP,100,90); gx = axg(f);
    geoplot(gx,laM,loM,'-','LineWidth',1.4,'Color',C(1,:));
    geoplot(gx,laB1,loB1,'-','LineWidth',1.2,'Color',C(4,:));
    geoplot(gx,laB2,loB2,'-','LineWidth',1.2,'Color',C(3,:));
    geoscatter(gx,22.66,113.24,60,C(2,:),'filled');
    geolimits(gx,[22.64 22.79],[113.20 113.52]);
    lgp(gx,PP,{'10 kV 主干线','支线 1','支线 2','变电站'},'northeast');
    save_png(f,FIGD,'10kV馈线路径');
catch ME, fprintf('[FAIL] 9-3 %s\n',ME.message); close all; end

% ---- 9-4 海底电缆登陆段走向 ----
try
    laSea = [22.36 22.40 22.44 22.48 22.52];
    loSea = [113.97 113.94 113.91 113.88 113.85];
    laLnd = [22.52 22.56 22.59 22.62 22.65];
    loLnd = [113.85 113.82 113.78 113.74 113.70];
    f = ml_fig(PP,100,90); gx = axg(f);
    geoplot(gx,laSea,loSea,'-o','LineWidth',1.4,'MarkerSize',3.2,'Color',C(1,:));
    geoplot(gx,laLnd,loLnd,'-s','LineWidth',1.4,'MarkerSize',3.2,'Color',C(3,:));
    geoscatter(gx,22.52,113.85,60,C(2,:),'filled');
    geolimits(gx,[22.33 22.68],[113.66 114.01]);
    lgp(gx,PP,{'海缆段 12.4 km','陆缆段 6.2 km','登陆点'},'northwest');
    save_png(f,FIGD,'海底电缆登陆段走向');
catch ME, fprintf('[FAIL] 9-4 %s\n',ME.message); close all; end

% ---- 9-5 抢修车辆行驶路径 ----
try
    la = [22.64 22.66 22.69 22.72 22.74 22.77 22.79];
    lo = [113.21 113.25 113.29 113.33 113.38 113.43 113.48];
    laAlt = [22.64 22.67 22.70 22.73 22.75 22.78 22.79];
    loAlt = [113.21 113.26 113.31 113.36 113.41 113.45 113.48];
    f = ml_fig(PP,100,90); gx = axg(f);
    geoplot(gx,la,lo,'-','LineWidth',1.4,'Color',C(2,:));
    geoplot(gx,laAlt,loAlt,'--','LineWidth',1.2,'Color',C(6,:));
    geoscatter(gx,22.64,113.21,55,C(1,:),'filled');
    geoscatter(gx,22.79,113.48,55,C(4,:),'filled');
    geolimits(gx,[22.62 22.82],[113.18 113.52]);
    lgp(gx,PP,{'最优路径 42 min','备选路径 55 min','抢修基地','故障点'},'northeast');
    save_png(f,FIGD,'抢修车辆行驶路径');
catch ME, fprintf('[FAIL] 9-5 %s\n',ME.message); close all; end