static void
smf__fits_copy_items( AstFitsChan * fromfits, AstFitsChan * tofits,
const char ** items, int * status ) {
size_t i = 0;
double value;
if (*status != SAI__OK) return;
while ( items[i] != NULL ) {
char card[ 81 ];
/* reset the position each time since we can not be sure that
"items" is in order */
astClear( fromfits, "Card" );
/* Look in fromfits for the card */
if ( astFindFits( fromfits, items[i], card, 0 ) ) {
/* now look in tofits for the card */
astClear( tofits, "Card" );
if ( astFindFits( tofits, items[i], NULL, 0 ) ) {
/* and if we find it replace it */
if (astGetI(fromfits, "CardType") == AST__FLOAT
&& astGetFitsF(fromfits, NULL, &value)) {
astSetFitsF(tofits, items[i], value, NULL, 1);
}
else {
astPutFits( tofits, card, 1 );
}
}
}
i++;
}
}
Hero * Hero::constructFromFitsFile(const QString &fname)
{
FitsParser parser;
bool parsedOk = parser.loadFile( FitsFileLocation::fromLocal( fname));
if( ! parsedOk) {
dbg(1) << "Parser failed to load " << fname;
Hero * heroPtr = new Hero;
heroPtr-> addError( "FitsParser failed to load the file");
return heroPtr;
}
// alias hdr
auto & hdr = parser.getHeaderInfo().headerLines;
Hero * heroPtr = new Hero;
Hero & hero = * heroPtr;
AstErrorGuard guard( heroPtr);
AstGCGuard gcGuard;
// set naxes in case AST fails to read this file
hero.m_ast.naxes = parser.getHeaderInfo().naxis;
// set up bunit
{
hero.m_bunit = parser.getHeaderInfo().bunit;
}
// and the nicer version of bunit
{
QString u = hero.m_bunit.simplified();
if( u.toLower() == "kelvin") {
hero.m_bunitNiceHtml = "K";
}
else {
hero.m_bunitNiceHtml = u;
}
}
// Create a FitsChan and feed it the fits header
AstFitsChan *fitschan;
#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wformat-zero-length"
fitschan = astFitsChan( NULL, NULL, "" );
#pragma GCC diagnostic pop
std::cout << "astOK = " << astOK << "\n";
// feed the header lines one by one and check for errors
for( const QString & s : hdr) {
std::string stdstr = s.toStdString();
astPutFits( fitschan, stdstr.c_str(), 1);
if( ! astOK) {
astClearStatus;
QString ss = s.trimmed();
std::cout << "Skipping bad card: " << ss << "\n";
hero.addError( "Skipping card: " + ss);
}
}
// reposition to the beginning of the channel (ast thing, it's required, hmmmkey)
astClear( fitschan, "Card" );
std::cout << "astOK = " << astOK << "\n";
std::cout << "Here\n";
auto encoding = AstWrappers::getC( fitschan, "Encoding" );
std::cout << "Encoding = " << encoding << "\n";
// do we have warnings?
AstKeyMap * warnings = static_cast<AstKeyMap *>( astWarnings( fitschan));
if( warnings && astOK ) {
std::cout << "Warnings:\n";
int iwarn = 1;
while( astOK ) {
std::string key = QString("Warning_%1").arg( iwarn).toStdString();
const char * message = nullptr;
if( astMapGet0C( warnings, key.c_str(), & message ) ) {
printf( "\n- %s\n", message );
hero.addError( QString( "Warning: %1").arg( message));
} else {
break;
}
}
}
else {
std::cout << "No warnings\n";
}
// create a frameset for this file
AstFrameSet * wcsinfo = static_cast<AstFrameSet *> ( astRead( fitschan ));
std::cout << "astOK = " << astOK << "\n";
if ( ! astOK ) {
std::cout << "astOK is not ok\n";
hero.addError( "astRead failed");
astClearStatus;
return heroPtr;
}
else if ( wcsinfo == AST__NULL ) {
hero.addError( "No WCS found in the fits file");
std::cout << "No WCS found\n";
return heroPtr;
}
else if ( AstWrappers::getC( wcsinfo, "Class" ) != "FrameSet") {
std::cout << "Some other weird error occured\n";
hero.addError( "AstLib returned non-frame-set");
return heroPtr;
}
// frame was read in OK, save it
hero.m_ast.origWcsInfo = wcsinfo;
astExempt( hero.m_ast.origWcsInfo);
hero.m_ast.currWcsInfo = astClone( hero.m_ast.origWcsInfo);
astExempt( hero.m_ast.currWcsInfo);
astShow( wcsinfo);
// extract the current sky system
// TODO: this assumes axis1 is a skycs
QString skysys = AstWrappers::getC( wcsinfo, "System(1)");
hero.m_currentSkyCs = string2skycs( skysys);
hero.m_originalSkyCs = hero.m_currentSkyCs;
// extract the labels/etc for axes
hero.m_ast.naxes = AstWrappers::getI( wcsinfo, "Naxes" );
hero.parseAxesInfo();
return heroPtr;
}
