def test_active_regions(ar_id, grav_long, grav_lat):
# Active regions position testing
if __name__ == '__main__':
ar_data = DataAccess('2003-09-26T00:00:00', '2003-09-26T01:00:00', 'AR', 'SOHO', 'MDI')
ar_chain_encoded = prep.decode_and_split(ar_data.get_chain_code())
ar_id = ar_data.get_ar_id()
ar_centers_lon = ar_data.get_grav_center_long()
ar_centers_lat = ar_data.get_grav_center_lat()
ar = ar_carr_synthesis[3]
prep.display_object([ar], "")
ar2 = np.array(ar)
x = ar2[0]
y = ar2[1]
length = len(ar2[0]) + len(ar2[1])
centroid = (sum(x) / len(ar2[0]), sum(y) / len(ar2[1]))
print(centroid)
# if the longitude difference is smaller than width of ar
# and latitude is smaller than height of ar then is success
# otherwise fail
if lon_diff < width and lat_diff < height:
print("SUCCESS", lon_diff)
success += 1
else:
print("FAIL", lon_diff)
print("long_center:", long)
print("calculated:", centroid[0])
fail += 1
# prep.display_object(carr_coords, [])
else:
ignored += 1
print("successes = ", success)
print("fail = ", fail)
print("ignored = ", ignored)
# Active regions position testing
if __name__ == '__main__':
ar_data = DataAccess('2003-09-28T00:00:00', '2003-10-23T01:00:00', 'AR',
'SOHO', 'MDI')
ar_chain_encoded = prep.decode_and_split(ar_data.get_chain_code())
ar_id = ar_data.get_ar_id()
ar_centers_lon = ar_data.get_grav_center_long()
ar_centers_lat = ar_data.get_grav_center_lat()
test_active_regions(ar_id, ar_centers_lon, ar_centers_lat)
lon_diff = centroid[0] - long
lat_diff = centroid[1] - lat
# if the longitude difference is smaller than 1.0
# and latitude is smaller than 1.0 then is success
# otherwise fail
if lon_diff < 1.0 and lat_diff < 1.0:
print("SUCCESS", lon_diff)
success += 1
else:
print("FAIL", lon_diff)
print("long_center:", long)
print("calculated:", centroid[0])
fail += 1
else:
ignored += 1
print("successes = ", success)
print("fail = ", fail)
print("ignored = ", ignored)
# Sunspots position testing
if __name__ == '__main__':
sp_data = DataAccess('2003-09-28T00:00:00', '2003-10-23T01:00:00', 'SP', 'SOHO', 'MDI')
sp_chain_encoded = prep.decode_and_split(sp_data.get_chain_code())
sp_id = sp_data.get_sp_id()
sp_centers_lon = sp_data.get_grav_center_long()
sp_centers_lat = sp_data.get_grav_center_lat()
test_sunspots(sp_id, sp_centers_lon, sp_centers_lat)
