Example #1
0
                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)
Example #2
0
    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,