void GeoCoordinateTransformer::transformRecursive_helper(CityObject& obj, GeoTransform& transformation) {
if (obj.getEnvelope().validBounds()) {
if (transformation.valid()) {
TVec3d lowerBound = obj.getEnvelope().getLowerBound();
TVec3d upperBound = obj.getEnvelope().getUpperBound();
transformation.transform(lowerBound);
transformation.transform(upperBound);
Envelope* newEnvelope = new Envelope(m_destinationSRS);
newEnvelope->setLowerBound(lowerBound);
newEnvelope->setUpperBound(upperBound);
obj.setEnvelope(newEnvelope);
} else {
CITYGML_LOG_WARN(m_logger, "No valid spatial reference system is given for CityObject with id '" << obj.getId() << "'. Envelope (Bounding Box) is not transformed.");
}
}
for (unsigned int i = 0; i < obj.getChildCityObjectsCount(); i++) {
transformRecursive(obj.getChildCityObject(i), transformation);
}
for (unsigned int i = 0; i < obj.getImplicitGeometryCount(); i++) {
transformRecursive(obj.getImplicitGeometry(i), transformation);
}
for (unsigned int i = 0; i < obj.getGeometriesCount(); i++) {
transform(obj.getGeometry(i), transformation);
}
}
void GeoCoordinateTransformer::transform(Geometry& obj, GeoTransform& transformation) {
if (!transformation.valid()) {
CITYGML_LOG_WARN(m_logger, "No valid spatial reference system is given for Geometry with id '" << obj.getId() << "'. Child Polygons are not transformed"
<< "(unless they are shared with another geometry for which a spatial reference system is defined)");
return;
}
for (unsigned int i = 0; i < obj.getPolygonsCount(); i++) {
const auto poly = obj.getPolygon(i);
auto it = m_transformedPolygonsSourceURNMap.find(poly.get());
if (it == m_transformedPolygonsSourceURNMap.end()) {
for (TVec3d& vertex : poly->getVertices()) {
transformation.transform(vertex);
}
m_transformedPolygonsSourceURNMap[poly.get()] = transformation.sourceURN();
} else if (it->second != transformation.sourceURN()) {
CITYGML_LOG_WARN(m_logger, "Polygon with id '" << poly->getId() << "' was already transformed from " << it->second << " to " << m_destinationSRS
<< ". But the spatial reference system of Geometry object with id '" << obj.getId() << "' that also contains the polygon is different "
<< "(" << transformation.sourceURN() << "). Ignoring new source SRS.");
}
}
for (unsigned int i = 0; i < obj.getLineStringCount(); i++) {
const auto lineString = obj.getLineString(i);
auto it = m_transformedLineStringsSourceURNMap.find(lineString.get());
if (it == m_transformedLineStringsSourceURNMap.end()) {
if (lineString->getDimensions() == 2) {
for (TVec2d& vertex : lineString->getVertices2D()) {
transformation.transform(vertex);
}
} else if (lineString->getDimensions() == 3) {
for (TVec3d& vertex : lineString->getVertices3D()) {
transformation.transform(vertex);
}
}
m_transformedLineStringsSourceURNMap[lineString.get()] = transformation.sourceURN();
} else if (it->second != transformation.sourceURN()) {
CITYGML_LOG_WARN(m_logger, "LineString with id '" << lineString->getId() << "' was already transformed from " << it->second << " to " << m_destinationSRS
<< ". But the spatial reference system of Geometry object with id '" << obj.getId() << "' that also contains the LineString is different "
<< "(" << transformation.sourceURN() << "). Ignoring new source SRS.");
}
}
for (unsigned int i = 0; i < obj.getGeometriesCount(); i++) {
transform(obj.getGeometry(i), transformation);
}
}
void update( float x, float y, osgViewer::View* view )
{
bool yes = false;
// look under the mouse:
osg::Vec3d world;
osgUtil::LineSegmentIntersector::Intersections hits;
if ( view->computeIntersections(x, y, hits) )
{
world = hits.begin()->getWorldIntersectPoint();
// convert to map coords:
GeoPoint mapPoint;
mapPoint.fromWorld( s_mapNode->getMapSRS(), world );
// Depending on the level of detail key you request, you will get a mesh that should line up exactly with the highest resolution mesh that the terrain engine will draw.
// At level 15 that is a 257x257 heightfield. If you select a higher lod, the mesh will be less dense.
TileKey key = s_mapNode->getMap()->getProfile()->createTileKey(mapPoint.x(), mapPoint.y(), 15);
OE_NOTICE << "Creating tile " << key.str() << std::endl;
osg::ref_ptr<osg::Node> node = s_mapNode->getTerrainEngine()->createTile(key);
if (node.valid())
{
// Extract the triangles from the node that was created and do our own rendering. Simulates what you would do when passing in the triangles to a physics engine.
OE_NOTICE << "Created tile for " << key.str() << std::endl;
CollectTrianglesVisitor v;
node->accept(v);
node = v.buildNode();
if (_node.valid())
{
s_root->removeChild( _node.get() );
}
osg::Group* group = new osg::Group;
// Show the actual mesh.
group->addChild( node.get() );
_node = group;
// Clamp the marker to the intersection of the triangles created by osgEarth. This should line up with the mesh that is actually rendered.
double z = 0.0;
s_mapNode->getTerrain()->getHeight( node, s_mapNode->getMapSRS(), mapPoint.x(), mapPoint.y(), &z);
GeoTransform* xform = new GeoTransform();
xform->setPosition( osgEarth::GeoPoint(s_mapNode->getMapSRS(),mapPoint.x(), mapPoint.y(), z, ALTMODE_ABSOLUTE) );
xform->addChild( marker.get() );
group->addChild( xform );
s_root->addChild( _node.get() );
}
else
{
OE_NOTICE << "Failed to create tile for " << key.str() << std::endl;
}
}
}
osg::Node*
installGeometry2(const SpatialReference* srs)
{
osg::Geode* geode = new osg::Geode();
geode->setComputeBoundingSphereCallback( new MyComputeBoundCallback(RADIUS2) );
geode->addDrawable( new osg::ShapeDrawable( new osg::Sphere(osg::Vec3f(0,0,0), RADIUS2) ) );
osg::Vec3f center = geode->getBound().center();
GeoTransform* xform = new GeoTransform();
xform->setPosition( GeoPoint(srs, 180.0, 0.0, 0.0, ALTMODE_ABSOLUTE) );
xform->addChild( geode );
return xform;
}
void GeoCoordinateTransformer::transformRecursive(ImplicitGeometry& obj, GeoTransform& transformation) {
if (!obj.getSRSName().empty() && !transformation.hasSourceSRS(obj.getSRSName())) {
GeoTransform childTransformation = transformation;
childTransformation.setSourceSRS(obj.getSRSName());
transformRecursive_helper(obj, childTransformation);
} else {
transformRecursive_helper(obj, transformation);
}
}
void GeoCoordinateTransformer::transformRecursive(CityObject& obj, GeoTransform& transformation) {
if (!obj.getEnvelope().srsName().empty() && !transformation.hasSourceSRS(obj.getEnvelope().srsName())) {
GeoTransform childTransformation = transformation;
childTransformation.setSourceSRS(obj.getEnvelope().srsName());
transformRecursive_helper(obj, childTransformation);
} else {
transformRecursive_helper(obj, transformation);
}
}
void GeoCoordinateTransformer::transformRecursive_helper(ImplicitGeometry& obj, GeoTransform& transformation) {
TVec3d referencePoint = obj.getReferencePoint();
if (transformation.valid()) {
transformation.transform(referencePoint);
obj.setReferencePoint(referencePoint);
obj.setSRSName(m_destinationSRS);
} else {
CITYGML_LOG_WARN(m_logger, "No valid spatial reference system is given for ImplicitGeometry with id '" << obj.getId() << "'. Reference Point is not transformed.");
}
// Do not transform the geometry of an ImplicitGeometry object. Implicit Geometries share Geometry objects but each implicit geometry
// defines its own transformation on the vertices of the shared geometry. Hence those vertices are in a local coordinate system. Read
// the citygml documentation for more details.
//for (unsigned int i = 0; i < obj.getGeometriesCount(); i++) {
// transform(obj.getGeometry(i), transformation);
//}
}
void Dataset::createGeoTransform ( GeoTransform &transfrom, const Extents& extents, unsigned int width, unsigned int height )
{
// Get the extents lower left and upper right.
Extents::Vertex ll ( extents.minimum() );
Extents::Vertex ur ( extents.maximum() );
// Get the length in x and y.
const double xLength ( ur[0] - ll[0] );
const double yLength ( ur[1] - ll[1] );
// Figure out the pixel resolution.
const double xResolution ( xLength / width );
const double yResolution ( yLength / height );
// Make sure there is enough room.
transfrom.resize ( 6 );
transfrom[0] = ll[0]; // top left x
transfrom[1] = xResolution; // w-e pixel resolution
transfrom[2] = 0; // rotation, 0 if image is "north up"
transfrom[3] = ur[1]; // top left y
transfrom[4] = 0; // rotation, 0 if image is "north up"
transfrom[5] = -yResolution; // n-s pixel resolution
}
