def get_currentDay_previusDay_noaa_activeRegions(new_feat_objects, noaa_data, get_df = False): ''' Output: nova_cents_labels_sameDay, tuple nova_cents_labels_prevDay, tuple ''' import pandas as pd # get noaa file from same date as image data same_day_df, year, month, day = get_current_day_noaa_df(new_feat_objects, noaa_data) current_day = year + month + day # these are NOAA active region centroid nova_cents_labels_sameDay, sameDay_df = extract.get_noaa_centroids(same_day_df, current_day) previous_day_df, previous_day = get_prevous_day_noaa_df(noaa_data, year, month, day) # these are NOAA active region centroid nova_cents_labels_prevDay, prevDay_df = extract.get_noaa_centroids(previous_day_df, previous_day) if get_df == False : return nova_cents_labels_sameDay, nova_cents_labels_prevDay else: return nova_cents_sameDay, nova_cents_prevDay, same_day_df, previous_day_df
def plot_sunspots_and_active_regions(df, scan_year, features, time_slice): ''' use scan_year to shift through noaa observations use time_slice to scan through image data ''' noaa_cents, _ = extract.get_noaa_centroids(df, scan_year) noaa_x, noaa_y = unpack_noaa_cents(noaa_cents) x_cents, y_cents= extract.get_image_active_region_centroids(features[time_slice]) plt.figure(figsize=(10,10)) noaa = plt.scatter(noaa_x, noaa_y, c='b', marker='o'); me = plt.scatter(x_cents, y_cents , c='r',marker='+'); plt.title("Sunspots & Active Regions " + scan_year); plt.legend((me, noaa), ('mydata-sunspots', 'noaa-AR'), scatterpoints=1, loc='lower right', ncol=2, fontsize=15);