Возможность отрисовки траекторий КА

Program.m

@@ -16,7 +16,7 @@ param.k=1.2; % коэффициент для расчёта доверитель
 % N = 4 - Охотское море [24.05.2013]
 % N = 5 - Залив Аляска  [23.01.2018]
 % N = 6 - Бейрут        [04.08.2020]
-N = 6;
+N = 1;
 prev=10;
 next=5;
 param.prev=prev;

analysis.m

@@ -1,4 +1,6 @@
 function param=analysis(param)
+fiC=[]
+LC=[]
 %% By L_DelOff
 % собираю данные по землетрясению
 %% Поехали
@@ -55,11 +57,18 @@ for i=1:length(param.dates)
                 fiB=DATA.DATA(j).GeoLat(k);        % 1. Географическая широта (градусы, север)(Geographic latitude (degrees, north))
                 LB=DATA.DATA(j).GeoLong(k);        % 2. Географическая долгота (градусы, восток)(Geographic longitude (degrees, east))
                 R=mydistance(fiA,LA,fiB,LB);
+                
+                fiC=[fiC fiB];
+                    LC=[LC LB];
+                
                 if R<=R_quake
                     HHeden(end+1)= mymean(DATA.DATA(j).RPA_HHeden((k-1)*5+1:(k-1)*5+5),1e10);
                     O2den(end+1) = mymean(DATA.DATA(j).RPA_O2den((k-1)*5+1:(k-1)*5+5),1e10);
                     Eden(end+1)  = mymean(DATA.DATA(j).EP_Eden((k-1)*5+1:(k-1)*5+5),1e10);
                     Etemp(end+1) = mymean(DATA.DATA(j).EP_Etemp((k-1)*5+1:(k-1)*5+5),1e10);
+                    
+                    
+                    
                 end
             end
         end

etopo1_ice_c_f4.hdr

@@ -0,0 +1,10 @@
+ncols         21600
+nrows         10800
+xllcorner     -180
+yllcorner     -90
+cellsize      0.01666666667
+NODATA_value  -99999
+byteorder     LSBFIRST
+NUMBERTYPE    4_BYTE_FLOAT
+MIN_VALUE     -10803.0
+MAX_VALUE     8333.0

map_draw.m

@@ -0,0 +1,35 @@
+%clear all
+close all
+
+% Загрузим береговые линии (векторные данные)
+load coast
+% загрузим топографические данные (растр)
+load topo
+% отобразим растровые данные с использованием соответствующей цветовой
+% карты
+
+
+[Z, refvec] = etopo('etopo1_ice_c_f4.flt',1, [30 50],[130 150]);
+
+geoshow(Z, refvec,'DisplayType','texturemap')
+demcmap(Z)
+
+
+%geoshow(topo,topolegend,'DisplayType','texturemap')
+%demcmap(topo)
+% нарисуем береговые линии поверх растровой карты
+%geoshow(lat,long,'Color','r')
+axis tight
+hold on
+
+plot(142.37,38.32,'*r','MarkerSize',10)
+
+% load matlab
+% a=1
+% b=12500
+% plot(LC(a:b),fiC(a:b),'r')
+% xlim([130 150])
+% ylim([30 50])
+
+xlabel('Долгота, град')
+ylabel('Широта, град')