blob: cd75445f42927d0523d82b9f4e5f1c10ac37098f (
plain) (
blame)
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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <mach.h>
# AHGM -- Accumulate the histogram of the input vector. The output vector
# HGM (the histogram) should be cleared prior to the first call.
procedure ahgmd (data, npix, hgm, nbins, z1, z2)
double data[ARB] # data vector
int npix # number of pixels
int hgm[ARB] # output histogram
int nbins # number of bins in histogram
double z1, z2 # greyscale values of first and last bins
double z
real dz
int bin, i
begin
dz = real (nbins - 1) / real (z2 - z1)
if (abs (dz - 1.0) < (EPSILONR * 2.0)) {
do i = 1, npix {
z = data[i]
if (z >= z1 && z <= z2) {
bin = int (z - z1) + 1
hgm[bin] = hgm[bin] + 1
}
}
} else {
do i = 1, npix {
z = data[i]
if (z >= z1 && z <= z2) {
bin = int ((z - z1) * dz) + 1
hgm[bin] = hgm[bin] + 1
}
}
}
end
|
