Indexing matrix to get corresponding values for condition,store in a new matrix each time?

Sorry for the perhaps confusing title...

Basically I have a 3x3 matrix containing elevation angle, azimuth angle and range. I want to generate new matrices each time elevation >5 deg. There are usually about 5 segments that have this data and I want to separate each one into a new matrix.

I know how to index but not sure how to put this condition in...

Thanks

sat_tcs=llh2tcsT(sat_llh,station_llh);
sat_elev=atan2(sat_tcs(3,:),sqrt(sat_tcs(1,:).^2+sat_tcs(2,:).^2)); sat_azim=atan2(-sat_tcs(2,:),sat_tcs(1,:));
range=sqrt(sat_tcs(1,:).^2+sat_tcs(2,:).^2+sat_tcs(3,:).^2);`                                  sat_elev(sat_elev < 5*deg2rad) = NaN;  
sat_look_tcs=[sat_elev;sat_azim;range];

Answer 1

[David Sanchez]

It is not very clear what you want, but I think it has to be something like this:

sat_look_tcs=[sat_elev;sat_azim;range];
elev_5 = find(sat_elev >5);
sat_look_tcs_5 = zeros(length(elev_5, 3, 3)); % initialize 3D matrix to hold your submatrices

for k=1:(length(elev_5)-1)
    sat_look_tcs_5(k+1,:,:) = sat_look_tcs((elev_5(k)+1):elev_5(k+1),:,:);
end

And you access your matrix with

sat_look_tcs_5(1,:,:)
sat_look_tcs_5(2,:,:)

etc...

[Buzz]

Thanks for your help. I tried your suggestion and it gives me an error for the 3d matrix saying there are too many inputs.

[Michael Haderlein]

There's a typo. It should read

sat_look_tcs_5 = zeros(length(elev_5),3,3);

I didn't check the rest of the code as I also didn't really understand the question. But maybe that's already what you need.

[Buzz]

Hi, yeah thanks. The problem I have is that I am getting a 3x3 with the elevation angle in each column, rather than the corresponding azimuth and range for each value.