fly/terragear
1
0
Fork 0
Code Activity

The resulting airport surfaces still had too much undulation from the noise

Browse source 
in the DEM data.  Coarsen the nurbs grid and average some surrounding data
to hide the noise.
This commit is contained in:
curt 2003-03-09 19:19:19 +00:00
parent 846473a424
commit c986802faa
1 changed files with 16 additions and 14 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

30
src/Airports/GenAirports/apt_surface.cxx
View file

@ -150,20 +150,20 @@ TGAptSurface::TGAptSurface( const string& path,
cout << "Area size = " << x_m << " x " << y_m << " (m)" << endl; cout << "Area size = " << x_m << " x " << y_m << " (m)" << endl;
int xdivs = (int)(x_m / 400.0) + 1; int xdivs = (int)(x_m / 800.0) + 1;
int ydivs = (int)(y_m / 400.0) + 1; int ydivs = (int)(y_m / 800.0) + 1;
cout << " M(" << xdivs << "," << ydivs << ")" << endl; cout << " M(" << xdivs << "," << ydivs << ")" << endl;
double dlon = x_deg / xdivs; double dlon = x_deg / xdivs;
double dlat = y_deg / ydivs; double dlat = y_deg / ydivs;
// Build the double res input grid // Build the double res input grid
int mult = 4;
Matrix_Point3Df dPts(xdivs*2 + 2, ydivs*2 + 2); Matrix_Point3Df dPts((xdivs+1) * mult, (ydivs + 1) * mult);
for ( int i = 0; i <= xdivs*2; ++i ) { for ( int i = 0; i < dPts.rows(); ++i ) {
for ( int j = 0; j <= ydivs*2; ++j ) { for ( int j = 0; j < dPts.cols(); ++j ) {
dPts(i,j) = Point3Df( min_deg.lon() + i * (dlon*0.5), dPts(i,j) = Point3Df( min_deg.lon() + i * (dlon/(double)mult),
min_deg.lat() + j * (dlat*0.5), min_deg.lat() + j * (dlat/(double)mult),
-9999 ); -9999 );
} }
} }
@ -171,18 +171,20 @@ TGAptSurface::TGAptSurface( const string& path,
// Determine elevation of the grid points // Determine elevation of the grid points
calc_elevations( path, dPts ); calc_elevations( path, dPts );
// Build the norma res input grid from the double res version // Build the normal res input grid from the double res version
Matrix_Point3Df Pts(xdivs + 1, ydivs + 1); Matrix_Point3Df Pts(xdivs + 1, ydivs + 1);
for ( int i = 0; i <= xdivs; ++i ) { for ( int i = 0; i <= xdivs; ++i ) {
for ( int j = 0; j <= ydivs; ++j ) { for ( int j = 0; j <= ydivs; ++j ) {
cout << i << "," << j << endl; cout << i << "," << j << endl;
double z1 = dPts(2*i,2*j).z(); double accum = 0.0;
double z2 = dPts(2*i+1,2*j).z(); for ( int ii = 0; ii < mult; ++ii ) {
double z3 = dPts(2*i,2*j+1).z(); for ( int jj = 0; jj < mult; ++jj ) {
double z4 = dPts(i*2+1,2*j+1).z(); accum += dPts(mult * i + ii, mult * j + jj).z();
}
}
Pts(i,j) = Point3Df( min_deg.lon() + i * dlon, Pts(i,j) = Point3Df( min_deg.lon() + i * dlon,
min_deg.lat() + j * dlat, min_deg.lat() + j * dlat,
(z1+z2+z3+z4) / 4.0 ); accum / (mult*mult) );
} }
} }