use_db = "mysql-Q" # 'sqlite' Use local sqlite file-based db # 'mysql-Q' Use the remote MySQL database on Q user = "******" # only needed for remote databases. password = "" # See example109 for setting-up an encrypted password. In this case leave password="", and # init_db_conn_old() will automatically find and use your saved password. Otherwise, enter your MySQL password here. # ------------ NO NEED TO UPDATE ANYTHING BELOW ------------- # start_time_tot = time.time() # connect to database if use_db == 'sqlite': # setup database connection to local sqlite file sqlite_path = os.path.join(database_dir, sqlite_filename) db_session = init_db_conn_old(db_name=use_db, chd_base=db_class.Base, sqlite_path=sqlite_path) elif use_db == 'mysql-Q': # setup database connection to MySQL database on Q db_session = init_db_conn_old(db_name=use_db, chd_base=db_class.Base, user=user, password=password) # create mu bin array mu_bin_array = np.array(range(n_mu_bins + 1), dtype="float") * 0.05 + 0.1 mu_bin_centers = (mu_bin_array[1:] + mu_bin_array[:-1]) / 2 # time arrays # returns array of moving averages center dates, based off start and end date moving_avg_centers, moving_width = lbcc.moving_averages(query_time_min,
# set movie frames-per-second fps = 12 # one day per second # set frame dots per inch dpi = 300 # set to None for automatic dpi estimation # designate which database to connect to use_db = "mysql-Q" # 'sqlite' Use local sqlite file-based db # 'mysql-Q' Use the remote MySQL database on Q # 'mysql-Q_test' Use the development database on Q user = "******" # only needed for remote databases. password = "" # See example109 for setting-up an encrypted password. In this case leave password="", and # init_db_conn_old() will automatically find and use your saved password. Otherwise, enter your MySQL password here. # Establish connection to database db_session = db_funs.init_db_conn_old(db_name=use_db, chd_base=DBClass.Base, user=user, password=password) # --- Begin execution ---------------------- # estimate dates for start/end of rotation CR_start = datetime.datetime(1853, 11, 9) CR_dur = datetime.timedelta(days=27.2753) estimated_center_date = CR_start + (cr_rot - 0.5) * CR_dur # get the decimal Carrington number for Earth at this time cr_earth = sunpy.coordinates.sun.carrington_rotation_number( estimated_center_date) # calc reasonable estimate of rotation start/end movie_start = estimated_center_date - (cr_earth - cr_rot) * CR_dur movie_end = estimated_center_date + (cr_rot + 1 - cr_earth) * CR_dur # query maps in time range
n_intensity_bins = 200 R0 = 1.01 log10 = True lat_band = [-np.pi / 2.4, np.pi / 2.4] # define database paths raw_data_dir = App.RAW_DATA_HOME hdf_data_dir = App.PROCESSED_DATA_HOME database_dir = App.DATABASE_HOME sqlite_filename = App.DATABASE_FNAME # setup database connection use_db = "sqlite" sqlite_path = os.path.join(database_dir, sqlite_filename) db_session = init_db_conn_old(db_name=use_db, chd_base=db_class.Base, sqlite_path=sqlite_path) # use_db = "mysql-Q" # user = "******" # password = "" # sqlite_path = os.path.join(database_dir, sqlite_filename) # db_session = db_funcs.init_db_conn_old(db_name=use_db, chd_base=db_class.Base, user=user, password=password) #### ----- GENERATE HISTOGRAMS ---- #### # start time start_time = time.time() #### GET REFERENCE INFO FOR LATER USE #### #### GET REFERENCE INFO FOR LATER USE #### # get index number of reference instrument ref_index = inst_list.index(ref_inst)
HISTORY: - This example was adapted by CD from commented code by TE in image_testing.py on 2021/01/04 """ import pandas as pd import datetime from chmap.database.db_funs import init_db_conn_old, update_image_val, query_euv_images from chmap.database.db_classes import Base, EUV_Images from chmap.utilities.file_io.io_helpers import print_full_dataframe # setup the mysql session use_db = "mysql-Q" user = "******" password = "" db_session = init_db_conn_old(db_name=use_db, chd_base=Base, user=user, password=password) # columns to print when running the script (so i can read it...) columns = [ 'image_id', 'date_obs', 'instrument', 'wavelength', 'flag', 'fname_raw', 'fname_hdf' ] # OPTIONAL query a specific date to get image_ids of UNFLAGGED IMAGES (right now query_euv_images only looks for flag=0) # set this to False if you want the script to continue do_check_ids = False if do_check_ids: period_start = datetime.datetime(2020, 2, 24, 1, 0, 0) period_end = datetime.datetime(2020, 2, 24, 3, 0, 0) query = query_euv_images(db_session=db_session,