bool
VerticalDatum::transform(const VerticalDatum* from,
const VerticalDatum* to,
const GeoExtent& extent,
osg::HeightField* hf )
{
if ( from == to )
return true;
unsigned cols = hf->getNumColumns();
unsigned rows = hf->getNumRows();
osg::Vec3d sw(extent.west(), extent.south(), 0.0);
osg::Vec3d ne(extent.east(), extent.north(), 0.0);
double xstep = abs(extent.east() - extent.west()) / double(cols-1);
double ystep = abs(extent.north() - extent.south()) / double(rows-1);
if ( !extent.getSRS()->isGeographic() )
{
const SpatialReference* geoSRS = extent.getSRS()->getGeographicSRS();
extent.getSRS()->transform(sw, geoSRS, sw);
extent.getSRS()->transform(ne, geoSRS, ne);
xstep = (ne.x()-sw.x()) / double(cols-1);
ystep = (ne.y()-sw.y()) / double(rows-1);
}
for( unsigned c=0; c<cols; ++c)
{
double lon = sw.x() + xstep*double(c);
for( unsigned r=0; r<rows; ++r)
{
double lat = sw.y() + ystep*double(r);
float& h = hf->getHeight(c, r);
if (h != NO_DATA_VALUE)
{
VerticalDatum::transform( from, to, lat, lon, h );
}
}
}
return true;
}
osg::Node*
SingleKeyNodeFactory::createTile(TileModel* model,
bool setupChildrenIfNecessary,
ProgressCallback* progress)
{
#ifdef EXPERIMENTAL_TILE_NODE_CACHE
osg::ref_ptr<TileNode> tileNode;
TileNodeCache::Record rec;
cache.get(model->_tileKey, rec);
if ( rec.valid() )
{
tileNode = rec.value().get();
}
else
{
tileNode = _modelCompiler->compile( model, _frame );
cache.insert(model->_tileKey, tileNode);
}
#else
// compile the model into a node:
TileNode* tileNode = _modelCompiler->compile(model, _frame, progress);
tileNode->setEngineUID( _engineUID );
#endif
// see if this tile might have children.
bool prepareForChildren =
setupChildrenIfNecessary &&
model->_tileKey.getLOD() < *_options.maxLOD();
osg::Node* result = 0L;
if ( prepareForChildren )
{
osg::BoundingSphere bs = tileNode->getBound();
TilePagedLOD* plod = new TilePagedLOD( _engineUID, _liveTiles, _deadTiles );
plod->setCenter ( bs.center() );
plod->addChild ( tileNode );
plod->setFileName( 1, Stringify() << tileNode->getKey().str() << "." << _engineUID << ".osgearth_engine_mp_tile" );
plod->setDebug ( _debug );
if ( _options.rangeMode().value() == osg::LOD::DISTANCE_FROM_EYE_POINT )
{
//Compute the min range based on the 2D size of the tile
GeoExtent extent = model->_tileKey.getExtent();
double radius = 0.0;
#if 0
// Test code to use the equitorial radius so that all of the tiles at the same level
// have the same range. This will make the poles page in more appropriately.
if (_frame.getMapInfo().isGeocentric())
{
GeoExtent equatorialExtent(
extent.getSRS(),
extent.west(),
-extent.height()/2.0,
extent.east(),
extent.height()/2.0 );
radius = equatorialExtent.getBoundingGeoCircle().getRadius();
}
else
#endif
{
GeoPoint lowerLeft(extent.getSRS(), extent.xMin(), extent.yMin(), 0.0, ALTMODE_ABSOLUTE);
GeoPoint upperRight(extent.getSRS(), extent.xMax(), extent.yMax(), 0.0, ALTMODE_ABSOLUTE);
osg::Vec3d ll, ur;
lowerLeft.toWorld( ll );
upperRight.toWorld( ur );
radius = (ur - ll).length() / 2.0;
}
float minRange = (float)(radius * _options.minTileRangeFactor().value());
plod->setRange( 0, minRange, FLT_MAX );
plod->setRange( 1, 0, minRange );
plod->setRangeMode( osg::LOD::DISTANCE_FROM_EYE_POINT );
}
else
{
// the *2 is because we page in 4-tile sets, not individual tiles.
float size = 2.0f * _options.tilePixelSize().value();
plod->setRange( 0, 0.0f, size );
plod->setRange( 1, size, FLT_MAX );
plod->setRangeMode( osg::LOD::PIXEL_SIZE_ON_SCREEN );
}
// Install a tile-aligned bounding box in the pager node itself so we can do
// visibility testing before paging in subtiles.
plod->setChildBoundingBoxAndMatrix(
1,
tileNode->getTerrainBoundingBox(),
tileNode->getMatrix() );
// DBPager will set a priority based on the ratio range/maxRange.
// This will offset that number with a full LOD #, giving LOD precedence.
// Experimental.
//plod->setPriorityScale( 1, model->_tileKey.getLOD()+1 );
#if USE_FILELOCATIONCALLBACK
osgDB::Options* options = plod->getOrCreateDBOptions();
options->setFileLocationCallback( new FileLocationCallback() );
#endif
result = plod;
// this one rejects back-facing tiles:
if ( _frame.getMapInfo().isGeocentric() && _options.clusterCulling() == true )
{
#if 1
osg::HeightField* hf =
model->_elevationData.getHeightField();
result->addCullCallback( HeightFieldUtils::createClusterCullingCallback(
hf,
tileNode->getKey().getProfile()->getSRS()->getEllipsoid(),
*_options.verticalScale() ) );
#else
// This works, but isn't quite as tight at the cluster culler.
// Re-evaluate down the road.
result->addCullCallback( new HorizonTileCuller(
_frame.getMapInfo().getSRS(),
tileNode->getTerrainBoundingBox(),
tileNode->getMatrix(), tileNode->getKey() ) );
#endif
}
}
else
{
result = tileNode;
}
return result;
}
