void qSlicerArielImplantPlanModuleWidgetPrivate::updateStemCutPlane() { Q_Q(qSlicerArielImplantPlanModuleWidget); if (curStemInfo.brand.empty() || curStemInfo.type.empty()) return; vtkNew<vtkMatrix4x4> stemPoseMatrix; stemPoseNode->GetMatrixTransformToParent(stemPoseMatrix); // 活动评估下 if (currentMode == qSlicerArielImplantPlanModuleWidgetPrivate::CollisionDetect) if (::isPoseValid(this->stemPoseMatrixRaw)) { stemPoseMatrix->DeepCopy(this->stemPoseMatrixRaw); } vtkNew<vtkMatrix4x4> identity; if (areMatricesEqual(stemPoseMatrix, identity)) return; double stemInfo[4]; stemInfo[0] = curStemInfo.neckLength; stemInfo[1] = curStemInfo.horizontalOffset; stemInfo[2] = curStemInfo.neckShaftAngle; stemInfo[3] = curStemInfo.cutAngle; // transform double origin[3], normal[3]; currentModuleLogic->getCutFemurPlaneOriginAndNormal(stemInfo, origin, normal); vtkNew<vtkPlane> stemCutPlane; currentModuleLogic->getTransformedPlane(origin, normal, stemPoseMatrix, stemCutPlane); stemCutPlane->GetOrigin(origin); stemCutPlane->GetNormal(normal); // Plane if (stemCutPlane && this->stemCutPlane) { std::cout << __FUNCTION__ << std::endl; //股骨柄切面高度 auto stemHeight = q->calculateCutFemurHeight(origin, normal); std::cout << u8"股骨柄切面高度:height=" << stemHeight << std::endl; //当前截骨面高度 auto distance = q->calculateCutFemurHeight(); std::cout << u8"矫正前截骨面高度:distance=" << distance << std::endl; //偏移=实际-理论 std::cout << u8"截骨偏移: cutPlaneOffset=" << cutPlaneOffset << std::endl; //高度or法向量 不相等:调节 if ((stemHeight + cutPlaneOffset) != distance) { //重置PlaneNode double newOrigin[3] = {0}; double newNormal[3] = {0}; memcpy(newNormal, normal, sizeof(double) * 3); if (cutPlaneOffset == 0.0) memcpy(newOrigin, origin, sizeof(double) * 3); else currentModuleLogic->calculatePositionOfCutPlaneByOffset(origin, normal, cutPlaneOffset, newOrigin); this->stemCutPlane->setPlaneNodeByPlaneType( PlaneType( {newNormal[0], newNormal[1], newNormal[2], newOrigin[0], newOrigin[1], newOrigin[2]}), this->cutPlaneRadius); //重置高度 q->resetCutFemurHeightSpinBoxValue(); //重新获取offset distance = q->calculateCutFemurHeight(); std::cout << u8"矫正后截骨面高度:distance=" << distance << ", stemHeight=" << stemHeight << std::endl; } } }
时间: 2023-02-09 09:52:46 浏览: 174
mpi_type.rar_If...
这段代码看起来是在更新一个叫做"stemCutPlane"的切割平面。首先,它检查当前模块是否处于"CollisionDetect"模式,如果是,它会把"stemPoseMatrix"的值设置为"stemPoseMatrixRaw"。然后,它使用当前的"stemInfo"数组来获取平面的原点和法向量,并使用"stemPoseMatrix"对平面进行变换。最后,它获取变换后平面的原点和法向量。
阅读全文