how to copy specific columns of a matrix to a table

I have a matrix of six columns each column for a variable (Altitude, Latitude, Longitude, Azimuth, TEC, Ne). and a table of 3 columns for the vraiables (altitude, Ne and date).the matrix we get it from reading some files. so we have a huge number of files (>1M) we read them into matrices then we check if the latittude and the longitude of a certain chosen row in that matrix are within a certain range then we extract only the columns of the altitude and Ne from the matrix an add them to our table. also each matrix is associated with a certain date so we also need to add the coressponding date to the table. I tried implementing what I want using the following lines but I am not sure if this is a correct implementations since I get no errors but also no output. and I am not sure how I can add the date of a certain matrix to the table with the oethr two colmns. if anyone can help..

attached are two samples of the table and the matrix.

I have modified this code to get rid of of the for loop that may be causing the issue but I am still facing the same problem.

  
             %[Height, Latitude, Longitude, Azimuth, TEC, Ne]

       
        EDPdate= datetime(Year, Month, Day, Hour, Min, 0,'TimeZone','UTC','Format', 'uuuuMMddHHmm');
        EDPdate.TimeZone='Asia/Dubai';
        
         EDPT= array2table(EDP);
         
         %(modification is done here) Searching for Profile Peak Parameters' Index 
       [pks,loc] = findpeaks(EDP(:,6), 'SortStr', 'descend', 'NPeaks', 1);
         
         %Searching for Profile Peak Parameters' Index
       %for ied=1:length(EDP)

             

         
          if EDP(loc,2)>(-13.5) && EDP(loc,2)<(-10.5) && EDP(loc,3)>(-79.8) && EDP(loc,3)<(-73.8) 
                  st = 1;
                 for ii = 1:length(EDP)
                   rows = size([tt;EDPT(:,[1,6])],1);%determine the number of rows in tt and EDPT(:,[1,6])
                   en = st+rows - 1;
                   tt3(st:en,ii) = [tt;EDPT(:,[1,6])];% assign found values to tt3
                   st = en +1;
                end
          end 

   
  

[Sijmen Duineveld]

First question: are the values [tt;EDPT(:,[1,6])] what you want?

You can simply add an extra line, with only [tt;EDPT(:,[1,6])] (no ";"), so it displays the values in your command window.

Answer 1

[Sargondjani]

It is not clear what format EDPT has. And also the variable tt3 will be assigned a different value during each loop. If you want store all values tt2 and EDPT(:,[1,6]) during your loop you need to do something like:

st = 1;
for ii = 1:i_lng
    rows = size([tt2;EDPT(:,[1,6])],1);%determine the number of rows in tt2 and EDPT(:,[1,6])
    en = st+rows - 1;
     tt3(st:en,ii) = [tt2;EDPT(:,[1,6])];% assign found values to tt3
     st = en +1;
end

This code is highly inefficient, but lets first try to get your code running. Once it runs you can work on efficiency if speed matters, and for example use concentenate (function cat)

[Salma fathi]

I have edited my code there was one line missing, EDPT is the matrix EDP conveted into table. I have tried the lines you provided but I am still getting no output as if the lines does not do anything in the code. I also added the loop that where we find the index imax from beacause in that loop I am runnning the lines.

[Sargondjani]

The problem is this part:

if EDP(imax,1) < 190

EDP(imax,:)=[];

elseif EDP(imax,1) > 190

break;

end

If the first statement is true, you assign EDP(imax,:) = []; which means the row EDP(imax,:) has only empty values.

In the second case you break the loop (ie. you go to the next iteration ied).

So what was the intention of this part of your code?

[Salma fathi]

this loop is for finding the index of what we call the peak of a profile, which is the max value between the values that has Altitude > 190. but i don't understand, whatever I am trying to do is being executed after this loop finishes all its iterations, so shouldn't it not affect whatever is executed next. please correct me if I am wrong and how e may fix that?

[Sargondjani]

Ok, so you have the location of the maximum value in column 6 which is 'loc'.

And then you check if the point is within a certain range of longitude and latitude:

if EDP(loc,2)>(-13.5) && EDP(loc,2)<(-10.5) && EDP(loc,3)>(-79.8) && EDP(loc,3)<(-73.8)

You should first check if this condition is ever true, because there has to be a peak which is within this range:

So after the 'if statement, put a line:

if EDP(loc,2)>(-13.5) && EDP(loc,2)<(-10.5) && EDP(loc,3)>(-79.8) && EDP(loc,3)<(-73.8)

display('If statement is passed. Lets have a look what happens next')

tt

EDPT(:,[1,6])]

%put a breakpoint here by right clicking in the left margin

end

Run the code, and see if this gives you what you expect. The variables tt and EDPT(:,[1,6])] are now displayed on screen.