file_name_20 = determine_file_name_e1('John', 'Doe', DATASET_CHAR, return_period_20) write_contour(iform_contour_20.coordinates[0][0], iform_contour_20.coordinates[0][1], folder_name + file_name_20, label_x=label_hs, label_y=label_tz) # Read the contours from the csv files. (contour_hs_1, contour_tz_1) = read_contour(folder_name + file_name_1) (contour_hs_20, contour_tz_20) = read_contour(folder_name + file_name_20) # Find datapoints that exceed the 20-yr contour. hs_outside, tz_outside, hs_inside, tz_inside = \ points_outside(contour_hs_20, contour_tz_20, np.asarray(sample_hs), np.asarray(sample_tz)) print('Number of points outside the contour: ' + str(len(hs_outside))) fig = plt.figure(figsize=(5, 5), dpi=150) ax = fig.add_subplot(111) # Plot the 1-year contour. plot_contour(x=contour_tz_1, y=contour_hs_1, ax=ax, contour_label=str(return_period_1) + '-yr contour', x_label=label_tz, y_label=label_hs, line_style='b--')
file_name_20 = determine_file_name_e1('Asta', 'Hannesdottir', DATASET_CHAR, T20) write_contour(x1_20, x0_20, folder_name + file_name_20, label_x=df.columns[1], label_y=df.columns[2]) # Read the contours from the csv files. (contour_hs_1, contour_tz_1) = read_contour(folder_name + file_name_1) (contour_hs_20, contour_tz_20) = read_contour(folder_name + file_name_20) # Find datapoints that exceed the 20-yr contour. hs_outside, tz_outside, hs_inside, tz_inside = \ points_outside(contour_hs_20, contour_tz_20, np.asarray(df[df.columns[1]].values), np.asarray(df[df.columns[2]].values)) print('Number of points outside the contour: ' + str(len(hs_outside))) #%% nan_mask = np.isnan(contour_tz_20) fig = plt.figure(figsize=(5, 5), dpi=150) ax = fig.add_subplot(111) plotted_sample = PlottedSample(x=np.asarray(sample_tz), y=np.asarray(sample_hs), ax=ax, x_inside=tz_inside, y_inside=hs_inside, x_outside=tz_outside, y_outside=hs_outside,
label_y=label_v) file_name_50 = determine_file_name_e1('John', 'Doe', DATASET_CHAR, return_period_50) write_contour(iform_contour_50.coordinates[0][0], iform_contour_50.coordinates[0][1], folder_name + file_name_50, label_x=label_hs, label_y=label_v) # Read the contour coordinates from the created csv files. (contour_hs_1, contour_v_1) = read_contour(folder_name + file_name_1) (contour_hs_50, contour_v_50) = read_contour(folder_name + file_name_50) # Find datapoints that exceed the contour. hs_outside, v_outside, hs_inside, v_inside = \ points_outside(contour_hs_50, contour_v_50, np.asarray(sample_hs), np.asarray(sample_v)) print('Number of points outside the contour: ' + str(len(hs_outside))) fig = plt.figure(figsize=(5, 5), dpi=150) ax = fig.add_subplot(111) # Plot the 1-year contour. plot_contour(x=contour_v_1, y=contour_hs_1, ax=ax, contour_label=str(return_period_1) + '-yr contour', x_label=label_v, y_label=label_hs, line_style='b--')
(contour_x_1, contour_y_1) = read_contour(folder_name + file_name_1) (contour_x_long, contour_y_long) = read_contour(folder_name + file_name_long_tr) # Switch the order of variables for plotting Hs over Tz. if dataset_char in ('A', 'B', 'C'): sample_x, sample_y = sample_y, sample_x label_x, label_y = label_y, label_x contour_x_1, contour_y_1 = contour_y_1, contour_x_1 contour_x_long, contour_y_long = contour_y_long, contour_x_long # Find datapoints that exceed the 20/50-yr contour. x_outside, y_outside, x_inside, y_inside = \ points_outside(contour_x_long, contour_y_long, np.asarray(sample_x), np.asarray(sample_y)) print('Number of points outside the contour: ' + str(len(x_outside))) fig = plt.figure(figsize=(5, 5), dpi=150) ax = fig.add_subplot(111) # Plot the 1-year contour. plot_contour(x=contour_x_1, y=contour_y_1, ax=ax, contour_label='1-yr contour', line_style='b--') # Plot the 20/50-year contour and the sample. plotted_sample = PlottedSample(x=np.asarray(sample_x),