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)