MATLAB绘制星下点轨迹消除经度±180°跳变的方法

以下我简单地给出解决思路，拒绝伸手党。（以后我有时间了会把整个代码发上来）

主要用到的函数是plot(longitude,latitude)，其中long为经度，在±180°之间；lat为纬度，在±90°之间。
我们所需要的星下点是连续的曲线，即使航天器星下点穿过了国际日界线，即±180°经线。如下图

但是可惜的是，如果简单的用下面这句函数，就无法画出上面的曲线

plot(Long*rad2deg, lat*rad2deg,'b-');
axis([-180,180,-90,90]);grid on;

可以发现，穿过日界线时，地图左右两端的点被直线连接起来了。解决方法是，判断跳变发生的时刻，分段画连续曲线。代码如下

la=length(Long);
diffLong=abs(diff(Long))>(npasses*180/la);
diffLong=find(diffLong);
if isempty(diffLong)
    plot(Long*rad2deg, lat*rad2deg,'b-'); 
else
    ka=1;
    totalSections=length(diffLong);
    for i=1:totalSections
        kb=diffLong(i);
        sec=ka:kb;
        ka=kb+1;
        plot(Long(sec)*rad2deg, lat(sec)*rad2deg,'b-'); 
    end
    sec=kb+1:la;
    plot(Long(sec)*rad2deg, lat(sec)*rad2deg,'b-'); 
end

hold off;
axis([-180,180,-90,90]);grid on;

其中diffLong 对经度数据取了差分；npasses是轨道圈数，这个变量可以大概估计一下；每次截取绘制的是向量Long的一部分sec=ka:kb；

如果用MATLAB的Mapping Toolbox，可以加载地球图形，就更好看了。也就是在星下点的原图底部，先有这样两句

hold on;
geoshow('landareas.shp','FaceColor','green'); % Plot map