function [La0,Lo0,La1,Lo1]=get_circumference3(lat,lon,yaw,altitude,range,rota)
    % get the circumference 
    La0=[];Lo0=[];
    La1=[];Lo1=[];
    for idx2=1:length(lat)
        if ~ismember(idx2,rota)
            [la0,lo0,d0]=getGeo(0.0,0.0,range,lat(idx2),lon(idx2),altitude(idx2),yaw(idx2),0.0);
            La0=[La0,la0];
            Lo0=[Lo0,lo0];
            [la1,lo1,d1]=getGeo(1.0,0.0,range,lat(idx2),lon(idx2),altitude(idx2),yaw(idx2),0.0);
            La1=[La1,la1];
            Lo1=[Lo1,lo1];
        end
    end
end