1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
|
;
; Josh Sixsmith 16 May 2012, under GNU GPL v3 or later
;
; A simple statistical test for distribution shape
;
; After HIST_EQUAL we check that:
; -- kurtosis of the image is roughly -1.2
; -- skewness of the image is roughly 0
; This should describe a Uniform distribution
;
pro TEST_HIST_EQUAL, samples=samples, test=test, verbose=verbose, help=help
;
if KEYWORD_SET(help) then begin
print, 'pro TEST_HIST_EQUAL, samples=samples, test=test, $'
print, ' verbose=verbose. help=help'
return
endif
;
if N_ELEMENTS(samples) eq 0 then samples=100
;
; We will calculate an image from a random normal distribution.
; This should (roughly) have a skewness of 0 and a kurtosis of 0.
img = RANDOMN(sd, samples,samples)
;
; compute the histogram equalisation on the image
scl = HIST_EQUAL(img)
;
; compute the moment of the raw image and the scaled image
m_raw = MOMENT(img)
m_scl = MOMENT(scl)
;
; The more samples used in the image creation, the closer the
; image is to being a normal distribution.
; As such, some variation will be expected, so an error of 0.1
; should suffice.
err = 0.1
total_pb = 0
;
if (ABS(m_scl[3] + 1.2) ge err) then begin
if KEYWORD_SET(verbose) then MESSAGE, 'pb with kurtosis ofimage', /continue
total_pb = total_pb + 1
endif
;
if (ABS(m_scl[2]) ge err) then begin
if KEYWORD_SET(verbose) then MESSAGE, 'pb with skewness of image', /continue
total_pb = total_pb + 1
endif
;
if KEYWORD_SET(verbose) then begin
print, 'Stats on Raw Image:'
print, 'Mean: ', m_raw[0] & print, 'Variance: ', m_raw[1]
print, 'Skewness: ', m_raw[2] & print, 'Kurtosis: ', m_raw[3]
print, 'Stats on Equalised Image:'
print, 'Mean: ', m_scl[0] & print, 'Variance: ', m_scl[1]
print, 'Skewness: ', m_scl[2] & print, 'Kurtosis: ', m_scl[3]
endif
;
if KEYWORD_SET(test) then stop
;
if (total_pb gt 0) then begin
if KEYWORD_SET(no_exit) then begin
MESSAGE, 'Tests failed. Increase sample SIZE?', /continue
endif else begin
EXIT, status=1
endelse
endif else begin
MESSAGE, 'Tests successful.', /continue
endelse
;
end
|
