39, 52, 35, 51, 63, 42]) max_value = arr.max() min_value = arr.min() start, stop, step = int(min_value), int(max_value + 2), 1 histogram, bin_edge = np.histogram(arr, xrange(start, stop, step)) compute_cfs(histogram) # >>> array([ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, # 1, 2, 3, 3, 3, 3, 4, 5, 5, 6, 6, 6, 6, # 6, 7, 7, 7, 7, 9, 9, 10, 10, 11, 12, 13, 14, # 15, 15, 16, 16, 16, 16, 17, 18, 20, 21, 21, 23, 24, # 25, 26, 26, 26, 27, 27, 28, 29, 29, 30]) calculate_pdf_for_pixel(arr, histogram, bin_edge, 54) # >>> 0.066666666666666666 calculate_pdf_for_pixel(arr, histogram, bin_edge, 20) # >>> 0.0 calculate_pdf_for_pixel(arr, histogram, bin_edge, 18) # >>> 0.033333333333333333 calculate_cdf_for_pixel(arr, histogram, bin_edge, 4) # >>> 0.033333333333333333 calculate_cdf_for_pixel(arr, histogram, bin_edge, 50) # >>> 0.59999999999999998 compute_cdfs(arr, histogram, bin_edge)
41, 50, 56, 4, 54, 39, 52, 35, 51, 63, 42 ]) max_value = arr.max() min_value = arr.min() start, stop, step = int(min_value), int(max_value + 2), 1 histogram, bin_edge = np.histogram(arr, range(start, stop, step)) compute_cfs(histogram) # >>> array([ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, # 1, 2, 3, 3, 3, 3, 4, 5, 5, 6, 6, 6, 6, # 6, 7, 7, 7, 7, 9, 9, 10, 10, 11, 12, 13, 14, # 15, 15, 16, 16, 16, 16, 17, 18, 20, 21, 21, 23, 24, # 25, 26, 26, 26, 27, 27, 28, 29, 29, 30]) calculate_pdf_for_pixel(arr, histogram, bin_edge, 54) # >>> 0.066666666666666666 calculate_pdf_for_pixel(arr, histogram, bin_edge, 20) # >>> 0.0 calculate_pdf_for_pixel(arr, histogram, bin_edge, 18) # >>> 0.033333333333333333 calculate_cdf_for_pixel(arr, histogram, bin_edge, 4) # >>> 0.033333333333333333 calculate_cdf_for_pixel(arr, histogram, bin_edge, 50) # >>> 0.59999999999999998 compute_cdfs(arr, histogram, bin_edge)
arr = np.array([31, 49, 19, 62, 24, 45, 23, 51, 55, 60, 40, 35,
54, 26, 57, 37, 43, 65, 18, 41, 50, 56, 4, 54,
39, 52, 35, 51, 63, 42])
max_value = arr.max()
min_value = arr.min()
start, stop, step = int(min_value), int(max_value + 2), 1
histogram, bin_edge = np.histogram(arr, xrange(start, stop, step)
compute_cfs(histogram)
>>> array([ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 2, 3, 3, 3, 3, 4, 5, 5, 6, 6, 6, 6,
6, 7, 7, 7, 7, 9, 9, 10, 10, 11, 12, 13, 14,
15, 15, 16, 16, 16, 16, 17, 18, 20, 21, 21, 23, 24,
25, 26, 26, 26, 27, 27, 28, 29, 29, 30])
calculate_pdf_for_pixel(arr, histogram, bin_edge, 54)
>>> 0.066666666666666666
calculate_pdf_for_pixel(arr, histogram, bin_edge, 20)
>>> 0.0
calculate_pdf_for_pixel(arr, histogram, bin_edge, 18)
>>> 0.033333333333333333
calculate_cdf_for_pixel(arr, histogram, bin_edge, 4)
>>> 0.033333333333333333
calculate_cdf_for_pixel(arr, histogram, bin_edge, 50)
>>> 0.59999999999999998
compute_cdfs(arr, histogram, bin_edge)
>>> array([ 0.03333333, 0.03333333, 0.03333333, 0.03333333, 0.03333333,
0.03333333, 0.03333333, 0.03333333, 0.03333333, 0.03333333,
0.03333333, 0.03333333, 0.03333333, 0.03333333, 0.06666667,
0.1 , 0.1 , 0.1 , 0.1 , 0.13333333,