Method GetImagePositionPatientFilenamePairs

Utilities/DICOMParser/DICOMAppHelper.cxx:1142–1203

}

void DICOMAppHelper::GetImagePositionPatientFilenamePairs(
  const std::string& seriesUID, std::vector<std::pair<float, std::string>>& v, bool ascending)
{
  v.clear();

  std::map<std::string, std::vector<std::string>, ltstdstr>::iterator miter =
    this->Implementation->SeriesUIDMap.find(seriesUID);

  if (miter == this->Implementation->SeriesUIDMap.end())
  {
    return;
  }

  // grab the filenames for the specified series
  std::vector<std::string> files = (*miter).second;

  for (std::vector<std::string>::iterator fileIter = files.begin(); fileIter != files.end();
       ++fileIter)
  {
    std::pair<float, std::string> p;
    p.second = std::string(*fileIter);

    float image_position;
    float normal[3];

    std::map<std::string, DICOMOrderingElements, ltstdstr>::iterator sn_iter =
      Implementation->SliceOrderingMap.find(*fileIter);

    if (sn_iter != Implementation->SliceOrderingMap.end())
    {
      // compute the image patient position wrt to the slice image
      // plane normal

      normal[0] = ((*sn_iter).second.ImageOrientationPatient[1] *
                    (*sn_iter).second.ImageOrientationPatient[5]) -
        ((*sn_iter).second.ImageOrientationPatient[2] *
          (*sn_iter).second.ImageOrientationPatient[4]);
      normal[1] = ((*sn_iter).second.ImageOrientationPatient[2] *
                    (*sn_iter).second.ImageOrientationPatient[3]) -
        ((*sn_iter).second.ImageOrientationPatient[0] *
          (*sn_iter).second.ImageOrientationPatient[5]);
      normal[2] = ((*sn_iter).second.ImageOrientationPatient[0] *
                    (*sn_iter).second.ImageOrientationPatient[4]) -
        ((*sn_iter).second.ImageOrientationPatient[1] *
          (*sn_iter).second.ImageOrientationPatient[3]);

      image_position = (normal[0] * (*sn_iter).second.ImagePositionPatient[0]) +
        (normal[1] * (*sn_iter).second.ImagePositionPatient[1]) +
        (normal[2] * (*sn_iter).second.ImagePositionPatient[2]);
      p.first = image_position;
      v.push_back(p);
    }
  }
  if (ascending)
  {
    std::sort(v.begin(), v.end(), lt_pair_float_string());
  }
  else