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);
