/**
    @class SolImageSignal.cpp
    @author Brian Magill
    @datecreated 3/8/2006

	$Date: 2006/04/06 18:25:15 $
    $Revision: 1.2 $

    @copyright
     (©) Copyright 2006 by GATS Inc.
     11864 Canon Blvd., Suite 101, Newport News, VA 23606

     All Rights Reserved. No part of this software or publication may be
     reproduced, stored in a retrieval system, or transmitted, in any form
     or by any means, electronic, mechanical, photocopying, recording, or
     otherwise without the prior written permission of GATS Inc.

 
      @brief class for a sequence of Solar Image data
*/
//----------------------------------------------------------------------
//

#include <iostream>
#include <stdexcept>
//#include "SolarImSignal.h"
#include "SolImageSignal.h"
//#include "SOFIE_Exception.h"
#include "SOFIE_namespace.h"

using namespace SOFIE;
using namespace std;


    SolImageSignal& SolImageSignal::operator = (const SolImageSignal& rhs ) {   

        if (this == &rhs)
            return *this;
                 
        SolImageInfo = rhs.SolImageInfo;            

        return *this;        
        }


    void SolImageSignal::dump() const
        {
    
       vector<SolImageSample>::const_iterator Ind;

        for (Ind = SolImageInfo.begin(); Ind != SolImageInfo.end(); Ind++)
            {
            (*Ind).dump();
            }            
        cout << endl;                
        }
        
    valarray<double> SolImageSignal::getTime() const
        {
        
        unsigned long count = SolImageInfo.size();
        valarray<double> time(count);
        unsigned long Ind;

        for (Ind = 0; Ind < count; Ind++)
            {
            time[Ind] = SolImageInfo[Ind].getTime();
            }            

        return time;
        }

    vector< valarray<int> > SolImageSignal::SolarData() const
        {
        
        valarray <int> temp(SolarImgSampleSize);                
        vector< valarray <int> > solar_avg;
        vector<SolImageSample>::const_iterator Ind;

        for (Ind = SolImageInfo.begin(); Ind != SolImageInfo.end(); Ind++)
            {
            temp = (*Ind).getSample();
            solar_avg.push_back(temp);
            }            

        return solar_avg;                
        }   


std::valarray<int> SolImageSignal::getSolarData(unsigned long indx) const
    {
        unsigned long j;
        unsigned long cnt = SolarImgSampleSize;
        valarray<int> image(cnt);
               
        for (j = 0; j < cnt; ++j)
            image[j] = SolImageInfo.at(indx).getImageData(j);

        return image;
    
    };