astapastrometryparser.cpp

/*
    SPDX-FileCopyrightText: 2016 Jasem Mutlaq <[email protected]>
 
    SPDX-License-Identifier: GPL-2.0-or-later
*/
 
#include "astapastrometryparser.h"
 
#include "align.h"
#include "ekos_align_debug.h"
#include "Options.h"
#include "indi/clientmanager.h"
#include "indi/driverinfo.h"
#include "indi/guimanager.h"
#include "indi/indidevice.h"
 
#include <indicom.h>
 
#include <KMessageBox>
 
namespace Ekos
{
ASTAPAstrometryParser::ASTAPAstrometryParser() : AstrometryParser()
{
}
 
bool ASTAPAstrometryParser::init()
{
    return QFile::exists(Options::aSTAPExecutable());
}
 
void ASTAPAstrometryParser::verifyIndexFiles(double, double)
{
}
 
bool ASTAPAstrometryParser::startSolver(const QString &filename, const QStringList &args, bool generated)
{
    INDI_UNUSED(generated);
 
    QStringList solverArgs = args;
 
    solverArgs << "-f" << filename;
 
    QString solutionFile = QDir::tempPath() + "/solution";
    solverArgs << "-o" << solutionFile;
 
    solver.clear();
    solver = new QProcess(this);
 
    connect(solver, static_cast<void (QProcess::*)(int exitCode, QProcess::ExitStatus exitStatus)>(&QProcess::finished),
            this, &ASTAPAstrometryParser::solverComplete);
    //    solver->setProcessChannelMode(QProcess::MergedChannels);
    //    connect(solver, SIGNAL(readyReadStandardOutput()), this, SLOT(logSolver()));
#if QT_VERSION > QT_VERSION_CHECK(5, 6, 0)
    connect(solver.data(), &QProcess::errorOccurred, this, [&]()
    {
        align->appendLogText(i18n("Error starting solver: %1", solver->errorString()));
        emit solverFailed();
    });
#else
    connect(solver, SIGNAL(error(QProcess::ProcessError)), this, SIGNAL(solverFailed()));
#endif
 
    solverTimer.start();
 
    QString solverPath = Options::aSTAPExecutable();
 
    solver->start(solverPath, solverArgs);
 
    align->appendLogText(i18n("Starting solver..."));
 
    if (Options::alignmentLogging())
    {
        QString command = solverPath + ' ' + solverArgs.join(' ');
        align->appendLogText(command);
    }
 
    return true;
}
 
void ASTAPAstrometryParser::solverComplete(int exitCode, QProcess::ExitStatus exitStatus)
{
    Q_UNUSED(exitCode)
    Q_UNUSED(exitStatus)
 
    QFile solution(QDir::tempPath() + "/solution.ini");
 
    if (!solution.open(QIODevice::ReadOnly))
    {
        qCritical(KSTARS_EKOS_ALIGN) << "Failed to open solution file" << QDir::tempPath() + "/solution.ini";
        emit solverFailed();
        return;
    }
 
    QTextStream in(&solution);
    QString line = in.readLine();
 
    QStringList ini = line.split("=");
    if (ini.count() <= 1 || ini[1] == "F")
    {
        align->appendLogText(i18n("Solver failed. Try again."));
        emit solverFailed();
        return;
    }
 
    double ra = 0, dec = 0, orientation = 0, pixscale = 0;
    double cd11 = 0;
    double cd22 = 0;
    double cd12 = 0;
    double cd21 = 0;
    bool ok[8] = {false};
 
    line = in.readLine();
    while (!line.isNull())
    {
        QStringList ini = line.split("=");
        if (ini[0] == "CRVAL1")
            ra = ini[1].trimmed().toDouble(&ok[0]);
        else if (ini[0] == "CRVAL2")
            dec = ini[1].trimmed().toDouble(&ok[1]);
        else if (ini[0] == "CDELT1")
            pixscale = ini[1].trimmed().toDouble(&ok[2]) * 3600.0;
        else if (ini[0] == "CROTA1")
            orientation = ini[1].trimmed().toDouble(&ok[3]);
        else if (ini[0] == "CD1_1")
            cd11 = ini[1].trimmed().toDouble(&ok[4]);
        else if (ini[0] == "CD1_2")
            cd12 = ini[1].trimmed().toDouble(&ok[5]);
        else if (ini[0] == "CD2_1")
            cd21 = ini[1].trimmed().toDouble(&ok[6]);
        else if (ini[0] == "CD2_2")
            cd22 = ini[1].trimmed().toDouble(&ok[7]);
 
        line = in.readLine();
    }
 
    if ( ok[0] && ok[1] && ok[2] && ok[3] )
    {
        Bool eastToTheRight = true;
        if ( ok[4] && ok[5] && ok[6] && ok[7] )
        {
            // Negative determinant = positive parity (i.e. east on the left).
            double det = cd11 * cd22 - cd12 * cd21;
            if(det < 0)
                eastToTheRight = false;
        }
        int elapsed = static_cast<int>(round(solverTimer.elapsed() / 1000.0));
        align->appendLogText(i18np("Solver completed in %1 second.", "Solver completed in %1 seconds.", elapsed));
        emit solverFinished(orientation, ra, dec, pixscale, eastToTheRight);
    }
    else
    {
        align->appendLogText(i18n("Solver failed. Try again."));
        emit solverFailed();
    }
}
 
bool ASTAPAstrometryParser::stopSolver()
{
    if (solver.isNull() == false)
    {
        solver->terminate();
        solver->disconnect();
    }
 
    return true;
}
}