Esempio n. 1
0
def load_file(filename):
    spec = s6fits.s6dataspec_t()
    spec.filename = filename
    spec.sortby_rfreq = 0
    spec.sortby_time = 0
    spec.sortby_bors = 0
    #spec.sortby_ifreq = 0
    # Setting up this mode will remove all those hits which shows change of RF due to
    # smart frequency switching experiment. It only reads those hits which occur maximum number of time.
    spec.filterby_rf_center_mode = 1
    #Get hits from file
    s6fits.get_s6data(spec)

    hitnum = int(s6fits.get_hits_over_file(spec.filename))

    m_rfreq = []
    m_unix_time = []
    m_meanpower = []

    base_filename = ntpath.basename(spec.filename)

    start_time = spec.s6hits[0].unix_time
    #MARK for i in range(0,hitnum-1):
    for i in range(hitnum):
        if spec.s6hits[i].unix_time == start_time:
            m_rfreq.append(spec.s6hits[i].rfreq)
            m_unix_time.append(spec.s6hits[i].unix_time)
            m_meanpower.append(spec.s6hits[i].mean_power)
        else:
            break
    sort_rfreq_index = np.argsort(m_rfreq)
    sort_rfreq = np.array(m_rfreq)[sort_rfreq_index]
    sort_meanpower = np.array(m_meanpower)[sort_rfreq_index]

    return sort_meanpower, sort_rfreq, base_filename
Esempio n. 2
0
def get_size_and_hits(filename, file_list):
    hits_size_file = []
    spec = s6fits.s6dataspec_t()
    spec.filename = filename
    s6fits.get_s6data(spec)
    total_hits = 0
    for hit in spec.s6hits:
        total_hits += 1
    hits_size_file = [
        filename, (os.stat(filename).st_size / (1024 * 1024)), total_hits
    ]
    file_list.append(hits_size_file)
    return file_list
Esempio n. 3
0
def get_size_and_hits(filename, file_list):
  hits_size_file = []
  spec = s6fits.s6dataspec_t()
  spec.filename = filename  
  s6fits.get_s6data(spec)
  total_hits = 0
  for hit in spec.s6hits:
    total_hits += 1 
  hits_size_file = [filename, 
                    (os.stat(filename).st_size/(1024*1024)),
                    total_hits]
  file_list.append(hits_size_file)  
  return file_list
Esempio n. 4
0
def main():
  if sys.argv is None:
    print "provide a valid fits file to run this"
    exit()
  """
  initialize your dataspec_t object. For now you must initialize it with no
  arguments and fill them in subsequent lines. If you only specify filename, the
  program will assume you want all the hits in the file you provide.
  """ 
  dataspec = s6fits.s6dataspec_t()
  dataspec.filename = sys.argv[1]
  dataspec.sortby_bors = 0
  dataspec.sortby_time = 0
  dataspec.sortby_ifreq = 0
  dataspec.sortby_rfreq = 0
  """
  get_s6data() will either create a new s6fits vector, or if your dataspec
  includes a vector, it will add hits to it. Right now we have no vector of hits
  we want to add to so we'll just leave it empty.
  """  
  s6fits.get_s6data(dataspec)
  s6fits.print_hits_table(dataspec.s6hits)
Esempio n. 5
0
def main():
  if sys.argv is None:
    print "provide a valid fits file to run this"
    exit()
  """
  initialize your dataspec_t object. For now you must initialize it with no
  arguments and fill them in subsequent lines. If you only specify filename, the
  program will assume you want all the hits in the file you provide.
  """ 
  dataspec = s6fits.s6dataspec_t()
  dataspec.filename = sys.argv[1]
  dataspec.sortby_bors = 0
  dataspec.sortby_time = 0
  dataspec.sortby_ifreq = 0
  dataspec.sortby_rfreq = 0
  dataspec.filterby_rf_reference_mode = 1
  """
  get_s6data() will either create a new s6fits vector, or if your dataspec
  includes a vector, it will add hits to it. Right now we have no vector of hits
  we want to add to so we'll just leave it empty.
  """  
  s6fits.get_s6data(dataspec)
  s6fits.print_hits_table(dataspec.s6hits)
Esempio n. 6
0
    plt.show()
    end_time = time.time()
    print("show plot was %g seconds" % (end_time - start_time))


if __name__ == "__main__":
    spec = s6fits.s6dataspec_t()

    spec.filename = sys.argv[1]

    #Set various data spec
    spec.sortby_rfreq = 0
    spec.sortby_time = 0
    spec.sortby_bors = 0
    spec.sortby_ifreq = 0

    # Setting up this mode will remove all those hits which shows change of RF due to
    # smart frequency switching experiment. It only reads those hits which occur maximum number of time.

    spec.filterby_rf_center_mode = 1

    #Get hits from file
    start_time = time.time()
    s6fits.get_s6data(spec)
    end_time = time.time()
    print("read time was %g seconds" % (end_time - start_time))
    #s6fits.get_s6hitsheaders(spec)

    #plot_data_with_matplotlib(spec)
    plot_data_with_gnuplot(spec)
Esempio n. 7
0
def s6fits_test(filename):
  spec = s6fits.s6dataspec_t()
  spec.filename = filename
  s6fits.get_s6data(spec)
Esempio n. 8
0
def run_s6(filename):
  spec = s6fits.s6dataspec_t()
  spec.filename = filename  
  s6fits.get_s6data(spec)
Esempio n. 9
0
import s6fits

a = s6fits.s6dataspec_t()
a.filename = "serendip6_eth2_AO_ALFA_2820_20150609_130436.fits"
s6fits.get_s6data(a)
vec = a.s6hits

for i in vec:
    print(i.ra)
Esempio n. 10
0
def s6fits_test(filename):
    spec = s6fits.s6dataspec_t()
    spec.filename = filename
    s6fits.get_s6data(spec)
Esempio n. 11
0
import s6fits

a = s6fits.s6dataspec_t()
a.filename = "serendip6_eth2_AO_ALFA_2820_20150609_130436.fits"
s6fits.get_s6data(a)
vec = a.s6hits

for i in vec:
  print (i.ra)