Skip to content
Snippets Groups Projects
Commit 95f2b03d authored by Oliver Lemke's avatar Oliver Lemke
Browse files

Add real maximum dsv counts to output file

parent 4439b9de
No related branches found
No related tags found
No related merge requests found
......@@ -12,6 +12,7 @@ ret.mhsfile = filename;
ret.x = -1;
ret.y = -1;
ret.maximum = -1;
ret.maxcounts = -1;
ret.lon=-1;
ret.lat=-1;
ret.height=-1;
......@@ -139,15 +140,15 @@ if intr_end < intr_start
disp('Not enough data in file')
return
end
y1 = double(squeeze(data.dsvcounts11));
y2 = double(squeeze(data.dsvcounts21));
y3 = double(squeeze(data.dsvcounts31));
y4 = double(squeeze(data.dsvcounts41));
y{1} = double(squeeze(data.dsvcounts11));
y{2} = double(squeeze(data.dsvcounts21));
y{3} = double(squeeze(data.dsvcounts31));
y{4} = double(squeeze(data.dsvcounts41));
for i = intr_start:intr_end
A(1,i) = median(y1(i-5:i+5))-(median(y1(i-70:i-60))+median(y1(i+60:i+70)))/2;
A(2,i) = median(y2(i-5:i+5))-(median(y2(i-70:i-60))+median(y2(i+60:i+70)))/2;
A(3,i) = median(y3(i-5:i+5))-(median(y3(i-70:i-60))+median(y3(i+60:i+70)))/2;
A(4,i) = median(y4(i-5:i+5))-(median(y4(i-70:i-60))+median(y4(i+60:i+70)))/2;
A(1,i) = median(y{1}(i-5:i+5))-(median(y{1}(i-70:i-60))+median(y{1}(i+60:i+70)))/2;
A(2,i) = median(y{2}(i-5:i+5))-(median(y{2}(i-70:i-60))+median(y{2}(i+60:i+70)))/2;
A(3,i) = median(y{3}(i-5:i+5))-(median(y{3}(i-70:i-60))+median(y{3}(i+60:i+70)))/2;
A(4,i) = median(y{4}(i-5:i+5))-(median(y{4}(i-70:i-60))+median(y{4}(i+60:i+70)))/2;
end
M(1) = max(A(1,:));
M(2) = max(A(2,:));
......@@ -160,6 +161,7 @@ if size(ret.x, 1)>1
ret.y = ret.y(1);
end
ret.maximum = repmat(maximum, size(ret.x));
ret.maxcounts = y{ret.x}(ret.y);
ret.lon=((double(data.lonfov45(ret.y))+double(data.lonfov46(ret.y)))/2E4)'; % longitude of nadir at maximum of light curve
ret.lat=((double(data.latfov45(ret.y))+double(data.latfov46(ret.y)))/2E4)'; % latitude of nadir at maximum of light curve
ret.height=(double(data.height(ret.y))/10)'; % altitude of satellite at maximum of light curve
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment