Method initialize

src/ArtAlgorithm.cpp:87–143

--------------------------------------------------------------------------------------- Initialize - Config

//---------------------------------------------------------------------------------------
// Initialize - Config
bool CArtAlgorithm::initialize(const Config& _cfg)
{
	assert(!m_bIsInitialized);

	ConfigReader<CAlgorithm> CR("ArtAlgorithm", this, _cfg);
	
	// initialization of parent class
	if (!CReconstructionAlgorithm2D::initialize(_cfg)) {
		return false;
	}

	// ray order
	std::string projOrder;
	if (!CR.getOptionString("RayOrder", projOrder, "sequential"))
		return false;
	m_iCurrentRay = 0;
	int iRayCount = m_pProjector->getProjectionGeometry().getProjectionAngleCount() *
		m_pProjector->getProjectionGeometry().getDetectorCount();
	if (projOrder == "sequential") {
		m_piProjectionOrder.resize(iRayCount);
		m_piDetectorOrder.resize(iRayCount);
		for (int i = 0; i < iRayCount; i++) {
			m_piProjectionOrder[i] = i / m_pProjector->getProjectionGeometry().getDetectorCount();
			m_piDetectorOrder[i] = i % m_pProjector->getProjectionGeometry().getDetectorCount();
		}
	} else if (projOrder == "custom") {
		std::vector<int> rayOrderList;
		if (!CR.hasOption("RayOrderList"))
			return false;
		if (!CR.getOptionIntArray("RayOrderList", rayOrderList))
			return false;
		iRayCount = rayOrderList.size() / 2;
		m_piProjectionOrder.resize(iRayCount);
		m_piDetectorOrder.resize(iRayCount);
		for (int i = 0; i < iRayCount; i++) {
			m_piProjectionOrder[i] = rayOrderList[2*i];
			m_piDetectorOrder[i] = rayOrderList[2*i+1];
		}
	} else {
		ASTRA_ERROR("Unknown RayOrder");
		return false;
	}

	bool ok = true;

	// "Lambda" is replaced by the more descriptive "Relaxation"
	if (CR.hasOption("Relaxation"))
		ok &= CR.getOptionNumerical("Relaxation", m_fLambda, 1.0f);
	else
		ok &= CR.getOptionNumerical("Lambda", m_fLambda, 1.0f);
	if (!ok)
		return false;

	// success
	m_bIsInitialized = _check();
	return m_bIsInitialized;
}