Method RequestData

Filters/Geometry/vtkStructuredGridGeometryFilter.cxx:29–411

}

int vtkStructuredGridGeometryFilter::RequestData(vtkInformation* vtkNotUsed(request),
  vtkInformationVector** inputVector, vtkInformationVector* outputVector)
{
  // get the info objects
  vtkInformation* inInfo = inputVector[0]->GetInformationObject(0);
  vtkInformation* outInfo = outputVector->GetInformationObject(0);

  // get the input and output
  vtkStructuredGrid* input =
    vtkStructuredGrid::SafeDownCast(inInfo->Get(vtkDataObject::DATA_OBJECT()));
  vtkPolyData* output = vtkPolyData::SafeDownCast(outInfo->Get(vtkDataObject::DATA_OBJECT()));

  int dimension, dir[3], diff[3];
  int i, j, k, extent[6], *inExt;
  vtkIdType ptIds[4], idx, startIdx, startCellIdx, cellId;
  vtkPoints* newPts = nullptr;
  vtkCellArray* newVerts = nullptr;
  vtkCellArray* newLines = nullptr;
  vtkCellArray* newPolys = nullptr;
  vtkIdType totPoints, pos, cellPos;
  int offset[3], cellOffset[3], numPolys;
  double x[3];
  vtkPointData *pd, *outPD;
  vtkCellData *cd, *outCD;

  vtkDebugMacro(<< "Extracting structured points geometry");

  if (input->GetPoints() == nullptr)
  {
    vtkDebugMacro(<< "No data to extract");
    return 1;
  }

  pd = input->GetPointData();
  outPD = output->GetPointData();
  outPD->CopyNormalsOff();
  cd = input->GetCellData();
  outCD = output->GetCellData();
  int dims[3];
  input->GetDimensions(dims);
  inExt = input->GetExtent();

  // Based on the dimensions of the structured data, and the extent of
  // the geometry, compute the combined extent plus the dimensionality
  // of the data.
  //
  dimension = 3;
  for (i = 0; i < 3; i++)
  {
    extent[2 * i] = this->Extent[2 * i];
    extent[2 * i] = std::max(extent[2 * i], inExt[2 * i]);
    extent[2 * i + 1] = this->Extent[2 * i + 1];
    extent[2 * i + 1] = std::min(extent[2 * i + 1], inExt[2 * i + 1]);

    // Handle empty extent.
    if (extent[2 * i] > extent[2 * i + 1])
    {
      return 1;