1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
|
#include <QApplication>
//Includes for this project
#include "structuredefinitions.h"
#include "stellarsolver.h"
#include "ssolverutils/fileio.h"
int main(int argc, char *argv[])
{
QApplication app(argc, argv);
#if defined(__linux__)
setlocale(LC_NUMERIC, "C");
#endif
fileio imageLoader;
imageLoader.logToSignal = false;
if(!imageLoader.loadImage("randomsky.fits"))
{
printf("Error in loading FITS file");
exit(1);
}
FITSImage::Statistic stats = imageLoader.getStats();
uint8_t *imageBuffer = imageLoader.getImageBuffer();
StellarSolver stellarSolver(SSolver::SOLVE, stats, imageBuffer);
stellarSolver.setIndexFolderPaths(QStringList() << "astrometry");
if(imageLoader.position_given)
{
printf("Using Position: %f hours, %f degrees\n", imageLoader.ra, imageLoader.dec);
stellarSolver.setSearchPositionRaDec(imageLoader.ra, imageLoader.dec);
}
if(imageLoader.scale_given)
{
stellarSolver.setSearchScale(imageLoader.scale_low, imageLoader.scale_high, imageLoader.scale_units);
printf("Using Scale: %f to %f, %s\n", imageLoader.scale_low, imageLoader.scale_high, SSolver::getScaleUnitString(imageLoader.scale_units).toUtf8().data());
}
printf("Starting to solve. . .\n");
fflush( stdout );
if(!stellarSolver.solve())
{
printf("Solver Failed");
exit(0);
}
FITSImage::Solution solution = stellarSolver.getSolution();
printf("+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++\n");
printf("Field center: (RA,Dec) = (%f, %f) deg.\n", solution.ra, solution.dec);
printf("Field size: %f x %f arcminutes\n", solution.fieldWidth, solution.fieldHeight);
printf("Pixel Scale: %f\"\n", solution.pixscale);
printf("Field rotation angle: up is %f degrees E of N\n", solution.orientation);
printf("Field parity: %s\n\n", FITSImage::getParityText(solution.parity).toUtf8().data());
return 0;
}
|
