fly/terragear
1
0
Fork 0
Code Activity

genapts: correctly calculate bucket center and radius for the airport btg files.

Browse source 
We already fixed this for construct in 2012 in commit eaec9905f7
This commit is contained in:
Christian Schmitt 2013-11-14 20:54:54 +01:00
parent c55eaeeb1d
commit 3191c6464d
1 changed files with 8 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

13
src/Airports/GenAirports850/airport.cxx
View file

@ -1111,15 +1111,18 @@ void Airport::BuildBtg(const std::string& root, const string_list& elev_src )
wgs84_nodes.push_back( SGVec3d::fromGeod(geod_nodes[i]) );
}
SGSphered d;
SGVec3d gbs_center = SGVec3d::fromGeod( b.get_center() );
double dist_squared, radius_squared = 0;
for ( unsigned int i = 0; i < wgs84_nodes.size(); ++i )
{
d.expandBy(wgs84_nodes[ i ]);
dist_squared = distSqr(gbs_center, wgs84_nodes[i]);
if ( dist_squared > radius_squared ) {
radius_squared = dist_squared;
}
}
SGVec3d gbs_center = d.getCenter();
double gbs_radius = d.getRadius();
TG_LOG(SG_GENERAL, SG_DEBUG, "gbs center = " << gbs_center);
double gbs_radius = sqrt(radius_squared);
TG_LOG(SG_GENERAL, SG_DEBUG, "Done with wgs84 node mapping");
TG_LOG(SG_GENERAL, SG_DEBUG, " center = " << gbs_center << " radius = " << gbs_radius );