def generate_image(filename, sourcename, band, size): f = open(sourcename, 'r') fo = open(filename, 'w') fo.write(f.readline()) #fo.write(f.readline()) for _ in range(size): z = f.readline().split() for ik in range(2): z[ik] = random.normal(z[ik], ERROR) fo.write('%s %s %s %s\n' % (z[0], z[1], z[2], 0.01)) parsetname = path.basename(filename) parsetname = parsetname[:parsetname.index('.')] + '.parset' write_parset(parsetname, filename, FREQUENCY[band], 180.0, 0.0, 5.0) f.close() fo.close()
def generate_image(filename, sourcename, band, size): f = open(sourcename, 'r') fo = open(filename, 'w') fo.write(f.readline()) #fo.write(f.readline()) for _ in xrange(size): z = f.readline().split() for ik in xrange(2): z[ik] = random.normal(z[ik], ERROR) fo.write('%s %s %s %s\n' % (z[0], z[1], z[2], 0.01)) parsetname = path.basename(filename) parsetname = parsetname[:parsetname.index('.')] + '.parset' write_parset(parsetname, filename, FREQUENCY[band], 180.0, 0.0, 5.0) f.close() fo.close()
def generate_snapshot(filename): conn = GSMConnectionManager(database='stress', use_monet=False).get_connection() decl_center = -50 while decl_center < -30: decl_center = degrees(acos(2 * random.random() - 1) - 0.5 * pi) ra_center = degrees(2 * pi * random.random()) band = random.random_integers(1, 3) sql = get_field(ra_center, decl_center, 3.0, 2) cur = conn.get_cursor(sql) f = open(filename, 'w') f.write('# RA DEC Total_flux e_Total_flux\n\n') for data in iter(cur.fetchone, None): f.write('%s %s %s %s\n' % (random.normal(data[0], ERROR), random.normal(data[1], ERROR), random.normal(data[2], FLUX_ERROR), data[3])) f.close() parsetname = path.basename(filename) parsetname = parsetname[:parsetname.index('.')] + '.parset' write_parset(parsetname, filename, FREQUENCY[band], ra_center, decl_center)
def generate_field_parset(filename, ra, decl, radius, size, frequency=8): generate_field_file(filename, ra, decl, radius, size) parsetname = path.basename(filename) parsetname = parsetname[:parsetname.index('.')] + '.parset' write_parset(parsetname, filename, FREQUENCY[frequency], ra, decl, radius)