Method CalculateStaticEdgeCost

Filters/Modeling/vtkDijkstraImageGeodesicPath.cxx:155–182

------------------------------------------------------------------------------

//------------------------------------------------------------------------------
double vtkDijkstraImageGeodesicPath::CalculateStaticEdgeCost(
  vtkDataSet* inData, vtkIdType u, vtkIdType v)
{
  vtkImageData* image = vtkImageData::SafeDownCast(inData);

  double p1[3];
  image->GetPoint(u, p1);
  double p2[3];
  image->GetPoint(v, p2);

  double pcoords[3];
  int ijk1[3];
  int ijk2[3];
  image->ComputeStructuredCoordinates(p1, ijk1, pcoords);
  image->ComputeStructuredCoordinates(p2, ijk2, pcoords);

  double cost = this->ImageWeight *
    (image->GetScalarComponentAsDouble(ijk1[0], ijk1[1], ijk1[2], 0) +
      image->GetScalarComponentAsDouble(ijk1[0], ijk1[1], ijk1[2], 0));

  if (this->EdgeLengthWeight != 0.0)
  {
    cost +=
      this->EdgeLengthWeight * (sqrt(vtkMath::Distance2BetweenPoints(p1, p2)) / this->PixelSize);
  }

  return cost;
}

//------------------------------------------------------------------------------
double vtkDijkstraImageGeodesicPath::CalculateDynamicEdgeCost(