p61743098
/
sagecal
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

You can now insert sources at random positions within some angle from 3C196

This commit is contained in:
Hanno Spreeuw 2017-10-20 15:45:48 +02:00
parent fd32097a60
commit 402efc6f17
1 changed files with 38 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

52
test/Generate_sources.py
View File

@ -21,16 +21,36 @@ except AssertionError:
# Best to center sources around 3C196 if using the sm.ms measurement set with a
# four degree tolerance
tol_degrees = 4.
tol_hours = tol_degrees * 24/360
tol_seconds_of_decl = 0.9
tol_seconds_of_RA = tol_seconds_of_decl * 24/360
ra_3C196 = 8.
ra_low = ra_3C196 - tol_hours
ra_high = ra_3C196 + tol_hours
RA_hours_3C196 = 8
RA_minutes_3C196 = 13
RA_seconds_3C196 = 35.981540
decl_3C196 = 48.
decl_low = decl_3C196 - tol_degrees
decl_high = decl_3C196 + tol_degrees
RA_seconds_3C196 = (RA_hours_3C196 * 60 + RA_minutes_3C196) *60 + RA_seconds_3C196
RA_seconds_low = RA_seconds_3C196 - tol_seconds_of_RA
RA_minutes_low, RA_seconds_low = divmod(RA_seconds_low, 60)
RA_hours_low, RA_minutes_low = divmod(RA_minutes_low, 60)
RA_seconds_high = RA_seconds_3C196 + tol_seconds_of_RA
RA_minutes_high, RA_seconds_high = divmod(RA_seconds_high, 60)
RA_hours_high, RA_minutes_high = divmod(RA_minutes_high, 60)
decl_degrees_3C196 = 48
decl_minutes_3C196 = 12
decl_seconds_3C196 = 59.174770
decl_seconds_3C196 = (decl_degrees_3C196 * 60 + decl_minutes_3C196) *60 + decl_seconds_3C196
decl_seconds_low = decl_seconds_3C196 - tol_seconds_of_decl
decl_minutes_low, decl_seconds_low = divmod(decl_seconds_low, 60)
decl_degrees_low, decl_minutes_low = divmod(decl_minutes_low, 60)
decl_seconds_high = decl_seconds_3C196 + tol_seconds_of_decl
decl_minutes_high, decl_seconds_high = divmod(decl_seconds_high, 60)
decl_degrees_high, decl_minutes_high = divmod(decl_minutes_high, 60)
I_low = 10
I_high = 100
@ -58,13 +78,17 @@ with warnings.catch_warnings():
sources_parameters.name[source_name] = 'P' + str(source_name).zfill(number_of_digits_for_sources)
# Right ascension can have all values.
sources_parameters.rah = np.random.randint(ra_low, ra_high, size=number_of_sources)
sources_parameters.ram = np.random.randint(0, 59, size=number_of_sources)
sources_parameters.ras = 60 * np.random.rand(number_of_sources)
# sources_parameters.rah = np.random.randint(ra_low, ra_high, size=number_of_sources)
# sources_parameters.ram = np.random.randint(0, 59, size=number_of_sources)
# sources_parameters.ras = 60 * np.random.rand(number_of_sources)
sources_parameters.dad = np.random.randint(decl_low, decl_high, size=number_of_sources)
sources_parameters.dam = np.random.randint(0, 59, size=number_of_sources)
sources_parameters.das = 60 * np.random.rand(number_of_sources)
sources_parameters.rah = (RA_hours_high - RA_hours_low) * np.random.random_sample(number_of_sources) + RA_hours_low
sources_parameters.ram = (RA_minutes_high - RA_minutes_low) * np.random.random_sample(number_of_sources) + RA_minutes_low
sources_parameters.ras = (RA_seconds_high - RA_seconds_low) * np.random.random_sample(number_of_sources) + RA_seconds_low
sources_parameters.dad = (decl_degrees_high - decl_degrees_low) * np.random.random_sample(number_of_sources) + decl_degrees_low
sources_parameters.dam = (decl_minutes_high - decl_minutes_low) * np.random.random_sample(number_of_sources) + decl_minutes_low
sources_parameters.das = (decl_seconds_high - decl_seconds_low) * np.random.random_sample(number_of_sources) + decl_seconds_low
sources_parameters.I = (I_high - I_low) * np.random.rand(number_of_sources) + I_low
sources_parameters.Q = 0