diff --git a/kstars/ekos/align/align.cpp b/kstars/ekos/align/align.cpp index 567fff31d..59ffa640a 100644 --- a/kstars/ekos/align/align.cpp +++ b/kstars/ekos/align/align.cpp @@ -1,4388 +1,4456 @@ /* Ekos Alignment Module Copyright (C) 2013 Jasem Mutlaq This application is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. */ #include #include "kstars.h" #include "kstarsdata.h" #include "align.h" #include "dms.h" #include "fov.h" #include "ekos/auxiliary/darklibrary.h" #include "Options.h" #include #include #include #include #include "auxiliary/QProgressIndicator.h" #include "indi/driverinfo.h" #include "indi/indicommon.h" #include "indi/clientmanager.h" #include "alignadaptor.h" #include "fitsviewer/fitsviewer.h" #include "fitsviewer/fitstab.h" #include "fitsviewer/fitsview.h" #include "ekos/ekosmanager.h" #include "onlineastrometryparser.h" #include "offlineastrometryparser.h" #include "remoteastrometryparser.h" #include "opsastrometry.h" #include "opsalign.h" #include "opsastrometrycfg.h" #include "opsastrometryindexfiles.h" #include "skymapcomposite.h" #include "dialogs/finddialog.h" #include "ui_mountmodel.h" #include "starobject.h" +#include "skymap.h" +#include "flagcomponent.h" #include #define PAH_CUTOFF_FOV 30 // Minimum FOV width in arcminutes for PAH to work #define MAXIMUM_SOLVER_ITERATIONS 10 #define AL_FORMAT_VERSION 1.0 namespace Ekos { // 30 arcmiutes RA movement const double Align::RAMotion = 0.5; // Sidereal rate, degrees/s const float Align::SIDRATE = 0.004178; Align::Align() { setupUi(this); new AlignAdaptor(this); QDBusConnection::sessionBus().registerObject("/KStars/Ekos/Align", this); dirPath = QDir::homePath(); state = ALIGN_IDLE; focusState = FOCUS_IDLE; pahStage = PAH_IDLE; currentCCD = NULL; currentTelescope = NULL; currentFilter = NULL; useGuideHead = false; canSync = false; //loadSlewMode = false; loadSlewState=IPS_IDLE; //m_isSolverComplete = false; //m_isSolverSuccessful = false; //m_slewToTargetSelected=false; m_wcsSynced=false; //isFocusBusy=false; ccd_hor_pixel = ccd_ver_pixel = focal_length = aperture = sOrientation = sRA = sDEC = -1; decDeviation = azDeviation = altDeviation = 0; rememberUploadMode = ISD::CCD::UPLOAD_CLIENT; currentFilter = NULL; filterPositionPending = false; lockedFilterIndex = currentFilterIndex = -1; retries=0; targetDiff=1e6; solverIterations=0; fov_x=fov_y=fov_pixscale=0; parser = NULL; solverFOV = new FOV(); solverFOV->setColor(KStars::Instance()->data()->colorScheme()->colorNamed( "SolverFOVColor" ).name()); onlineParser = NULL; offlineParser = NULL; remoteParser = NULL; showFITSViewerB->setIcon(QIcon::fromTheme("kstars_fitsviewer", QIcon(":/icons/breeze/default/kstars_fitsviewer.svg"))); showFITSViewerB->setAttribute(Qt::WA_LayoutUsesWidgetRect); connect(showFITSViewerB, SIGNAL(clicked()), this, SLOT(showFITSViewer())); toggleFullScreenB->setIcon(QIcon::fromTheme("view-fullscreen", QIcon(":/icons/breeze/default/view-fullscreen.svg"))); toggleFullScreenB->setAttribute(Qt::WA_LayoutUsesWidgetRect); connect(toggleFullScreenB, SIGNAL(clicked()), this, SLOT(toggleAlignWidgetFullScreen())); alignView = new AlignView(alignWidget, FITS_ALIGN); alignView->setSizePolicy(QSizePolicy::Expanding, QSizePolicy::Expanding); alignView->setBaseSize(alignWidget->size()); alignView->createFloatingToolBar(); QVBoxLayout *vlayout = new QVBoxLayout(); vlayout->addWidget(alignView); alignWidget->setLayout(vlayout); connect(solveB, SIGNAL(clicked()), this, SLOT(captureAndSolve())); connect(stopB, SIGNAL(clicked()), this, SLOT(abort())); connect(measureAltB, SIGNAL(clicked()), this, SLOT(measureAltError())); connect(measureAzB, SIGNAL(clicked()), this, SLOT(measureAzError())); connect(CCDCaptureCombo, SIGNAL(activated(QString)), this, SLOT(setDefaultCCD(QString))); connect(CCDCaptureCombo, SIGNAL(currentIndexChanged(int)), this, SLOT(checkCCD(int))); connect(correctAltB, SIGNAL(clicked()), this, SLOT(correctAltError())); connect(correctAzB, SIGNAL(clicked()), this, SLOT(correctAzError())); connect(loadSlewB, SIGNAL(clicked()), this, SLOT(loadAndSlew())); gotoModeButtonGroup->setId(syncR, GOTO_SYNC); gotoModeButtonGroup->setId(slewR, GOTO_SLEW); gotoModeButtonGroup->setId(nothingR, GOTO_NOTHING); connect(gotoModeButtonGroup, static_cast(&QButtonGroup::buttonClicked), this, [=](int id){ this->currentGotoMode = static_cast(id); }); alignTimer.setInterval(Options::astrometryTimeout()*1000); connect(&alignTimer, SIGNAL(timeout()), this, SLOT(checkAlignmentTimeout())); currentGotoMode = static_cast(Options::solverGotoOption()); gotoModeButtonGroup->button(currentGotoMode)->setChecked(true); editOptionsB->setIcon(QIcon::fromTheme("document-edit", QIcon(":/icons/breeze/default/document-edit.svg"))); editOptionsB->setAttribute(Qt::WA_LayoutUsesWidgetRect); KConfigDialog* dialog = new KConfigDialog(this, "alignsettings", Options::self()); #ifdef Q_OS_OSX dialog->setWindowFlags(Qt::Tool| Qt::WindowStaysOnTopHint); #endif opsAlign = new OpsAlign(this); dialog->addPage(opsAlign, i18n("Astrometry.net")); opsAstrometry = new OpsAstrometry(this); dialog->addPage(opsAstrometry, i18n("Solver Options")); #ifdef Q_OS_OSX opsAstrometryCfg = new OpsAstrometryCfg(this); dialog->addPage(opsAstrometryCfg, i18n("Astrometry.cfg")); #endif #ifndef Q_OS_WIN opsAstrometryIndexFiles = new OpsAstrometryIndexFiles(this); dialog->addPage(opsAstrometryIndexFiles, i18n("Astrometry Index Files")); #endif connect(editOptionsB, SIGNAL(clicked()), dialog, SLOT(show())); appendLogText(i18n("Idle.")); pi = new QProgressIndicator(this); stopLayout->addWidget(pi); exposureIN->setValue(Options::alignExposure()); altStage = ALT_INIT; azStage = AZ_INIT; // Online/Offline/Remote solver check solverTypeGroup->setId(onlineSolverR, SOLVER_ONLINE); solverTypeGroup->setId(offlineSolverR, SOLVER_OFFLINE); solverTypeGroup->setId(remoteSolverR, SOLVER_REMOTE); #ifdef Q_OS_WIN offlineSolverR->setEnabled(false); offlineSolverR->setToolTip(i18n("Offline solver is not supported under Windows. Please use either the Online or Remote solvers.")); #endif solverTypeGroup->button(Options::solverType())->setChecked(true); connect(solverTypeGroup, SIGNAL(buttonClicked(int)), SLOT(setSolverType(int))); switch(solverTypeGroup->checkedId()) { case SOLVER_ONLINE: onlineParser = new Ekos::OnlineAstrometryParser(); parser = onlineParser; break; case SOLVER_OFFLINE: offlineParser = new OfflineAstrometryParser(); parser = offlineParser; break; case SOLVER_REMOTE: remoteParser = new RemoteAstrometryParser(); parser = remoteParser; break; } parser->setAlign(this); if (parser->init() == false) setEnabled(false); else { connect(parser, SIGNAL(solverFinished(double,double,double, double)), this, SLOT(solverFinished(double,double,double, double)), Qt::UniqueConnection); connect(parser, SIGNAL(solverFailed()), this, SLOT(solverFailed()), Qt::UniqueConnection); } //solverOptions->setText(Options::solverOptions()); // Which telescope info to use for FOV calculations //kcfg_solverOTA->setChecked(Options::solverOTA()); guideScopeCCDs = Options::guideScopeCCDs(); connect(FOVScopeCombo, SIGNAL(currentIndexChanged(int)), this, SLOT(updateGuideScopeCCDs(int))); accuracySpin->setValue(Options::solverAccuracyThreshold()); alignDarkFrameCheck->setChecked(Options::alignDarkFrame()); delaySpin->setValue(Options::settlingTime()); connect(delaySpin, SIGNAL(editingFinished()), this, SLOT(saveSettleTime())); connect(binningCombo, SIGNAL(currentIndexChanged(int)), this, SLOT(setBinningIndex(int))); // PAH Connections connect(PAHRestartB, SIGNAL(clicked()), this, SLOT(restartPAHProcess())); connect(PAHStartB, SIGNAL(clicked()), this, SLOT(startPAHProcess())); connect(PAHFirstCaptureB, &QPushButton::clicked, this, [this]() { PAHFirstCaptureB->setEnabled(false); captureAndSolve();}); connect(PAHSecondCaptureB, &QPushButton::clicked, this, [this]() { PAHSecondCaptureB->setEnabled(false); captureAndSolve();}); connect(PAHThirdCaptureB, &QPushButton::clicked, this, [this]() { PAHThirdCaptureB->setEnabled(false); captureAndSolve();}); connect(PAHFirstRotateB, &QPushButton::clicked, this, [this]() { PAHFirstRotateB->setEnabled(false); rotatePAH();}); connect(PAHSecondRotateB, &QPushButton::clicked, this, [this]() { PAHSecondRotateB->setEnabled(false); rotatePAH();}); connect(PAHCorrectionsNextB, SIGNAL(clicked()), this, SLOT(setPAHCorrectionSelectionComplete())); connect(PAHRefreshB, SIGNAL(clicked()), this, SLOT(startPAHRefreshProcess())); connect(PAHDoneB, SIGNAL(clicked()), this, SLOT(setPAHRefreshComplete())); if (solverOptions->text().contains("no-fits2fits")) appendLogText(i18n("Warning: If using astrometry.net v0.68 or above, remove the --no-fits2fits from the astrometry options.")); hemisphere = KStarsData::Instance()->geo()->lat()->Degrees() > 0 ? NORTH_HEMISPHERE : SOUTH_HEMISPHERE; double accuracyRadius=accuracySpin->value(); alignPlot->setBackground(QBrush(Qt::black)); alignPlot->setSelectionTolerance(10); alignPlot->xAxis->setBasePen(QPen(Qt::white, 1)); alignPlot->yAxis->setBasePen(QPen(Qt::white, 1)); alignPlot->xAxis->setTickPen(QPen(Qt::white, 1)); alignPlot->yAxis->setTickPen(QPen(Qt::white, 1)); alignPlot->xAxis->setSubTickPen(QPen(Qt::white, 1)); alignPlot->yAxis->setSubTickPen(QPen(Qt::white, 1)); alignPlot->xAxis->setTickLabelColor(Qt::white); alignPlot->yAxis->setTickLabelColor(Qt::white); alignPlot->xAxis->setLabelColor(Qt::white); alignPlot->yAxis->setLabelColor(Qt::white); alignPlot->xAxis->setLabelFont(QFont(font().family(), 10)); alignPlot->yAxis->setLabelFont(QFont(font().family(), 10)); alignPlot->xAxis->setLabelPadding(2); alignPlot->yAxis->setLabelPadding(2); alignPlot->xAxis->grid()->setPen(QPen(QColor(140, 140, 140), 1, Qt::DotLine)); alignPlot->yAxis->grid()->setPen(QPen(QColor(140, 140, 140), 1, Qt::DotLine)); alignPlot->xAxis->grid()->setSubGridPen(QPen(QColor(80, 80, 80), 1, Qt::DotLine)); alignPlot->yAxis->grid()->setSubGridPen(QPen(QColor(80, 80, 80), 1, Qt::DotLine)); alignPlot->xAxis->grid()->setZeroLinePen(QPen(Qt::yellow)); alignPlot->yAxis->grid()->setZeroLinePen(QPen(Qt::yellow)); alignPlot->xAxis->setLabel(i18n("dRA (arcsec)")); alignPlot->yAxis->setLabel(i18n("dDE (arcsec)")); alignPlot->xAxis->setRange(-accuracyRadius*3, accuracyRadius*3); alignPlot->yAxis->setRange(-accuracyRadius*3, accuracyRadius*3); alignPlot->setInteractions(QCP::iRangeZoom); alignPlot->setInteraction(QCP::iRangeDrag, true); alignPlot->addGraph(); alignPlot->graph(0)->setLineStyle(QCPGraph::lsNone); alignPlot->graph(0)->setScatterStyle(QCPScatterStyle(QCPScatterStyle::ssDisc, Qt::white, 15)); buildTarget(); connect(alignPlot, SIGNAL(mouseMove(QMouseEvent*)), this, SLOT(handlePointTooltip(QMouseEvent*))); connect(rightLayout, SIGNAL(splitterMoved(int,int)), this, SLOT(handleVerticalPlotSizeChange())); connect(alignSplitter, SIGNAL(splitterMoved(int,int)), this, SLOT(handleHorizontalPlotSizeChange())); connect(accuracySpin, SIGNAL(valueChanged(int)), this, SLOT(buildTarget())); alignPlot->resize(190, 190); alignPlot->replot(); solutionTable->setColumnWidth(0,70); solutionTable->setColumnWidth(1,75); solutionTable->setColumnWidth(2,80); solutionTable->setColumnWidth(3,30); solutionTable->setColumnWidth(4,100); solutionTable->setColumnWidth(5,100); clearAllSolutionsB->setIcon(QIcon::fromTheme("application-exit", QIcon(":/icons/breeze/default/application-exit.svg") )); clearAllSolutionsB->setAttribute(Qt::WA_LayoutUsesWidgetRect); removeSolutionB->setIcon(QIcon::fromTheme("list-remove", QIcon(":/icons/breeze/default/list-remove.svg") )); removeSolutionB->setAttribute(Qt::WA_LayoutUsesWidgetRect); exportSolutionsCSV->setIcon(QIcon::fromTheme("document-save-as", QIcon(":/icons/breeze/default/document-save-as.svg") )); exportSolutionsCSV->setAttribute(Qt::WA_LayoutUsesWidgetRect); mountModel.setupUi(&mountModelDialog); mountModelDialog.setWindowTitle("Mount Model Tool"); mountModelDialog.setWindowFlags(Qt::Tool| Qt::WindowStaysOnTopHint); mountModel.alignTable->setColumnWidth(0,70); mountModel.alignTable->setColumnWidth(1,75); mountModel.alignTable->setColumnWidth(2,130); mountModel.alignTable->setColumnWidth(3,30); mountModel.wizardAlignB->setIcon(QIcon::fromTheme("tools-wizard", QIcon(":/icons/breeze/default/tools-wizard.svg") )); mountModel.wizardAlignB->setAttribute(Qt::WA_LayoutUsesWidgetRect); mountModel.clearAllAlignB->setIcon(QIcon::fromTheme("application-exit", QIcon(":/icons/breeze/default/application-exit.svg") )); mountModel.clearAllAlignB->setAttribute(Qt::WA_LayoutUsesWidgetRect); mountModel.removeAlignB->setIcon(QIcon::fromTheme("list-remove", QIcon(":/icons/breeze/default/list-remove.svg") )); mountModel.removeAlignB->setAttribute(Qt::WA_LayoutUsesWidgetRect); mountModel.addAlignB->setIcon(QIcon::fromTheme("list-add", QIcon(":/icons/breeze/default/list-add.svg") )); mountModel.addAlignB->setAttribute(Qt::WA_LayoutUsesWidgetRect); mountModel.findAlignB->setIcon(QIcon::fromTheme("edit-find", QIcon(":/icons/breeze/default/edit-find.svg") )); mountModel.findAlignB->setAttribute(Qt::WA_LayoutUsesWidgetRect); mountModel.alignTable->verticalHeader()->setDragDropOverwriteMode(false); mountModel.alignTable->verticalHeader()->setSectionsMovable(true); mountModel.alignTable->verticalHeader()->setDragEnabled(true); mountModel.alignTable->verticalHeader()->setDragDropMode(QAbstractItemView::InternalMove); connect(mountModel.alignTable->verticalHeader(), SIGNAL(sectionMoved(int, int , int)), this, SLOT(moveAlignPoint(int, int , int))); mountModel.loadAlignB->setIcon(QIcon::fromTheme("document-open", QIcon(":/icons/breeze/default/document-open.svg") )); mountModel.loadAlignB->setAttribute(Qt::WA_LayoutUsesWidgetRect); mountModel.saveAsAlignB->setIcon(QIcon::fromTheme("document-save-as", QIcon(":/icons/breeze/default/document-save-as.svg") )); mountModel.saveAsAlignB->setAttribute(Qt::WA_LayoutUsesWidgetRect); mountModel.saveAlignB->setIcon(QIcon::fromTheme("document-save", QIcon(":/icons/breeze/default/document-save.svg") )); mountModel.saveAlignB->setAttribute(Qt::WA_LayoutUsesWidgetRect); + mountModel.previewB->setIcon(QIcon::fromTheme("kstars_grid", QIcon(":/icons/breeze/default/kstars_grid.svg") )); + mountModel.previewB->setAttribute(Qt::WA_LayoutUsesWidgetRect); + mountModel.previewB->setCheckable(true); + mountModel.stopAlignB->setIcon(QIcon::fromTheme("media-playback-stop", QIcon(":/icons/breeze/default/media-playback-stop.svg") )); mountModel.stopAlignB->setAttribute(Qt::WA_LayoutUsesWidgetRect); mountModel.startAlignB->setIcon(QIcon::fromTheme("media-playback-start", QIcon(":/icons/breeze/default/media-playback-start.svg") )); mountModel.startAlignB->setAttribute(Qt::WA_LayoutUsesWidgetRect); connect(clearAllSolutionsB, SIGNAL(clicked()), this, SLOT(slotClearAllSolutionPoints())); connect(removeSolutionB, SIGNAL(clicked()), this, SLOT(slotRemoveSolutionPoint())); connect(exportSolutionsCSV, SIGNAL(clicked()), this, SLOT(exportSolutionPoints())); connect(mountModelB, SIGNAL(clicked()), this, SLOT(slotMountModel())); connect(solutionTable, SIGNAL(cellClicked(int, int)), this, SLOT(selectSolutionTableRow(int, int))); connect(mountModel.wizardAlignB, SIGNAL(clicked()), this, SLOT(slotWizardAlignmentPoints())); connect(mountModel.alignTypeBox, SIGNAL(currentIndexChanged(const QString)), this, SLOT(alignTypeChanged(const QString))); connect(mountModel.starListBox, SIGNAL(currentIndexChanged(const QString)), this, SLOT(slotStarSelected(const QString))); connect(mountModel.greekStarListBox, SIGNAL(currentIndexChanged(const QString)), this, SLOT(slotStarSelected(const QString))); connect(mountModel.loadAlignB, SIGNAL(clicked()), this, SLOT(slotLoadAlignmentPoints())); connect(mountModel.saveAsAlignB, SIGNAL(clicked()), this, SLOT(slotSaveAsAlignmentPoints())); connect(mountModel.saveAlignB, SIGNAL(clicked()), this, SLOT(slotSaveAlignmentPoints())); connect(mountModel.clearAllAlignB, SIGNAL(clicked()), this, SLOT(slotClearAllAlignPoints())); connect(mountModel.removeAlignB, SIGNAL(clicked()), this, SLOT(slotRemoveAlignPoint())); connect(mountModel.addAlignB, SIGNAL(clicked()), this, SLOT(slotAddAlignPoint())); connect(mountModel.findAlignB, SIGNAL(clicked()), this, SLOT(slotFindAlignObject())); + connect(mountModel.previewB, SIGNAL(clicked()), this, SLOT(togglePreviewAlignPoints())); connect(mountModel.stopAlignB, SIGNAL(clicked()), this, SLOT(resetAlignmentProcedure())); connect(mountModel.startAlignB, SIGNAL(clicked()), this, SLOT(startStopAlignmentProcedure())); } Align::~Align() { delete(pi); delete(solverFOV); delete(parser); if (alignWidget->parent() == NULL) toggleAlignWidgetFullScreen(); // Remove temporary FITS files left before by the solver QDir dir(QDir::tempPath()); dir.setNameFilters(QStringList() << "fits*" << "tmp.*"); dir.setFilter(QDir::Files); foreach(QString dirFile, dir.entryList()) dir.remove(dirFile); } void Align::selectSolutionTableRow(int row, int column){ solutionTable->selectRow(row); for(int i=0;iitemCount();i++){ QCPAbstractItem *abstractItem=alignPlot->item(i); if(abstractItem){ QCPItemText *item = qobject_cast(abstractItem); if(item){ if(i==row){ item->setColor(Qt::black); item->setBrush(Qt::yellow); }else{ item->setColor(Qt::red); item->setBrush(Qt::white); } } } } alignPlot->replot(); } void Align::handleHorizontalPlotSizeChange(){ alignPlot->xAxis->setScaleRatio(alignPlot->yAxis,1.0); alignPlot->replot(); } void Align::handleVerticalPlotSizeChange(){ alignPlot->yAxis->setScaleRatio(alignPlot->xAxis,1.0); alignPlot->replot(); } void Align::resizeEvent(QResizeEvent * event){ if(event->oldSize().width()!=-1){ if(event->oldSize().width()!=size().width()) handleHorizontalPlotSizeChange(); else if(event->oldSize().height()!=size().height()) handleVerticalPlotSizeChange(); } else{ QTimer::singleShot(10, this, SLOT(handleHorizontalPlotSizeChange())); } } void Align::handlePointTooltip(QMouseEvent* event){ QCPAbstractItem *item =alignPlot->itemAt(event->localPos()); if(item){ QCPItemText *label = qobject_cast(item); if(label){ QString labelText=label->text(); int point=labelText.toInt()-1; if (point < 0) return; QToolTip::showText(event->globalPos(), tr("" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "" "
Object %L1: %L2
RA:%L3
DE:%L4
dRA:%L5
dDE:%L6
"). arg(point+1). arg(solutionTable->item(point, 2)->text()). arg(solutionTable->item(point, 0)->text()). arg(solutionTable->item(point, 1)->text()). arg(solutionTable->item(point, 4)->text()). arg(solutionTable->item(point, 5)->text()), alignPlot, alignPlot->rect()); } } } void Align::buildTarget(){ double accuracyRadius=accuracySpin->value(); if(centralTarget){ concentricRings->data()->clear(); redTarget->data()->clear(); yellowTarget->data()->clear(); centralTarget->data()->clear(); }else{ concentricRings = new QCPCurve(alignPlot->xAxis, alignPlot->yAxis); redTarget = new QCPCurve(alignPlot->xAxis, alignPlot->yAxis); yellowTarget = new QCPCurve(alignPlot->xAxis, alignPlot->yAxis); centralTarget = new QCPCurve(alignPlot->xAxis, alignPlot->yAxis); } const int pointCount = 200; QVector circleRings(pointCount*(5));//Have to multiply by the number of rings, Rings at : 25%, 50%, 75%, 125%, 175% QVector circleCentral(pointCount); QVector circleYellow(pointCount); QVector circleRed(pointCount); int circleRingPt=0; for (int i=0; isetLineStyle(QCPCurve::lsNone); concentricRings->setScatterStyle(QCPScatterStyle(QCPScatterStyle::ssDisc, QColor(255, 255, 255, 150), 1)); concentricRings->data()->set(circleRings, true); redTarget->data()->set(circleRed, true); yellowTarget->data()->set(circleYellow, true); centralTarget->data()->set(circleCentral, true); concentricRings->setPen(QPen(Qt::white)); redTarget->setPen(QPen(Qt::red)); yellowTarget->setPen(QPen(Qt::yellow)); centralTarget->setPen(QPen(Qt::green)); concentricRings->setBrush(Qt::NoBrush); redTarget->setBrush(QBrush(QColor(255, 0, 0, 50))); yellowTarget->setBrush(QBrush(QColor(0, 255, 0, 50))); //Note this is actually yellow. It is green on top of red with equal opacity. centralTarget->setBrush(QBrush(QColor(0, 255, 0, 50))); if(alignPlot->size().width()>0) alignPlot->replot(); } void Align::slotWizardAlignmentPoints(){ int points=mountModel.alignPtNum->value(); if(points<2) //The minimum is 2 because the wizard calculations require the calculation of an angle between points. return; //It should not be less than 2 because the minimum in the spin box is 2. int minAlt=mountModel.minAltBox->value(); KStarsData *data = KStarsData::Instance(); GeoLocation *geo=data->geo(); double lat=geo->lat()->Degrees(); if(mountModel.alignTypeBox->currentText()=="Fixed DEC"){ double decAngle=mountModel.alignDec->value(); //Dec that never rises. if(lat>0){ if(decAngle < lat - 90 + minAlt){ //Min altitude possible at minAlt deg above horizon KMessageBox::sorry( 0, i18n( "DEC is below the altitude limit" ) ); return; } } else{ if(decAngle > lat + 90 - minAlt){ //Max altitude possible at minAlt deg above horizon KMessageBox::sorry( 0, i18n( "DEC is below the altitude limit" ) ); return; } } } + //If there are less than 6 points, keep them all in the same DEC, + //any more, set the num per row to be the sqrt of the points to evenly distribute in RA and DEC + int numRAperDEC=5; + if(points>5) + numRAperDEC=qSqrt(points); + //These calculations rely on modulus and int division counting beginning at 0, but the #s start at 1. - int decPoints=(points-1)/5+1; - int lastSetRAPoints=(points-1)%5+1; + int decPoints=(points-1)/numRAperDEC+1; + int lastSetRAPoints=(points-1)%numRAperDEC+1; double decIncrement = -1; double initDEC = -1; SkyPoint spTest; if(mountModel.alignTypeBox->currentText()=="Fixed DEC"){ decPoints=1; initDEC=mountModel.alignDec->value(); decIncrement=0; } else if(decPoints==1) { decIncrement=0; spTest.setAlt(minAlt); //The goal here is to get the point exactly West at the minAlt so that we can use that DEC spTest.setAz(270); spTest.HorizontalToEquatorial(KStars::Instance()->data()->lst(), KStars::Instance()->data()->geo()->lat()); initDEC=spTest.dec().Degrees(); } else{ spTest.setAlt(minAlt+10); //We don't want to be right at the minAlt because there would be only 1 point on the dec circle above the alt. spTest.setAz(180); spTest.HorizontalToEquatorial(KStars::Instance()->data()->lst(), KStars::Instance()->data()->geo()->lat()); initDEC=spTest.dec().Degrees(); if(lat>0) decIncrement=(80-initDEC)/(decPoints); //Don't quite want to reach NCP else decIncrement=(initDEC-80)/(decPoints); //Don't quite want to reach SCP } qDebug()<0) dec = initDEC + d * decIncrement; else dec = initDEC - d * decIncrement; if(mountModel.alignTypeBox->currentText()=="Fixed DEC"){ raPoints=points; }else if(d==decPoints-1){ raPoints=lastSetRAPoints; }else{ - raPoints=5; + raPoints=numRAperDEC; } //This computes both the initRA and the raIncrement. calculateAngleForRALine(raIncrement, initRA, dec, lat, raPoints, minAlt); if(raIncrement==-1||decIncrement==-1){ KMessageBox::sorry( 0, i18n( "Point Calculation Error." ) ); return; } for(int i=0;ira().Hours(),o->dec().Degrees(),ra_report, dec_report); name=o->longname(); } else{ getFormattedCoords(dms(ra).Hours(),dec,ra_report, dec_report); name="None"; } int currentRow = mountModel.alignTable->rowCount(); mountModel.alignTable->insertRow(currentRow); QTableWidgetItem *RAReport = new QTableWidgetItem(); RAReport->setText(ra_report); RAReport->setTextAlignment(Qt::AlignHCenter); mountModel.alignTable->setItem(currentRow, 0, RAReport); QTableWidgetItem *DECReport = new QTableWidgetItem(); DECReport->setText(dec_report); DECReport->setTextAlignment(Qt::AlignHCenter); mountModel.alignTable->setItem(currentRow, 1, DECReport); QTableWidgetItem *ObjNameReport = new QTableWidgetItem(); ObjNameReport->setText(name); ObjNameReport->setTextAlignment(Qt::AlignHCenter); mountModel.alignTable->setItem(currentRow, 2, ObjNameReport); QTableWidgetItem *disabledBox= new QTableWidgetItem(); disabledBox->setFlags(Qt::ItemIsSelectable); mountModel.alignTable->setItem(currentRow, 3, disabledBox); } } + if(previewShowing) + updatePreviewAlignPoints(); } void Align::calculateAngleForRALine(double &raIncrement,double &initRA, double initDEC, double lat, double raPoints, double minAlt){ SkyPoint spEast; SkyPoint spWest; //Circumpolar dec if(fabs(initDEC) > (90 - fabs(lat) + minAlt)){ raIncrement = 360 / (raPoints-1); initRA=0; }else{ dms AZEast,AZWest; calculateAZPointsForDEC(dms(initDEC), dms(minAlt), AZEast,AZWest); spEast.setAlt(minAlt); spEast.setAz(AZEast.Degrees()); spEast.HorizontalToEquatorial(KStars::Instance()->data()->lst(), KStars::Instance()->data()->geo()->lat()); spWest.setAlt(minAlt); spWest.setAz(AZWest.Degrees()); spWest.HorizontalToEquatorial(KStars::Instance()->data()->lst(), KStars::Instance()->data()->geo()->lat()); dms angleSep; if(spEast.ra().Degrees()>spWest.ra().Degrees()) angleSep=spEast.ra()-spWest.ra(); else angleSep=spEast.ra()+dms(360)-spWest.ra(); initRA=spWest.ra().Degrees(); raIncrement=angleSep.Degrees()/(raPoints-1); } } void Align::calculateAZPointsForDEC(dms dec, dms alt, dms &AZEast, dms &AZWest){ KStarsData *data = KStarsData::Instance(); GeoLocation *geo=data->geo(); double AZRad; double sindec, cosdec, sinlat, coslat; double sinAlt, cosAlt; geo->lat()->SinCos( sinlat, coslat ); dec.SinCos( sindec, cosdec ); alt.SinCos( sinAlt, cosAlt ); double arg = ( sindec - sinlat*sinAlt )/( coslat*cosAlt ); AZRad = acos( arg ); AZEast.setRadians(AZRad); AZWest.setRadians(2.0*dms::PI - AZRad); } const SkyObject* Align::getWizardAlignObject(double ra, double dec, double angle){ double maxSearch=5.0; if(mountModel.alignTypeBox->currentText()=="Any Object") return KStarsData::Instance()->skyComposite()->objectNearest(new SkyPoint(dms(ra),dms(dec)), maxSearch ); else if(mountModel.alignTypeBox->currentText()=="Fixed DEC"||mountModel.alignTypeBox->currentText()=="Fixed Grid") return NULL; else if(mountModel.alignTypeBox->currentText()=="Any Stars") return KStarsData::Instance()->skyComposite()->starNearest(new SkyPoint(dms(ra),dms(dec)), maxSearch ); //If they want named stars, then try to search for and return the closest Align Star to the requested location dms bestDiff=dms(360); double index=-1; for(int i=0;ihasName()){ SkyPoint thisPt(ra,dec); dms thisDiff=thisPt.angularDistanceTo(star); if(thisDiff.Degrees() < bestDiff.Degrees()){ index=i; bestDiff=thisDiff; } } } } if(index==-1) return KStarsData::Instance()->skyComposite()->starNearest(new SkyPoint(dms(ra),dms(dec)), maxSearch ); return alignStars.value(index); } void Align::alignTypeChanged(const QString alignType){ if(alignType=="Fixed DEC") mountModel.alignDec->setEnabled(true); else mountModel.alignDec->setEnabled(false); } void Align::slotStarSelected(const QString selectedStar){ for(int i=0;iname()==selectedStar||star->gname().simplified()==selectedStar){ int currentRow = mountModel.alignTable->rowCount(); mountModel.alignTable->insertRow(currentRow); QString ra_report, dec_report; getFormattedCoords(star->ra().Hours(),star->dec().Degrees(),ra_report, dec_report); QTableWidgetItem *RAReport = new QTableWidgetItem(); RAReport->setText(ra_report); RAReport->setTextAlignment(Qt::AlignHCenter); mountModel.alignTable->setItem(currentRow, 0, RAReport); QTableWidgetItem *DECReport = new QTableWidgetItem(); DECReport->setText(dec_report); DECReport->setTextAlignment(Qt::AlignHCenter); mountModel.alignTable->setItem(currentRow, 1, DECReport); QTableWidgetItem *ObjNameReport = new QTableWidgetItem(); ObjNameReport->setText(star->longname()); ObjNameReport->setTextAlignment(Qt::AlignHCenter); mountModel.alignTable->setItem(currentRow, 2, ObjNameReport); QTableWidgetItem *disabledBox= new QTableWidgetItem(); disabledBox->setFlags(Qt::ItemIsSelectable); mountModel.alignTable->setItem(currentRow, 3, disabledBox); mountModel.starListBox->setCurrentIndex(0); mountModel.greekStarListBox->setCurrentIndex(0); return; } } } + if(previewShowing) + updatePreviewAlignPoints(); } void Align::generateAlignStarList(){ alignStars.clear(); mountModel.starListBox->clear(); mountModel.greekStarListBox->clear(); KStarsData *data = KStarsData::Instance(); QVector> listStars; listStars.append(data->skyComposite()->objectLists(SkyObject::STAR)); for(int i=0;i pair=listStars.value(i); const StarObject *star = dynamic_cast(pair.second); if(star) alignStars.append(star); } QStringList boxNames; QStringList greekBoxNames; for(int i=0;ihasLatinName()) boxNames << star->name(); else{ if(!star->gname().isEmpty()) greekBoxNames << star->gname().simplified(); } } } } boxNames.sort(Qt::CaseInsensitive); boxNames.removeDuplicates(); greekBoxNames.removeDuplicates(); qSort(greekBoxNames.begin(), greekBoxNames.end(), [](const QString &a, const QString &b) { QStringList aParts=a.split(" "); QStringList bParts=b.split(" "); if(aParts.length()<2||bParts.length()<2) return a < b; //This should not happen, they should all have 2 words in the string. if( aParts[1]==bParts[1]){ return aParts[0] < bParts[0]; //This compares the greek letter when the constellation is the same } else return aParts[1] < bParts[1]; //This compares the constellation names } ); mountModel.starListBox->addItem("Select one:"); mountModel.greekStarListBox->addItem("Select one:"); for(int i=0;iaddItem(boxNames.at(i)); for(int i=0;iaddItem(greekBoxNames.at(i)); } bool Align::isVisible(const SkyObject *so){ return (getAltitude(so) > 30); } double Align::getAltitude(const SkyObject *so){ KStarsData *data = KStarsData::Instance(); GeoLocation *geo=data->geo(); CachingDms *lst=data->lst(); KStarsDateTime ut = geo->LTtoUT(KStarsDateTime(QDateTime::currentDateTime().toLocalTime())); SkyPoint sp = so->recomputeCoords(ut, geo); //check altitude of object at this time. sp.EquatorialToHorizontal(lst, geo->lat()); return sp.alt().Degrees(); } +void Align::togglePreviewAlignPoints(){ + previewShowing=!previewShowing; + mountModel.previewB->setChecked(previewShowing); + updatePreviewAlignPoints(); +} + +void Align::updatePreviewAlignPoints(){ + + FlagComponent *flags = KStarsData::Instance()->skyComposite()->flags(); + for(int i=0;isize();i++){ + if(flags->label(i).startsWith("Align")){ + flags->remove(i); + i--; + } + } + if(previewShowing){ + for (int i=0; i < mountModel.alignTable->rowCount(); i++) + { + QTableWidgetItem *raCell=mountModel.alignTable->item(i,0); + QTableWidgetItem *deCell=mountModel.alignTable->item(i,1); + QTableWidgetItem *objNameCell=mountModel.alignTable->item(i,2); + + if(raCell&&deCell&&objNameCell){ + QString raString=raCell->text(); + QString deString=deCell->text(); + dms raDMS = dms::fromString(raString,false); + dms decDMS = dms::fromString(deString,true); + + QString objString=objNameCell->text(); + + + SkyPoint flagPoint( raDMS, decDMS ); + flags->add( flagPoint, QString::number(KStarsDateTime::currentDateTime().epoch()), "Default", "Align " + QString::number(i+1)+ " " + objString, "white" ); + } + } + } + KStars::Instance()->map()->forceUpdate(true); + +} + void Align::slotLoadAlignmentPoints(){ QUrl fileURL = QFileDialog::getOpenFileUrl(KStars::Instance(), i18n("Open Ekos Alignment List"), alignURLPath, "Ekos AlignmentList (*.eal)"); if (fileURL.isEmpty()) return; if (fileURL.isValid() == false) { QString message = i18n( "Invalid URL: %1", fileURL.toLocalFile() ); KMessageBox::sorry( 0, message, i18n( "Invalid URL" ) ); return; } alignURLPath = QUrl(fileURL.url(QUrl::RemoveFilename)); loadAlignmentPoints(fileURL.toLocalFile()); + if(previewShowing) + updatePreviewAlignPoints(); } bool Align::loadAlignmentPoints(const QString &fileURL){ QFile sFile; sFile.setFileName(fileURL); if ( !sFile.open( QIODevice::ReadOnly)) { QString message = i18n( "Unable to open file %1", fileURL); KMessageBox::sorry( 0, message, i18n( "Could Not Open File" ) ); return false; } mountModel.alignTable->setRowCount(0); LilXML *xmlParser = newLilXML(); char errmsg[MAXRBUF]; XMLEle *root = NULL; char c; while ( sFile.getChar(&c)) { root = readXMLEle(xmlParser, c, errmsg); if (root) { double sqVersion= atof(findXMLAttValu(root, "version")); if (sqVersion < AL_FORMAT_VERSION) { appendLogText(i18n("Deprecated sequence file format version %1. Please construct a new sequence file.", sqVersion)); return false; } XMLEle *ep = NULL; XMLEle *subEP = NULL; int currentRow=0; for (ep = nextXMLEle(root, 1) ; ep != NULL ; ep = nextXMLEle(root, 0)) { if (!strcmp(tagXMLEle(ep), "AlignmentPoint")) { mountModel.alignTable->insertRow(currentRow); subEP = findXMLEle(ep, "RA"); if (subEP){ QTableWidgetItem *RAReport = new QTableWidgetItem(); RAReport->setText(pcdataXMLEle(subEP)); RAReport->setTextAlignment(Qt::AlignHCenter); mountModel.alignTable->setItem(currentRow, 0, RAReport); } else return false; subEP = findXMLEle(ep, "DE"); if (subEP){ QTableWidgetItem *DEReport = new QTableWidgetItem(); DEReport->setText(pcdataXMLEle(subEP)); DEReport->setTextAlignment(Qt::AlignHCenter); mountModel.alignTable->setItem(currentRow, 1, DEReport); } else return false; subEP = findXMLEle(ep, "NAME"); if (subEP){ QTableWidgetItem *ObjReport = new QTableWidgetItem(); ObjReport->setText(pcdataXMLEle(subEP)); ObjReport->setTextAlignment(Qt::AlignHCenter); mountModel.alignTable->setItem(currentRow, 2, ObjReport); } else return false; } currentRow++; } return true; } } return false; } void Align::slotSaveAsAlignmentPoints(){ alignURL.clear(); slotSaveAlignmentPoints(); } void Align::slotSaveAlignmentPoints() { QUrl backupCurrent = alignURL; if (alignURL.toLocalFile().startsWith("/tmp/") || alignURL.toLocalFile().contains("/Temp")) alignURL.clear(); if (alignURL.isEmpty()) { alignURL = QFileDialog::getSaveFileUrl(KStars::Instance(), i18n("Save Ekos Alignment List"), alignURLPath, "Ekos Alignment List (*.eal)"); // if user presses cancel if (alignURL.isEmpty()) { alignURL = backupCurrent; return; } alignURLPath = QUrl(alignURL.url(QUrl::RemoveFilename)); if (alignURL.toLocalFile().endsWith(".eal") == false) alignURL.setPath(alignURL.toLocalFile() + ".eal"); if (QFile::exists(alignURL.toLocalFile())) { int r = KMessageBox::warningContinueCancel(0, i18n( "A file named \"%1\" already exists. " "Overwrite it?", alignURL.fileName() ), i18n( "Overwrite File?" ), KStandardGuiItem::overwrite() ); if(r==KMessageBox::Cancel) return; } } if ( alignURL.isValid() ) { if ( (saveAlignmentPoints(alignURL.toLocalFile())) == false) { KMessageBox::error(KStars::Instance(), i18n("Failed to save alignment list"), i18n("Save")); return; } } else { QString message = i18n( "Invalid URL: %1", alignURL.url() ); KMessageBox::sorry(KStars::Instance(), message, i18n( "Invalid URL" ) ); } } bool Align::saveAlignmentPoints(const QString &path){ QFile file; file.setFileName(path); if ( !file.open( QIODevice::WriteOnly)) { QString message = i18n( "Unable to write to file %1", path); KMessageBox::sorry( 0, message, i18n( "Could Not Open File" ) ); return false; } QTextStream outstream(&file); outstream << "" << endl; outstream << "" << endl; for (int i=0; i < mountModel.alignTable->rowCount(); i++) { QTableWidgetItem *raCell=mountModel.alignTable->item(i,0); QTableWidgetItem *deCell=mountModel.alignTable->item(i,1); QTableWidgetItem *objNameCell=mountModel.alignTable->item(i,2); if(!raCell||!deCell||!objNameCell) return false; QString raString=raCell->text(); QString deString=deCell->text(); QString objString=objNameCell->text(); outstream << "" << endl; outstream << "" << raString << "" << endl; outstream << "" << deString << "" << endl; outstream << "" << objString << "" << endl; outstream << "" << endl; } outstream << "" << endl; appendLogText(i18n("Alignment List saved to %1", path)); file.close(); return true; } void Align::exportSolutionPoints() { if (solutionTable->rowCount() == 0) return; QUrl exportFile = QFileDialog::getSaveFileUrl(KStars::Instance(), i18n("Export Solution Points"), alignURLPath, "CSV File (*.csv)"); if (exportFile.isEmpty()) // if user presses cancel return; if (exportFile.toLocalFile().endsWith(".csv") == false) exportFile.setPath(exportFile.toLocalFile() + ".csv"); QString path=exportFile.toLocalFile(); if (QFile::exists(path)) { int r = KMessageBox::warningContinueCancel(0, i18n( "A file named \"%1\" already exists. " "Overwrite it?", exportFile.fileName() ), i18n( "Overwrite File?" ), KStandardGuiItem::overwrite() ); if(r==KMessageBox::Cancel) return; } if ( !exportFile.isValid() ) { QString message = i18n( "Invalid URL: %1", exportFile.url() ); KMessageBox::sorry(KStars::Instance(), message, i18n( "Invalid URL" ) ); return; } QFile file; file.setFileName(path); if ( !file.open( QIODevice::WriteOnly)) { QString message = i18n( "Unable to write to file %1", path); KMessageBox::sorry( 0, message, i18n( "Could Not Open File" ) ); return; } QTextStream outstream(&file); outstream << "RA,DE,Name,RA Error (arcsec),DE Error (arcsec)"<< endl; for (int i=0; i < solutionTable->rowCount(); i++) { QTableWidgetItem *raCell=solutionTable->item(i,0); QTableWidgetItem *deCell=solutionTable->item(i,1); QTableWidgetItem *objNameCell=solutionTable->item(i,2); QTableWidgetItem *raErrorCell=solutionTable->item(i,4); QTableWidgetItem *deErrorCell=solutionTable->item(i,5); if(!raCell||!deCell||!objNameCell||!raErrorCell||!deErrorCell){ KMessageBox::sorry( 0, i18n( "Error in table structure." ) ); return; } outstream << raCell->text()<<"," << deCell->text()<<"," << objNameCell->text()<<"," << raErrorCell->text().remove("\"")<<"," << deErrorCell->text().remove("\"") << endl; } appendLogText(i18n("Solution Points Saved as: %1", path)); file.close(); } void Align::slotClearAllSolutionPoints() { if (solutionTable->rowCount() == 0) return; if (KMessageBox::questionYesNo(KStars::Instance(), i18n("Are you sure you want to clear all of the solution points?"), i18n("Clear Solution Points"), KStandardGuiItem::yes(), KStandardGuiItem::no()) == KMessageBox::Yes) { solutionTable->setRowCount(0); alignPlot->graph(0)->data()->clear(); alignPlot->clearItems(); buildTarget(); double accuracyRadius=accuracySpin->value(); alignPlot->xAxis->setRange(-accuracyRadius*3, accuracyRadius*3); alignPlot->yAxis->setRange(-accuracyRadius*3, accuracyRadius*3); alignPlot->xAxis->setScaleRatio(alignPlot->yAxis,1.0); alignPlot->replot(); } } void Align::slotClearAllAlignPoints() { if (mountModel.alignTable->rowCount() == 0) return; if (KMessageBox::questionYesNo(KStars::Instance(), i18n("Are you sure you want to clear all the alignment points?"), i18n("Clear Align Points")) == KMessageBox::Yes) mountModel.alignTable->setRowCount(0); + + if(previewShowing) + updatePreviewAlignPoints(); } void Align::slotRemoveSolutionPoint(){ QCPAbstractItem *abstractItem=alignPlot->item(solutionTable->currentRow()); if(abstractItem){ QCPItemText *item = qobject_cast(abstractItem); if(item){ double point=item->position->key(); alignPlot->graph(0)->data()->remove(point); } } alignPlot->removeItem(solutionTable->currentRow()); for(int i=0;iitemCount();i++){ QCPAbstractItem *abstractItem=alignPlot->item(i); if(abstractItem){ QCPItemText *item = qobject_cast(abstractItem); if(item) item->setText(QString::number(i+1)); } } solutionTable->removeRow(solutionTable->currentRow()); alignPlot->replot(); } void Align::slotRemoveAlignPoint(){ mountModel.alignTable->removeRow(mountModel.alignTable->currentRow()); + if(previewShowing) + updatePreviewAlignPoints(); } void Align::moveAlignPoint(int logicalIndex, int oldVisualIndex, int newVisualIndex) { for (int i=0; i < mountModel.alignTable->columnCount(); i++) { QTableWidgetItem *oldItem = mountModel.alignTable->takeItem(oldVisualIndex, i); QTableWidgetItem *newItem = mountModel.alignTable->takeItem(newVisualIndex, i); mountModel.alignTable->setItem(newVisualIndex, i, oldItem); mountModel.alignTable->setItem(oldVisualIndex, i, newItem); } disconnect(mountModel.alignTable->verticalHeader(), SIGNAL(sectionMoved(int, int , int)), this, SLOT(moveAlignPoint(int, int , int))); mountModel.alignTable->verticalHeader()->moveSection(newVisualIndex,oldVisualIndex); connect(mountModel.alignTable->verticalHeader(), SIGNAL(sectionMoved(int, int , int)), this, SLOT(moveAlignPoint(int, int , int))); + if(previewShowing) + updatePreviewAlignPoints(); } void Align::slotMountModel(){ generateAlignStarList(); SkyPoint spWest; spWest.setAlt(30); spWest.setAz(270); spWest.HorizontalToEquatorial(KStars::Instance()->data()->lst(), KStars::Instance()->data()->geo()->lat()); mountModel.alignDec->setValue((int)spWest.dec().Degrees()); mountModelDialog.show(); } void Align::slotAddAlignPoint(){ int currentRow = mountModel.alignTable->rowCount(); mountModel.alignTable->insertRow(currentRow); QTableWidgetItem *disabledBox= new QTableWidgetItem(); disabledBox->setFlags(Qt::ItemIsSelectable); mountModel.alignTable->setItem(currentRow, 3, disabledBox); } void Align::slotFindAlignObject() { QPointer fd = new FindDialog( KStars::Instance() ); if ( fd->exec() == QDialog::Accepted ) { SkyObject *o = fd->targetObject(); if( o != 0 ) { int currentRow = mountModel.alignTable->rowCount(); mountModel.alignTable->insertRow(currentRow); QString ra_report, dec_report; getFormattedCoords(o->ra().Hours(),o->dec().Degrees(),ra_report, dec_report); QTableWidgetItem *RAReport = new QTableWidgetItem(); RAReport->setText(ra_report); RAReport->setTextAlignment(Qt::AlignHCenter); mountModel.alignTable->setItem(currentRow, 0, RAReport); QTableWidgetItem *DECReport = new QTableWidgetItem(); DECReport->setText(dec_report); DECReport->setTextAlignment(Qt::AlignHCenter); mountModel.alignTable->setItem(currentRow, 1, DECReport); QTableWidgetItem *ObjNameReport = new QTableWidgetItem(); ObjNameReport->setText(o->longname()); ObjNameReport->setTextAlignment(Qt::AlignHCenter); mountModel.alignTable->setItem(currentRow, 2, ObjNameReport); QTableWidgetItem *disabledBox= new QTableWidgetItem(); disabledBox->setFlags(Qt::ItemIsSelectable); mountModel.alignTable->setItem(currentRow, 3, disabledBox); } } delete fd; + if(previewShowing) + updatePreviewAlignPoints(); } void Align::resetAlignmentProcedure(){ mountModel.alignTable->setCellWidget(currentAlignmentPoint, 3,new QWidget()); QTableWidgetItem *statusReport = new QTableWidgetItem(); statusReport->setFlags(Qt::ItemIsSelectable); statusReport->setIcon(QIcon(":/icons/AlignWarning.svg")); mountModel.alignTable->setItem(currentAlignmentPoint, 3, statusReport); appendLogText(i18n("The Mount Model Tool is Reset.")); mountModel.startAlignB->setIcon(QIcon::fromTheme("media-playback-start", QIcon(":/icons/breeze/default/media-playback-start.svg") )); mountModelRunning=false; currentAlignmentPoint=0; abort(); } bool Align::alignmentPointsAreBad(){ for (int i=0; i < mountModel.alignTable->rowCount(); i++) { QTableWidgetItem *raCell=mountModel.alignTable->item(i,0); if(!raCell) return true; QString raString=raCell->text(); if(dms().setFromString(raString,false)==false) return true; QTableWidgetItem *decCell=mountModel.alignTable->item(i,1); if(!decCell) return true; QString decString=decCell->text(); if(dms().setFromString(decString,true)==false) return true; } return false; } void Align::startStopAlignmentProcedure(){ if(!mountModelRunning){ if(mountModel.alignTable->rowCount()>0){ if(alignmentPointsAreBad()){ KMessageBox::error(0, i18n("Please Check the Alignment Points.")); return; } if(currentAlignmentPoint==0){ for(int row=0;rowrowCount();row++){ QTableWidgetItem *statusReport = new QTableWidgetItem(); statusReport->setIcon(QIcon()); mountModel.alignTable->setItem(row, 3, statusReport); } } mountModel.startAlignB->setIcon(QIcon::fromTheme("media-playback-pause", QIcon(":/icons/breeze/default/media-playback-pause.svg") )); mountModelRunning=true; appendLogText(i18n("The Mount Model Tool is Starting.")); startAlignmentPoint(); } }else{ mountModel.startAlignB->setIcon(QIcon::fromTheme("media-playback-start", QIcon(":/icons/breeze/default/media-playback-start.svg") )); mountModel.alignTable->setCellWidget(currentAlignmentPoint, 3,new QWidget()); appendLogText(i18n("The Mount Model Tool is Paused.")); abort(); mountModelRunning=false; QTableWidgetItem *statusReport = new QTableWidgetItem(); statusReport->setFlags(Qt::ItemIsSelectable); statusReport->setIcon(QIcon(":/icons/AlignWarning.svg")); mountModel.alignTable->setItem(currentAlignmentPoint, 3, statusReport); } } void Align::startAlignmentPoint(){ if(mountModelRunning&¤tAlignmentPoint>=0&¤tAlignmentPointrowCount()){ QTableWidgetItem *raCell=mountModel.alignTable->item(currentAlignmentPoint,0); QString raString=raCell->text(); dms raDMS = dms::fromString(raString,false); double ra = raDMS.Hours(); QTableWidgetItem *decCell=mountModel.alignTable->item(currentAlignmentPoint,1); QString decString=decCell->text(); dms decDMS = dms::fromString(decString,true); double dec= decDMS.Degrees(); QProgressIndicator *alignIndicator = new QProgressIndicator(this); mountModel.alignTable->setCellWidget(currentAlignmentPoint, 3, alignIndicator); alignIndicator->startAnimation(); targetCoord.setRA(ra); targetCoord.setDec(dec); Slew(); } } void Align::finishAlignmentPoint(bool solverSucceeded){ if(mountModelRunning&¤tAlignmentPoint>=0&¤tAlignmentPointrowCount()){ mountModel.alignTable->setCellWidget(currentAlignmentPoint, 3,new QWidget()); QTableWidgetItem *statusReport = new QTableWidgetItem(); statusReport->setFlags(Qt::ItemIsSelectable); if(solverSucceeded) statusReport->setIcon(QIcon(":/icons/AlignSuccess.svg")); else statusReport->setIcon(QIcon(":/icons/AlignFailure.svg")); mountModel.alignTable->setItem(currentAlignmentPoint, 3, statusReport); currentAlignmentPoint++; if(currentAlignmentPointrowCount()){ startAlignmentPoint(); } else{ mountModelRunning=false; mountModel.startAlignB->setIcon(QIcon::fromTheme("media-playback-start", QIcon(":/icons/breeze/default/media-playback-start.svg") )); appendLogText(i18n("The Mount Model Tool is Finished.")); currentAlignmentPoint=0; } } } bool Align::isParserOK() { bool rc = parser->init(); if (rc) { connect(parser, SIGNAL(solverFinished(double,double,double,double)), this, SLOT(solverFinished(double,double,double,double)), Qt::UniqueConnection); connect(parser, SIGNAL(solverFailed()), this, SLOT(solverFailed()), Qt::UniqueConnection); } return rc; } void Align::checkAlignmentTimeout() { if (loadSlewState != IPS_IDLE || ++solverIterations == MAXIMUM_SOLVER_ITERATIONS) abort(); else if (loadSlewState == IPS_IDLE) { appendLogText(i18n("Solver timed out")); captureAndSolve(); } // TODO must also account for loadAndSlew. Retain file name } void Align::setSolverType(int type) { switch(type) { case SOLVER_ONLINE: loadSlewB->setEnabled(true); if (onlineParser != NULL) { parser = onlineParser; return; } onlineParser = new Ekos::OnlineAstrometryParser(); parser = onlineParser; break; case SOLVER_OFFLINE: loadSlewB->setEnabled(true); if (offlineParser != NULL) { parser = offlineParser; return; } offlineParser = new Ekos::OfflineAstrometryParser(); parser = offlineParser; break; case SOLVER_REMOTE: loadSlewB->setEnabled(true); if (remoteParser != NULL) { parser = remoteParser; (dynamic_cast(parser))->setAstrometryDevice(remoteParserDevice); return; } remoteParser = new Ekos::RemoteAstrometryParser(); parser = remoteParser; (dynamic_cast(parser))->setAstrometryDevice(remoteParserDevice); break; } parser->setAlign(this); if (parser->init()) { connect(parser, SIGNAL(solverFinished(double,double,double, double)), this, SLOT(solverFinished(double,double,double, double)), Qt::UniqueConnection); connect(parser, SIGNAL(solverFailed()), this, SLOT(solverFailed()), Qt::UniqueConnection); } else parser->disconnect(); } bool Align::setCCD(QString device) { for (int i=0; i < CCDCaptureCombo->count(); i++) if (device == CCDCaptureCombo->itemText(i)) { CCDCaptureCombo->setCurrentIndex(i); return true; } return false; } void Align::setDefaultCCD(QString ccd) { Options::setDefaultAlignCCD(ccd); } void Align::checkCCD(int ccdNum) { if (ccdNum == -1) { ccdNum = CCDCaptureCombo->currentIndex(); if (ccdNum == -1) return; } if (ccdNum <= CCDs.count()) { currentCCD = CCDs.at(ccdNum); //if (solverTypeGroup->checkedId() == SOLVER_REMOTE) //(dynamic_cast(parser))->setCCD(currentCCD); } FOVScopeCombo->setCurrentIndex(guideScopeCCDs.contains(currentCCD->getDeviceName()) ? ISD::CCD::TELESCOPE_GUIDE : ISD::CCD::TELESCOPE_PRIMARY); syncCCDInfo(); } void Align::addCCD(ISD::GDInterface *newCCD) { if (CCDs.contains(static_cast(newCCD))) { syncCCDInfo(); return; } CCDs.append(static_cast(newCCD)); CCDCaptureCombo->addItem(newCCD->getDeviceName()); } void Align::setTelescope(ISD::GDInterface *newTelescope) { currentTelescope = static_cast (newTelescope); connect(currentTelescope, SIGNAL(numberUpdated(INumberVectorProperty*)), this, SLOT(processTelescopeNumber(INumberVectorProperty*))); syncTelescopeInfo(); } void Align::syncTelescopeInfo() { if (currentTelescope == NULL) return; canSync = currentTelescope->canSync(); if (canSync == false && syncR->isEnabled()) { slewR->setChecked(true); appendLogText(i18n("Mount does not support syncing.")); } syncR->setEnabled(canSync); INumberVectorProperty * nvp = currentTelescope->getBaseDevice()->getNumber("TELESCOPE_INFO"); double primaryFL=-1, primaryAperture=-1, guideFL=-1, guideAperture=-1; if (nvp) { INumber *np = NULL; np = IUFindNumber(nvp, "TELESCOPE_APERTURE"); if (np && np->value > 0) primaryAperture = np->value; np = IUFindNumber(nvp, "GUIDER_APERTURE"); if (np && np->value > 0) guideAperture = np->value; aperture = primaryAperture; if (currentCCD && Options::guideScopeCCDs().contains(currentCCD->getDeviceName())) aperture = guideAperture; np = IUFindNumber(nvp, "TELESCOPE_FOCAL_LENGTH"); if (np && np->value > 0) primaryFL = np->value; np = IUFindNumber(nvp, "GUIDER_FOCAL_LENGTH"); if (np && np->value > 0) guideFL = np->value; focal_length = primaryFL; if (currentCCD && Options::guideScopeCCDs().contains(currentCCD->getDeviceName())) focal_length = guideFL; } if (focal_length == -1 || aperture == -1) return; if (ccd_hor_pixel != -1 && ccd_ver_pixel != -1 && focal_length != -1 && aperture != -1) { FOVScopeCombo->setItemData(ISD::CCD::TELESCOPE_PRIMARY, i18nc("F-Number, Focal Length, Aperture", "F%1 Focal Length: %2 mm Aperture: %3 mm2", QString::number(primaryFL/primaryAperture, 'f', 1), QString::number(primaryFL, 'f', 2), QString::number(primaryAperture, 'f', 2)), Qt::ToolTipRole); FOVScopeCombo->setItemData(ISD::CCD::TELESCOPE_GUIDE, i18nc("F-Number, Focal Length, Aperture", "F%1 Focal Length: %2 mm Aperture: %3 mm2", QString::number(guideFL/guideAperture, 'f', 1), QString::number(guideFL, 'f', 2), QString::number(guideAperture, 'f', 2)), Qt::ToolTipRole); calculateFOV(); generateArgs(); } } void Align::syncCCDInfo() { INumberVectorProperty * nvp = NULL; if (currentCCD == NULL) return; if (useGuideHead) nvp = currentCCD->getBaseDevice()->getNumber("GUIDER_INFO"); else nvp = currentCCD->getBaseDevice()->getNumber("CCD_INFO"); if (nvp) { INumber *np = IUFindNumber(nvp, "CCD_PIXEL_SIZE_X"); if (np && np->value >0) ccd_hor_pixel = ccd_ver_pixel = np->value; np = IUFindNumber(nvp, "CCD_PIXEL_SIZE_Y"); if (np && np->value >0) ccd_ver_pixel = np->value; np = IUFindNumber(nvp, "CCD_PIXEL_SIZE_Y"); if (np && np->value >0) ccd_ver_pixel = np->value; } ISD::CCDChip *targetChip = currentCCD->getChip(useGuideHead ? ISD::CCDChip::GUIDE_CCD : ISD::CCDChip::PRIMARY_CCD); ISwitchVectorProperty *svp = currentCCD->getBaseDevice()->getSwitch("WCS_CONTROL"); if (svp) setWCSEnabled(Options::astrometrySolverWCS()); targetChip->setImageView(alignView, FITS_ALIGN); targetChip->getFrameMinMax(NULL, NULL, NULL, NULL, NULL, &ccd_width, NULL, &ccd_height); //targetChip->getFrame(&x,&y,&ccd_width,&ccd_height); binningCombo->setEnabled(targetChip->canBin()); if (targetChip->canBin()) { binningCombo->blockSignals(true); int binx=1,biny=1; targetChip->getMaxBin(&binx, &biny); binningCombo->clear(); for (int i=0; i < binx; i++) binningCombo->addItem(QString("%1x%2").arg(i+1).arg(i+1)); binningCombo->setCurrentIndex(Options::solverBinningIndex()); binningCombo->blockSignals(false); } if (ccd_hor_pixel == -1 || ccd_ver_pixel == -1) return; if (ccd_hor_pixel != -1 && ccd_ver_pixel != -1 && focal_length != -1 && aperture != -1) { calculateFOV(); generateArgs(); } } void Align::getFOVScale(double &fov_w, double & fov_h, double &fov_scale) { fov_w = fov_x; fov_h = fov_y; fov_scale = fov_pixscale; } void Align::calculateFOV() { // Calculate FOV // FOV in arcsecs fov_x = 206264.8062470963552 * ccd_width * ccd_hor_pixel / 1000.0 / focal_length; fov_y = 206264.8062470963552 * ccd_height * ccd_ver_pixel / 1000.0 / focal_length; // Pix Scale fov_pixscale = (fov_x * (Options::solverBinningIndex()+1)) / ccd_width; // FOV in arcmins fov_x /= 60.0; fov_y /= 60.0; solverFOV->setSize(fov_x, fov_y); FOVOut->setText(QString("%1' x %2'").arg(fov_x, 0, 'g', 3).arg(fov_y, 0, 'g', 3)); if ( ((fov_x + fov_y) / 2.0) > PAH_CUTOFF_FOV) { PAHWidgets->setEnabled(true); PAHWidgets->setToolTip(QString()); FOVDisabledLabel->hide(); } else { PAHWidgets->setEnabled(false); PAHWidgets->setToolTip(i18n("

Polar Alignment Helper tool requires the following:

1. German Equatorial Mount

2. FOV >" " 0.5 degrees

For small FOVs, use the Legacy Polar Alignment Tool.

")); FOVDisabledLabel->show(); } if (opsAstrometry->kcfg_AstrometryUseImageScale->isChecked()) { int unitType = opsAstrometry->kcfg_AstrometryImageScaleUnits->currentIndex(); // Degrees if (unitType == 0) { double fov_low = qMin(fov_x/60, fov_y/60); double fov_high= qMax(fov_x/60, fov_y/60); opsAstrometry->kcfg_AstrometryImageScaleLow->setValue(fov_low); opsAstrometry->kcfg_AstrometryImageScaleHigh->setValue(fov_high); Options::setAstrometryImageScaleLow(fov_low); Options::setAstrometryImageScaleHigh(fov_high); } // Arcmins else if (unitType == 1) { double fov_low = qMin(fov_x, fov_y); double fov_high= qMax(fov_x, fov_y); opsAstrometry->kcfg_AstrometryImageScaleLow->setValue(fov_low); opsAstrometry->kcfg_AstrometryImageScaleHigh->setValue(fov_high); Options::setAstrometryImageScaleLow(fov_low); Options::setAstrometryImageScaleHigh(fov_high); } // Arcsec per pixel else { opsAstrometry->kcfg_AstrometryImageScaleLow->setValue(fov_pixscale*0.9); opsAstrometry->kcfg_AstrometryImageScaleHigh->setValue(fov_pixscale*1.1); // 10% boundary Options::setAstrometryImageScaleLow(fov_pixscale*0.9); Options::setAstrometryImageScaleHigh(fov_pixscale*1.1); } } } QStringList Align::generateOptions(const QVariantMap & optionsMap) { // -O overwrite // -3 Expected RA // -4 Expected DEC // -5 Radius (deg) // -L lower scale of image in arcminutes // -H upper scale of image in arcmiutes // -u aw set scale to be in arcminutes // -W solution.wcs name of solution file // apog1.jpg name of target file to analyze //solve-field -O -3 06:40:51 -4 +09:49:53 -5 1 -L 40 -H 100 -u aw -W solution.wcs apod1.jpg QStringList solver_args; // Start with always-used arguments solver_args << "-O" << "--no-plots"; // Now go over boolean options // noverify if (optionsMap.contains("noverify")) solver_args << "--no-verify"; // noresort if (optionsMap.contains("resort")) solver_args << "--resort"; // fits2fits if (optionsMap.contains("nofits2fits")) solver_args << "--no-fits2fits"; // downsample if (optionsMap.contains("downsample")) solver_args << "--downsample" << QString::number(optionsMap.value("downsample", 2).toInt()); // image scale low if (optionsMap.contains("scaleL")) solver_args << "-L" << QString::number(optionsMap.value("scaleL").toDouble()); // image scale high if (optionsMap.contains("scaleH")) solver_args << "-H" << QString::number(optionsMap.value("scaleH").toDouble()); // image scale units if (optionsMap.contains("scaleUnits")) solver_args << "-u" << optionsMap.value("scaleUnits").toString(); // RA if (optionsMap.contains("ra")) solver_args << "-3" << QString::number(optionsMap.value("ra").toDouble()); // DE if (optionsMap.contains("de")) solver_args << "-4" << QString::number(optionsMap.value("de").toDouble()); // Radius if (optionsMap.contains("radius")) solver_args << "-5" << QString::number(optionsMap.value("radius").toDouble()); // Custom if (optionsMap.contains("custom")) solver_args << optionsMap.value("custom").toString(); return solver_args; } //This will generate the high and low scale of the imager field size based on the stated units. void Align::generateFOVBounds(double fov_h, double fov_v, QString &fov_low, QString &fov_high) { double fov_lower, fov_upper; // let's stretch the boundaries by 5% fov_lower = ((fov_h < fov_v) ? (fov_h * 0.95) : (fov_v * 0.95)); fov_upper = ((fov_h > fov_v) ? (fov_h * 1.05) : (fov_v * 1.05)); //No need to do anything if they are aw, since that is the default fov_low = QString::number(fov_lower); fov_high = QString::number(fov_upper); } void Align::generateArgs() { // -O overwrite // -3 Expected RA // -4 Expected DEC // -5 Radius (deg) // -L lower scale of image in arcminutes // -H upper scale of image in arcmiutes // -u aw set scale to be in arcminutes // -W solution.wcs name of solution file // apog1.jpg name of target file to analyze //solve-field -O -3 06:40:51 -4 +09:49:53 -5 1 -L 40 -H 100 -u aw -W solution.wcs apod1.jpg QVariantMap optionsMap; if (Options::astrometryUseNoVerify()) optionsMap["noverify"] = true; if (Options::astrometryUseResort()) optionsMap["resort"] = true; if (Options::astrometryUseNoFITS2FITS()) optionsMap["nofits2fits"] = true; if (Options::astrometryUseDownsample()) optionsMap["downsample"] = Options::astrometryDownsample(); if (Options::astrometryUseImageScale()) { QString units = ImageScales[Options::astrometryImageScaleUnits()]; if (Options::astrometryAutoUpdateImageScale()) { QString fov_low,fov_high; double fov_w = fov_x; double fov_h = fov_y; if (units == "dw") { fov_w /= 60; fov_h /= 60; } else if (units == "app") { fov_w = fov_pixscale; fov_h = fov_pixscale; } generateFOVBounds(fov_w, fov_h, fov_low, fov_high); optionsMap["scaleL"] = fov_low; optionsMap["scaleH"] = fov_high; optionsMap["scaleUnits"] = units; } else { optionsMap["scaleL"] = Options::astrometryImageScaleLow(); optionsMap["scaleH"] = Options::astrometryImageScaleHigh(); optionsMap["scaleUnits"] = units; } } if (Options::astrometryUsePosition()) { double ra=0,dec=0; currentTelescope->getEqCoords(&ra, &dec); optionsMap["ra"] = ra*15.0; optionsMap["de"] = dec; optionsMap["radius"] = Options::astrometryRadius(); } if (Options::astrometryCustomOptions().isEmpty() == false) optionsMap["custom"] = Options::astrometryCustomOptions(); QStringList solverArgs = generateOptions(optionsMap); QString options = solverArgs.join(" "); solverOptions->setText(options); solverOptions->setToolTip(options); } bool Align::captureAndSolve() { //m_isSolverComplete = false; if (currentCCD == NULL) return false; if (currentCCD->isConnected() == false) { appendLogText(i18n("Error: Lost connection to CCD.")); KNotification::event( QLatin1String( "AlignFailed"), i18n("Astrometry alignment failed") ); return false; } if (parser->init() == false) return false; if (focal_length == -1 || aperture == -1) { KMessageBox::error(0, i18n("Telescope aperture and focal length are missing. Please check your driver settings and try again.")); return false; } if (ccd_hor_pixel == -1 || ccd_ver_pixel == -1) { KMessageBox::error(0, i18n("CCD pixel size is missing. Please check your driver settings and try again.")); return false; } if (currentFilter != NULL && lockedFilterIndex != -1) { if (lockedFilterIndex != currentFilterIndex) { int lockedFilterPosition = lockedFilterIndex + 1; filterPositionPending = true; currentFilter->runCommand(INDI_SET_FILTER, &lockedFilterPosition); return true; } } if (currentCCD->getDriverInfo()->getClientManager()->getBLOBMode(currentCCD->getDeviceName(), "CCD1") == B_NEVER) { if (KMessageBox::questionYesNo(0, i18n("Image transfer is disabled for this camera. Would you like to enable it?")) == KMessageBox::Yes) { currentCCD->getDriverInfo()->getClientManager()->setBLOBMode(B_ALSO, currentCCD->getDeviceName(), "CCD1"); currentCCD->getDriverInfo()->getClientManager()->setBLOBMode(B_ALSO, currentCCD->getDeviceName(), "CCD2"); } else { return false; } } double seqExpose = exposureIN->value(); ISD::CCDChip *targetChip = currentCCD->getChip(useGuideHead ? ISD::CCDChip::GUIDE_CCD : ISD::CCDChip::PRIMARY_CCD); if (focusState >= FOCUS_PROGRESS) { appendLogText(i18n("Cannot capture while focus module is busy! Retrying...")); QTimer::singleShot(5000, this, SLOT(captureAndSolve())); return false; } if (targetChip->isCapturing()) { appendLogText(i18n("Cannot capture while CCD exposure is in progress! Retrying...")); QTimer::singleShot(1000, this, SLOT(captureAndSolve())); return false; } alignView->setBaseSize(alignWidget->size()); connect(currentCCD, SIGNAL(BLOBUpdated(IBLOB*)), this, SLOT(newFITS(IBLOB*))); connect(currentCCD, SIGNAL(newExposureValue(ISD::CCDChip*,double,IPState)), this, SLOT(checkCCDExposureProgress(ISD::CCDChip*,double,IPState))); // In case of remote solver, we set mode to UPLOAD_BOTH if (solverTypeGroup->checkedId() == SOLVER_REMOTE && remoteParser) { //rememberUploadMode = currentCCD->getUploadMode(); //currentCCD->setUploadMode(ISD::CCD::UPLOAD_BOTH); // For solver remote we need to start solver BEFORE capture //startSolving(QString()); // Enable remote parse dynamic_cast(remoteParser)->setEnabled(true); QString options = solverOptions->text().simplified(); QStringList solverArgs = options.split(" "); dynamic_cast(remoteParser)->sendArgs(solverArgs); if (solverIterations == 0) { double ra,dec; currentTelescope->getEqCoords(&ra, &dec); targetCoord.setRA(ra); targetCoord.setDec(dec); } } //else //{ if (currentCCD->getUploadMode() == ISD::CCD::UPLOAD_LOCAL) { rememberUploadMode = ISD::CCD::UPLOAD_LOCAL; currentCCD->setUploadMode(ISD::CCD::UPLOAD_CLIENT); } // Remove temporary FITS files left before by the solver QDir dir(QDir::tempPath()); dir.setNameFilters(QStringList() << "fits*" << "tmp.*"); dir.setFilter(QDir::Files); foreach(QString dirFile, dir.entryList()) dir.remove(dirFile); //} currentCCD->setTransformFormat(ISD::CCD::FORMAT_FITS); targetChip->resetFrame(); targetChip->setBatchMode(false); targetChip->setCaptureMode(FITS_ALIGN); targetChip->setFrameType(FRAME_LIGHT); int bin = Options::solverBinningIndex()+1; targetChip->setBinning(bin, bin); // In case we're in refresh phase of the polar alignment helper then we use capture value from there if (pahStage == PAH_REFRESH) targetChip->capture(PAHExposure->value()); else targetChip->capture(seqExpose); Options::setAlignExposure(seqExpose); solveB->setEnabled(false); stopB->setEnabled(true); pi->startAnimation(); state = ALIGN_PROGRESS; emit newStatus(state); appendLogText(i18n("Capturing image...")); return true; } void Align::newFITS(IBLOB *bp) { // Ignore guide head if there is any. if (!strcmp(bp->name, "CCD2")) return; disconnect(currentCCD, SIGNAL(BLOBUpdated(IBLOB*)), this, SLOT(newFITS(IBLOB*))); disconnect(currentCCD, SIGNAL(newExposureValue(ISD::CCDChip*,double,IPState)), this, SLOT(checkCCDExposureProgress(ISD::CCDChip*,double,IPState))); appendLogText(i18n("Image received.")); blobType = *(static_cast(bp->aux1)); blobFileName = QString(static_cast(bp->aux2)); if (solverTypeGroup->checkedId() != SOLVER_REMOTE) { if (blobType == ISD::CCD::BLOB_FITS) { ISD::CCDChip *targetChip = currentCCD->getChip(useGuideHead ? ISD::CCDChip::GUIDE_CCD : ISD::CCDChip::PRIMARY_CCD); if (alignDarkFrameCheck->isChecked()) { int x,y,w,h,binx=1,biny=1; targetChip->getFrame(&x,&y,&w,&h); targetChip->getBinning(&binx, &biny); FITSData *darkData = NULL; uint16_t offsetX = x / binx; uint16_t offsetY = y / biny; darkData = DarkLibrary::Instance()->getDarkFrame(targetChip, exposureIN->value()); connect(DarkLibrary::Instance(), SIGNAL(darkFrameCompleted(bool)), this, SLOT(setCaptureComplete())); connect(DarkLibrary::Instance(), SIGNAL(newLog(QString)), this, SLOT(appendLogText(QString))); if (darkData) DarkLibrary::Instance()->subtract(darkData, alignView, FITS_NONE, offsetX, offsetY); else { bool rc = DarkLibrary::Instance()->captureAndSubtract(targetChip, alignView, exposureIN->value(), offsetX, offsetY); alignDarkFrameCheck->setChecked(rc); } return; } } setCaptureComplete(); } } void Align::setCaptureComplete() { DarkLibrary::Instance()->disconnect(this); if (pahStage == PAH_REFRESH) { captureAndSolve(); return; } if (solverTypeGroup->checkedId() == SOLVER_ONLINE && Options::astrometryUseJPEG()) { ISD::CCDChip *targetChip = currentCCD->getChip(useGuideHead ? ISD::CCDChip::GUIDE_CCD : ISD::CCDChip::PRIMARY_CCD); if (targetChip) { QString jpegFile = blobFileName + ".jpg"; bool rc = alignView->getDisplayImage()->save(jpegFile, "JPG"); if (rc) blobFileName = jpegFile; } } if (fov()) fov()->setImageDisplay(alignView->getDisplayImage()); startSolving(blobFileName); } void Align::setGOTOMode(int mode) { gotoModeButtonGroup->button(mode)->setChecked(true); currentGotoMode = static_cast(mode); } void Align::startSolving(const QString &filename, bool isGenerated) { QStringList solverArgs; double ra,dec; QString options = solverOptions->text().simplified(); if (isGenerated) { solverArgs = options.split(" "); // Replace RA and DE with LST & 90/-90 pole if (pahStage == PAH_FIRST_CAPTURE) { for (int i=0; i lst()->Degrees()); // DE. +90 for Northern hemisphere. -90 for southern hemisphere else if (solverArgs[i] == "-4") solverArgs[i+1] = QString::number(hemisphere == NORTH_HEMISPHERE ? 90 : -90); } } } else if (filename.endsWith("fits") || filename.endsWith("fit")) { solverArgs = getSolverOptionsFromFITS(filename); appendLogText(i18n("Using solver options: %1", solverArgs.join(" "))); } else { KGuiItem blindItem(i18n("Blind solver"), QString(), i18n("Blind solver takes a very long time to solve but can reliably solve any image any where in the sky given enough time.")); KGuiItem existingItem(i18n("Use existing settings"), QString(), i18n("Mount must be pointing close to the target location and current field of view must match the image's field of view.")); int rc = KMessageBox::questionYesNoCancel(0, i18n("No metadata is available in this image. Do you want to use the blind solver or the existing solver settings?"), i18n("Astrometry solver"), blindItem, existingItem, KStandardGuiItem::cancel(), "blind_solver_or_existing_solver_option"); if (rc == KMessageBox::Yes) { QVariantMap optionsMap; if (Options::astrometryUseNoVerify()) optionsMap["noverify"] = true; if (Options::astrometryUseResort()) optionsMap["resort"] = true; if (Options::astrometryUseNoFITS2FITS()) optionsMap["nofits2fits"] = true; if (Options::astrometryUseDownsample()) optionsMap["downsample"] = Options::astrometryDownsample(); solverArgs = generateOptions(optionsMap); } else if (rc == KMessageBox::No) solverArgs = options.split(" "); else { abort(); return; } } currentTelescope->getEqCoords(&ra, &dec); if (solverIterations == 0) { targetCoord.setRA(ra); targetCoord.setDec(dec); } //This block of code will create the row in the solution table and populate RA, DE, and object name. //It also starts the progress indicator. if (loadSlewState == IPS_IDLE){ int currentRow = solutionTable->rowCount(); solutionTable->insertRow(currentRow); for(int i=4;i<6;i++){ QTableWidgetItem *disabledBox= new QTableWidgetItem(); disabledBox->setFlags(Qt::ItemIsSelectable); solutionTable->setItem(currentRow, i, disabledBox); } QTableWidgetItem *RAReport = new QTableWidgetItem(); RAReport->setText(ScopeRAOut->text()); RAReport->setTextAlignment(Qt::AlignHCenter); RAReport->setFlags(Qt::ItemIsSelectable); solutionTable->setItem(currentRow, 0, RAReport); QTableWidgetItem *DECReport = new QTableWidgetItem(); DECReport->setText(ScopeDecOut->text()); DECReport->setTextAlignment(Qt::AlignHCenter); DECReport->setFlags(Qt::ItemIsSelectable); solutionTable->setItem(currentRow, 1, DECReport); double maxrad = 1.0; SkyObject *so = KStarsData::Instance()->skyComposite()->objectNearest(new SkyPoint(dms(ra*15),dms(dec)), maxrad ); QString name; if (so){ name = so->longname(); } else{ name = "None"; } QTableWidgetItem *ObjNameReport = new QTableWidgetItem(); ObjNameReport->setText(name); ObjNameReport->setTextAlignment(Qt::AlignHCenter); ObjNameReport->setFlags(Qt::ItemIsSelectable); solutionTable->setItem(currentRow, 2, ObjNameReport); QProgressIndicator *alignIndicator = new QProgressIndicator(this); solutionTable->setCellWidget(currentRow, 3, alignIndicator); alignIndicator->startAnimation(); } Options::setSolverType(solverTypeGroup->checkedId()); //Options::setSolverOptions(solverOptions->text()); Options::setGuideScopeCCDs(guideScopeCCDs); Options::setSolverAccuracyThreshold(accuracySpin->value()); Options::setAlignDarkFrame(alignDarkFrameCheck->isChecked()); Options::setSolverGotoOption(currentGotoMode); //m_isSolverComplete = false; //m_isSolverSuccessful = false; parser->verifyIndexFiles(fov_x, fov_y); solverTimer.start(); alignTimer.start(); if (currentGotoMode == GOTO_SLEW) appendLogText(i18n("Solver iteration #%1", solverIterations+1)); state = ALIGN_PROGRESS; emit newStatus(state); parser->startSovler(filename, solverArgs, isGenerated); } void Align::solverFinished(double orientation, double ra, double dec, double pixscale) { pi->stopAnimation(); stopB->setEnabled(false); solveB->setEnabled(true); sOrientation = orientation; sRA = ra; sDEC = dec; alignTimer.stop(); if (solverTypeGroup->checkedId() == SOLVER_REMOTE && remoteParser) { // Disable remote parse dynamic_cast(remoteParser)->setEnabled(false); } int binx, biny; ISD::CCDChip *targetChip = currentCCD->getChip(useGuideHead ? ISD::CCDChip::GUIDE_CCD : ISD::CCDChip::PRIMARY_CCD); targetChip->getBinning(&binx, &biny); if (Options::astrometrySolverVerbose()) appendLogText(i18n("Solver RA (%1) DEC (%2) Orientation (%3) Pixel Scale (%4)", QString::number(ra, 'g' , 5), QString::number(dec, 'g' , 5), QString::number(orientation, 'g' , 5), QString::number(pixscale, 'g' , 5))); if (pixscale > 0 && loadSlewState == IPS_IDLE) { double solver_focal_length = (206.264 * ccd_hor_pixel) / pixscale * binx; if (fabs(focal_length - solver_focal_length) > 1) appendLogText(i18n("Current focal length is %1 mm while computed focal length from the solver is %2 mm. Please update the mount focal length to obtain accurate results.", QString::number(focal_length, 'g' , 5), QString::number(solver_focal_length, 'g' , 5))); } alignCoord.setRA0(ra/15.0); alignCoord.setDec0(dec); RotOut->setText(QString::number(orientation, 'g', 5)); // Convert to JNow alignCoord.apparentCoord((long double) J2000, KStars::Instance()->data()->ut().djd()); // Get horizontal coords alignCoord.EquatorialToHorizontal(KStarsData::Instance()->lst(), KStarsData::Instance()->geo()->lat()); double raDiff = (alignCoord.ra().Degrees()-targetCoord.ra().Degrees()) * 3600; double deDiff = (alignCoord.dec().Degrees()-targetCoord.dec().Degrees()) * 3600; dms RADiff(fabs(raDiff)/3600.0), DEDiff(deDiff/3600.0); dRAOut->setText(QString("%1%2").arg(raDiff > 0 ? "+" : "-").arg(RADiff.toHMSString())); dDEOut->setText(DEDiff.toDMSString(true)); pixScaleOut->setText(QString::number(pixscale, 'f', 2)); emit newSolutionDeviation(raDiff, deDiff); targetDiff = sqrt(raDiff*raDiff + deDiff*deDiff); solverFOV->setCenter(alignCoord); solverFOV->setRotation(sOrientation); solverFOV->setImageDisplay(Options::astrometrySolverOverlay()); QString ra_dms, dec_dms; getFormattedCoords(alignCoord.ra().Hours(), alignCoord.dec().Degrees(), ra_dms, dec_dms); SolverRAOut->setText(ra_dms); SolverDecOut->setText(dec_dms); //This block of code will write the result into the solution table and plot it on the graph. int currentRow = solutionTable->rowCount()-1; if (loadSlewState == IPS_IDLE){ QTableWidgetItem *dRAReport = new QTableWidgetItem(); if(dRAReport){ dRAReport->setText(QString::number(raDiff,'f',3)+"\""); dRAReport->setTextAlignment(Qt::AlignHCenter); dRAReport->setFlags(Qt::ItemIsSelectable); solutionTable->setItem(currentRow, 4, dRAReport); } QTableWidgetItem *dDECReport = new QTableWidgetItem(); if(dDECReport){ dDECReport->setText(QString::number(deDiff,'f',3)+"\""); dDECReport->setTextAlignment(Qt::AlignHCenter); dDECReport->setFlags(Qt::ItemIsSelectable); solutionTable->setItem(currentRow, 5, dDECReport); } double raPlot=raDiff; double decPlot=deDiff; alignPlot->graph(0)->addData(raPlot, decPlot); QCPItemText *textLabel = new QCPItemText(alignPlot); textLabel->setPositionAlignment(Qt::AlignVCenter|Qt::AlignHCenter); textLabel->position->setType(QCPItemPosition::ptPlotCoords); textLabel->position->setCoords(raPlot,decPlot); textLabel->setColor(Qt::red); textLabel->setPadding(QMargins(0,0,0,0)); textLabel->setBrush(Qt::white); //textLabel->setBrush(Qt::NoBrush); textLabel->setPen(Qt::NoPen); textLabel->setText(" " + QString::number(solutionTable->rowCount()) + " "); textLabel->setFont(QFont(font().family(), 8)); if(!alignPlot->xAxis->range().contains(raDiff)){ alignPlot->graph(0)->rescaleKeyAxis(true); alignPlot->yAxis->setScaleRatio(alignPlot->xAxis,1.0); } if(!alignPlot->yAxis->range().contains(deDiff)){ alignPlot->graph(0)->rescaleValueAxis(true); alignPlot->xAxis->setScaleRatio(alignPlot->yAxis,1.0); } alignPlot->replot(); } if (Options::astrometrySolverWCS()) { INumberVectorProperty *ccdRotation = currentCCD->getBaseDevice()->getNumber("CCD_ROTATION"); if (ccdRotation) { INumber *rotation = IUFindNumber(ccdRotation, "CCD_ROTATION_VALUE"); if (rotation) { ClientManager *clientManager = currentCCD->getDriverInfo()->getClientManager(); rotation->value = orientation; clientManager->sendNewNumber(ccdRotation); if (m_wcsSynced == false) { appendLogText(i18n("WCS information updated. Images captured from this point forward shall have valid WCS.")); // Just send telescope info in case the CCD driver did not pick up before. INumberVectorProperty *telescopeInfo = currentTelescope->getBaseDevice()->getNumber("TELESCOPE_INFO"); if (telescopeInfo) clientManager->sendNewNumber(telescopeInfo); m_wcsSynced=true; } } } } retries=0; appendLogText(i18n("Solution coordinates: RA (%1) DEC (%2) Telescope Coordinates: RA (%3) DEC (%4)", alignCoord.ra().toHMSString(), alignCoord.dec().toDMSString(), telescopeCoord.ra().toHMSString(), telescopeCoord.dec().toDMSString())); if (loadSlewState == IPS_IDLE && currentGotoMode == GOTO_SLEW) { dms diffDeg(targetDiff/3600.0); appendLogText(i18n("Target is within %1 degrees of solution coordinates.", diffDeg.toDMSString())); } if (rememberUploadMode != currentCCD->getUploadMode()) currentCCD->setUploadMode(rememberUploadMode); //if (syncR->isChecked() || nothingR->isChecked() || targetDiff <= accuracySpin->value()) // CONTINUE HERE //This block of code along with some sections in the switch below will set the status report in the solution table for this item. QTableWidgetItem *statusReport = new QTableWidgetItem(); if (loadSlewState == IPS_IDLE){ solutionTable->setCellWidget(currentRow, 3,new QWidget()); statusReport->setFlags(Qt::ItemIsSelectable); } switch (currentGotoMode) { case GOTO_SYNC: executeGOTO(); if (loadSlewState == IPS_IDLE){ statusReport->setIcon(QIcon(":/icons/AlignSuccess.svg")); solutionTable->setItem(currentRow, 3, statusReport); } if(mountModelRunning) finishAlignmentPoint(true); return; break; case GOTO_SLEW: if (loadSlewState == IPS_BUSY || targetDiff > accuracySpin->value()) { if (loadSlewState == IPS_IDLE && ++solverIterations == MAXIMUM_SOLVER_ITERATIONS) { appendLogText(i18n("Maximum number of iterations reached. Solver failed.")); if (loadSlewState == IPS_IDLE){ statusReport->setIcon(QIcon(":/icons/AlignFailure.svg")); solutionTable->setItem(currentRow, 3, statusReport); } solverFailed(); if(mountModelRunning) finishAlignmentPoint(false); return; } targetAccuracyNotMet=true; if (loadSlewState == IPS_IDLE){ statusReport->setIcon(QIcon(":/icons/AlignWarning.svg")); solutionTable->setItem(currentRow, 3, statusReport); } executeGOTO(); return; } if (loadSlewState == IPS_IDLE){ statusReport->setIcon(QIcon(":/icons/AlignSuccess.svg")); solutionTable->setItem(currentRow, 3, statusReport); } appendLogText(i18n("Target is within acceptable range. Astrometric solver is successful.")); if(mountModelRunning){ finishAlignmentPoint(true); if(mountModelRunning) return; } break; case GOTO_NOTHING: if (loadSlewState == IPS_IDLE){ statusReport->setIcon(QIcon(":/icons/AlignSuccess.svg")); solutionTable->setItem(currentRow, 3, statusReport); } if(mountModelRunning){ finishAlignmentPoint(true); if(mountModelRunning) return; } break; } KNotification::event( QLatin1String( "AlignSuccessful"), i18n("Astrometry alignment completed successfully") ); state = ALIGN_COMPLETE; emit newStatus(state); solverIterations=0; if (pahStage != PAH_IDLE) processPAHStage(orientation, ra, dec, pixscale); else if (azStage > AZ_INIT || altStage > ALT_INIT) executePolarAlign(); else { solveB->setEnabled(true); loadSlewB->setEnabled(true); } } void Align::solverFailed() { KNotification::event( QLatin1String( "AlignFailed"), i18n("Astrometry alignment failed with errors") ); pi->stopAnimation(); stopB->setEnabled(false); solveB->setEnabled(true); alignTimer.stop(); azStage = AZ_INIT; altStage = ALT_INIT; //loadSlewMode = false; loadSlewState=IPS_IDLE; PAHFirstCaptureB->setEnabled(true); PAHSecondCaptureB->setEnabled(true); PAHThirdCaptureB->setEnabled(true); solverIterations=0; retries=0; //emit solverComplete(false); state = ALIGN_FAILED; emit newStatus(state); int currentRow = solutionTable->rowCount()-1; solutionTable->setCellWidget(currentRow, 3,new QWidget()); QTableWidgetItem *statusReport = new QTableWidgetItem(); statusReport->setIcon(QIcon(":/icons/AlignFailure.svg")); statusReport->setFlags(Qt::ItemIsSelectable); solutionTable->setItem(currentRow, 3, statusReport); } void Align::abort() { parser->stopSolver(); pi->stopAnimation(); stopB->setEnabled(false); solveB->setEnabled(true); loadSlewB->setEnabled(true); azStage = AZ_INIT; altStage = ALT_INIT; //loadSlewMode = false; loadSlewState=IPS_IDLE; PAHFirstCaptureB->setEnabled(true); PAHSecondCaptureB->setEnabled(true); PAHThirdCaptureB->setEnabled(true); solverIterations=0; retries=0; alignTimer.stop(); //currentCCD->disconnect(this); disconnect(currentCCD, SIGNAL(BLOBUpdated(IBLOB*)), this, SLOT(newFITS(IBLOB*))); disconnect(currentCCD, SIGNAL(newExposureValue(ISD::CCDChip*,double,IPState)), this, SLOT(checkCCDExposureProgress(ISD::CCDChip*,double,IPState))); if (rememberUploadMode != currentCCD->getUploadMode()) currentCCD->setUploadMode(rememberUploadMode); ISD::CCDChip *targetChip = currentCCD->getChip(useGuideHead ? ISD::CCDChip::GUIDE_CCD : ISD::CCDChip::PRIMARY_CCD); // If capture is still in progress, let's stop that. if (pahStage == PAH_REFRESH) { if (targetChip->isCapturing()) targetChip->abortExposure(); appendLogText(i18n("Refresh is complete.")); } else { if (targetChip->isCapturing()) { targetChip->abortExposure(); appendLogText(i18n("Capture aborted.")); } else { int elapsed = (int) round(solverTimer.elapsed()/1000.0); appendLogText(i18np("Solver aborted after %1 second.", "Solver aborted after %1 seconds", elapsed)); } } state = ALIGN_ABORTED; emit newStatus(state); int currentRow = solutionTable->rowCount()-1; solutionTable->setCellWidget(currentRow, 3,new QWidget()); QTableWidgetItem *statusReport = new QTableWidgetItem(); statusReport->setIcon(QIcon(":/icons/AlignFailure.svg")); statusReport->setFlags(Qt::ItemIsSelectable); solutionTable->setItem(currentRow, 3, statusReport); } QList Align::getSolutionResult() { QList result; result << sOrientation << sRA << sDEC; return result; } void Align::appendLogText(const QString &text) { logText.insert(0, i18nc("log entry; %1 is the date, %2 is the text", "%1 %2", QDateTime::currentDateTime().toString("yyyy-MM-ddThh:mm:ss"), text)); if (Options::alignmentLogging()) qDebug() << "Alignment: " << text; emit newLog(); } void Align::clearLog() { logText.clear(); emit newLog(); } void Align::processTelescopeNumber(INumberVectorProperty *coord) { QString ra_dms, dec_dms; if (!strcmp(coord->name, "EQUATORIAL_EOD_COORD")) { getFormattedCoords(coord->np[0].value, coord->np[1].value, ra_dms, dec_dms); telescopeCoord.setRA(coord->np[0].value); telescopeCoord.setDec(coord->np[1].value); telescopeCoord.EquatorialToHorizontal(KStarsData::Instance()->lst(), KStarsData::Instance()->geo()->lat()); ScopeRAOut->setText(ra_dms); ScopeDecOut->setText(dec_dms); switch (coord->s) { case IPS_OK: { // Update the boxes as the mount just finished slewing if (isSlewDirty && Options::astrometryAutoUpdatePosition()) { opsAstrometry->estRA->setText(ra_dms); opsAstrometry->estDec->setText(dec_dms); Options::setAstrometryPositionRA(coord->np[0].value*15); Options::setAstrometryPositionDE(coord->np[1].value); generateArgs(); } if (isSlewDirty && pahStage == PAH_FIND_CP) { isSlewDirty = false; appendLogText(i18n("Mount completed slewing near celestial pole. Capture again to verify.")); PAHFirstCaptureB->setEnabled(true); setGOTOMode(GOTO_NOTHING); pahStage = PAH_FIRST_CAPTURE; return; } if (isSlewDirty && pahStage == PAH_FIRST_ROTATE) { isSlewDirty = false; appendLogText(i18n("Mount first rotation is complete.")); pahStage = PAH_SECOND_CAPTURE; PAHWidgets->setCurrentWidget(PAHSecondCapturePage); if (PAHAutoModeCheck->isChecked()) { PAHSecondCaptureB->setEnabled(true); PAHSecondCaptureB->click(); } } else if (isSlewDirty && pahStage == PAH_SECOND_ROTATE) { isSlewDirty = false; appendLogText(i18n("Mount second rotation is complete.")); pahStage = PAH_THIRD_CAPTURE; PAHWidgets->setCurrentWidget(PAHThirdCapturePage); if (PAHAutoModeCheck->isChecked()) { PAHThirdCaptureB->setEnabled(true); PAHThirdCaptureB->click(); } } switch (state) { case ALIGN_PROGRESS: break; case ALIGN_SYNCING: { isSlewDirty = false; if (currentGotoMode == GOTO_SLEW) { Slew(); return; } else { appendLogText(i18n("Mount is synced to solution coordinates. Astrometric solver is successful.")); KNotification::event( QLatin1String( "AlignSuccessful"), i18n("Astrometry alignment completed successfully") ); state = ALIGN_COMPLETE; emit newStatus(state); } } break; case ALIGN_SLEWING: if (isSlewDirty == false) break; isSlewDirty = false; if (loadSlewState == IPS_BUSY) { loadSlewState = IPS_IDLE; if (Options::alignmentLogging()) qDebug() << "Alignment: loadSlewState is IDLE."; state = ALIGN_PROGRESS; emit newStatus(state); QTimer::singleShot(delaySpin->value(), this, SLOT(captureAndSolve())); return; } else if(currentGotoMode == GOTO_SLEW||mountModelRunning) { if(targetAccuracyNotMet) appendLogText(i18n("Slew complete. Target accuracy is not met, running solver again...")); else appendLogText(i18n("Slew complete. Solving Alignment Point. . .")); targetAccuracyNotMet=false; state = ALIGN_PROGRESS; emit newStatus(state); QTimer::singleShot(delaySpin->value(), this, SLOT(captureAndSolve())); return; } break; default: { isSlewDirty = false; } break; } } break; case IPS_BUSY: { isSlewDirty = true; } break; case IPS_ALERT: { if (state == ALIGN_SYNCING || state == ALIGN_SLEWING) { if (state == ALIGN_SYNCING) appendLogText(i18n("Syncing failed!")); else appendLogText(i18n("Slewing failed!")); if (++retries == 3) { abort(); return; } else { if (currentGotoMode == GOTO_SLEW) Slew(); else Sync(); } } return; } break; default: break; } /*if (Options::alignmentLogging()) qDebug() << "Alignment: State is " << Ekos::getAlignStatusString(state) << " isSlewing? " << currentTelescope->isSlewing() << " slew Dirty? " << slew_dirty << " Current GOTO Mode? " << currentGotoMode << " LoadSlewState? " << pstateStr(loadSlewState);*/ switch (azStage) { case AZ_SYNCING: if (currentTelescope->isSlewing()) azStage=AZ_SLEWING; break; case AZ_SLEWING: if (currentTelescope->isSlewing() == false) { azStage = AZ_SECOND_TARGET; measureAzError(); } break; case AZ_CORRECTING: if (currentTelescope->isSlewing() == false) { appendLogText(i18n("Slew complete. Please adjust azimuth knob until the target is in the center of the view.")); azStage = AZ_INIT; } break; default: break; } switch (altStage) { case ALT_SYNCING: if (currentTelescope->isSlewing()) altStage = ALT_SLEWING; break; case ALT_SLEWING: if (currentTelescope->isSlewing() == false) { altStage = ALT_SECOND_TARGET; measureAltError(); } break; case ALT_CORRECTING: if (currentTelescope->isSlewing() == false) { appendLogText(i18n("Slew complete. Please adjust altitude knob until the target is in the center of the view.")); altStage = ALT_INIT; } break; default: break; } } // N.B. EkosManager already mananges TELESCOPE_INFO, why here again? //if (!strcmp(coord->name, "TELESCOPE_INFO")) //syncTelescopeInfo(); } void Align::executeGOTO() { if (loadSlewState == IPS_BUSY) { //if (loadSlewIterations == loadSlewIterationsSpin->value()) //loadSlewCoord = alignCoord; //targetCoord = loadSlewCoord; targetCoord = alignCoord; SlewToTarget(); } else if (currentGotoMode == GOTO_SYNC) Sync(); else if (currentGotoMode == GOTO_SLEW) SlewToTarget(); } void Align::Sync() { state = ALIGN_SYNCING; if (currentTelescope->Sync(&alignCoord)) { emit newStatus(state); appendLogText(i18n("Syncing to RA (%1) DEC (%2)", alignCoord.ra().toHMSString(), alignCoord.dec().toDMSString())); } else { state = ALIGN_IDLE; emit newStatus(state); appendLogText(i18n("Syncing failed.")); } } void Align::Slew() { state = ALIGN_SLEWING; emit newStatus(state); isSlewDirty = currentTelescope->Slew(&targetCoord); appendLogText(i18n("Slewing to target coordinates: RA (%1) DEC (%2).", targetCoord.ra().toHMSString(), targetCoord.dec().toDMSString())); } void Align::SlewToTarget() { if (canSync && loadSlewState == IPS_IDLE) { Sync(); return; } Slew(); } void Align::executePolarAlign() { appendLogText(i18n("Processing solution for polar alignment...")); switch (azStage) { case AZ_FIRST_TARGET: case AZ_FINISHED: measureAzError(); break; default: break; } switch (altStage) { case ALT_FIRST_TARGET: case ALT_FINISHED: measureAltError(); break; default: break; } } void Align::measureAzError() { static double initRA=0, initDEC=0, finalRA=0, finalDEC=0, initAz=0; if (pahStage != PAH_IDLE && (KMessageBox::warningContinueCancel(KStars::Instance(), i18n("Polar Alignment Helper is still active. Do you want to continue using legacy polar alignment tool?"))!=KMessageBox::Continue)) return; pahStage = PAH_IDLE; if (Options::alignmentLogging()) qDebug() << "Polar Alignment: Measureing Azimuth Error..."; switch (azStage) { case AZ_INIT: // Display message box confirming user point scope near meridian and south if (KMessageBox::warningContinueCancel( 0, hemisphere == NORTH_HEMISPHERE ? i18n("Point the telescope at the southern meridian. Press continue when ready.") : i18n("Point the telescope at the northern meridian. Press continue when ready.") , i18n("Polar Alignment Measurement"), KStandardGuiItem::cont(), KStandardGuiItem::cancel(), "ekos_measure_az_error")!=KMessageBox::Continue) return; appendLogText(i18n("Solving first frame near the meridian.")); azStage = AZ_FIRST_TARGET; captureAndSolve(); break; case AZ_FIRST_TARGET: // start solving there, find RA/DEC initRA = alignCoord.ra().Degrees(); initDEC = alignCoord.dec().Degrees(); initAz = alignCoord.az().Degrees(); if (Options::alignmentLogging()) qDebug() << "Polar Alignment: initRA " << alignCoord.ra().toHMSString() << " initDEC " << alignCoord.dec().toDMSString() << " initlAz " << alignCoord.az().toDMSString() << " initAlt " << alignCoord.alt().toDMSString(); // Now move 30 arcminutes in RA if (canSync) { azStage = AZ_SYNCING; currentTelescope->Sync(initRA/15.0, initDEC); currentTelescope->Slew((initRA - RAMotion)/15.0, initDEC); } // If telescope doesn't sync, we slew relative to its current coordinates else { azStage = AZ_SLEWING; currentTelescope->Slew(telescopeCoord.ra().Hours() - RAMotion/15.0, telescopeCoord.dec().Degrees()); } appendLogText(i18n("Slewing 30 arcminutes in RA...")); break; case AZ_SECOND_TARGET: // We reached second target now // Let now solver for RA/DEC appendLogText(i18n("Solving second frame near the meridian.")); azStage = AZ_FINISHED; captureAndSolve(); break; case AZ_FINISHED: // Measure deviation in DEC // Call function to report error // set stage to AZ_FIRST_TARGET again appendLogText(i18n("Calculating azimuth alignment error...")); finalRA = alignCoord.ra().Degrees(); finalDEC = alignCoord.dec().Degrees(); if (Options::alignmentLogging()) qDebug() << "Polar Alignment: finalRA " << alignCoord.ra().toHMSString() << " finalDEC " << alignCoord.dec().toDMSString() << " finalAz " << alignCoord.az().toDMSString() << " finalAlt " << alignCoord.alt().toDMSString(); // Slew back to original position if (canSync) currentTelescope->Slew(initRA/15.0, initDEC); else { currentTelescope->Slew(telescopeCoord.ra().Hours() + RAMotion/15.0, telescopeCoord.dec().Degrees()); } appendLogText(i18n("Slewing back to original position...")); calculatePolarError(initRA, initDEC, finalRA, finalDEC, initAz); azStage = AZ_INIT; break; default: break; } } void Align::measureAltError() { static double initRA=0, initDEC=0, finalRA=0, finalDEC=0, initAz=0; if (pahStage != PAH_IDLE && (KMessageBox::warningContinueCancel(KStars::Instance(), i18n("Polar Alignment Helper is still active. Do you want to continue using legacy polar alignment tool?"))!=KMessageBox::Continue)) return; pahStage = PAH_IDLE; if (Options::alignmentLogging()) qDebug() << "Polar Alignment: Measureing Altitude Error..."; switch (altStage) { case ALT_INIT: // Display message box confirming user point scope near meridian and south if (KMessageBox::warningContinueCancel( 0, i18n("Point the telescope to the eastern or western horizon with a minimum altitude of 20 degrees. Press continue when ready.") , i18n("Polar Alignment Measurement"), KStandardGuiItem::cont(), KStandardGuiItem::cancel(), "ekos_measure_alt_error")!=KMessageBox::Continue) return; appendLogText(i18n("Solving first frame.")); altStage = ALT_FIRST_TARGET; captureAndSolve(); break; case ALT_FIRST_TARGET: // start solving there, find RA/DEC initRA = alignCoord.ra().Degrees(); initDEC = alignCoord.dec().Degrees(); initAz = alignCoord.az().Degrees(); if (Options::alignmentLogging()) qDebug() << "Polar Alignment: initRA " << alignCoord.ra().toHMSString() << " initDEC " << alignCoord.dec().toDMSString() << " initlAz " << alignCoord.az().toDMSString() << " initAlt " << alignCoord.alt().toDMSString(); // Now move 30 arcminutes in RA if (canSync) { altStage = ALT_SYNCING; currentTelescope->Sync(initRA/15.0, initDEC); currentTelescope->Slew((initRA - RAMotion)/15.0, initDEC); } // If telescope doesn't sync, we slew relative to its current coordinates else { altStage = ALT_SLEWING; currentTelescope->Slew(telescopeCoord.ra().Hours() - RAMotion/15.0, telescopeCoord.dec().Degrees()); } appendLogText(i18n("Slewing 30 arcminutes in RA...")); break; case ALT_SECOND_TARGET: // We reached second target now // Let now solver for RA/DEC appendLogText(i18n("Solving second frame.")); altStage = ALT_FINISHED; captureAndSolve(); break; case ALT_FINISHED: // Measure deviation in DEC // Call function to report error appendLogText(i18n("Calculating altitude alignment error...")); finalRA = alignCoord.ra().Degrees(); finalDEC = alignCoord.dec().Degrees(); if (Options::alignmentLogging()) qDebug() << "Polar Alignment: finalRA " << alignCoord.ra().toHMSString() << " finalDEC " << alignCoord.dec().toDMSString() << " finalAz " << alignCoord.az().toDMSString() << " finalAlt " << alignCoord.alt().toDMSString(); // Slew back to original position if (canSync) currentTelescope->Slew(initRA/15.0, initDEC); // If telescope doesn't sync, we slew relative to its current coordinates else { currentTelescope->Slew(telescopeCoord.ra().Hours() + RAMotion/15.0, telescopeCoord.dec().Degrees()); } appendLogText(i18n("Slewing back to original position...")); calculatePolarError(initRA, initDEC, finalRA, finalDEC, initAz); altStage = ALT_INIT; break; default: break; } } void Align::calculatePolarError(double initRA, double initDEC, double finalRA, double finalDEC, double initAz) { double raMotion = finalRA - initRA; decDeviation = finalDEC - initDEC; // East/West of meridian int horizon = (initAz > 0 && initAz <= 180) ? 0 : 1; // How much time passed siderrally form initRA to finalRA? //double RATime = fabs(raMotion / SIDRATE) / 60.0; // 2016-03-30: Diff in RA is sufficient for time difference // raMotion in degrees. RATime in minutes. double RATime = fabs(raMotion) * 60.0; // Equation by Frank Berret (Measuring Polar Axis Alignment Error, page 4) // In degrees double deviation = (3.81 * (decDeviation * 3600) ) / ( RATime * cos(initDEC * dms::DegToRad)) / 60.0; dms devDMS(fabs(deviation)); KLocalizedString deviationDirection; switch (hemisphere) { // Northern hemisphere case NORTH_HEMISPHERE: if (azStage == AZ_FINISHED) { if (decDeviation > 0) deviationDirection = ki18n("%1 too far east"); else deviationDirection = ki18n("%1 too far west"); } else if (altStage == ALT_FINISHED) { switch (horizon) { // East case 0: if (decDeviation > 0) deviationDirection = ki18n("%1 too far high"); else deviationDirection = ki18n("%1 too far low"); break; // West case 1: if (decDeviation > 0) deviationDirection = ki18n("%1 too far low"); else deviationDirection = ki18n("%1 too far high"); break; default: break; } } break; // Southern hemisphere case SOUTH_HEMISPHERE: if (azStage == AZ_FINISHED) { if (decDeviation > 0) deviationDirection = ki18n("%1 too far west"); else deviationDirection = ki18n("%1 too far east"); } else if (altStage == ALT_FINISHED) { switch (horizon) { // East case 0: if (decDeviation > 0) deviationDirection = ki18n("%1 too far low"); else deviationDirection = ki18n("%1 too far high"); break; // West case 1: if (decDeviation > 0) deviationDirection = ki18n("%1 too far high"); else deviationDirection = ki18n("%1 too far low"); break; default: break; } } break; default: break; } if (Options::alignmentLogging()) { qDebug() << "Polar Alignment: Hemisphere is " << ((hemisphere == NORTH_HEMISPHERE) ? "North" : "South") << " --- initAz " << initAz; qDebug() << "Polar Alignment: initRA " << initRA << " initDEC " << initDEC << " finalRA " << finalRA << " finalDEC " << finalDEC; qDebug() << "Polar Alignment: decDeviation " << decDeviation*3600 << " arcsec " << " RATime " << RATime << " minutes"; qDebug() << "Polar Alignment: Raw Deviaiton " << deviation << " degrees."; } if (azStage == AZ_FINISHED) { azError->setText(deviationDirection.subs(QString("%1").arg(devDMS.toDMSString())).toString()); //azError->setText(deviationDirection.subs(QString("%1")azDMS.toDMSString()); azDeviation = deviation * (decDeviation > 0 ? 1 : -1); if (Options::alignmentLogging()) qDebug() << "Polar Alignment: Azimuth Deviation " << azDeviation << " degrees."; correctAzB->setEnabled(true); } if (altStage == ALT_FINISHED) { //altError->setText(deviationDirection.subs(QString("%1").arg(fabs(deviation), 0, 'g', 3)).toString()); altError->setText(deviationDirection.subs(QString("%1").arg(devDMS.toDMSString())).toString()); altDeviation = deviation * (decDeviation > 0 ? 1 : -1); if (Options::alignmentLogging()) qDebug() << "Polar Alignment: Altitude Deviation " << altDeviation << " degrees."; correctAltB->setEnabled(true); } } void Align::correctAltError() { double newRA, newDEC; SkyPoint currentCoord (telescopeCoord); dms targetLat; if (Options::alignmentLogging()) { qDebug() << "Polar Alignment: Correcting Altitude Error..."; qDebug() << "Polar Alignment: Current Mount RA " << currentCoord.ra().toHMSString() << " DEC " << currentCoord.dec().toDMSString() << "Az " << currentCoord.az().toDMSString() << " Alt " << currentCoord.alt().toDMSString(); } // An error in polar alignment altitude reflects a deviation in the latitude of the mount from actual latitude of the site // Calculating the latitude accounting for the altitude deviation. This is the latitude at which the altitude deviation should be zero. targetLat.setD(KStars::Instance()->data()->geo()->lat()->Degrees() + altDeviation); // Calculate the Az/Alt of the mount if it were located at the corrected latitude currentCoord.EquatorialToHorizontal(KStars::Instance()->data()->lst(), &targetLat ); // Convert corrected Az/Alt to RA/DEC given the local sideral time and current (not corrected) latitude currentCoord.HorizontalToEquatorial(KStars::Instance()->data()->lst(), KStars::Instance()->data()->geo()->lat()); // New RA/DEC should reflect the position in the sky at which the polar alignment altitude error is minimal. newRA = currentCoord.ra().Hours(); newDEC = currentCoord.dec().Degrees(); altStage = ALT_CORRECTING; if (Options::alignmentLogging()) { qDebug() << "Polar Alignment: Target Latitude = Latitude " << KStars::Instance()->data()->geo()->lat()->Degrees() << " + Altitude Deviation " << altDeviation << " = " << targetLat.Degrees(); qDebug() << "Polar Alignment: Slewing to calibration position..."; } currentTelescope->Slew(newRA, newDEC); appendLogText(i18n("Slewing to calibration position, please wait until telescope completes slewing.")); } void Align::correctAzError() { double newRA, newDEC, currentAlt, currentAz; SkyPoint currentCoord (telescopeCoord); if (Options::alignmentLogging()) { qDebug() << "Polar Alignment: Correcting Azimuth Error..."; qDebug() << "Polar Alignment: Current Mount RA " << currentCoord.ra().toHMSString() << " DEC " << currentCoord.dec().toDMSString() << "Az " << currentCoord.az().toDMSString() << " Alt " << currentCoord.alt().toDMSString(); qDebug() << "Polar Alignment: Target Azimuth = Current Azimuth " << currentCoord.az().Degrees() << " + Azimuth Deviation " << azDeviation << " = " << currentCoord.az().Degrees() + azDeviation; } // Get current horizontal coordinates of the mount currentCoord.EquatorialToHorizontal(KStars::Instance()->data()->lst(), KStars::Instance()->data()->geo()->lat()); // Keep Altitude as it is and change Azimuth to account for the azimuth deviation // The new sky position should be where the polar alignment azimuth error is minimal currentAlt = currentCoord.alt().Degrees(); currentAz = currentCoord.az().Degrees() + azDeviation; // Update current Alt and Azimuth to new values currentCoord.setAlt(currentAlt); currentCoord.setAz(currentAz); // Convert Alt/Az back to equatorial coordinates currentCoord.HorizontalToEquatorial(KStars::Instance()->data()->lst(), KStars::Instance()->data()->geo()->lat()); // Get new RA and DEC newRA = currentCoord.ra().Hours(); newDEC = currentCoord.dec().Degrees(); azStage = AZ_CORRECTING; if (Options::alignmentLogging()) qDebug() << "Polar Alignment: Slewing to calibration position..."; currentTelescope->Slew(newRA, newDEC); appendLogText(i18n("Slewing to calibration position, please wait until telescope completes slewing.")); } void Align::getFormattedCoords(double ra, double dec, QString &ra_str, QString &dec_str) { dms ra_s,dec_s; ra_s.setH(ra); dec_s.setD(dec); ra_str = QString("%1:%2:%3").arg(ra_s.hour(), 2, 10, QChar('0')).arg(ra_s.minute(), 2, 10, QChar('0')).arg(ra_s.second(), 2, 10, QChar('0')); if (dec_s.Degrees() < 0) dec_str = QString("-%1:%2:%3").arg(abs(dec_s.degree()), 2, 10, QChar('0')).arg(abs(dec_s.arcmin()), 2, 10, QChar('0')).arg(dec_s.arcsec(), 2, 10, QChar('0')); else dec_str = QString("%1:%2:%3").arg(dec_s.degree(), 2, 10, QChar('0')).arg(dec_s.arcmin(), 2, 10, QChar('0')).arg(dec_s.arcsec(), 2, 10, QChar('0')); } void Align::loadAndSlew(QString fileURL) { /*if (solverTypeGroup->checkedId() == SOLVER_REMOTE) { appendLogText(i18n("Load and Slew is not supported in remote solver mode.")); loadSlewB->setEnabled(false); return; }*/ if (fileURL.isEmpty()) fileURL = QFileDialog::getOpenFileName(KStars::Instance(), i18n("Load Image"), dirPath, "Images (*.fits *.fit *.jpg *.jpeg)"); if (fileURL.isEmpty()) return; QFileInfo fileInfo(fileURL); dirPath = fileInfo.absolutePath(); //loadSlewMode = true; loadSlewState=IPS_BUSY; restartPAHProcess(); slewR->setChecked(true); currentGotoMode = GOTO_SLEW; solveB->setEnabled(false); stopB->setEnabled(true); pi->startAnimation(); startSolving(fileURL, false); } void Align::setExposure(double value) { exposureIN->setValue(value); } void Align::setBinningIndex(int binIndex) { Options::setSolverBinningIndex(binIndex); // If sender is not our combo box, then we need to update the combobox itself if ( dynamic_cast(sender()) != binningCombo) { binningCombo->blockSignals(true); binningCombo->setCurrentIndex(binIndex); binningCombo->blockSignals(false); } // Need to calculate FOV and args for APP if (Options::astrometryImageScaleUnits() == OpsAstrometry::SCALE_ARCSECPERPIX) { calculateFOV(); generateArgs(); } } void Align::setSolverArguments(const QString & value) { solverOptions->setText(value); } void Align::setSolverSearchCoords(double ra, double dec) { dms RA, DEC; RA.setH(ra); DEC.setD(dec); //TODO //raBox->setText(RA.toHMSString()); //decBox->setText(DEC.toDMSString()); } void Align::setUseOAGT(bool enabled) { FOVScopeCombo->setCurrentIndex(enabled ? ISD::CCD::TELESCOPE_GUIDE : ISD::CCD::TELESCOPE_PRIMARY); } FOV* Align::fov() { if (sOrientation == -1) return NULL; else return solverFOV; } void Align::setLockedFilter(ISD::GDInterface *filter, int lockedPosition) { currentFilter = filter; if (currentFilter) { lockedFilterIndex = lockedPosition; INumberVectorProperty *filterSlot = filter->getBaseDevice()->getNumber("FILTER_SLOT"); if (filterSlot) currentFilterIndex = filterSlot->np[0].value-1; connect(currentFilter, SIGNAL(numberUpdated(INumberVectorProperty*)), this, SLOT(processFilterNumber(INumberVectorProperty*)), Qt::UniqueConnection); } } void Align::processFilterNumber(INumberVectorProperty *nvp) { if (currentFilter && !strcmp(nvp->name, "FILTER_SLOT") && !strcmp(nvp->device, currentFilter->getDeviceName())) { currentFilterIndex = nvp->np[0].value - 1; if (filterPositionPending) { if (currentFilterIndex == lockedFilterIndex) { filterPositionPending = false; captureAndSolve(); } } } } void Align::setWCSEnabled(bool enable) { if (currentCCD == NULL) return; ISwitchVectorProperty *wcsControl = currentCCD->getBaseDevice()->getSwitch("WCS_CONTROL"); ISwitch *wcs_enable = IUFindSwitch(wcsControl, "WCS_ENABLE"); ISwitch *wcs_disable = IUFindSwitch(wcsControl, "WCS_DISABLE"); if (!wcs_enable || !wcs_disable) return; if ( (wcs_enable->s == ISS_ON && enable) || (wcs_disable->s == ISS_ON && !enable) ) return; IUResetSwitch(wcsControl); if (enable) { appendLogText(i18n("World Coordinate System (WCS) is enabled. CCD rotation must be set either manually in the CCD driver or by solving an image before proceeding to capture any further images, otherwise the WCS information may be invalid.")); wcs_enable->s = ISS_ON; } else { wcs_disable->s = ISS_ON; m_wcsSynced=false; appendLogText(i18n("World Coordinate System (WCS) is disabled.")); } ClientManager *clientManager = currentCCD->getDriverInfo()->getClientManager(); clientManager->sendNewSwitch(wcsControl); } void Align::checkCCDExposureProgress(ISD::CCDChip *targetChip, double remaining, IPState state) { INDI_UNUSED(targetChip); INDI_UNUSED(remaining); if (state == IPS_ALERT) { if (++retries == 3 && pahStage != PAH_REFRESH) { appendLogText(i18n("Capture error! Aborting...")); abort(); return; } appendLogText(i18n("Restarting capture attempt #%1", retries)); captureAndSolve(); } } void Align::setFocusStatus(Ekos::FocusState state) { focusState = state; } QStringList Align::getSolverOptionsFromFITS(const QString &filename) { int status=0, fits_ccd_width, fits_ccd_height, fits_focal_length=-1, fits_binx=1, fits_biny=1; char comment[128], error_status[512]; fitsfile * fptr = NULL; double ra=0,dec=0, fits_fov_x, fits_fov_y, fov_lower, fov_upper, fits_ccd_hor_pixel=-1, fits_ccd_ver_pixel=-1; QString fov_low,fov_high; QStringList solver_args; QVariantMap optionsMap; if (Options::astrometryUseNoVerify()) optionsMap["noverify"] = true; if (Options::astrometryUseResort()) optionsMap["resort"] = true; if (Options::astrometryUseNoFITS2FITS()) optionsMap["nofits2fits"] = true; if (Options::astrometryUseDownsample()) optionsMap["downsample"] = Options::astrometryDownsample(); if (Options::astrometryCustomOptions().isEmpty() == false) optionsMap["custom"] = Options::astrometryCustomOptions(); solver_args = generateOptions(optionsMap); if (fits_open_image(&fptr, filename.toLatin1(), READONLY, &status)) { fits_report_error(stderr, status); fits_get_errstatus(status, error_status); qWarning() << "Could not open file " << filename << " Error: " << QString::fromUtf8(error_status); return solver_args; } if (fits_read_key(fptr, TINT, "NAXIS1", &fits_ccd_width, comment, &status )) { fits_report_error(stderr, status); fits_get_errstatus(status, error_status); appendLogText(i18n("FITS header: Cannot find NAXIS1.")); return solver_args; } if (fits_read_key(fptr, TINT, "NAXIS2", &fits_ccd_height, comment, &status )) { fits_report_error(stderr, status); fits_get_errstatus(status, error_status); appendLogText(i18n("FITS header: Cannot find NAXIS2.")); return solver_args; } bool coord_ok = true; if (fits_read_key(fptr, TDOUBLE, "OBJCTRA", &ra, comment, &status )) { char objectra_str[32]; if (fits_read_key(fptr, TSTRING, "OBJCTRA", objectra_str, comment, &status )) { fits_report_error(stderr, status); fits_get_errstatus(status, error_status); coord_ok=false; appendLogText(i18n("FITS header: Cannot find OBJCTRA.")); } else { dms raDMS = dms::fromString(objectra_str, false); ra = raDMS.Hours(); } } if (coord_ok && fits_read_key(fptr, TDOUBLE, "OBJCTDEC", &dec, comment, &status )) { char objectde_str[32]; if (fits_read_key(fptr, TSTRING, "OBJCTDEC", objectde_str, comment, &status )) { fits_report_error(stderr, status); fits_get_errstatus(status, error_status); coord_ok=false; appendLogText(i18n("FITS header: Cannot find OBJCTDEC.")); } else { dms deDMS = dms::fromString(objectde_str, true); dec = deDMS.Degrees(); } } /*if (coord_ok == false) { ra = telescopeCoord.ra0().Hours(); dec = telescopeCoord.dec0().Degrees(); }*/ if (coord_ok && Options::astrometryUsePosition()) solver_args << "-3" << QString::number(ra*15.0) << "-4" << QString::number(dec) << "-5 15"; if (fits_read_key(fptr, TINT, "FOCALLEN", &fits_focal_length, comment, &status )) { fits_report_error(stderr, status); fits_get_errstatus(status, error_status); appendLogText(i18n("FITS header: Cannot find FOCALLEN.")); return solver_args; } if (fits_read_key(fptr, TDOUBLE, "PIXSIZE1", &fits_ccd_hor_pixel, comment, &status )) { fits_report_error(stderr, status); fits_get_errstatus(status, error_status); appendLogText(i18n("FITS header: Cannot find PIXSIZE1.")); return solver_args; } if (fits_read_key(fptr, TDOUBLE, "PIXSIZE2", &fits_ccd_ver_pixel, comment, &status )) { fits_report_error(stderr, status); fits_get_errstatus(status, error_status); appendLogText(i18n("FITS header: Cannot find PIXSIZE2.")); return solver_args; } fits_read_key(fptr, TINT, "XBINNING", &fits_binx, comment, &status ); fits_read_key(fptr, TINT, "YBINNING", &fits_biny, comment, &status ); // Calculate FOV fits_fov_x = 206264.8062470963552 * fits_ccd_width * fits_ccd_hor_pixel / 1000.0 / fits_focal_length * fits_binx; fits_fov_y = 206264.8062470963552 * fits_ccd_height * fits_ccd_ver_pixel / 1000.0 / fits_focal_length* fits_biny; fits_fov_x /= 60.0; fits_fov_y /= 60.0; // let's stretch the boundaries by 10% fov_lower = qMin(fits_fov_x, fits_fov_y); fov_upper = qMax(fits_fov_x, fits_fov_y); fov_lower *= 0.90; fov_upper *= 1.10; fov_low = QString::number(fov_lower); fov_high = QString::number(fov_upper); if (Options::astrometryUseImageScale()) solver_args << "-L" << fov_low << "-H" << fov_high << "-u" << "aw"; return solver_args; } void Align::saveSettleTime() { Options::setSettlingTime(delaySpin->value()); } void Align::setCaptureStatus(CaptureState newState) { switch (newState) { case CAPTURE_ALIGNING: QTimer::singleShot(Options::settlingTime(), this, SLOT(captureAndSolve())); break; default: break; } } void Align::showFITSViewer() { FITSData *data = alignView->getImageData(); if (data) { QUrl url = QUrl::fromLocalFile(data->getFilename()); if (fv.isNull()) { if (Options::singleWindowCapturedFITS()) fv = KStars::Instance()->genericFITSViewer(); else { fv = new FITSViewer(Options::independentWindowFITS() ? NULL : KStars::Instance()); KStars::Instance()->getFITSViewersList().append(fv); } fv->addFITS(&url); FITSView *currentView = fv->getCurrentView(); if (currentView) currentView->getImageData()->setAutoRemoveTemporaryFITS(false); } else fv->updateFITS(&url, 0); fv->show(); } } void Align::toggleAlignWidgetFullScreen() { if (alignWidget->parent() == NULL) { alignWidget->setParent(this); rightLayout->insertWidget(0, alignWidget); //rightLayout->setStretch(0, 2); // rightLayout->setStretch(1, 1); alignWidget->showNormal(); } else { alignWidget->setParent(0); alignWidget->setWindowTitle(i18n("Align Frame")); alignWidget->setWindowFlags(Qt::Window | Qt::WindowTitleHint | Qt::CustomizeWindowHint); alignWidget->showMaximized(); alignWidget->show(); } } void Align::startPAHProcess() { pahStage = PAH_FIRST_CAPTURE; nothingR->setChecked(true); currentGotoMode = GOTO_NOTHING; loadSlewB->setEnabled(false); if (Options::limitedResourcesMode()) appendLogText(i18n("Warning: Equatorial Grid Lines will not be drawn due to limited resources mode.")); appendLogText(i18n("Clearing mount Alignment Model...")); currentTelescope->clearAlignmentModel(); PAHWidgets->setCurrentWidget(PAHFirstCapturePage); } void Align::restartPAHProcess() { if (pahStage == PAH_IDLE) return; // Only display dialog if user explicity restarts if ( (static_cast(sender()) == PAHRestartB) && KMessageBox::questionYesNo(KStars::Instance(), i18n("Are you sure you want to restart the polar alignment process?"), i18n("Polar Alignment Assistant"), KStandardGuiItem::yes(), KStandardGuiItem::no(), "restart_PAA_process_dialog") == KMessageBox::No) return; pahStage = PAH_IDLE; PAHFirstCaptureB->setEnabled(true); PAHSecondCaptureB->setEnabled(true); PAHThirdCaptureB->setEnabled(true); PAHFirstRotateB->setEnabled(true); PAHSecondRotateB->setEnabled(true); PAHRefreshB->setEnabled(true); PAHWidgets->setCurrentWidget(PAHIntroPage); qDeleteAll(pahImageInfos); pahImageInfos.clear(); correctionVector = QLineF(); correctionOffset = QPointF(); alignView->setCorrectionParams(correctionVector); alignView->setCorrectionOffset(correctionOffset); alignView->setRACircle(QVector3D()); disconnect(alignView, SIGNAL(trackingStarSelected(int,int)), this, SLOT(setPAHCorrectionOffset(int,int))); } void Align::rotatePAH() { double raDiff = (pahStage == PAH_FIRST_ROTATE) ? PAHFirstRotationSpin->value() : PAHSecondRotationSpin->value(); bool westMeridian = (pahStage == PAH_FIRST_ROTATE) ? PAHFirstWestMeridianR->isChecked() : PAHSecondWestMeridianR->isChecked(); // North if (hemisphere == NORTH_HEMISPHERE) { // West if (westMeridian) raDiff *= -1; // East else raDiff *= 1; } // South else { // West if (westMeridian) raDiff *= 1; // East else raDiff *= -1; } SkyPoint targetPAH; double newTelescopeRA = (telescopeCoord.ra().Degrees() + raDiff) / 15.0; targetPAH.setRA(newTelescopeRA); targetPAH.setDec(telescopeCoord.dec()); // Convert to JNow //targetPAH.apparentCoord((long double) J2000, KStars::Instance()->data()->ut().djd()); // Get horizontal coords //targetPAH.EquatorialToHorizontal(KStarsData::Instance()->lst(), KStarsData::Instance()->geo()->lat()); currentTelescope->Slew(&targetPAH); appendLogText(i18n("Please wait until mount completes rotating to RA (%1) DE (%2)", targetPAH.ra().toHMSString(), targetPAH.dec().toDMSString())); } void Align::calculatePAHError() { QVector3D RACircle; bool rc = findRACircle(RACircle); if (rc == false) { appendLogText(i18n("Failed to find a solution. Try again.")); restartPAHProcess(); return; } if (alignView->isEQGridShown() == false) alignView->toggleEQGrid(); alignView->setRACircle(RACircle); FITSData *imageData = alignView->getImageData(); QPointF RACenterPoint(RACircle.x(), RACircle.y()); SkyPoint RACenter; rc = imageData->pixelToWCS(RACenterPoint, RACenter); if (rc == false) { appendLogText(i18n("Failed to find RA Axis center: %1.", imageData->getLastError())); return; } SkyPoint CP(0, (hemisphere == NORTH_HEMISPHERE) ? 90 : -90); RACenter.setRA(RACenter.ra0()); RACenter.setDec(RACenter.dec0()); dms polarError = RACenter.angularDistanceTo(&CP); if (Options::alignmentLogging()) { qDebug() << "Alignment: RA Axis Circle X: " << RACircle.x() << " Y: " << RACircle.y() << " Radius: " << RACircle.z(); qDebug() << "Alignment: RA Axis Location RA: " << RACenter.ra0().toHMSString() << "DE: " << RACenter.dec0().toDMSString(); qDebug() << "Alignment: RA Axis Offset: " << polarError.toDMSString(); } PAHErrorLabel->setText(polarError.toDMSString()); correctionVector.setP1(celestialPolePoint); correctionVector.setP2(RACenterPoint); /* bool RAAxisInside = imageData->contains(RACenterPoint); bool CPPointInside= imageData->contains(celestialPolePoint); if (RAAxisInside == false && CPPointInside == false) appendLogText(i18n("Warning: Mount axis and celestial pole are outside the field of view. Correction vector may be inaccurate.")); */ alignView->setCorrectionParams(correctionVector); connect(alignView, SIGNAL(trackingStarSelected(int,int)), this, SLOT(setPAHCorrectionOffset(int,int))); } void Align::setPAHCorrectionOffset(int x, int y) { correctionOffset.setX(x); correctionOffset.setY(y); alignView->setCorrectionOffset(correctionOffset); } void Align::setPAHCorrectionSelectionComplete() { pahStage = PAH_PRE_REFRESH; PAHWidgets->setCurrentWidget(PAHRefreshPage); } void Align::startPAHRefreshProcess() { pahStage = PAH_REFRESH; PAHRefreshB->setEnabled(false); // Hide EQ Grids if shown if (alignView->isEQGridShown()) alignView->toggleEQGrid(); // We for refresh, just capture really captureAndSolve(); } void Align::setPAHRefreshComplete() { pahStage = PAH_REFRESH; abort(); restartPAHProcess(); } void Align::processPAHStage(double orientation, double ra, double dec, double pixscale) { // Create temporary file to hold all WCS data QTemporaryFile tmpFile(QDir::tempPath() + "/fitswcsXXXXXX"); tmpFile.setAutoRemove(false); tmpFile.open(); QString newWCSFile = tmpFile.fileName(); tmpFile.close(); if (pahStage == PAH_FIND_CP) { setGOTOMode(GOTO_NOTHING); appendLogText(i18n("Mount is synced to celestial pole. You can now continue Polar Alignment Assistant procedure.")); pahStage = PAH_FIRST_CAPTURE; return; } bool rc = false; if (pahStage == PAH_FIRST_CAPTURE) { // Set First PAH Center PAHImageInfo *solution = new PAHImageInfo(); solution->skyCenter.setRA0(alignCoord.ra0()); solution->skyCenter.setDec0(alignCoord.dec0()); solution->orientation = orientation; solution->pixelScale = pixscale; pahImageInfos.append(solution); // Only invoke this if limited resource mode is false since we want to use CPU heavy WCS if (Options::limitedResourcesMode() == false) { rc = alignView->createWCSFile(newWCSFile, orientation, ra, dec, pixscale); if (rc == false) { appendLogText(i18n("Error creating WCS file: %1", alignView->getImageData()->getLastError())); return; } // Find Celestial pole location SkyPoint CP(0, (hemisphere == NORTH_HEMISPHERE) ? 90 : -90); FITSData *imageData = alignView->getImageData(); QPointF pixelPoint, imagePoint; rc = imageData->wcsToPixel(CP, pixelPoint, imagePoint); // TODO check if pixelPoint is located TOO far from the current position as well // i.e. if X > Width * 2..etc if (rc == false) { appendLogText(i18n("Failed to process World Coordinate System: %1. Try again.", imageData->getLastError())); return; } // If celestial pole out of range, ask the user if they want to move to it if (pixelPoint.x() < (-1*imageData->getWidth()) || pixelPoint.x() > (imageData->getWidth()*2) || pixelPoint.y() < (-1*imageData->getHeight())|| pixelPoint.y() > (imageData->getHeight()*2)) { if (currentTelescope->canSync() && KMessageBox::questionYesNo(0, i18n("Celestial pole is located outside of the field of view. Would you like to sync and slew the telescope to the celestial pole? WARNING: Slewing near poles may cause your mount to end up in unsafe position. Proceed with caution.")) == KMessageBox::Yes) { pahStage = PAH_FIND_CP; targetCoord.setRA(KStarsData::Instance()->lst()->Hours()); targetCoord.setDec(CP.dec().Degrees() > 0 ? 89.5 : -89.5); qDeleteAll(pahImageInfos); pahImageInfos.clear(); setGOTOMode(GOTO_SLEW); Sync(); return; } else appendLogText(i18n("Warning: Celestial pole is located outside the field of view. Move the mount closer to the celestial pole.")); } } pahStage = PAH_FIRST_ROTATE; PAHWidgets->setCurrentWidget(PAHFirstRotatePage); } else if (pahStage == PAH_SECOND_CAPTURE) { rc= alignView->createWCSFile(newWCSFile, orientation, ra, dec, pixscale); if (rc == false) { appendLogText(i18n("Error creating WCS file: %1", alignView->getImageData()->getLastError())); // Not critical error //return; } // Set 2nd PAH Center PAHImageInfo *solution = new PAHImageInfo(); solution->skyCenter.setRA0(alignCoord.ra0()); solution->skyCenter.setDec0(alignCoord.dec0()); solution->orientation = orientation; solution->pixelScale = pixscale; pahImageInfos.append(solution); pahStage = PAH_SECOND_ROTATE; PAHWidgets->setCurrentWidget(PAHSecondRotatePage); if (PAHAutoModeCheck->isChecked()) { PAHSecondRotateB->setEnabled(true); PAHSecondRotateB->click(); } } else if (pahStage == PAH_THIRD_CAPTURE) { rc = alignView->createWCSFile(newWCSFile, orientation, ra, dec, pixscale); if (rc == false) { appendLogText(i18n("Error creating WCS file: %1", alignView->getImageData()->getLastError())); return; } // Set Third PAH Center PAHImageInfo *solution = new PAHImageInfo(); solution->skyCenter.setRA0(alignCoord.ra0()); solution->skyCenter.setDec0(alignCoord.dec0()); solution->orientation = orientation; solution->pixelScale = pixscale; pahImageInfos.append(solution); // Find Celstial pole location SkyPoint CP(0, (hemisphere == NORTH_HEMISPHERE) ? 90 : -90); FITSData *imageData = alignView->getImageData(); QPointF imagePoint; imageData->wcsToPixel(CP, celestialPolePoint, imagePoint); // Now find pixel locations for all recorded center coordinates in the 3rd frame reference imageData->wcsToPixel(pahImageInfos[0]->skyCenter, pahImageInfos[0]->pixelCenter, imagePoint); imageData->wcsToPixel(pahImageInfos[1]->skyCenter, pahImageInfos[1]->pixelCenter, imagePoint); imageData->wcsToPixel(pahImageInfos[2]->skyCenter, pahImageInfos[2]->pixelCenter, imagePoint); if (Options::alignmentLogging()) { qDebug() << "Alignment: P1 RA: " << pahImageInfos[0]->skyCenter.ra0().toHMSString() << "DE: " << pahImageInfos[0]->skyCenter.dec0().toDMSString(); qDebug() << "Alignment: P2 RA: " << pahImageInfos[1]->skyCenter.ra0().toHMSString() << "DE: " << pahImageInfos[1]->skyCenter.dec0().toDMSString(); qDebug() << "Alignment: P3 RA: " << pahImageInfos[2]->skyCenter.ra0().toHMSString() << "DE: " << pahImageInfos[2]->skyCenter.dec0().toDMSString(); qDebug() << "Alignment: P1 X: " << pahImageInfos[0]->pixelCenter.x() << "Y: " << pahImageInfos[0]->pixelCenter.y(); qDebug() << "Alignment: P2 X: " << pahImageInfos[1]->pixelCenter.x() << "Y: " << pahImageInfos[1]->pixelCenter.y(); qDebug() << "Alignment: P3 X: " << pahImageInfos[2]->pixelCenter.x() << "Y: " << pahImageInfos[2]->pixelCenter.y(); } // We have 3 points which uniquely defines a circle with its center representing the RA Axis // We have celestial pole location. So correction vector is just the vector between these two points calculatePAHError(); pahStage = PAH_STAR_SELECT; PAHWidgets->setCurrentWidget(PAHCorrectionPage); } } void Align::syncGuideScopeCCDs() { updateGuideScopeCCDs(FOVScopeCombo->currentIndex()); } void Align::updateGuideScopeCCDs(int index) { if (currentCCD == NULL) return; if (index == ISD::CCD::TELESCOPE_GUIDE && guideScopeCCDs.contains(currentCCD->getDeviceName()) == false) guideScopeCCDs.append(currentCCD->getDeviceName()); else if (index == ISD::CCD::TELESCOPE_PRIMARY) guideScopeCCDs.removeOne(currentCCD->getDeviceName()); if (guideScopeCCDs.contains(currentCCD->getDeviceName())) currentCCD->setTelescopeType(ISD::CCD::TELESCOPE_GUIDE); else currentCCD->setTelescopeType(ISD::CCD::TELESCOPE_PRIMARY); Options::setGuideScopeCCDs(guideScopeCCDs); syncTelescopeInfo(); } // Function adapted from https://rosettacode.org/wiki/Circles_of_given_radius_through_two_points Align::CircleSolution Align::findCircleSolutions(const QPointF & p1, const QPointF p2, double angle, QPair &circleSolutions) { QPointF solutionOne(1,1), solutionTwo(1,1); double radius = distance(p1, p2) / (dms::DegToRad * angle); if (p1 == p2) { if (angle == 0) { circleSolutions = qMakePair(p1, p2); appendLogText(i18n("Only one solution is found.")); return ONE_CIRCLE_SOLUTION; } else { circleSolutions = qMakePair(solutionOne, solutionTwo); appendLogText(i18n("Infinite number of solutions found.")); return INFINITE_CIRCLE_SOLUTION; } } QPointF center( p1.x()/2 + p2.x()/2, p1.y()/2 + p2.y()/2); double halfDistance = distance(center, p1); if (halfDistance > radius) { circleSolutions = qMakePair(solutionOne, solutionTwo); appendLogText(i18n("No solution is found. Points are too far away")); return NO_CIRCLE_SOLUTION; } if (halfDistance - radius == 0) { circleSolutions = qMakePair(center, solutionTwo); appendLogText(i18n("Only one solution is found.")); return ONE_CIRCLE_SOLUTION; } double root = std::hypotf(radius, halfDistance) / distance (p1, p2); solutionOne.setX(center.x() + root * (p1.y() - p2.y())); solutionOne.setY(center.y() + root * (p2.x() - p1.x())); solutionTwo.setX(center.x() - root * (p1.y() - p2.y())); solutionTwo.setY(center.y() - root * (p2.x() - p1.x())); circleSolutions = qMakePair(solutionOne, solutionTwo); return TWO_CIRCLE_SOLUTION; } double Align::distance(const QPointF & p1, const QPointF & p2) { return std::hypotf(p2.x()-p1.x(), p2.y()-p1.y()); } bool Align::findRACircle(QVector3D & RACircle) { bool rc = false; QPointF p1 = pahImageInfos[0]->pixelCenter; QPointF p2 = pahImageInfos[1]->pixelCenter; QPointF p3 = pahImageInfos[2]->pixelCenter; if (!isPerpendicular(p1, p2, p3)) rc = calcCircle(p1, p2, p3, RACircle); else if (!isPerpendicular(p1, p3, p2) ) rc = calcCircle(p1, p3, p2, RACircle); else if (!isPerpendicular(p2, p1, p3) ) rc = calcCircle(p2, p1, p3, RACircle); else if (!isPerpendicular(p2, p3, p1) ) rc = calcCircle(p2, p3, p1, RACircle); else if (!isPerpendicular(p3, p2, p1) ) rc = calcCircle(p3, p2, p1, RACircle); else if (!isPerpendicular(p3, p1, p2) ) rc = calcCircle(p3, p1, p2, RACircle); else { //TRACE("\nThe three pts are perpendicular to axis\n"); return false; } return rc; } bool Align::isPerpendicular(const QPointF &p1, const QPointF &p2, const QPointF &p3) // Check the given point are perpendicular to x or y axis { double yDelta_a= p2.y() - p1.y(); double xDelta_a= p2.x() - p1.x(); double yDelta_b= p3.y() - p2.y(); double xDelta_b= p3.x() - p2.x(); // checking whether the line of the two pts are vertical if (fabs(xDelta_a) <= 0.000000001 && fabs(yDelta_b) <= 0.000000001) { //TRACE("The points are pependicular and parallel to x-y axis\n"); return false; } if (fabs(yDelta_a) <= 0.0000001) { //TRACE(" A line of two point are perpendicular to x-axis 1\n"); return true; } else if (fabs(yDelta_b) <= 0.0000001) { //TRACE(" A line of two point are perpendicular to x-axis 2\n"); return true; } else if (fabs(xDelta_a)<= 0.000000001) { //TRACE(" A line of two point are perpendicular to y-axis 1\n"); return true; } else if (fabs(xDelta_b)<= 0.000000001) { //TRACE(" A line of two point are perpendicular to y-axis 2\n"); return true; } else return false; } bool Align::calcCircle(const QPointF &p1, const QPointF &p2, const QPointF &p3, QVector3D & RACircle) { double yDelta_a= p2.y() - p1.y(); double xDelta_a= p2.x() - p1.x(); double yDelta_b= p3.y() - p2.y(); double xDelta_b= p3.x() - p2.x(); if (fabs(xDelta_a) <= 0.000000001 && fabs(yDelta_b) <= 0.000000001) { RACircle.setX(0.5*(p2.x() + p3.x())); RACircle.setY(0.5*(p1.y() + p2.y())); QPointF center(RACircle.x(), RACircle.y()); RACircle.setZ(distance(center,p1)); return true; } // IsPerpendicular() assure that xDelta(s) are not zero double aSlope=yDelta_a/xDelta_a; // double bSlope=yDelta_b/xDelta_b; if (fabs(aSlope-bSlope) <= 0.000000001) { // checking whether the given points are colinear. //TRACE("The three ps are colinear\n"); return false; } // calc center RACircle.setX((aSlope*bSlope*(p1.y() - p3.y()) + bSlope*(p1.x() + p2.x()) - aSlope*(p2.x()+p3.x()) )/(2* (bSlope-aSlope) ) ); RACircle.setY(-1*(RACircle.x() - (p1.x()+p2.x())/2)/aSlope + (p1.y()+p2.y())/2 ); QPointF center(RACircle.x(), RACircle.y()); RACircle.setZ(distance(center,p1)); return true; } void Align::setMountStatus(ISD::Telescope::TelescopeStatus newState) { switch (newState) { case ISD::Telescope::MOUNT_PARKING: case ISD::Telescope::MOUNT_SLEWING: solveB->setEnabled(false); loadSlewB->setEnabled(false); PAHFirstCaptureB->setEnabled(false); PAHSecondCaptureB->setEnabled(false); PAHThirdCaptureB->setEnabled(false); break; default: if (state != ALIGN_PROGRESS) { solveB->setEnabled(true); PAHFirstCaptureB->setEnabled(true); PAHSecondCaptureB->setEnabled(true); PAHThirdCaptureB->setEnabled(true); if (pahStage == PAH_IDLE) loadSlewB->setEnabled(true); } break; } } void Align::setAstrometryDevice(ISD::GDInterface *newAstrometry) { remoteParserDevice = newAstrometry; remoteSolverR->setEnabled(true); if (remoteParser) remoteParser->setAstrometryDevice(remoteParserDevice); } } diff --git a/kstars/ekos/align/align.h b/kstars/ekos/align/align.h index b21775846..a11fd886e 100644 --- a/kstars/ekos/align/align.h +++ b/kstars/ekos/align/align.h @@ -1,640 +1,643 @@ /* Ekos Polar Alignment Tool Copyright (C) 2013 Jasem Mutlaq This application is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. */ #ifndef ALIGN_H #define ALIGN_H #include #include #include #include #include "ui_align.h" #include "alignview.h" #include "ekos/ekos.h" #include "indi/inditelescope.h" #include "indi/indiccd.h" #include "indi/indistd.h" #include "starobject.h" #include "ui_mountmodel.h" class FOV; class QProgressIndicator; namespace Ekos { class AstrometryParser; class OnlineAstrometryParser; class OfflineAstrometryParser; class RemoteAstrometryParser; class OpsAstrometry; class OpsAlign; class OpsAstrometryCfg; class OpsAstrometryIndexFiles; /** *@class Align *@short Align class handles plate-solving and polar alignment measurement and correction using astrometry.net * The align class can capture images from the CCD and use either online or offline astrometry.net engine to solve the plate constants and find the center RA/DEC coordinates. The user selects the action * to perform when the solver completes successfully. * Align module provide Polar Align Helper tool which enables easy-to-follow polar alignment procedure given wide FOVs (> 1.5 degrees) * For small FOVs, the Legacy polar alignment measurement should be used. * LEGACY: Measurement of polar alignment errors is performed by capturing two images on selected points in the sky and measuring the declination drift to calculate * the error in the mount's azimutn and altitude displacement from optimal. Correction is carried by asking the user to re-center a star by adjusting the telescope's azimuth and/or altitude knobs. *@author Jasem Mutlaq *@version 1.3 */ class Align : public QWidget, public Ui::Align { Q_OBJECT Q_CLASSINFO("D-Bus Interface", "org.kde.kstars.Ekos.Align") public: Align(); ~Align(); typedef enum { AZ_INIT, AZ_FIRST_TARGET, AZ_SYNCING, AZ_SLEWING, AZ_SECOND_TARGET, AZ_CORRECTING, AZ_FINISHED } AZStage; typedef enum { ALT_INIT, ALT_FIRST_TARGET, ALT_SYNCING, ALT_SLEWING, ALT_SECOND_TARGET, ALT_CORRECTING, ALT_FINISHED } ALTStage; typedef enum { GOTO_SYNC, GOTO_SLEW, GOTO_NOTHING } GotoMode; typedef enum { SOLVER_ONLINE, SOLVER_OFFLINE, SOLVER_REMOTE} SolverType; typedef enum { PAH_IDLE, PAH_FIRST_CAPTURE, PAH_FIND_CP, PAH_FIRST_ROTATE, PAH_SECOND_CAPTURE, PAH_SECOND_ROTATE, PAH_THIRD_CAPTURE, PAH_STAR_SELECT, PAH_PRE_REFRESH, PAH_REFRESH, PAH_ERROR } PAHStage; typedef enum { NORTH_HEMISPHERE, SOUTH_HEMISPHERE } HemisphereType; // Image Scales const QStringList ImageScales = { "dw", "aw", "app" }; enum CircleSolution { NO_CIRCLE_SOLUTION, ONE_CIRCLE_SOLUTION, TWO_CIRCLE_SOLUTION, INFINITE_CIRCLE_SOLUTION }; /** @defgroup AlignDBusInterface Ekos DBus Interface - Align Module * Ekos::Align interface provides advanced scripting capabilities to solve images using online or offline astrometry.net */ /*@{*/ /** DBUS interface function. * Select CCD * @param device CCD device name * @return Returns true if device if found and selected, false otherwise. */ Q_SCRIPTABLE bool setCCD(QString device); /** DBUS interface function. * Start the plate-solving process given the passed image file. * @param filename Name of image file to solve. FITS and JPG/JPG/TIFF formats are accepted. * @param isGenerated Set to true if filename is generated from a CCD capture operation. If the file is loaded from any storage or network media, pass false. * @return Returns true if device if found and selected, false otherwise. */ Q_SCRIPTABLE Q_NOREPLY void startSolving(const QString &filename, bool isGenerated=true); /** DBUS interface function. * Select Goto Mode of Solver. The solver mode is the action the solver performs upon successful solving. * @param mode 0 for Sync, 1 for SlewToTarget, 2 for Nothing */ Q_SCRIPTABLE Q_NOREPLY void setGOTOMode(int mode); /** DBUS interface function. * Returns the solver's solution results * @return Returns array of doubles. First item is RA in degrees. Second item is DEC in degrees. */ Q_SCRIPTABLE QList getSolutionResult(); /** DBUS interface function. * Returns the solver's current status * @return Returns solver status (Ekos::AlignState) */ Q_SCRIPTABLE int getStatus() { return state; } #if 0 /** DBUS interface function. * @return Returns true if the solver process completed or aborted, false otherwise. */ Q_SCRIPTABLE bool isSolverComplete() { return m_isSolverComplete; } /** DBUS interface function. * @return Returns true if the solver process completed successfully, false otherwise. */ Q_SCRIPTABLE bool isSolverSuccessful() { return m_isSolverSuccessful; } #endif /** DBUS interface function. * @return Returns State of load slew procedure. Idle if not started. Busy if in progress. Ok if complete. Alert if procedure failed. */ Q_SCRIPTABLE int getLoadAndSlewStatus() { return loadSlewState; } /** DBUS interface function. * Sets the exposure of the selected CCD device. * @param value Exposure value in seconds */ Q_SCRIPTABLE Q_NOREPLY void setExposure(double value); /** DBUS interface function. * Sets the arguments that gets passed to the astrometry.net offline solver. * @param value space-separated arguments. */ Q_SCRIPTABLE Q_NOREPLY void setSolverArguments(const QString & value); /** DBUS interface function. * Sets the solver search area options * @param ra center RA for search pattern, in hours. * @param dec center DEC for search pattern, in degrees. * @param radius radius of search pattern, in degrees. */ Q_SCRIPTABLE Q_NOREPLY void setSolverSearchCoords(double ra, double dec); /** DBUS interface function. * Sets the solver's option * @param useOAGT if true, use the Off-Axis Guide Telescope focal length and aperture for FOV calculations. Otherwise, the main telescope's focal length and aperture are used for the calculation. */ Q_SCRIPTABLE Q_NOREPLY void setUseOAGT(bool enabled); /** @}*/ /** * @brief Add CCD to the list of available CCD. * @param newCCD pointer to CCD device. */ void addCCD(ISD::GDInterface *newCCD); /** * @brief Set the current telescope * @newTelescope pointer to telescope device. */ void setTelescope(ISD::GDInterface *newTelescope); /** * @brief setAstrometryDevice * @param newAstrometry */ void setAstrometryDevice(ISD::GDInterface *newAstrometry); /** * @brief CCD information is updated, sync them. */ void syncCCDInfo(); /** * @brief Generate arguments we pass to the online and offline solvers. Keep user own arguments in place. */ void generateArgs(); /** * @brief Does our parser exist in the system? */ bool isParserOK(); /** * @brief syncGuideScopeCCDs When TELESCOPE_TYPE is first defined by the CCD driver, update its value given the current state of useGuideScope checkbox * If checked, we set TELESCOPE_TYPE to GUIDE, otherwise it is set to PRIMARY */ void syncGuideScopeCCDs(); // Log QString getLogText() { return logText.join("\n"); } void clearLog(); /** * @brief Return FOV object used to represent the solved image orientation on the sky map. */ FOV *fov(); /** * @brief getFOVScale Returns calculated FOV values * @param fov_w FOV width in arcmins * @param fov_h FOV height in arcmins * @param fov_scale FOV scale in arcsec per pixel */ void getFOVScale(double &fov_w, double & fov_h, double &fov_scale); /** * @brief generateOptions Generate astrometry.net option given the supplied map * @param optionsMap List of key=value pairs for all astrometry.net options * @return String List of valid astrometry.net options */ static QStringList generateOptions(const QVariantMap & optionsMap); static void generateFOVBounds(double fov_h, double fov_v, QString &fov_low, QString &fov_high); public slots: /** * @brief Process updated telescope coordinates. * @coord pointer to telescope coordinate property. */ void processTelescopeNumber(INumberVectorProperty *coord); /** * @brief Check CCD and make sure information is updated and FOV is re-calculated. * @param CCDNum By default, we check the already selected CCD in the dropdown menu. If CCDNum is specified, the check is made against this specific CCD in the dropdown menu. CCDNum is the index of the CCD in the dropdown menu. */ void checkCCD(int CCDNum=-1); /** * @brief checkCCDExposureProgress Track the progress of CCD exposure * @param targeChip Target chip under exposure * @param remaining how many seconds remaining * @param state status of exposure */ void checkCCDExposureProgress(ISD::CCDChip* targetChip, double remaining, IPState state); /** * @brief Process new FITS received from CCD. * @param bp pointer to blob property */ void newFITS(IBLOB *bp); /** \addtogroup AlignDBusInterface * @{ */ /** DBUS interface function. * Aborts the solving operation. */ Q_SCRIPTABLE Q_NOREPLY void abort(); /** DBUS interface function. * Select the solver type * @param type Set solver type. 0 online, 1 offline, 2 remote */ Q_SCRIPTABLE Q_NOREPLY void setSolverType(int type); /** DBUS interface function. * Capture and solve an image using the astrometry.net engine * @return Returns true if the procedure started successful, false otherwise. */ Q_SCRIPTABLE bool captureAndSolve(); /** DBUS interface function. * Loads an image (FITS or JPG/TIFF) and solve its coordinates, then it slews to the solved coordinates and an image is captured and solved to ensure * the telescope is pointing to the same coordinates of the image. * @param fileURL URL to the image to solve */ Q_SCRIPTABLE Q_NOREPLY void loadAndSlew(QString fileURL = QString()); /** DBUS interface function. * Sets the binning of the selected CCD device. * @param binIndex Index of binning value. Default values range from 0 (binning 1x1) to 3 (binning 4x4) */ Q_SCRIPTABLE Q_NOREPLY void setBinningIndex(int binIndex); /** @}*/ /** * @brief Solver finished successfully, process the data and execute the required actions depending on the mode. * @param orientation Orientation of image in degrees (East of North) * @param ra Center RA in solved image, degrees. * @param dec Center DEC in solved image, degrees. * @param pixscale Image scale is arcsec/pixel */ void solverFinished(double orientation, double ra, double dec, double pixscale); /** * @brief Process solver failure. */ void solverFailed(); /** * @brief We received new telescope info, process them and update FOV. */ void syncTelescopeInfo(); /** * @brief setWCSEnabled enables/disables World Coordinate System settings in the CCD driver. * @param enable true to enable WCS, false to disable. */ void setWCSEnabled(bool enable); void setLockedFilter(ISD::GDInterface *filter, int lockedPosition); void setFocusStatus(Ekos::FocusState state); // Log void appendLogText(const QString &); // Capture void setCaptureComplete(); // Update Capture Module status void setCaptureStatus(Ekos::CaptureState newState); // Update Mount module status void setMountStatus(ISD::Telescope::TelescopeStatus newState); private slots: /* Polar Alignment */ void measureAltError(); void measureAzError(); void correctAzError(); void correctAltError(); void processFilterNumber(INumberVectorProperty *nvp); void setDefaultCCD(QString ccd); void saveSettleTime(); // Solver timeout void checkAlignmentTimeout(); void updateGuideScopeCCDs(int index); // External View void showFITSViewer(); void toggleAlignWidgetFullScreen(); // Polar Alignment Helper slots void startPAHProcess(); void restartPAHProcess(); void rotatePAH(); void setPAHCorrectionOffset(int x, int y); void setPAHCorrectionSelectionComplete(); void startPAHRefreshProcess(); void setPAHRefreshComplete(); //Solutions Display slots void buildTarget(); void handlePointTooltip(QMouseEvent * event); void handleVerticalPlotSizeChange(); void handleHorizontalPlotSizeChange(); void selectSolutionTableRow(int row, int column); void slotClearAllSolutionPoints(); void slotRemoveSolutionPoint(); void slotMountModel(); //Mount Model Slots void slotWizardAlignmentPoints(); void slotStarSelected(const QString selectedStar); void slotLoadAlignmentPoints(); void slotSaveAsAlignmentPoints(); void slotSaveAlignmentPoints(); void slotClearAllAlignPoints(); void slotRemoveAlignPoint(); void slotAddAlignPoint(); void slotFindAlignObject(); void resetAlignmentProcedure(); void startStopAlignmentProcedure(); void startAlignmentPoint(); void finishAlignmentPoint(bool solverSucceeded); void moveAlignPoint(int logicalIndex, int oldVisualIndex, int newVisualIndex); void exportSolutionPoints(); void alignTypeChanged(const QString alignType); + void togglePreviewAlignPoints(); signals: void newLog(); void solverComplete(bool); void solverSlewComplete(); void newStatus(Ekos::AlignState state); void newSolutionDeviation(double ra_arcsecs, double de_arcsecs); private: /** * @brief Calculate Field of View of CCD+Telescope combination that we need to pass to astrometry.net solver. */ void calculateFOV(); /** * @brief After a solver process is completed successfully, sync, slew to target, or do nothing as set by the user. */ void executeGOTO(); /** * @brief After a solver process is completed successfully, measure Azimuth or Altitude error as requested by the user. */ void executePolarAlign(); /** * @brief Sync the telescope to the solved alignment coordinate. */ void Sync(); /** * @brief Slew the telescope to the solved alignment coordinate. */ void Slew(); /** * @brief Sync the telescope to the solved alignment coordinate, and then slew to the target coordinate. */ void SlewToTarget(); /** * @brief Calculate polar alignment error magnitude and direction. * The calculation is performed by first capturing and solving a frame, then slewing 30 arcminutes and solving another frame to find the exact coordinates, then computing the error. * @param initRA RA of first frame. * @param initDEC DEC of first frame * @param finalRA RA of second frame * @param finalDEC DEC of second frame * @param initAz Azimuth of first frame */ void calculatePolarError(double initRA, double initDEC, double finalRA, double finalDEC, double initAz); /** * @brief Get formatted RA & DEC coordinates compatible with astrometry.net format. * @param ra Right ascension * @param dec Declination * @param ra_str will contain the formatted RA string * @param dec_str will contain the formatted DEC string */ void getFormattedCoords(double ra, double dec, QString &ra_str, QString &dec_str); /** * @brief getSolverOptionsFromFITS Generates a set of solver options given the supplied FITS image. The function reads FITS keyword headers and build the argument list accordingly. In case of a missing header keyword, it falls back to * the Alignment module existing values. * @param filename FITS path * @return List of Solver options */ QStringList getSolverOptionsFromFITS(const QString &filename); /** * @brief calculatePAHError Calculate polar alignment error in the Polar Alignment Helper (PAH) method */ void calculatePAHError(); /** * @brief processPAHStage After solver is complete, handle PAH Stage processing */ void processPAHStage(double orientation, double ra, double dec, double pixscale); CircleSolution findCircleSolutions(const QPointF & p1, const QPointF p2, double angle, QPair & circleSolutions); double distance(const QPointF & p1, const QPointF & p2); bool findRACircle(QVector3D & RACircle); bool isPerpendicular(const QPointF &p1, const QPointF &p2, const QPointF &p3); bool calcCircle(const QPointF &p1, const QPointF &p2, const QPointF &p3, QVector3D &RACircle); // Which chip should we invoke in the current CCD? bool useGuideHead; // Can the mount sync its coordinates to those set by Ekos? bool canSync; // LoadSlew mode is when we load an image and solve it, no capture is done. //bool loadSlewMode; // If load and slew is solved successfully, coordinates obtained, slewed to target, and then captured, solved, and re-slewed to target again. IPState loadSlewState; // Solver iterations count uint8_t solverIterations; // FOV double ccd_hor_pixel, ccd_ver_pixel, focal_length, aperture, fov_x, fov_y, fov_pixscale; int ccd_width, ccd_height; // Keep track of solver results double sOrientation, sRA, sDEC; // Solver alignment coordinates SkyPoint alignCoord; // Target coordinates we need to slew to SkyPoint targetCoord; // Actual current telescope coordinates SkyPoint telescopeCoord; // Coord from Load & Slew SkyPoint loadSlewCoord; // Difference between solution and target coordinate double targetDiff; // Progress icon if the solver is running QProgressIndicator *pi; // Keep track of how long the solver is running QTime solverTimer; // Polar Alignment AZStage azStage; ALTStage altStage; double azDeviation, altDeviation; double decDeviation; static const double RAMotion; static const float SIDRATE; // Have we slewed? bool isSlewDirty=false; // Online and Offline parsers AstrometryParser *parser; OnlineAstrometryParser *onlineParser; OfflineAstrometryParser *offlineParser; RemoteAstrometryParser *remoteParser; ISD::GDInterface *remoteParserDevice = NULL; // Pointers to our devices ISD::Telescope *currentTelescope; ISD::CCD *currentCCD; QList CCDs; // Optional device filter ISD::GDInterface *currentFilter; int lockedFilterIndex; int currentFilterIndex; // True if we need to change filter position and wait for result before continuing capture bool filterPositionPending; // Keep track of solver FOV to be plotted in the skymap after each successful solve operation FOV *solverFOV; // WCS bool m_wcsSynced; // Log QStringList logText; // Capture retries int retries; // State AlignState state; FocusState focusState; // Track which upload mode the CCD is set to. If set to UPLOAD_LOCAL, then we need to switch it to UPLOAD_CLIENT in order to do focusing, and then switch it back to UPLOAD_LOCAL ISD::CCD::UploadMode rememberUploadMode; GotoMode currentGotoMode; QString dirPath; // Timer QTimer alignTimer; // BLOB Type ISD::CCD::BlobType blobType; QString blobFileName; // Align Frame AlignView *alignView; // FITS Viewer in case user want to display in it instead of internal view QPointer fv; // Polar Alignment Helper PAHStage pahStage; // Sky centers typedef struct { SkyPoint skyCenter; QPointF pixelCenter; double pixelScale; double orientation; } PAHImageInfo; QVector pahImageInfos; // User desired offset when selecting a bright star in the image QPointF celestialPolePoint, correctionOffset; // Correction vector line between mount RA Axis and celestial pole QLineF correctionVector; // CCDs using Guide Scope for parameters QStringList guideScopeCCDs; // Which hemisphere are we located on? HemisphereType hemisphere; // Astrometry Options OpsAstrometry *opsAstrometry; OpsAlign *opsAlign; OpsAstrometryCfg *opsAstrometryCfg; OpsAstrometryIndexFiles *opsAstrometryIndexFiles; void resizeEvent(QResizeEvent * event); QCPCurve *centralTarget = NULL; QCPCurve *yellowTarget = NULL; QCPCurve *redTarget = NULL; QCPCurve *concentricRings = NULL; QDialog mountModelDialog; Ui_mountModel mountModel; int currentAlignmentPoint=0; bool mountModelRunning=false; bool targetAccuracyNotMet=false; bool alignmentPointsAreBad(); bool loadAlignmentPoints(const QString &fileURL); bool saveAlignmentPoints(const QString &path); QUrl alignURL; QUrl alignURLPath; QVector alignStars; void generateAlignStarList(); bool isVisible(const SkyObject *so); double getAltitude(const SkyObject *so); const SkyObject* getWizardAlignObject(double ra, double de, double angle); void calculateAngleForRALine(double &raIncrement,double &initRA, double initDEC, double lat, double raPoints, double minAlt); void calculateAZPointsForDEC(dms dec, dms alt, dms &AZEast, dms &AZWest); + bool previewShowing=false; + void updatePreviewAlignPoints(); }; } #endif // ALIGN_H diff --git a/kstars/ekos/align/mountmodel.ui b/kstars/ekos/align/mountmodel.ui index 65c0681ef..3ee2734b1 100644 --- a/kstars/ekos/align/mountmodel.ui +++ b/kstars/ekos/align/mountmodel.ui @@ -1,587 +1,609 @@ mountModel 0 0 - 360 + 438 535 Dialog 3 3 3 3 3 <html><head/><body><p>This tool can help your mount build a better pointing model by moving the mount to different points in the sky and solving captured images. This can increase the accuracy of your GOTOs. To qualify, your mount must be capable of improving its internal pointing model after each <b>Sync</b> commanded by Ekos. If this is not the case, then this tool is not suitable for your mount.</p></body></html> Qt::AlignJustify|Qt::AlignVCenter true Mount Model Wizard Alignment Points: - <html><head/><body><p>The Number of points the wizard should generate. If this number is less than 5, it will generate them at the same DEC. After that, every 5 points will be at a new dec. Note, if Fixed DEC is selected, all of the points will be at the same DEC.</p></body></html> + <html><head/><body><p>The Number of points the wizard should generate. If this number is less than 5, it will generate them at the same DEC. After that they will be evenly distributed over RA and DEC. Note, if Fixed DEC is selected, all of the points will be at the same DEC.</p></body></html> 2 Qt::Horizontal 171 20 0 0 <html><head/><body><p>Automatically generate the specified number of alignment points in the table below. They will be generated based on the options selected in this wizard.</p></body></html> Generate - Minimum Alt; + Minimum Alt: <html><head/><body><p>This is the minimum altitude above the horizon to use in generating points with the wizard.</p></body></html> 80 30 Decination: false <html><head/><body><p>The Declination of the points that will be generated. This option only applies to the Fixed DEC Option.</p></body></html> -89 89 Qt::Horizontal 77 20 Object Type: <html><head/><body><p>This allows you to select the type of objects/points added by the wizard:</p><p>Note: all of the options except Fixed DEC start with a grid of RA/DEC points.</p><p>Any Stars--The wizard searches for the nearest star.</p><p>Any object --The wizard searches for the nearest object of any type.</p><p>Named Stars--The wizard searches for the nearest star in the currently visible star list. Note that the first named star might be fairly far from the intended point and also sometimes the same star could be the closest one for multiple points.</p><p>Fixed DEC -- The wizard generates all points at the chosen DEC.</p><p>Fixed Grid -- The wizard jut uses the original grid without trying to pair it with objects.</p></body></html> Any Stars Named Stars Any Object Fixed DEC Fixed Grid Qt::Horizontal 77 20 Add Currently Visible Stars 3 3 3 3 Common Names: 0 0 <html><head/><body><p>These are the common names for the named stars that are currently up at your location. If you select a star, it will be added to the table below.</p></body></html> Qt::Horizontal 125 20 Greek: 0 0 <html><head/><body><p>These are the Greek names for the brighter stars that are currently up at your location. If you select a star, it will be added to the table below.</p></body></html> Mount Model Alignment Points 3 3 3 3 0 0 32 32 32 32 <html><head/><body><p>Clear the entire table of alignment points below. Note that this will not affect your mount in any way. It will not clear your pointing model. It just clears the list of points so that if you run the mount model routine again, these points will not be included.</p></body></html> 32 32 32 32 <html><head/><body><p>Remove the selected row from the alignment point table below. Note that it does <b>not</b> command your mount to clear that point from the model, it just means that if you run the mount model routine, that point will not be included in the list.</p></body></html> 32 32 32 32 <html><head/><body><p>Add an empty row to the alignment list so that you can manually add an RA and DE for a position you would like the telescope to slew to when doing the mount model routine.</p></body></html> 32 32 32 32 <html><head/><body><p>Open the Find Object dialog so that you can select an object to get an RA and DEC automatically that you would like the telescope to slew to when doing the mount model routine.</p></body></html> Qt::Horizontal 40 20 <html><head/><body><p>This is a table of alignment points that you would like the mount model routine to use in aligning your telescope. If you would like to add or remove points, you may use the buttons above. If you would like to run, pause, or stop the mount model routine, you may use the buttons below. You can directly edit an RA or DE for a particular point directly in the cells. Just be sure to follow the format hh:mm:ss. The object name is just provided for convenience and is not used. If you want to change the order of the alignment points, you can drag the row up or down using the vertical header for that row on the left. </p></body></html> RA AlignCenter DEC AlignCenter Object AlignCenter ~~ AlignCenter 0 32 32 32 32 <html><head/><body><p>Load a saved alignment point list so that you can run the same alignment you ran on a previous occasion.</p></body></html> 32 32 32 32 <html><head/><body><p>Save the alignment list so that you can run this alignment again in the future.</p></body></html> 32 32 32 32 <html><head/><body><p>Save the alignment list so that you can run this alignment again in the future.</p></body></html> + + + + + 32 + 32 + + + + + 32 + 32 + + + + <html><head/><body><p>Display the alignment points on the skymap as flags. Clicking it again will remove the flags. The flags will not be saved.</p></body></html> + + + + + + Qt::Horizontal 40 20 32 32 32 32 <html><head/><body><p>Start or pause the mount model routine. It will slew to and astrometrically solve the list of points in the table above using the settings in the align module. If the routine was previously paused, it will pick up where it left off. If it was stopped or it had finished it will start the routine over again.</p></body></html> 32 32 32 32 <html><head/><body><p>Stop the mount model routine. It will <b>not</b> clear any points from your telescope's pointing model. It will stop the routine and any points currently being solved. If you run the mount model again after hitting stop, it will start the routine over again with the first point.</p></body></html>