void mitk::BaseGeometry::ChangeImageGeometryConsideringOriginOffset(const bool isAnImageGeometry)
{
// If Geometry is switched to ImageGeometry, you have to put an offset to the origin, because
// imageGeometries origins are pixel-center-based
// ... and remove the offset, if you switch an imageGeometry back to a normal geometry
// For more information please see the Geometry documentation page
if (m_ImageGeometry == isAnImageGeometry)
return;
const BoundingBox::BoundsArrayType &boundsarray = this->GetBoundingBox()->GetBounds();
Point3D originIndex;
FillVector3D(originIndex, boundsarray[0], boundsarray[2], boundsarray[4]);
if (isAnImageGeometry == true)
FillVector3D(originIndex, originIndex[0] + 0.5, originIndex[1] + 0.5, originIndex[2] + 0.5);
else
FillVector3D(originIndex, originIndex[0] - 0.5, originIndex[1] - 0.5, originIndex[2] - 0.5);
Point3D originWorld;
originWorld = GetIndexToWorldTransform()->TransformPoint(originIndex);
// instead could as well call IndexToWorld(originIndex,originWorld);
SetOrigin(originWorld);
this->SetImageGeometry(isAnImageGeometry);
}
bool
PlaneGeometry::Project(const mitk::Vector3D &vec3d_mm, mitk::Vector3D &projectedVec3d_mm) const
{
assert(this->IsBoundingBoxNull()==false);
Vector3D vec3d_units;
Superclass::WorldToIndex(vec3d_mm, vec3d_units);
vec3d_units[2] = 0;
projectedVec3d_mm = GetIndexToWorldTransform()->TransformVector(vec3d_units);
return true;
}
bool
PlaneGeometry::Project(
const mitk::Point3D &pt3d_mm, mitk::Point3D &projectedPt3d_mm) const
{
assert(this->IsBoundingBoxNull()==false);
Point3D pt3d_units;
Superclass::WorldToIndex(pt3d_mm, pt3d_units);
pt3d_units[2] = 0;
projectedPt3d_mm = GetIndexToWorldTransform()->TransformPoint(pt3d_units);
return const_cast<BoundingBox*>(this->GetBoundingBox())->IsInside(pt3d_units);
}
void
PlaneGeometry::Map(const mitk::Point2D &pt2d_mm, mitk::Point3D &pt3d_mm) const
{
//pt2d_mm is measured from the origin of the world geometry (at leats it called form BaseRendere::Mouse...Event)
Point3D pt3d_units;
pt3d_units[0] = pt2d_mm[0] / (GetExtentInMM(0) / GetExtent(0));
pt3d_units[1] = pt2d_mm[1] / (GetExtentInMM(1) / GetExtent(1));
pt3d_units[2]=0;
//pt3d_units is a continuos index. We divided it with the Scale Factor (= spacing in x and y) to convert it from mm to index units.
//
pt3d_mm = GetIndexToWorldTransform()->TransformPoint(pt3d_units);
//now we convert the 3d index to a 3D world point in mm. We could have used IndexToWorld as well as GetITW->Transform...
}
void mitk::GizmoInteractor::ApplyTranslationToManipulatedObject(const Vector3D& translation)
{
assert(m_ManipulatedObjectGeometry.IsNotNull());
auto manipulatedGeometry = m_InitialManipulatedObjectGeometry->Clone();
m_FinalDoOperation.reset(new PointOperation(OpMOVE, translation));
if (m_UndoEnabled)
{
m_FinalUndoOperation.reset(new PointOperation(OpMOVE, -translation));
}
manipulatedGeometry->ExecuteOperation(m_FinalDoOperation.get());
m_ManipulatedObjectGeometry->SetIdentity();
m_ManipulatedObjectGeometry->Compose( manipulatedGeometry->GetIndexToWorldTransform() );
}
bool
PlaneGeometry::Project(const mitk::Point3D & atPt3d_mm,
const mitk::Vector3D &vec3d_mm, mitk::Vector3D &projectedVec3d_mm) const
{
MITK_WARN << "Deprecated function! Call Project(vec3D,vec3D) instead.";
assert(this->IsBoundingBoxNull()==false);
Vector3D vec3d_units;
Superclass::WorldToIndex(atPt3d_mm, vec3d_mm, vec3d_units);
vec3d_units[2] = 0;
projectedVec3d_mm = GetIndexToWorldTransform()->TransformVector(vec3d_units);
Point3D pt3d_units;
Superclass::WorldToIndex(atPt3d_mm, pt3d_units);
return const_cast<BoundingBox*>(this->GetBoundingBox())->IsInside(pt3d_units);
}
void mitk::GizmoInteractor::ApplyEqualScalingToManipulatedObject(double scalingFactor)
{
assert(m_ManipulatedObjectGeometry.IsNotNull());
auto manipulatedGeometry = m_InitialManipulatedObjectGeometry->Clone();
m_FinalDoOperation.reset(new ScaleOperation(OpSCALE,
scalingFactor - 1.0,
m_InitialGizmoCenter3D));
if (m_UndoEnabled)
{
m_FinalUndoOperation.reset(new ScaleOperation(OpSCALE,
-(scalingFactor - 1.0),
m_InitialGizmoCenter3D));
}
manipulatedGeometry->ExecuteOperation(m_FinalDoOperation.get());
m_ManipulatedObjectGeometry->SetIdentity();
m_ManipulatedObjectGeometry->Compose( manipulatedGeometry->GetIndexToWorldTransform() );
}
void mitk::GizmoInteractor::ApplyRotationToManipulatedObject(double angle_deg)
{
assert(m_ManipulatedObjectGeometry.IsNotNull());
auto manipulatedGeometry = m_InitialManipulatedObjectGeometry->Clone();
m_FinalDoOperation.reset(new RotationOperation(OpROTATE,
m_InitialGizmoCenter3D,
m_AxisOfRotation,
angle_deg));
if (m_UndoEnabled)
{
m_FinalUndoOperation.reset(new RotationOperation(OpROTATE,
m_InitialGizmoCenter3D,
m_AxisOfRotation,
-angle_deg));
}
manipulatedGeometry->ExecuteOperation(m_FinalDoOperation.get());
m_ManipulatedObjectGeometry->SetIdentity();
m_ManipulatedObjectGeometry->Compose( manipulatedGeometry->GetIndexToWorldTransform() );
}