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
|
include "identify.h"
define NBIN 10 # Bin parameter for mode determination
# ID_SHIFT -- Determine a shift by correlating feature user positions
# with peaks in the image data.
double procedure id_shift (id)
pointer id # ID pointer
int i, j, npeaks, ndiff, find_peaks()
real d, dmin
double pix, id_center(), id_fitpt()
pointer x, y, diff
errchk malloc, find_peaks
begin
# Find the peaks in the image data and center.
call malloc (x, ID_NPTS(id), TY_REAL)
npeaks = find_peaks (IMDATA(id,1), Memr[x], ID_NPTS(id), 0.,
int (ID_MINSEP(id)), 0, ID_MAXFEATURES(id), 0., false)
# Center the peaks and convert to user coordinates.
call malloc (y, npeaks, TY_DOUBLE)
j = 0
do i = 1, npeaks {
pix = id_center (id, double(Memr[x+i-1]), 1, ID_FWIDTH(id),
ID_FTYPE(id), NO)
if (!IS_INDEFD (pix)) {
Memd[y+j] = id_fitpt (id, pix)
j = j + 1
}
}
npeaks = j
# Compute differences with feature list.
ndiff = npeaks * ID_NFEATURES(id)
call malloc (diff, ndiff, TY_REAL)
ndiff = 0
do i = 1, ID_NFEATURES(id) {
do j = 1, npeaks {
Memr[diff+ndiff] = Memd[y+j-1] - FIT(id,i)
ndiff = ndiff + 1
}
}
call mfree (x, TY_REAL)
call mfree (y, TY_DOUBLE)
# Sort the differences and find the mode.
call asrtr (Memr[diff], Memr[diff], ndiff)
dmin = Memr[diff+ndiff-1] - Memr[diff]
do i = 0, ndiff-NBIN-1 {
j = i + NBIN
d = Memr[diff+j] - Memr[diff+i]
if (d < dmin) {
dmin = d
pix = Memr[diff+i] + d / 2.
}
}
call mfree (diff, TY_REAL)
return (pix)
end
|
