function getFlightFromIGC($filename)
{
if (!is_file($filename)) {
DEBUG("IGC", 1, "getFlightFromIGC: File was not found:{$filename}<br>");
return 0;
}
set_time_limit(100);
if ($this->forceBounds) {
$startTime = $this->START_TIME;
$endTime = $this->END_TIME;
}
$this->resetData();
if ($this->forceBounds) {
$this->START_TIME = $startTime;
$this->END_TIME = $endTime;
}
$this->setAllowedParams();
$this->filename = basename($filename);
$this->checkForOLCfile($filename);
$done = 0;
$try_no_takeoff_detection = 0;
$tm1 = time() + 3600 * 24 * 30;
$tm2 = 0;
while (!$done) {
$lines = file($filename);
$linesNum = count($lines);
DEBUG("IGC", 1, "File has total {$linesNum} lines<br>");
$points = 0;
$outputBuffer = "";
// process points
// filter bad ones
$p = 0;
$Brecords = 0;
$getPointsNum = 5;
// + 4 points + this one
for ($i = 0; $i < count($lines) - 10; $i++) {
$pointOK = 1;
$line = trim($lines[$i]);
if (strlen($line) == 0) {
continue;
}
if ($line[0] != 'B') {
continue;
}
/* if ($i>2){
$prevline=trim($lines[$i-1]);
if (!substr_compare($line,$prevline,1,6)) continue;
//throw out Brecords with same time signatures P. Wild 4.11.2012
} */
$Brecords++;
// also check for bad points
// 012345678901234567890123456789
// B1522144902558N00848090EV0058400000
if (strlen($line) < 23 || $line[24] == 'V') {
$lines[$i][1] = 'X';
continue;
}
$neighboors = array();
$nextPointPos = $i;
for ($t1 = 0; $t1 < $getPointsNum; $t1++) {
$thisPoint = new gpsPoint(trim($lines[$nextPointPos]), $this->timezone);
$neighboors[$t1] = $thisPoint;
$nextPointPos = $this->getNextPointPos($lines, $nextPointPos);
}
// got all next points
// find mean values
$mean_speed = 0;
$mean_vario = 0;
for ($t1 = 1; $t1 < $getPointsNum; $t1++) {
// for 4 (5-1) points in a row
// create arrays
$T_distance[$t1] = $neighboors[$t1]->calcDistance($neighboors[$t1 - 1]);
$T_alt[$t1] = $neighboors[$t1]->getAlt();
$T_deltaseconds[$t1] = $neighboors[$t1]->getTime() - $neighboors[$t1 - 1]->getTime();
$T_speed[$t1] = $T_deltaseconds[$t1] ? $T_distance[$t1] * 3.6 / $T_deltaseconds[$t1] : 0.0;
/* in km/h */
if ($T_deltaseconds[$t1]) {
$T_vario[$t1] = ($T_alt[$t1] - $neighboors[$t1 - 1]->getAlt()) / $T_deltaseconds[$t1];
}
$mean_speed += $T_speed[$t1];
$mean_vario += $T_vario[$t1];
}
$mean_speed = $mean_speed / ($getPointsNum - 1);
$mean_vario = ($neighboors[$getPointsNum - 1]->getAlt() - $neighboors[0]->getAlt()) / ($neighboors[$getPointsNum - 1]->getTime() - $neighboors[0]->getTime());
//if ($mean_vario<0) {
// DEBUG("IGC",8,"[$Brecords-$p] mean_vario :$mean_vario <0 <br>");
//}
$data_speed[$i] = $mean_speed;
$data_vario[$i] = $mean_vario;
// $mean_vario = $mean_vario/($getPointsNum-1); // mean vario is wrong
if ($neighboors[0]->getTime() < $tm1) {
$tm1 = $neighboors[0]->getTime();
}
if ($neighboors[0]->getTime() > $tm2) {
$tm2 = $neighboors[0]->getTime();
}
if ($neighboors[0]->getTime() - $neighboors[1]->getTime() > 0) {
// the next point is more than one hour in the past
// echo "#"; $pointOK=0;
}
if ($T_distance[1] < 0.5 && $T_deltaseconds[1] > 2) {
// less than 0.5 m distance
$pointOK = 0;
DEBUG("IGC", 8, "[{$Brecords}-{$p}] Distance <0.5 m <br>");
$REJ_T_distance++;
// we dont go through the other tests
//echo "{$lines[$i]} =>";
$lines[$i][1] = 'X';
//echo "{$lines[$i]}<br>";
continue;
//echo $T_distance[1]."*<br>";
}
if (abs($mean_speed - $T_speed[1]) > 40) {
// diff more than 40 km/h
$pointOK = 0;
$REJ_T_mean_speed++;
DEBUG("IGC", 8, "[{$Brecords}-{$p}] Mean speed > 40 km/h <br>");
//echo "@";
}
//if ( abs ($mean_vario - $T_vario[1] ) > 6 ) { // diff more than 6 m/sec
// $pointOK=0;
//echo "#";
//}
if ($T_deltaseconds[1] == 0) {
$pointOK = 0;
$REJ_T_zero_time_diff++;
DEBUG("IGC", 8, "[{$Brecords}-{$p}] No time Diff<br>");
}
if ($T_alt[1] > $this->maxAllowedHeight) {
$pointOK = 0;
$REJ_max_alt++;
}
if (abs($T_speed[1]) > $this->maxAllowedSpeed) {
$pointOK = 0;
$REJ_max_speed++;
DEBUG("IGC", 8, "[{$Brecords}-{$p}] > " . abs($T_speed[1]) . "km/h max allowed speed<br>");
// echo "S";
}
if (abs($T_vario[1]) > $this->maxAllowedVario) {
$pointOK = 0;
$REJ_max_vario++;
// echo "V";
DEBUG("IGC", 8, "[{$Brecords}-{$p}] > " . abs($T_vario[1]) . "m/sec > max allowed vario<br>");
}
if ($p < 5 && !$try_no_takeoff_detection && !$this->forceBounds) {
// first 5 points need special care
$takeoffMaxSpeed = $this->maxAllowedSpeed * 0.5;
DEBUG("IGC", 8, "[{$Brecords}-{$p}] TAKEOFF sequence SPEED: " . abs($T_speed[1]) . " max:{$takeoffMaxSpeed}<br>");
if (abs($T_speed[1]) > $takeoffMaxSpeed) {
$pointOK = 0;
$REJ_max_speed_start++;
// echo "s";
}
if (abs($T_vario[1]) > $this->maxAllowedVario * 0.4) {
$pointOK = 0;
$REJ_max_vario_start++;
//echo "v";
}
}
if (!$pointOK) {
$lines[$i][1] = 'X';
} else {
$p++;
if ($p == 5) {
DEBUG("IGC", 1, "Passed the strict testing (p=5)<br>");
}
}
}
DEBUG("IGC", 1, "REJ: [{$REJ_T_distance}] <0.5 distance<br>");
DEBUG("IGC", 1, "REJ: [{$REJ_T_zero_time_diff}] zero_time_diff<br>");
DEBUG("IGC", 1, "REJ: [{$REJ_T_mean_speed}] mean_speed diff >40km/h<br>");
DEBUG("IGC", 1, "REJ: [{$REJ_max_alt}] >max_alt<br>");
DEBUG("IGC", 1, "REJ: [{$REJ_max_speed}] >max_speed<br>");
DEBUG("IGC", 1, "REJ: [{$REJ_max_vario}] >max_vario<br>");
DEBUG("IGC", 1, "REJ: [{$REJ_max_speed_start}] >max_speed_start<br>");
DEBUG("IGC", 1, "REJ: [{$REJ_max_vario_start}] >max_vario_start<br>");
DEBUG("IGC", 1, "Found {$p} valid B records out of {$Brecords} total<br>");
if ($p > 0) {
$done = 1;
} else {
if ($Brecords > 0 && $REJ_T_zero_time_diff / $Brecords > 0.9) {
// more than 90% stopped points
$lines = file($filename);
$done = 1;
$garminSpecialCase = 1;
$p = $Brecords;
DEBUG("IGC", 1, "Many Stopped points, it is a Garmin Special Case<br>");
} else {
if (!$try_no_takeoff_detection) {
$try_no_takeoff_detection = 1;
DEBUG("IGC", 1, "Will try no_takeoff_detection<br>");
} else {
$done = 1;
}
}
}
}
// while not done
//
if ($p == 0) {
DEBUG("IGC", 1, "NO VALID POINTS FOUND");
return 0;
// no valid points found
}
$mod = 0;
if ($p > $this->maxPointNum) {
$reduceArray = getReduceArray($p, $this->maxPointNum);
// print_r($recudeArray);
$mod = count($reduceArray);
// $mod= ceil( $p / $this->maxPointNum );
}
DEBUG("IGC", 1, "will use a reduce array of length {$mod}<br>");
$duration = $tm2 - $tm1;
$intervalSecs = round($duration / $p);
// echo "<hr>good points: $p duration:$duration, intervalSecs:$intervalSecs<hr>";
$pointsNeededForTakeoff = 5;
if ($intervalSecs >= 8) {
$pointsNeededForTakeoff = 2;
}
$alreadyInPoints = 0;
$stopReadingPoints = 0;
$this->timezone = 1000;
$day_offset = 0;
$foundNewTrack = 0;
$slow_points = 0;
$slow_points_dt = 0;
$stillOnGround = 1;
$tmpDate = 0;
foreach ($lines as $iii => $line) {
if ($foundNewTrack) {
break;
}
$outputLine = $line;
$line = trim($line);
if (strlen($line) == 0) {
continue;
}
// if (strtoupper(substr($line,0,1)) !="B" ) echo "@$line<BR>";
if (strtoupper(substr($line, 0, 3)) == "OLC") {
continue;
}
// it is an olc file , dont put the OLC... line in to the saned file
if (strtoupper(substr($line, 0, 5)) == "HFDTE" || strtoupper(substr($line, 0, 5)) == "HPDTE") {
// HFDTE170104 OR HPDTE310805
if ($alreadyInPoints && $points > 0) {
if ($prevPoint->gpsTime < 86200) {
// if last good point is > 86200 (200 secs before day change at 86400) we dont treat this as a new track
$stopReadingPoints = 1;
DEBUG("IGC", 1, "[{$points}] {$line}<br>");
DEBUG("IGC", 1, "[{$points}] Found a new track (NEW HFDTE)<br>");
}
} else {
$this->DATE = substr($line, 5, 6);
$yr_last = substr($this->DATE, 4, 2);
// case of YY=0 (1 digit) HFDTE08070
if ($yr_last == "0") {
$yr_last = "00";
}
if ($yr_last > 80) {
$yr = "19" . $yr_last;
} else {
$yr = "20" . $yr_last;
}
$this->DATE = $yr . "-" . substr($this->DATE, 2, 2) . "-" . substr($this->DATE, 0, 2);
$alreadyInPoints = 1;
}
} else {
if (strtoupper(substr($line, 0, 13)) == "HFTZOTIMEZONE") {
// HFTZOTimezone:3 OR HFTZOTimezone:-8
$this->timezone = substr($line, 14) + 0;
// echo $this->timezone."#^^";
} else {
if (strtoupper(substr($line, 2, 13)) == "GTYGLIDERTYPE") {
// HOGTYGLIDERTYPE: Gradient Bliss 26 OR HPGTYGliderType:Gradient Nevada
if (!$this->glider) {
$this->glider = trim(substr($line, 16));
}
// HFGTYGLIDERTYPE
// HOGTYGLIDERTYPE
} else {
if ($line[0] == 'B') {
if ($stopReadingPoints) {
continue;
}
if ($line[1] == 'X') {
DEBUG("IGC", 1, "[{$points}] BAD : {$line}<br>");
continue;
// MARKED BAD from BEFORE
}
if (strlen($line) < 23 || strlen($line) > 100) {
continue;
}
if ($points == 0) {
// first point
// echo "######## first point <br>";
$firstPoint = new gpsPoint($line, $this->timezone);
if ($this->timezone == 1000) {
// no timezone in the file
// echo "calc timezone<br>";
$this->timezone = getUTMtimeOffset($firstPoint->lat, $firstPoint->lon, $this->DATE);
$this->timezone = getTZ($firstPoint->lat, $firstPoint->lon, $this->DATE);
// echo $this->timezone;
$firstPoint->timezone = $this->timezone;
}
$tmpTime = $firstPoint->gpsTime + $this->timezone * 3600;
// Now also check if we are one day minus (for US flights )
if ($tmpTime < 0 && $tmpDate == 0) {
// one day before!
$this->DATE = dates::moveDaysFromDate($this->DATE, -1);
$tmpDate = 1;
}
// take care of day change for timezones in australia/nz etc
if ($tmpTime > 86400 && $tmpDate == 0) {
// UTC date in the igc file needs to be +1
$this->DATE = dates::moveDaysFromDate($this->DATE, 1);
$tmpDate = 1;
}
// sanity checks
if ($firstPoint->getAlt() > $this->maxAllowedHeight) {
continue;
}
// echo "times: ".$firstPoint->gpsTime.", ".$firstPoint->getTime()." start_time: ".$this->START_TIME ."<BR> ";
if ($this->forceBounds && !$this->checkBound($firstPoint->getTime())) {
continue;
}
// not inside time window
//$this->FIRST_POINT=$line;
$this->firstPointTM = $firstPoint->gpsTime;
$this->firstLat = $firstPoint->lat();
$this->firstLon = $firstPoint->lon();
$this->TAKEOFF_ALT = $firstPoint->getAlt();
$this->MIN_ALT = $firstPoint->getAlt();
if (!$this->forceBounds) {
$this->START_TIME = $firstPoint->getTime();
}
$prevPoint = new gpsPoint($line, $this->timezone);
$prevPoint2 = new gpsPoint($line, $this->timezone);
} else {
$lastPoint = new gpsPoint($line, $this->timezone);
$lastPoint->gpsTime += $day_offset;
if ($this->forceBounds && !$this->checkBound($lastPoint->getTime())) {
$lastPoint = $prevPoint;
continue;
// not inside time window
}
$time_diff = $lastPoint->getTime() - $prevPoint->getTime();
$time_diff2 = $lastPoint->getTime() - $prevPoint2->getTime();
// echo "time diff: $time_diff # $line<br>";
if ($time_diff < 0 && $time_diff > -36000) {
// if the time is less than 10 hours in the past we just ignore it
// $day_offset = 86400; // if time seems to have gone backwards, add a day
DEBUG("IGC", 1, "[{$points}] {$line}<br>");
DEBUG("IGC", 1, "[{$points}] Point in the past<br>");
continue;
} else {
if ($time_diff < 0 && $time_diff > -86000) {
// CHANGING DAY , means the flight is at night
$lastPoint = $prevPoint;
$foundNewTrack = 1;
DEBUG("IGC", 1, "[{$points}] {$line}<br>");
DEBUG("IGC", 1, "[{$points}] Flight at night ????<br>");
continue;
} else {
if ($time_diff > $this->max_allowed_time_gap) {
// found time gap
// if we are forceBounds check to see if inside window
if ($this->forceBounds && !$this->checkBound($lastPoint->getTime()) || !$this->forceBounds) {
// if we are still on the ground we dont care about time gap
if (!$stillOnGround) {
// not inside time window OR not checking go ahead
$lastPoint = $prevPoint;
$foundNewTrack = 1;
DEBUG("IGC", 1, "[{$points}] {$line}<br>");
DEBUG("IGC", 1, "[{$points}] Found a new track (Time diff of {$time_diff} secs)<br>");
continue;
}
} else {
// inside the window, forced to continue
}
}
}
}
$this->LAST_POINT = $line;
// compute some things
$tmp = $lastPoint->calcDistance($prevPoint);
$alt = $lastPoint->getAlt();
// GUS begin
$deltaseconds = $lastPoint->getTime() - $prevPoint->getTime();
$speedDeltaSecs = $deltaseconds + ($deltaseconds < -85000 ? 86400 : 0);
$speed = $speedDeltaSecs ? $tmp * 3.6 / $speedDeltaSecs : 0.0;
/* in km/h */
if ($speedDeltaSecs) {
$vario = ($alt - $prevPoint->getAlt()) / $speedDeltaSecs;
}
// GUS end
if (!$garminSpecialCase && !$this->forceBounds) {
if (($fast_points >= $pointsNeededForTakeoff || $fast_points_dt > 30) && $stillOnGround) {
// found 5 flying points or 30 secs
$stillOnGround = 0;
DEBUG("IGC", 1, "[{$points}] {$line}<br>");
DEBUG("IGC", 1, "[{$points}] Found Takeoff <br>");
}
if ($stillOnGround) {
//takeoff scan
// either speed >= 15 or if we already have 2 fast points settle with speed>=10
if ($speed >= 15 || $speed >= 10 && $fast_points >= 2) {
$fast_points++;
$fast_points_dt += $speedDeltaSecs;
DEBUG("IGC", 1, "[{$points}] {$line}<br>");
DEBUG("IGC", 1, "[{$points}] Found a fast speed point <br>");
} else {
// reset takeoff scan
DEBUG("IGC", 1, "[{$points}] {$line}<br>");
DEBUG("IGC", 1, "[{$points}] takeoff scan: speed: {$speed} time_diff: {$time_diff}<br>");
$fast_points = 0;
$fast_points_dt = 0;
}
$points = 0;
continue;
} else {
//landing scan
if ($speed < 5) {
$slow_points++;
$slow_points_dt += $speedDeltaSecs;
DEBUG("IGC", 1, "[{$points}] {$line}<br>");
DEBUG("IGC", 1, "[{$points}] Found a slow speed point (speed, dt)=({$speed},{$speedDeltaSecs})<br>");
} else {
$slow_points = 0;
$slow_points_dt = 0;
}
}
// found landing (5 stopped points and >2mins) or 5 mins (300 secs)
if ($slow_points > $pointsNeededForTakeoff && $slow_points_dt > 180 || $slow_points_dt > 300) {
$foundNewTrack = 1;
DEBUG("IGC", 1, "[{$points}] {$line}<br>");
DEBUG("IGC", 1, "[{$points}] Found a new track /landing <br>");
}
}
// sanity checks
if ($deltaseconds == 0 && !$garminSpecialCase) {
continue;
}
if ($alt > $this->maxAllowedHeight) {
continue;
}
if (abs($speed) > $this->maxAllowedSpeed) {
continue;
}
if (abs($vario) > $this->maxAllowedVario) {
continue;
}
$takeoffDistance = $lastPoint->calcDistance($firstPoint);
if ($takeoffDistance > $this->LINEAR_DISTANCE) {
$this->LINEAR_DISTANCE = $takeoffDistance;
}
if ($time_diff2 > 10) {
$tmp = $lastPoint->calcDistance($prevPoint2);
$alt = $lastPoint->getAlt();
$deltaseconds = $time_diff2;
// GUS begin
$speedDeltaSecs = $deltaseconds + ($deltaseconds < -85000 ? 86400 : 0);
$speed = $speedDeltaSecs ? $tmp * 3.6 / $speedDeltaSecs : 0.0;
/* in km/h */
if ($speedDeltaSecs) {
$vario = ($alt - $prevPoint2->getAlt()) / $speedDeltaSecs;
}
// GUS end
$prevPoint2 = new gpsPoint($line, $this->timezone);
$prevPoint2->gpsTime += $day_offset;
// update maximum speed
if ($speed > $this->MAX_SPEED) {
$this->MAX_SPEED = $speed;
}
$this->MEAN_SPEED += $speed;
$MEAN_SPEED_POINTS++;
// UPDATE MIN-MAX VARIO
if ($vario > $this->MAX_VARIO) {
$this->MAX_VARIO = $vario;
}
if ($vario < $this->MIN_VARIO) {
$this->MIN_VARIO = $vario;
}
}
//$speed=$data_speed[$iii-1]+0;
// $vario=$data_vario[$iii-1]+0;
// UPDATE MIN-MAX ALT
if ($alt > $this->MAX_ALT) {
$this->MAX_ALT = $alt;
}
if ($alt < $this->MIN_ALT) {
$this->MIN_ALT = $alt;
}
// end computing
$prevPoint = new gpsPoint($line, $this->timezone);
$prevPoint->gpsTime += $day_offset;
if ($mod >= 1) {
if ($reduceArray[$points % $mod] == 0) {
$outputLine = "";
}
}
}
$points++;
}
}
}
}
// end else
$outputBuffer .= $outputLine;
}
// end main loop
// echo "<HR>MIN VARIO".$this->MIN_VARIO."<HR>";
//
if ($stillOnGround && $this->LINEAR_DISTANCE < 50) {
DEBUG("IGC", 1, "NO TAKEOFF FOUND: ");
return 0;
// no valid points found
}
$path_igc = dirname($this->getIGCFilename(1));
if (!is_dir($path_igc)) {
makeDir($path_igc, 0755);
}
/*write the full saned file */
$fullSanedFile = '';
foreach ($lines as $line) {
$line = trim($line);
if (strlen($line) == 0) {
continue;
}
if ($line[0] == 'B' && $line[1] == 'X') {
continue;
}
// MARKED BAD from BEFORE
// if ( strlen($line) < 23 || strlen($line) > 100 ) continue;
$fullSanedFile .= $line . "\n";
}
if (!writeFile($this->getIGCFilename(2), $fullSanedFile)) {
echo "Problem writing to file (" . $this->getIGCFilename(2) . ")";
}
DEBUG("IGC", 1, "<HR>" . $this->getIGCFilename(2) . ' size: ' . strlen($fullSanedFile) . "<HR>");
// echo "<HR><HR>". $this->getIGCFilename(2) .strlen($fullSanedFile)."<HR><HR>";
/* done wrting the full saned file */
// write saned IGC file
if (!writeFile($this->getIGCFilename(1), $outputBuffer)) {
echo "Problem writing to file (" . $this->getIGCFilename(1) . ")";
}
// done write saned IGC file
DEBUG("IGC", 1, "<HR>" . $this->getIGCFilename(1) . ' size: ' . strlen($outputBuffer) . "<HR>");
if ($lastPoint) {
$this->lastPointTM = $lastPoint->gpsTime;
$this->lastLon = $lastPoint->lon();
$this->lastLat = $lastPoint->lat();
$this->LANDING_ALT = $lastPoint->getAlt();
$this->END_TIME = $lastPoint->getTime();
} else {
$this->lastPointTM = 0;
$this->lastLon = 0;
$this->lastLat = 0;
$this->LANDING_ALT = 0;
$this->END_TIME = 0;
}
$this->DURATION = $this->END_TIME - $this->START_TIME;
if ($this->DURATION < 0) {
$this->DURATION += 86400;
}
$this->MEAN_SPEED = $this->MEAN_SPEED / $MEAN_SPEED_POINTS;
return 1;
}
function checkAirspace($filename)
{
$lines = file($filename);
if (!$lines) {
echo "Cant read file";
return;
}
$i = 0;
// find bounding box of flight
$min_lat = 1000;
$max_lat = -1000;
$min_lon = 1000;
$max_lon = -1000;
foreach ($lines as $line) {
$line = trim($line);
if (strlen($line) == 0) {
continue;
}
if ($line[0] == 'B') {
if (strlen($line) < 23) {
continue;
}
$thisPoint = new gpsPoint($line, 0);
if ($thisPoint->lat > $max_lat) {
$max_lat = $thisPoint->lat;
}
if ($thisPoint->lat < $min_lat) {
$min_lat = $thisPoint->lat;
}
if (-$thisPoint->lon > $max_lon) {
$max_lon = -$thisPoint->lon;
}
if (-$thisPoint->lon < $min_lon) {
$min_lon = -$thisPoint->lon;
}
}
}
// echo "$min_lat, $max_lat,$min_lon, $max_lon<BR>";
// now find the bounding boxes that have common points
// !( A1<X0 || A0>X1 ) && !( B1<Y0 || B0>Y1 )
// X,A -> lon
// Y,B -> lat
// X0 -> $min_lon A0-> $area->minx
// X1 -> $max_lon A1-> $area->maxx
// Y0 -> $min_lat B0-> $area->miny
// Y1 -> $max_lat B1-> $area->maxy
// !( $area->maxx<$min_lon || $area->minx>$max_lon ) && !( $area->maxx<$min_lat || $area->miny>$max_lat )
global $AirspaceArea, $NumberOfAirspaceAreas;
getAirspaceFromDB($min_lon, $max_lon, $min_lat, $max_lat);
$NumberOfAirspaceAreas = count($AirspaceArea);
/*
foreach($AirspaceArea as $i=>$area) {
if ( !( $area->maxx<$min_lon || $area->minx>$max_lon ) &&
!( $area->maxy<$min_lat || $area->miny>$max_lat )
) {
if ($area->Shape==1) $shape="Area "; else $shape="Circle";
DEBUG("checkAirspace",1, "Found $shape [$i] => ".$area->Name.'<BR>');
// print_r($area);
if ($area->Shape==1) {
$area->Points[]=$area->Points[0];
// $area->NumPoints=count($area->Points);
}
$selAirspaceArea[]=$area;
}
}
$AirspaceArea=$selAirspaceArea;
$NumberOfAirspaceAreas=count($AirspaceArea);
*/
// echo '<HR>';
// print_r($AirspaceArea);
// echo '<HR>';
// echo '<HR>';
$violations = array();
$i = 0;
foreach ($lines as $line) {
$line = trim($line);
if (strlen($line) == 0) {
continue;
}
if ($line[0] == 'B') {
if (strlen($line) < 23) {
continue;
}
$thisPoint = new gpsPoint($line, 0);
$alt = $thisPoint->getAlt(1);
// prefer vario alt
// $insideArea=-1;
$insideAreas = array();
$insideAreas = FindAirspaceArea(-$thisPoint->lon, $thisPoint->lat, $alt);
if (count($insideAreas) > 0) {
//echo "point [$i] INSIDE AIRSPACE areas: ";
//foreach($insideAreas as $areaInfo) echo $AirspaceArea[$areaInfo[0]]->Name." areaID[$areaInfo[0]] disInside[$areaInfo[1]] altInside[$areaInfo[2]] ";
foreach ($insideAreas as $areaInfo) {
$areaID = $areaInfo[0];
$disInside = $areaInfo[1];
$altInside = $areaInfo[2];
if ($disInside > $violations[$areaID]['maxDistance']) {
$violations[$areaID]['maxDistance'] = $disInside;
}
if ($altInside > $violations[$areaID]['maxAlt']) {
$violations[$areaID]['maxAlt'] = $altInside;
}
}
//echo "<BR>";
} else {
// echo "OK<BR>";
}
$i++;
}
}
$resStr = '';
$resStr1 = '';
if (count($violations) > 0) {
foreach ($violations as $i => $violatedArea) {
$resStr1 .= $AirspaceArea[$i]->id . ",";
$resStr .= "HorDist: " . floor($violatedArea['maxDistance']) . "m, VertDist:" . floor($violatedArea['maxAlt']) . "m, ";
$resStr .= 'Airspace: ' . $AirspaceArea[$i]->Name . ' [' . $AirspaceArea[$i]->Type . '] ' . floor($AirspaceArea[$i]->Base->Altitude) . '-' . floor($AirspaceArea[$i]->Top->Altitude) . "m\n";
// COMMENT: '.$AirspaceArea[$i]->Comment."\n";
}
if ($resStr1) {
$resStr1 = substr($resStr1, 0, -1) . "\n";
$resStr = trim($resStr);
}
}
//$m1=memory_get_usage();
//echo "ReadAltitude: mem usage: $m1 <BR>";
return $resStr1 . $resStr;
}
function drawFlight()
{
global $gradientMap, $gradient_colors_num;
global $minAlt, $maxAlt;
$size = 1;
$startSize = 10;
$pointColor = imagecolorallocate($this->img, 255, 0, 0);
$startColor = imagecolorallocate($this->img, 0, 250, 0);
$endColor = imagecolorallocate($this->img, 252, 250, 66);
// yellow
$shadowColor = imagecolorallocate($this->img, 0, 0, 0);
$lines = file($this->IGCfilename);
$points = 0;
$minAlt = 10000;
$maxAlt = 0;
//first pass to find min/max
foreach ($lines as $line) {
$line = trim($line);
if (strlen($line) == 0) {
continue;
}
if ($line[0] == 'B') {
$thisPoint = new gpsPoint($line);
$thisAlt = $thisPoint->getAlt();
if ($thisAlt < $minAlt) {
$minAlt = $thisAlt;
}
if ($thisAlt > $maxAlt) {
$maxAlt = $thisAlt;
}
}
}
// Make shadow first
foreach ($lines as $line) {
// $line=trim($line);
if (strlen($line) == 0) {
continue;
}
if ($line[0] == 'B') {
$thisPoint = new gpsPoint($line);
$thisColor = $gradientMap[floor(($gradient_colors_num - 1) * $altRelPos)];
list($ptX, $ptY) = $this->lonlat2xy(-$thisPoint->lon, $thisPoint->lat);
//DEBUG("MAP",8,"lon:".$thisPoint->lon." lat:".$thisPoint->lat." -> X:$ptX, Y:$ptY <br>");
$ptY = $this->pixelHeight - $ptY;
if ($points == 0) {
$firstPoint = new gpsPoint($line);
} else {
list($lastptX, $lastptY) = $this->lonlat2xy(-$lastPoint->lon, $lastPoint->lat);
$lastptY = $this->pixelHeight - $lastptY;
imageline($this->img, $lastptX - 1, $lastptY - 1, $ptX - 1, $ptY - 1, $shadowColor);
imageline($this->img, $lastptX, $lastptY - 1, $ptX, $ptY - 1, $shadowColor);
imageline($this->img, $lastptX + 1, $lastptY - 1, $ptX + 1, $ptY - 1, $shadowColor);
imageline($this->img, $lastptX - 1, $lastptY, $ptX - 1, $ptY, $shadowColor);
// imageline($this->img,$lastptX,$lastptY,$ptX,$ptY,$shadowColor);
imageline($this->img, $lastptX + 1, $lastptY, $ptX + 1, $ptY, $shadowColor);
imageline($this->img, $lastptX - 1, $lastptY + 1, $ptX - 1, $ptY + 1, $shadowColor);
imageline($this->img, $lastptX, $lastptY + 1, $ptX, $ptY + 1, $shadowColor);
imageline($this->img, $lastptX + 1, $lastptY + 1, $ptX + 1, $ptY + 1, $shadowColor);
}
$lastPoint = new gpsPoint($line);
$points++;
}
}
DEBUG("MAP", 8, "<hr>");
$points = 0;
foreach ($lines as $line) {
// $line=trim($line);
if (strlen($line) == 0) {
continue;
}
if ($line[0] == 'B') {
$thisPoint = new gpsPoint($line);
if ($maxAlt - $minAlt) {
$altRelPos = ($thisPoint->getAlt() - $minAlt) / ($maxAlt - $minAlt);
} else {
$altRelPos = 0;
}
$thisColor = $gradientMap[floor(($gradient_colors_num - 1) * $altRelPos)];
list($ptX, $ptY) = $this->lonlat2xy(-$thisPoint->lon, $thisPoint->lat);
$ptY = $this->pixelHeight - $ptY;
if ($points == 0) {
$firstPoint = new gpsPoint($line);
} else {
list($lastptX, $lastptY) = $this->lonlat2xy(-$lastPoint->lon, $lastPoint->lat);
$lastptY = $this->pixelHeight - $lastptY;
imageline($this->img, $lastptX, $lastptY, $ptX, $ptY, $thisColor);
imageline($this->img, $lastptX, $lastptY + 1, $ptX, $ptY + 1, $thisColor);
// leave this for extra line width + shadow
// imageline($this->img,$lastptX+1,$lastptY,$ptX+1,$ptY,$shadowColor);
}
$lastPoint = new gpsPoint($line);
$points++;
}
}
list($ptX, $ptY) = $this->lonlat2xy(-$firstPoint->lon, $firstPoint->lat);
$ptY = $this->pixelHeight - $ptY;
imagearc($this->img, $ptX, $ptY, $startSize, $startSize, 0, 360, $startColor);
imagearc($this->img, $ptX, $ptY, $startSize - 1, $startSize - 1, 0, 360, $startColor);
list($ptX, $ptY) = $this->lonlat2xy(-$lastPoint->lon, $lastPoint->lat);
$ptY = $this->pixelHeight - $ptY;
imagearc($this->img, $ptX + 1, $ptY + 1, $startSize, $startSize, 0, 360, $shadowColor);
imagearc($this->img, $ptX + 1, $ptY + 1, $startSize - 1, $startSize - 1, 0, 360, $shadowColor);
imagearc($this->img, $ptX, $ptY, $startSize, $startSize, 0, 360, $endColor);
imagearc($this->img, $ptX, $ptY, $startSize - 1, $startSize - 1, 0, 360, $endColor);
}
