1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
|
#include "poseidon_reader.h"
#include "common/ccsds/ccsds_time.h"
namespace jason3
{
namespace poseidon
{
PoseidonReader::PoseidonReader()
{
frames = 0;
}
PoseidonReader::~PoseidonReader()
{
}
void PoseidonReader::work(ccsds::CCSDSPacket &packet)
{
if (packet.payload.size() < 930)
return;
frames++;
// We need to know where the satellite was when that packet was created
double currentTime = ccsds::parseCCSDSTimeFull(packet, 16743, 1);
timestamps.push_back(currentTime);
/*predict_orbit(jason3_object, &jason3_orbit, predict_to_julian(currentTime));
// Scale to the map
int map_height = map_image_height.height();
int map_width = map_image_height.width();
int imageLat = map_height - ((90.0f + (jason3_orbit.latitude * 180.0f / M_PI)) / 180.0f) * map_height;
int imageLon = ((jason3_orbit.longitude * 180.0f / M_PI) / 360.0f) * map_width + (map_width / 2);
if (imageLon >= map_width)
imageLon -= map_width;
// Here is a very naive way of detecting when we got a radar echo...
// And how stable it was...
//
// The first just checks the position of the first value over a thresold,
// while the second counts all values over that thresold.
// All 8 samples are averaged to smooth the data out.
//
// My first time dealing with this kind of data so well...
// Hopefulyl not too bad :-)
int thresold = 10;
int firstOverThresoldAvg = 0, countOverThresoldAvg = 0;
for (int y = 0; y < 8; y++)
{
int firstOverThresold = 0;
for (int i = 0; i < 104; i++)
{
if (packet.payload[96 + (104 * y) + i] > thresold)
{
firstOverThresold = i;
break;
}
}
int countOverThresold = 0;
for (int i = 0; i < 104; i++)
{
if (packet.payload[96 + (104 * y) + i] > thresold)
countOverThresold++;
}
firstOverThresoldAvg += firstOverThresold;
countOverThresoldAvg += countOverThresold;
}
// Average
firstOverThresoldAvg /= 8;
countOverThresoldAvg /= 8;
// "Scatter" samples
int sampleScatter = 65 + ((float)countOverThresoldAvg / 104.0f) * 255;
if (sampleScatter < 0)
sampleScatter = 0;
if (sampleScatter > 255)
sampleScatter = 255;
// Height samples
int sampleHeight = 200 - ((float)firstOverThresoldAvg / 104.0f) * 255;
if (sampleHeight < 0)
sampleHeight = 0;
if (sampleHeight > 255)
sampleHeight = 255;
// Write on the map
unsigned char colorHeight[] = {(unsigned char)sampleHeight, (unsigned char)std::max(0, 255 - sampleHeight), 0};
map_image_height.draw_circle(imageLon, imageLat, 2, colorHeight);
// Write on the map
unsigned char colorScatter[] = {(unsigned char)sampleScatter, (unsigned char)std::max(0, 255 - sampleScatter), 0};
map_image_scatter.draw_circle(imageLon, imageLat, 2, colorScatter);*/
}
} // namespace modis
} // namespace eos
|
