bool
mitk::PlanarFigureInteractor::IsMousePositionAcceptableAsNewControlPoint(
const mitk::InteractionPositionEvent* positionEvent,
const PlanarFigure* planarFigure )
{
assert(positionEvent && planarFigure);
BaseRenderer* renderer = positionEvent->GetSender();
assert(renderer);
// Get the timestep to support 3D+t
int timeStep( renderer->GetTimeStep( planarFigure ) );
bool tooClose(false);
const PlaneGeometry *renderingPlane = renderer->GetCurrentWorldPlaneGeometry();
mitk::PlaneGeometry *planarFigureGeometry =
dynamic_cast< mitk::PlaneGeometry * >( planarFigure->GetGeometry( timeStep ) );
mitk::AbstractTransformGeometry *abstractTransformGeometry =
dynamic_cast< mitk::AbstractTransformGeometry * >( planarFigure->GetGeometry( timeStep ) );
if ( abstractTransformGeometry != NULL )
return false;
Point2D point2D, correctedPoint;
// Get the point2D from the positionEvent
if ( !this->TransformPositionEventToPoint2D( positionEvent, planarFigureGeometry, point2D ) )
{
return false;
}
// apply the controlPoint constraints of the planarFigure to get the
// coordinates that would actually be used.
correctedPoint = const_cast<PlanarFigure*>( planarFigure )->ApplyControlPointConstraints( 0, point2D );
// map the 2D coordinates of the new point to world-coordinates
// and transform those to display-coordinates
mitk::Point3D newPoint3D;
planarFigureGeometry->Map( correctedPoint, newPoint3D );
mitk::Point2D newDisplayPosition;
renderingPlane->Map( newPoint3D, newDisplayPosition );
renderer->GetDisplayGeometry()->WorldToDisplay( newDisplayPosition, newDisplayPosition );
for( int i=0; i < (int)planarFigure->GetNumberOfControlPoints(); i++ )
{
if ( i != planarFigure->GetSelectedControlPoint() )
{
// Try to convert previous point to current display coordinates
mitk::Point3D previousPoint3D;
// map the 2D coordinates of the control-point to world-coordinates
planarFigureGeometry->Map( planarFigure->GetControlPoint( i ), previousPoint3D );
if ( renderer->GetDisplayGeometry()->Distance( previousPoint3D ) < 0.1 ) // ugly, but assert makes this work
{
mitk::Point2D previousDisplayPosition;
// transform the world-coordinates into display-coordinates
renderingPlane->Map( previousPoint3D, previousDisplayPosition );
renderer->GetDisplayGeometry()->WorldToDisplay( previousDisplayPosition, previousDisplayPosition );
//Calculate the distance. We use display-coordinates here to make
// the check independent of the zoom-level of the rendering scene.
double a = newDisplayPosition[0] - previousDisplayPosition[0];
double b = newDisplayPosition[1] - previousDisplayPosition[1];
// If point is to close, do not set a new point
tooClose = (a * a + b * b < m_MinimumPointDistance );
}
if ( tooClose )
return false; // abort loop early
}
}
return !tooClose; // default
}
bool SlicesRotator::DoDecideBetweenRotationAndSliceSelection(Action*, const StateEvent* e)
{
// Decide between moving and rotation slices.
// For basic decision logic see class documentation.
/*
Detail logic:
1. Find the SliceNavigationController that has sent the event: this one defines our rendering plane and will NOT be rotated. Needs not even be counted or checked.
2. Inspect every other SliceNavigationController
- calculate the line intersection of this SliceNavigationController's plane with our rendering plane
- if there is NO interesection, ignore and continue
- IF there is an intersection
- check the mouse cursor's distance from that line.
0. if the line is NOT near the cursor, remember the plane as "one of the other planes" (which can be rotated in "locked" mode)
1. on first line near the cursor, just remember this intersection line as THE other plane that we want to rotate
2. on every consecutive line near the cursor, check if the line is geometrically identical to the line that we want to rotate
- if yes, we just push this line to the "other" lines and rotate it along
- if no, then we have a situation where the mouse is near two other lines (e.g. crossing point) and don't want to rotate
*/
const DisplayPositionEvent* posEvent = dynamic_cast<const DisplayPositionEvent*>(e->GetEvent());
if (!posEvent) return false;
BaseRenderer* clickedRenderer = e->GetEvent()->GetSender();
const PlaneGeometry* ourViewportGeometry = dynamic_cast<const PlaneGeometry*>( clickedRenderer->GetCurrentWorldPlaneGeometry() );
// These sanity checks were introduced with bug 17877, since plane geometries are now a shared base class of several geometries
// They may ultimately be unecessary
const mitk::AbstractTransformGeometry* abstractGeometry = dynamic_cast< const AbstractTransformGeometry * > (ourViewportGeometry);
if (abstractGeometry != NULL) MITK_WARN << "SliceRotator recieved an AbstractTransformGeometry, expecting a simple PlainGeometry, behaviour should be verified.";
const mitk::DisplayGeometry* displayGeometry = dynamic_cast< const DisplayGeometry * > (ourViewportGeometry);
if (displayGeometry != NULL) MITK_WARN << "SliceRotator recieved a DisplayGeometry, expecting a simple PlainGeometry, behaviour should be verified.";
// End sanity checks
if (!ourViewportGeometry) return false;
DisplayGeometry* clickedDisplayGeometry = clickedRenderer->GetDisplayGeometry();
if (!clickedDisplayGeometry) return false;
MITK_DEBUG << "=============================================";
MITK_DEBUG << "Renderer under cursor is " << clickedRenderer->GetName();
Point3D cursorPosition = posEvent->GetWorldPosition();
const PlaneGeometry* geometryToBeRotated = NULL; // this one is under the mouse cursor
const PlaneGeometry* anyOtherGeometry = NULL; // this is also visible (for calculation of intersection ONLY)
Line3D intersectionLineWithGeometryToBeRotated;
bool hitMultipleLines(false);
m_SNCsToBeRotated.clear();
const double threshholdDistancePixels = 12.0;
for (SNCVector::iterator iter = m_RotatableSNCs.begin(); iter != m_RotatableSNCs.end(); ++iter)
{
// If the mouse cursor is in 3D Renderwindow, do not check for intersecting planes.
if (clickedRenderer->GetMapperID() == BaseRenderer::Standard3D)
break;
const PlaneGeometry* otherRenderersRenderPlane = (*iter)->GetCurrentPlaneGeometry();
if (otherRenderersRenderPlane == NULL) continue; // ignore, we don't see a plane
MITK_DEBUG << " Checking plane of renderer " << (*iter)->GetRenderer()->GetName();
// check if there is an intersection
Line3D intersectionLine; // between rendered/clicked geometry and the one being analyzed
if (!ourViewportGeometry->IntersectionLine( otherRenderersRenderPlane, intersectionLine ))
{
continue; // we ignore this plane, it's parallel to our plane
}
// check distance from intersection line
double distanceFromIntersectionLine = intersectionLine.Distance( cursorPosition );
ScalarType distancePixels = distanceFromIntersectionLine / clickedDisplayGeometry->GetScaleFactorMMPerDisplayUnit();
MITK_DEBUG << " Distance of plane from cursor " << distanceFromIntersectionLine << " mm, which is around " << distancePixels << " px" ;
// far away line, only remember for linked rotation if necessary
if (distanceFromIntersectionLine > threshholdDistancePixels)
{
MITK_DEBUG << " Plane is too far away --> remember as otherRenderersRenderPlane";
anyOtherGeometry = otherRenderersRenderPlane; // we just take the last one, so overwrite each iteration (we just need some crossing point)
// TODO what about multiple crossings? NOW we have undefined behavior / random crossing point is used
if (m_LinkPlanes)
{
m_SNCsToBeRotated.push_back(*iter);
}
}
else // close to cursor
{
MITK_DEBUG << " Plane is close enough to cursor...";
if ( geometryToBeRotated == NULL ) // first one close to the cursor
{
MITK_DEBUG << " It is the first close enough geometry, remember as geometryToBeRotated";
geometryToBeRotated = otherRenderersRenderPlane;
intersectionLineWithGeometryToBeRotated = intersectionLine;
m_SNCsToBeRotated.push_back(*iter);
}
else
{
MITK_DEBUG << " Second or later close enough geometry";
// compare to the line defined by geometryToBeRotated: if identical, just rotate this otherRenderersRenderPlane together with the primary one
// if different, DON'T rotate
if ( intersectionLine.IsParallel( intersectionLineWithGeometryToBeRotated )
&& intersectionLine.Distance( intersectionLineWithGeometryToBeRotated.GetPoint1() ) < mitk::eps )
{
MITK_DEBUG << " This line is the same as intersectionLineWithGeometryToBeRotated which we already know";
m_SNCsToBeRotated.push_back(*iter);
}
else
{
MITK_DEBUG << " This line is NOT the same as intersectionLineWithGeometryToBeRotated which we already know";
hitMultipleLines = true;
}
}
}
}
bool moveSlices(true);
if ( geometryToBeRotated && anyOtherGeometry && ourViewportGeometry && !hitMultipleLines )
{
// assure all three are valid, so calculation of center of rotation can be done
moveSlices = false;
}
MITK_DEBUG << "geometryToBeRotated: " << (void*)geometryToBeRotated;
MITK_DEBUG << "anyOtherGeometry: " << (void*)anyOtherGeometry;
MITK_DEBUG << "ourViewportGeometry: " << (void*)ourViewportGeometry;
MITK_DEBUG << "hitMultipleLines? " << hitMultipleLines;
MITK_DEBUG << "moveSlices? " << moveSlices;
std::auto_ptr<StateEvent> decidedEvent;
// question in state machine is: "rotate?"
if (moveSlices) // i.e. NOT rotate
{
// move all planes to posEvent->GetWorldPosition()
decidedEvent.reset( new StateEvent(EIDNO, e->GetEvent()) );
MITK_DEBUG << "Rotation not possible, not enough information (other planes crossing rendering plane) ";
}
else
{ // we DO have enough information for rotation
m_LastCursorPosition = intersectionLineWithGeometryToBeRotated.Project(cursorPosition); // remember where the last cursor position ON THE LINE has been observed
if (anyOtherGeometry->IntersectionPoint(intersectionLineWithGeometryToBeRotated, m_CenterOfRotation)) // find center of rotation by intersection with any of the OTHER lines
{
decidedEvent.reset( new StateEvent(EIDYES, e->GetEvent()) );
MITK_DEBUG << "Rotation possible";
}
else
{
MITK_DEBUG << "Rotation not possible, cannot determine the center of rotation!?";
decidedEvent.reset( new StateEvent(EIDNO, e->GetEvent()) );
}
}
this->HandleEvent( decidedEvent.get() );
return true;
}