Practical Algorithms for Image Analysis

home *** CD-ROM | disk | FTP | other *** search
/ Practical Algorithms for Image Analysis / Practical Algorithms for Image Analysis.iso / TARFILE.GZ / tarfile / ch_3.5 / bcd / edge < prev next >
SHell self-extracting ARchive | 1999-09-11 | 36.4 KB
open in:
MacOS 8.1 extracted |
Win98 extracted |
DOS extracted
view JSON data |
view as text

This file was processed as: SHell self-extracting ARchive (archive/shar).

You can browse this item here: edge
Confidence	Program	Detection	Match Type	Support
`100%`	dexvert	SHell self-extracting ARchive `(archive/shar)`	magic	Supported
`1%`	dexvert	Text File `(text/txt)`	fallback	Supported
`100%`	file	data	default
`100%`	TrID	shar SHell self-extracting aRchive (var.2)	default
`100%`	perlTextCheck	Likely Text (Perl)	default
`100%`	detectItEasy	Format: plain text[CRLF]	default (weak)
hex view
+--------+-------------------------+-------------------------+--------+--------+
|00000000| 23 20 54 6f 20 75 6e 62 | 75 6e 64 6c 65 2c 20 73 |# To unb|undle, s|
|00000010| 68 20 74 68 69 73 20 66 | 69 6c 65 0d 0a 65 63 68 |h this f|ile..ech|
|00000020| 6f 20 4d 61 6b 65 66 69 | 6c 65 20 31 3e 26 32 0d |o Makefi|le 1>&2.|
|00000030| 0a 73 65 64 20 27 73 2f | 2e 2f 2f 27 20 3e 4d 61 |.sed 's/|.//' >Ma|
|00000040| 6b 65 66 69 6c 65 20 3c | 3c 27 2f 2f 47 4f 2e 53 |kefile <|<'//GO.S|
|00000050| 59 53 49 4e 20 44 44 20 | 4d 61 6b 65 66 69 6c 65 |YSIN DD |Makefile|
|00000060| 27 0d 0a 2d 43 43 20 3d | 20 63 63 0d 0a 2d 23 43 |'..-CC =| cc..-#C|
|00000070| 46 4c 41 47 53 20 3d 20 | 2d 67 20 2d 44 44 45 42 |FLAGS = |-g -DDEB|
|00000080| 55 47 20 2d 66 36 38 38 | 38 31 20 23 66 6f 72 20 |UG -f688|81 #for |
|00000090| 64 65 62 75 67 67 69 6e | 67 0d 0a 2d 23 43 46 4c |debuggin|g..-#CFL|
|000000a0| 41 47 53 20 3d 20 2d 4f | 34 20 23 20 66 6f 72 20 |AGS = -O|4 # for |
|000000b0| 73 75 6e 34 0d 0a 2d 23 | 43 46 4c 41 47 53 20 3d |sun4..-#|CFLAGS =|
|000000c0| 20 2d 4f 20 2d 66 36 38 | 38 38 31 20 23 66 6f 72 | -O -f68|881 #for|
|000000d0| 20 73 75 6e 33 0d 0a 2d | 43 46 4c 41 47 53 20 3d | sun3..-|CFLAGS =|
|000000e0| 20 2d 4f 33 20 23 20 53 | 69 6c 69 63 6f 6e 20 47 | -O3 # S|ilicon G|
|000000f0| 72 61 70 68 69 63 73 0d | 0a 2d 43 4c 49 42 53 20 |raphics.|.-CLIBS |
|00000100| 3d 20 2d 6c 6d 0d 0a 2d | 4f 42 4a 45 43 54 53 20 |= -lm..-|OBJECTS |
|00000110| 3d 20 63 6f 6e 76 2e 6f | 20 63 74 6f 72 2e 6f 20 |= conv.o| ctor.o |
|00000120| 64 65 74 2e 6f 20 64 74 | 6f 72 2e 6f 20 65 64 67 |det.o dt|or.o edg|
|00000130| 65 5f 6d 61 70 73 2e 6f | 20 66 69 6c 74 65 72 2e |e_maps.o| filter.|
|00000140| 6f 20 5c 0d 0a 2d 09 09 | 09 66 69 6e 64 5f 65 64 |o \..-..|.find_ed|
|00000150| 67 65 73 2e 6f 20 20 68 | 79 73 74 2e 6f 20 6c 6f |ges.o  h|yst.o lo|
|00000160| 63 2e 6f 20 74 68 69 6e | 2e 6f 20 74 68 6c 64 2e |c.o thin|.o thld.|
|00000170| 6f 20 77 72 69 74 65 2e | 6f 0d 0a 2d 0d 0a 2d 62 |o write.|o..-..-b|
|00000180| 63 64 3a 20 62 63 64 2e | 6f 20 24 28 4f 42 4a 45 |cd: bcd.|o $(OBJE|
|00000190| 43 54 53 29 0d 0a 2d 09 | 24 28 43 43 29 20 2d 6f |CTS)..-.|$(CC) -o|
|000001a0| 20 62 63 64 20 24 28 43 | 46 4c 41 47 53 29 20 62 | bcd $(C|FLAGS) b|
|000001b0| 63 64 2e 6f 20 24 28 4f | 42 4a 45 43 54 53 29 20 |cd.o $(O|BJECTS) |
|000001c0| 24 28 43 4c 49 42 53 29 | 0d 0a 2d 0d 0a 2d 62 63 |$(CLIBS)|..-..-bc|
|000001d0| 64 2e 6f 3a 20 65 64 67 | 65 5f 66 69 6e 64 65 72 |d.o: edg|e_finder|
|000001e0| 2e 68 0d 0a 2d 0d 0a 2d | 63 62 2e 6f 3a 20 65 64 |.h..-..-|cb.o: ed|
|000001f0| 67 65 5f 66 69 6e 64 65 | 72 2e 68 0d 0a 2d 0d 0a |ge_finde|r.h..-..|
|00000200| 2d 63 6f 6e 76 2e 6f 3a | 20 65 64 67 65 5f 66 69 |-conv.o:| edge_fi|
|00000210| 6e 64 65 72 2e 68 0d 0a | 2d 0d 0a 2d 63 74 6f 72 |nder.h..|-..-ctor|
|00000220| 2e 6f 3a 20 65 64 67 65 | 5f 66 69 6e 64 65 72 2e |.o: edge|_finder.|
|00000230| 68 0d 0a 2d 0d 0a 2d 64 | 65 74 2e 6f 3a 20 65 64 |h..-..-d|et.o: ed|
|00000240| 67 65 5f 66 69 6e 64 65 | 72 2e 68 0d 0a 2d 0d 0a |ge_finde|r.h..-..|
|00000250| 2d 64 74 6f 72 2e 6f 3a | 20 65 64 67 65 5f 66 69 |-dtor.o:| edge_fi|
|00000260| 6e 64 65 72 2e 68 0d 0a | 2d 0d 0a 2d 65 64 67 65 |nder.h..|-..-edge|
|00000270| 5f 6d 61 70 73 2e 6f 3a | 20 65 64 67 65 5f 66 69 |_maps.o:| edge_fi|
|00000280| 6e 64 65 72 2e 68 0d 0a | 2d 0d 0a 2d 66 69 6c 74 |nder.h..|-..-filt|
|00000290| 65 72 2e 6f 3a 20 65 64 | 67 65 5f 66 69 6e 64 65 |er.o: ed|ge_finde|
|000002a0| 72 2e 68 0d 0a 2d 0d 0a | 2d 66 69 6e 64 5f 65 64 |r.h..-..|-find_ed|
|000002b0| 67 65 73 2e 6f 3a 20 65 | 64 67 65 5f 66 69 6e 64 |ges.o: e|dge_find|
|000002c0| 65 72 2e 68 0d 0a 2d 0d | 0a 2d 68 79 73 74 2e 6f |er.h..-.|.-hyst.o|
|000002d0| 3a 20 65 64 67 65 5f 66 | 69 6e 64 65 72 2e 68 0d |: edge_f|inder.h.|
|000002e0| 0a 2d 0d 0a 2d 6c 6f 63 | 2e 6f 3a 20 65 64 67 65 |.-..-loc|.o: edge|
|000002f0| 5f 66 69 6e 64 65 72 2e | 68 0d 0a 2d 0d 0a 2d 74 |_finder.|h..-..-t|
|00000300| 68 69 6e 2e 6f 3a 20 65 | 64 67 65 5f 66 69 6e 64 |hin.o: e|dge_find|
|00000310| 65 72 2e 68 0d 0a 2d 0d | 0a 2d 74 68 6c 64 2e 6f |er.h..-.|.-thld.o|
|00000320| 3a 20 65 64 67 65 5f 66 | 69 6e 64 65 72 2e 68 0d |: edge_f|inder.h.|
|00000330| 0a 2d 0d 0a 2d 77 72 69 | 74 65 2e 6f 3a 20 65 64 |.-..-wri|te.o: ed|
|00000340| 67 65 5f 66 69 6e 64 65 | 72 2e 68 0d 0a 2d 0d 0a |ge_finde|r.h..-..|
|00000350| 2d 63 6c 65 61 6e 3a 0d | 0a 2d 09 72 6d 20 2a 2e |-clean:.|.-.rm *.|
|00000360| 6f 0d 0a 2d 0d 0a 2d 63 | 6c 65 61 6e 69 6d 61 67 |o..-..-c|leanimag|
|00000370| 65 73 3a 0d 0a 2d 09 72 | 6d 20 64 65 74 31 33 35 |es:..-.r|m det135|
|00000380| 20 64 65 74 34 35 20 64 | 65 74 78 20 64 65 74 79 | det45 d|etx dety|
|00000390| 20 6d 61 70 20 74 68 69 | 6e 5f 6d 61 70 20 68 79 | map thi|n_map hy|
|000003a0| 73 74 5f 6d 61 70 20 68 | 69 5f 6d 61 70 20 67 61 |st_map h|i_map ga|
|000003b0| 75 73 73 0d 0a 2d 09 0d | 0a 2f 2f 47 4f 2e 53 59 |uss..-..|.//GO.SY|
|000003c0| 53 49 4e 20 44 44 20 4d | 61 6b 65 66 69 6c 65 0d |SIN DD M|akefile.|
|000003d0| 0a 65 63 68 6f 20 52 45 | 41 44 4d 45 20 31 3e 26 |.echo RE|ADME 1>&|
|000003e0| 32 0d 0a 73 65 64 20 27 | 73 2f 2e 2f 2f 27 20 3e |2..sed '|s/.//' >|
|000003f0| 52 45 41 44 4d 45 20 3c | 3c 27 2f 2f 47 4f 2e 53 |README <|<'//GO.S|
|00000400| 59 53 49 4e 20 44 44 20 | 52 45 41 44 4d 45 27 0d |YSIN DD |README'.|
|00000410| 0a 2d 2f 2a 0d 0a 2d 20 | 2a 20 28 63 29 20 43 6f |.-/*..- |* (c) Co|
|00000420| 70 79 72 69 67 68 74 20 | 31 39 38 38 20 62 79 0d |pyright |1988 by.|
|00000430| 0a 2d 20 2a 20 52 6f 62 | 6f 74 69 63 73 20 50 72 |.- * Rob|otics Pr|
|00000440| 69 6e 63 69 70 6c 65 73 | 20 52 65 73 65 61 72 63 |inciples| Researc|
|00000450| 68 20 44 65 70 61 72 74 | 6d 65 6e 74 2c 20 41 54 |h Depart|ment, AT|
|00000460| 54 20 42 65 6c 6c 20 4c | 61 62 6f 72 61 74 6f 72 |T Bell L|aborator|
|00000470| 69 65 73 2e 0d 0a 2d 20 | 2a 20 41 6c 6c 20 72 69 |ies...- |* All ri|
|00000480| 67 68 74 73 20 72 65 73 | 65 72 76 65 64 2e 0d 0a |ghts res|erved...|
|00000490| 2d 20 2a 20 4c 61 73 74 | 20 6d 6f 64 69 66 69 65 |- * Last| modifie|
|000004a0| 64 20 36 2f 32 31 2f 38 | 38 20 49 6e 67 65 6d 61 |d 6/21/8|8 Ingema|
|000004b0| 72 20 4a 2e 20 43 6f 78 | 0d 0a 2d 20 2a 20 52 65 |r J. Cox|..- * Re|
|000004c0| 76 69 73 65 64 20 31 31 | 2f 31 31 2f 38 38 20 49 |vised 11|/11/88 I|
|000004d0| 6e 67 65 6d 61 72 20 4a | 2e 20 63 6f 78 0d 0a 2d |ngemar J|. cox..-|
|000004e0| 20 2a 2f 0d 0a 2d 20 20 | 20 20 20 20 20 20 20 20 | */..-  |        |
|000004f0| 20 20 20 20 20 20 20 20 | 20 20 20 42 4f 49 45 2d |        |   BOIE-|
|00000500| 43 4f 58 20 45 44 47 45 | 20 44 45 54 45 43 54 4f |COX EDGE| DETECTO|
|00000510| 52 0d 0a 2d 0d 0a 2d 54 | 77 6f 20 6d 61 69 6e 20 |R..-..-T|wo main |
|00000520| 70 72 6f 67 72 61 6d 73 | 2c 20 62 63 20 61 6e 64 |programs|, bc and|
|00000530| 20 63 62 2e 0d 0a 2d 0d | 0a 2d 62 63 64 20 2d 68 | cb...-.|.-bcd -h|
|00000540| 20 2d 74 20 46 69 6c 74 | 65 72 5f 57 69 64 74 68 | -t Filt|er_Width|
|00000550| 20 4f 70 74 69 6f 6e 61 | 6c 5f 54 68 72 65 73 68 | Optiona|l_Thresh|
|00000560| 6f 6c 64 20 3c 20 49 6e | 70 75 74 5f 46 69 6c 65 |old < In|put_File|
|00000570| 20 3e 20 4f 75 74 70 75 | 74 5f 46 69 6c 65 0d 0a | > Outpu|t_File..|
|00000580| 2d 0d 0a 2d 62 63 64 20 | 74 61 6b 65 73 20 69 74 |-..-bcd |takes it|
|00000590| 73 20 49 6e 70 75 74 5f | 46 69 6c 65 20 66 72 6f |s Input_|File fro|
|000005a0| 6d 20 73 74 61 6e 64 61 | 72 64 20 69 6e 70 75 74 |m standa|rd input|
|000005b0| 2c 20 61 6e 64 20 77 72 | 69 74 65 73 20 74 6f 20 |, and wr|ites to |
|000005c0| 73 74 61 6e 64 61 72 64 | 0d 0a 2d 6f 75 74 2e 0d |standard|..-out..|
|000005d0| 0a 2d 54 68 65 20 6f 70 | 74 69 6f 6e 20 22 2d 68 |.-The op|tion "-h|
|000005e0| 22 20 44 49 53 41 42 4c | 45 53 20 68 79 73 74 65 |" DISABL|ES hyste|
|000005f0| 72 65 73 69 73 20 74 68 | 72 65 73 68 6f 6c 64 69 |resis th|resholdi|
|00000600| 6e 67 2e 0d 0a 2d 54 68 | 65 20 6f 70 74 69 6f 6e |ng...-Th|e option|
|00000610| 20 22 2d 74 22 20 44 49 | 53 41 42 4c 45 53 20 74 | "-t" DI|SABLES t|
|00000620| 68 69 6e 6e 69 6e 67 0d | 0a 2d 0d 0a 2d 42 6f 74 |hinning.|.-..-Bot|
|00000630| 68 20 74 68 65 20 49 6e | 70 75 74 5f 46 69 6c 65 |h the In|put_File|
|00000640| 20 61 6e 64 20 74 68 65 | 20 4f 75 74 70 75 74 5f | and the| Output_|
|00000650| 46 69 6c 65 20 61 72 65 | 20 69 6e 20 62 69 6e 61 |File are| in bina|
|00000660| 72 79 20 66 6f 72 6d 2e | 20 49 74 20 69 73 0d 0a |ry form.| It is..|
|00000670| 2d 61 73 73 75 6d 65 64 | 20 74 68 61 74 20 74 68 |-assumed| that th|
|00000680| 65 20 69 6d 61 67 65 20 | 69 73 20 73 74 6f 72 65 |e image |is store|
|00000690| 64 20 72 6f 77 20 61 74 | 20 61 20 74 69 6d 65 20 |d row at| a time |
|000006a0| 77 69 74 68 20 61 20 31 | 30 20 62 79 74 65 20 68 |with a 1|0 byte h|
|000006b0| 65 61 64 65 72 0d 0a 2d | 61 74 20 74 68 65 20 66 |eader..-|at the f|
|000006c0| 72 6f 6e 74 2e 20 54 68 | 65 20 68 65 61 64 65 72 |ront. Th|e header|
|000006d0| 20 63 6f 6e 73 69 73 74 | 73 20 6f 66 20 61 20 74 | consist|s of a t|
|000006e0| 77 6f 20 62 79 74 65 20 | 74 79 70 65 20 66 69 65 |wo byte |type fie|
|000006f0| 6c 64 2e 0d 0a 2d 46 6f | 6c 6c 6f 77 65 64 20 62 |ld...-Fo|llowed b|
|00000700| 79 20 6e 78 20 28 74 77 | 6f 20 62 79 74 65 73 29 |y nx (tw|o bytes)|
|00000710| 20 61 6e 64 20 6e 79 20 | 28 74 77 6f 20 62 79 74 | and ny |(two byt|
|00000720| 65 73 29 2c 0d 0a 2d 74 | 68 65 20 72 6f 77 20 61 |es),..-t|he row a|
|00000730| 6e 64 20 63 6f 6c 75 6d | 6e 20 64 69 6d 65 6e 73 |nd colum|n dimens|
|00000740| 69 6f 6e 73 2e 20 54 68 | 65 20 72 65 6d 61 69 6e |ions. Th|e remain|
|00000750| 69 6e 67 20 66 6f 75 72 | 20 62 79 74 65 73 20 61 |ing four| bytes a|
|00000760| 72 65 20 75 6e 75 73 65 | 64 2e 0d 0a 2d 2a 2a 2a |re unuse|d...-***|
|00000770| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000780| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000790| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000007a0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000007b0| 2a 0d 0a 2d 57 41 52 4e | 49 4e 47 3a 20 54 68 65 |*..-WARN|ING: The|
|000007c0| 20 66 69 6c 65 73 20 62 | 63 64 2e 63 20 61 6e 64 | files b|cd.c and|
|000007d0| 20 77 72 69 74 65 2e 63 | 20 73 68 6f 75 6c 64 20 | write.c| should |
|000007e0| 62 65 20 6d 6f 64 69 66 | 69 65 64 20 74 6f 20 66 |be modif|ied to f|
|000007f0| 6f 6c 6c 6f 77 0d 0a 2d | 74 68 65 20 69 6d 61 67 |ollow..-|the imag|
|00000800| 65 20 66 6f 72 6d 61 74 | 20 75 73 65 64 20 61 74 |e format| used at|
|00000810| 20 74 68 69 73 20 73 69 | 74 65 2e 0d 0a 2d 2a 2a | this si|te...-**|
|00000820| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000830| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000840| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000850| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000860| 2a 2a 0d 0a 2d 0d 0a 2d | 54 68 65 20 65 64 67 65 |**..-..-|The edge|
|00000870| 20 64 65 74 65 63 74 6f | 72 20 72 65 71 75 69 72 | detecto|r requir|
|00000880| 65 73 20 74 68 65 20 75 | 73 65 72 20 74 6f 20 70 |es the u|ser to p|
|00000890| 72 6f 76 69 64 65 20 62 | 6f 74 68 20 74 68 65 20 |rovide b|oth the |
|000008a0| 73 69 67 6d 61 20 6f 66 | 20 74 68 65 0d 0a 2d 66 |sigma of| the..-f|
|000008b0| 69 6c 74 65 72 20 28 6d | 65 61 73 75 72 65 64 20 |ilter (m|easured |
|000008c0| 69 6e 20 70 69 78 65 6c | 73 29 20 61 6e 64 20 61 |in pixel|s) and a|
|000008d0| 20 74 68 72 65 73 68 6f | 6c 64 2e 0d 0a 2d 46 6f | thresho|ld...-Fo|
|000008e0| 72 20 6d 6f 73 74 20 69 | 6d 61 67 65 73 2c 20 31 |r most i|mages, 1|
|000008f0| 2e 35 20 61 70 70 65 61 | 72 73 0d 0a 2d 74 6f 20 |.5 appea|rs..-to |
|00000900| 62 65 20 61 20 63 6f 6d | 6d 6f 6e 6c 79 20 75 73 |be a com|monly us|
|00000910| 65 64 20 76 61 6c 75 65 | 20 66 6f 72 20 73 69 67 |ed value| for sig|
|00000920| 6d 61 2e 0d 0a 2d 0d 0a | 2d 49 74 20 69 73 20 72 |ma...-..|-It is r|
|00000930| 65 63 6f 6d 6d 65 6e 64 | 65 64 20 74 68 61 74 20 |ecommend|ed that |
|00000940| 6e 6f 20 74 68 72 65 73 | 68 6f 6c 64 20 62 65 20 |no thres|hold be |
|00000950| 73 75 70 70 6c 69 65 64 | 2e 0d 0a 2d 48 6f 77 65 |supplied|...-Howe|
|00000960| 76 65 72 2c 20 74 68 65 | 20 61 75 74 6f 6d 61 74 |ver, the| automat|
|00000970| 69 63 20 74 68 72 65 73 | 68 6f 6c 64 69 6e 67 20 |ic thres|holding |
|00000980| 63 61 6e 20 62 65 20 6f | 76 65 72 72 69 64 65 64 |can be o|verrided|
|00000990| 20 62 65 20 73 70 65 63 | 69 66 79 69 6e 67 0d 0a | be spec|ifying..|
|000009a0| 2d 74 68 65 20 64 65 73 | 69 72 65 64 20 6e 75 6d |-the des|ired num|
|000009b0| 65 72 69 63 61 6c 20 74 | 68 72 65 73 68 6f 6c 64 |erical t|hreshold|
|000009c0| 2e 0d 0a 2d 0d 0a 2f 2f | 47 4f 2e 53 59 53 49 4e |...-..//|GO.SYSIN|
|000009d0| 20 44 44 20 52 45 41 44 | 4d 45 0d 0a 65 63 68 6f | DD READ|ME..echo|
|000009e0| 20 62 63 64 2e 63 20 31 | 3e 26 32 0d 0a 73 65 64 | bcd.c 1|>&2..sed|
|000009f0| 20 27 73 2f 2e 2f 2f 27 | 20 3e 62 63 64 2e 63 20 | 's/.//'| >bcd.c |
|00000a00| 3c 3c 27 2f 2f 47 4f 2e | 53 59 53 49 4e 20 44 44 |<<'//GO.|SYSIN DD|
|00000a10| 20 62 63 64 2e 63 27 0d | 0a 2d 2f 2a 0d 0a 2d 20 | bcd.c'.|.-/*..- |
|00000a20| 2a 20 28 63 29 20 43 6f | 70 79 72 69 67 68 74 20 |* (c) Co|pyright |
|00000a30| 31 39 38 38 20 62 79 0d | 0a 2d 20 2a 20 52 6f 62 |1988 by.|.- * Rob|
|00000a40| 6f 74 69 63 73 20 50 72 | 69 6e 63 69 70 6c 65 73 |otics Pr|inciples|
|00000a50| 20 52 65 73 65 61 72 63 | 68 20 44 65 70 61 72 74 | Researc|h Depart|
|00000a60| 6d 65 6e 74 2c 20 41 54 | 54 20 42 65 6c 6c 20 4c |ment, AT|T Bell L|
|00000a70| 61 62 6f 72 61 74 6f 72 | 69 65 73 2e 0d 0a 2d 20 |aborator|ies...- |
|00000a80| 2a 20 41 6c 6c 20 72 69 | 67 68 74 73 20 72 65 73 |* All ri|ghts res|
|00000a90| 65 72 76 65 64 2e 0d 0a | 2d 20 2a 20 4c 61 73 74 |erved...|- * Last|
|00000aa0| 20 6d 6f 64 69 66 69 65 | 64 20 32 2f 38 2f 38 38 | modifie|d 2/8/88|
|00000ab0| 20 49 6e 67 65 6d 61 72 | 20 4a 2e 20 43 6f 78 0d | Ingemar| J. Cox.|
|00000ac0| 0a 2d 20 2a 20 52 65 76 | 69 73 65 64 20 34 2f 32 |.- * Rev|ised 4/2|
|00000ad0| 32 2f 38 38 20 49 6e 67 | 65 6d 61 72 20 4a 2e 20 |2/88 Ing|emar J. |
|00000ae0| 43 6f 78 0d 0a 2d 20 2a | 20 52 65 76 69 73 65 64 |Cox..- *| Revised|
|00000af0| 20 31 31 2f 31 31 2f 38 | 38 20 49 6e 67 65 6d 61 | 11/11/8|8 Ingema|
|00000b00| 72 20 4a 2e 20 43 6f 78 | 0d 0a 2d 20 2a 2f 0d 0a |r J. Cox|..- */..|
|00000b10| 2d 23 69 6e 63 6c 75 64 | 65 20 3c 73 74 64 69 6f |-#includ|e <stdio|
|00000b20| 2e 68 3e 0d 0a 2d 23 69 | 6e 63 6c 75 64 65 20 22 |.h>..-#i|nclude "|
|00000b30| 65 64 67 65 5f 66 69 6e | 64 65 72 2e 68 22 0d 0a |edge_fin|der.h"..|
|00000b40| 2d 0d 0a 2d 23 64 65 66 | 69 6e 65 20 55 42 59 54 |-..-#def|ine UBYT|
|00000b50| 45 5f 48 44 20 30 78 34 | 35 30 30 0d 0a 2d 0d 0a |E_HD 0x4|500..-..|
|00000b60| 2d 6d 61 69 6e 28 61 72 | 67 63 2c 20 61 72 67 76 |-main(ar|gc, argv|
|00000b70| 29 0d 0a 2d 69 6e 74 20 | 61 72 67 63 3b 0d 0a 2d |)..-int |argc;..-|
|00000b80| 63 68 61 72 20 2a 61 72 | 67 76 5b 5d 3b 0d 0a 2d |char *ar|gv[];..-|
|00000b90| 7b 0d 0a 2d 09 73 68 6f | 72 74 20 68 65 61 64 65 |{..-.sho|rt heade|
|00000ba0| 72 5b 35 5d 3b 0d 0a 2d | 09 69 6e 74 20 6e 78 2c |r[5];..-|.int nx,|
|00000bb0| 20 6e 79 3b 0d 0a 2d 09 | 73 74 61 74 69 63 20 75 | ny;..-.|static u|
|00000bc0| 6e 73 69 67 6e 65 64 20 | 63 68 61 72 20 2a 70 69 |nsigned |char *pi|
|00000bd0| 63 3b 0d 0a 2d 09 66 6c | 6f 61 74 20 73 69 67 6d |c;..-.fl|oat sigm|
|00000be0| 61 3b 0d 0a 2d 09 69 6e | 74 20 74 68 72 65 73 68 |a;..-.in|t thresh|
|00000bf0| 6f 6c 64 3b 0d 0a 2d 09 | 69 6e 74 20 68 79 73 74 |old;..-.|int hyst|
|00000c00| 65 72 65 73 69 73 20 3d | 20 31 3b 0d 0a 2d 09 69 |eresis =| 1;..-.i|
|00000c10| 6e 74 20 74 68 69 6e 6e | 69 6e 67 20 3d 20 31 3b |nt thinn|ing = 1;|
|00000c20| 0d 0a 2d 09 69 6e 74 20 | 63 6e 74 2c 20 69 3b 0d |..-.int |cnt, i;.|
|00000c30| 0a 2d 09 0d 0a 2d 09 75 | 6e 73 69 67 6e 65 64 20 |.-...-.u|nsigned |
|00000c40| 63 68 61 72 20 2a 65 64 | 67 65 5f 6d 61 70 3b 0d |char *ed|ge_map;.|
|00000c50| 0a 2d 09 75 6e 73 69 67 | 6e 65 64 20 63 68 61 72 |.-.unsig|ned char|
|00000c60| 20 2a 69 6d 61 67 65 5f | 66 69 6e 64 5f 65 64 67 | *image_|find_edg|
|00000c70| 65 73 28 29 3b 0d 0a 2d | 09 75 6e 73 69 67 6e 65 |es();..-|.unsigne|
|00000c80| 64 20 63 68 61 72 20 2a | 69 6d 61 67 65 5f 68 79 |d char *|image_hy|
|00000c90| 73 74 65 72 65 73 69 73 | 28 29 3b 0d 0a 2d 09 75 |steresis|();..-.u|
|00000ca0| 6e 73 69 67 6e 65 64 20 | 63 68 61 72 20 2a 69 6d |nsigned |char *im|
|00000cb0| 61 67 65 5f 74 68 69 6e | 28 29 3b 0d 0a 2d 09 66 |age_thin|();..-.f|
|00000cc0| 6c 6f 61 74 20 61 74 6f | 66 28 29 3b 0d 0a 2d 09 |loat ato|f();..-.|
|00000cd0| 0d 0a 2d 09 69 66 28 61 | 72 67 63 3c 32 29 0d 0a |..-.if(a|rgc<2)..|
|00000ce0| 2d 09 7b 0d 0a 2d 09 09 | 66 70 72 69 6e 74 66 28 |-.{..-..|fprintf(|
|00000cf0| 73 74 64 65 72 72 2c 20 | 22 5c 6e 2a 2a 2a 65 72 |stderr, |"\n***er|
|00000d00| 72 6f 72 3a 20 77 72 6f | 6e 67 20 6e 75 6d 62 65 |ror: wro|ng numbe|
|00000d10| 72 20 6f 66 20 61 72 67 | 75 6d 65 6e 74 73 2a 2a |r of arg|uments**|
|00000d20| 2a 5c 6e 22 29 3b 0d 0a | 2d 09 09 66 70 72 69 6e |*\n");..|-..fprin|
|00000d30| 74 66 28 73 74 64 65 72 | 72 2c 20 0d 0a 2d 09 09 |tf(stder|r, ..-..|
|00000d40| 09 22 5c 6e 25 73 20 2d | 68 20 2d 74 20 73 69 67 |."\n%s -|h -t sig|
|00000d50| 6d 61 20 6f 70 74 5f 74 | 68 6c 64 20 3c 20 69 6e |ma opt_t|hld < in|
|00000d60| 70 75 74 20 3e 20 6f 75 | 74 70 75 74 5c 6e 22 2c |put > ou|tput\n",|
|00000d70| 20 61 72 67 76 5b 30 5d | 29 3b 0d 0a 2d 09 09 65 | argv[0]|);..-..e|
|00000d80| 78 69 74 28 31 29 3b 0d | 0a 2d 09 7d 0d 0a 2d 09 |xit(1);.|.-.}..-.|
|00000d90| 63 6e 74 20 3d 20 31 3b | 0d 0a 2d 09 77 68 69 6c |cnt = 1;|..-.whil|
|00000da0| 65 28 61 72 67 76 5b 63 | 6e 74 5d 5b 30 5d 3d 3d |e(argv[c|nt][0]==|
|00000db0| 27 2d 27 29 0d 0a 2d 09 | 7b 0d 0a 2d 09 09 73 77 |'-')..-.|{..-..sw|
|00000dc0| 69 74 63 68 28 61 72 67 | 76 5b 63 6e 74 5d 5b 31 |itch(arg|v[cnt][1|
|00000dd0| 5d 29 0d 0a 2d 09 09 7b | 0d 0a 2d 09 09 09 63 61 |])..-..{|..-...ca|
|00000de0| 73 65 20 27 68 27 3a 0d | 0a 2d 09 09 09 09 68 79 |se 'h':.|.-....hy|
|00000df0| 73 74 65 72 65 73 69 73 | 20 3d 20 30 3b 0d 0a 2d |steresis| = 0;..-|
|00000e00| 09 09 09 09 63 6e 74 2b | 2b 3b 0d 0a 2d 09 09 09 |....cnt+|+;..-...|
|00000e10| 09 62 72 65 61 6b 3b 0d | 0a 2d 09 09 09 63 61 73 |.break;.|.-...cas|
|00000e20| 65 20 27 74 27 3a 0d 0a | 2d 09 09 09 09 74 68 69 |e 't':..|-....thi|
|00000e30| 6e 6e 69 6e 67 20 3d 20 | 30 3b 0d 0a 2d 09 09 09 |nning = |0;..-...|
|00000e40| 09 63 6e 74 2b 2b 3b 0d | 0a 2d 09 09 09 09 62 72 |.cnt++;.|.-....br|
|00000e50| 65 61 6b 3b 0d 0a 2d 09 | 09 09 64 65 66 61 75 6c |eak;..-.|..defaul|
|00000e60| 74 3a 0d 0a 2d 09 09 09 | 09 66 70 72 69 6e 74 66 |t:..-...|.fprintf|
|00000e70| 28 73 74 64 65 72 72 2c | 20 22 25 73 3a 20 73 79 |(stderr,| "%s: sy|
|00000e80| 6e 74 61 78 20 65 72 72 | 6f 72 5c 6e 22 2c 20 61 |ntax err|or\n", a|
|00000e90| 72 67 76 5b 30 5d 29 3b | 0d 0a 2d 09 09 09 09 66 |rgv[0]);|..-....f|
|00000ea0| 70 72 69 6e 74 66 28 73 | 74 64 65 72 72 2c 20 0d |printf(s|tderr, .|
|00000eb0| 0a 2d 09 09 09 09 09 22 | 25 73 20 2d 68 20 2d 74 |.-....."|%s -h -t|
|00000ec0| 20 73 69 67 6d 61 20 6f | 70 74 5f 74 68 6c 64 20 | sigma o|pt_thld |
|00000ed0| 3c 20 69 6e 70 75 74 20 | 3e 20 6f 75 74 70 75 74 |< input |> output|
|00000ee0| 5c 6e 22 2c 20 61 72 67 | 76 5b 30 5d 29 3b 0d 0a |\n", arg|v[0]);..|
|00000ef0| 2d 09 09 09 09 65 78 69 | 74 28 31 29 3b 0d 0a 2d |-....exi|t(1);..-|
|00000f00| 09 09 7d 0d 0a 2d 09 7d | 0d 0a 2d 09 73 69 67 6d |..}..-.}|..-.sigm|
|00000f10| 61 3d 61 74 6f 66 28 61 | 72 67 76 5b 63 6e 74 2b |a=atof(a|rgv[cnt+|
|00000f20| 2b 5d 29 3b 0d 0a 2d 09 | 69 66 28 61 72 67 63 3e |+]);..-.|if(argc>|
|00000f30| 63 6e 74 29 0d 0a 2d 09 | 09 74 68 72 65 73 68 6f |cnt)..-.|.thresho|
|00000f40| 6c 64 3d 61 74 6f 69 28 | 61 72 67 76 5b 63 6e 74 |ld=atoi(|argv[cnt|
|00000f50| 5d 29 3b 0d 0a 2d 09 65 | 6c 73 65 0d 0a 2d 09 09 |]);..-.e|lse..-..|
|00000f60| 74 68 72 65 73 68 6f 6c | 64 3d 30 3b 0d 0a 2d 2f |threshol|d=0;..-/|
|00000f70| 2a 0d 0a 2d 20 2a 20 6e | 6f 77 20 72 65 61 64 20 |*..- * n|ow read |
|00000f80| 69 6e 20 69 6d 61 67 65 | 0d 0a 2d 20 2a 2f 0d 0a |in image|..- */..|
|00000f90| 2d 09 72 65 61 64 28 30 | 2c 20 28 63 68 61 72 2a |-.read(0|, (char*|
|00000fa0| 29 68 65 61 64 65 72 2c | 20 31 30 29 3b 0d 0a 2d |)header,| 10);..-|
|00000fb0| 09 69 66 28 68 65 61 64 | 65 72 5b 30 5d 20 21 3d |.if(head|er[0] !=|
|00000fc0| 20 55 42 59 54 45 5f 48 | 44 29 0d 0a 2d 09 7b 0d | UBYTE_H|D)..-.{.|
|00000fd0| 0a 2d 09 09 66 70 72 69 | 6e 74 66 28 73 74 64 65 |.-..fpri|ntf(stde|
|00000fe0| 72 72 2c 20 22 25 73 3a | 20 45 52 52 4f 52 20 69 |rr, "%s:| ERROR i|
|00000ff0| 6d 61 67 65 20 6e 6f 74 | 20 6f 66 20 75 6e 73 69 |mage not| of unsi|
|00001000| 67 6e 65 64 20 63 68 61 | 72 5c 6e 22 2c 20 61 72 |gned cha|r\n", ar|
|00001010| 67 76 5b 30 5d 29 3b 0d | 0a 2d 09 09 65 78 69 74 |gv[0]);.|.-..exit|
|00001020| 28 31 29 3b 0d 0a 2d 09 | 7d 0d 0a 2d 09 6e 78 20 |(1);..-.|}..-.nx |
|00001030| 3d 20 68 65 61 64 65 72 | 5b 31 5d 3b 0d 0a 2d 09 |= header|[1];..-.|
|00001040| 6e 79 20 3d 20 68 65 61 | 64 65 72 5b 32 5d 3b 0d |ny = hea|der[2];.|
|00001050| 0a 2d 09 66 70 72 69 6e | 74 66 28 73 74 64 65 72 |.-.fprin|tf(stder|
|00001060| 72 2c 20 22 5c 6e 69 6d | 61 67 65 20 73 69 7a 65 |r, "\nim|age size|
|00001070| 20 69 73 20 25 64 20 78 | 20 25 64 5c 6e 22 2c 20 | is %d x| %d\n", |
|00001080| 6e 78 2c 20 6e 79 29 3b | 0d 0a 2d 09 0d 0a 2d 09 |nx, ny);|..-...-.|
|00001090| 69 66 28 28 70 69 63 20 | 3d 20 28 75 6e 73 69 67 |if((pic |= (unsig|
|000010a0| 6e 65 64 20 63 68 61 72 | 20 2a 29 6d 61 6c 6c 6f |ned char| *)mallo|
|000010b0| 63 28 6e 78 2a 6e 79 2a | 73 69 7a 65 6f 66 28 75 |c(nx*ny*|sizeof(u|
|000010c0| 6e 73 69 67 6e 65 64 20 | 63 68 61 72 29 29 29 3d |nsigned |char)))=|
|000010d0| 3d 30 29 0d 0a 2d 09 7b | 0d 0a 2d 09 09 66 70 72 |=0)..-.{|..-..fpr|
|000010e0| 69 6e 74 66 28 73 74 64 | 65 72 72 2c 20 0d 0a 2d |intf(std|err, ..-|
|000010f0| 09 09 09 22 25 73 3a 20 | 45 52 52 4f 52 20 63 61 |..."%s: |ERROR ca|
|00001100| 6e 6e 6f 74 20 61 6c 6c | 6f 63 61 74 65 20 69 6d |nnot all|ocate im|
|00001110| 61 67 65 20 61 72 72 61 | 79 5c 6e 22 2c 20 61 72 |age arra|y\n", ar|
|00001120| 67 76 5b 30 5d 29 3b 0d | 0a 2d 09 09 65 78 69 74 |gv[0]);.|.-..exit|
|00001130| 28 31 29 3b 0d 0a 2d 09 | 7d 0d 0a 2d 0d 0a 2d 09 |(1);..-.|}..-..-.|
|00001140| 66 6f 72 28 69 3d 30 3b | 20 69 3c 6e 79 3b 20 69 |for(i=0;| i<ny; i|
|00001150| 2b 2b 29 0d 0a 2d 09 7b | 0d 0a 2d 09 09 72 65 61 |++)..-.{|..-..rea|
|00001160| 64 28 30 2c 20 28 63 68 | 61 72 2a 29 70 69 63 2b |d(0, (ch|ar*)pic+|
|00001170| 69 2a 6e 78 2c 20 6e 78 | 29 3b 0d 0a 2d 09 7d 0d |i*nx, nx|);..-.}.|
|00001180| 0a 2d 09 69 6d 61 67 65 | 5f 69 6e 69 74 28 70 69 |.-.image|_init(pi|
|00001190| 63 2c 20 6e 78 2c 20 6e | 79 29 3b 0d 0a 2d 0d 0a |c, nx, n|y);..-..|
|000011a0| 2d 09 65 64 67 65 5f 6d | 61 70 20 3d 20 69 6d 61 |-.edge_m|ap = ima|
|000011b0| 67 65 5f 66 69 6e 64 5f | 65 64 67 65 73 28 73 69 |ge_find_|edges(si|
|000011c0| 67 6d 61 2c 20 26 74 68 | 72 65 73 68 6f 6c 64 29 |gma, &th|reshold)|
|000011d0| 3b 0d 0a 2d 09 69 66 28 | 68 79 73 74 65 72 65 73 |;..-.if(|hysteres|
|000011e0| 69 73 29 0d 0a 2d 09 7b | 0d 0a 2d 09 09 66 70 72 |is)..-.{|..-..fpr|
|000011f0| 69 6e 74 66 28 73 74 64 | 65 72 72 2c 20 22 5c 6e |intf(std|err, "\n|
|00001200| 61 70 70 6c 79 69 6e 67 | 20 68 79 73 74 65 72 65 |applying| hystere|
|00001210| 73 69 73 5c 6e 22 29 3b | 0d 0a 2d 09 09 65 64 67 |sis\n");|..-..edg|
|00001220| 65 5f 6d 61 70 20 3d 20 | 69 6d 61 67 65 5f 68 79 |e_map = |image_hy|
|00001230| 73 74 65 72 65 73 69 73 | 28 74 68 72 65 73 68 6f |steresis|(thresho|
|00001240| 6c 64 2f 33 29 3b 0d 0a | 2d 09 7d 0d 0a 2d 09 69 |ld/3);..|-.}..-.i|
|00001250| 66 28 74 68 69 6e 6e 69 | 6e 67 29 0d 0a 2d 09 7b |f(thinni|ng)..-.{|
|00001260| 0d 0a 2d 09 09 66 70 72 | 69 6e 74 66 28 73 74 64 |..-..fpr|intf(std|
|00001270| 65 72 72 2c 20 22 5c 6e | 74 68 69 6e 6e 69 6e 67 |err, "\n|thinning|
|00001280| 20 69 6d 61 67 65 5c 6e | 22 29 3b 0d 0a 2d 09 09 | image\n|");..-..|
|00001290| 65 64 67 65 5f 6d 61 70 | 20 3d 20 69 6d 61 67 65 |edge_map| = image|
|000012a0| 5f 74 68 69 6e 28 29 3b | 0d 0a 2d 09 7d 0d 0a 2d |_thin();|..-.}..-|
|000012b0| 09 77 72 69 74 65 28 31 | 2c 20 28 63 68 61 72 2a |.write(1|, (char*|
|000012c0| 29 68 65 61 64 65 72 2c | 20 31 30 29 3b 0d 0a 2d |)header,| 10);..-|
|000012d0| 09 77 72 69 74 65 28 31 | 2c 20 28 63 68 61 72 2a |.write(1|, (char*|
|000012e0| 29 65 64 67 65 5f 6d 61 | 70 2c 20 6e 78 2a 6e 79 |)edge_ma|p, nx*ny|
|000012f0| 29 3b 0d 0a 2d 0d 0a 2d | 09 69 6d 61 67 65 5f 63 |);..-..-|.image_c|
|00001300| 6c 65 61 6e 75 70 28 29 | 3b 0d 0a 2d 7d 0d 0a 2d |leanup()|;..-}..-|
|00001310| 0d 0a 2f 2f 47 4f 2e 53 | 59 53 49 4e 20 44 44 20 |..//GO.S|YSIN DD |
|00001320| 62 63 64 2e 63 0d 0a 65 | 63 68 6f 20 63 6f 6e 76 |bcd.c..e|cho conv|
|00001330| 2e 63 20 31 3e 26 32 0d | 0a 73 65 64 20 27 73 2f |.c 1>&2.|.sed 's/|
|00001340| 2e 2f 2f 27 20 3e 63 6f | 6e 76 2e 63 20 3c 3c 27 |.//' >co|nv.c <<'|
|00001350| 2f 2f 47 4f 2e 53 59 53 | 49 4e 20 44 44 20 63 6f |//GO.SYS|IN DD co|
|00001360| 6e 76 2e 63 27 0d 0a 2d | 2f 2a 0d 0a 2d 20 2a 20 |nv.c'..-|/*..- * |
|00001370| 28 63 29 20 43 6f 70 79 | 72 69 67 68 74 20 31 39 |(c) Copy|right 19|
|00001380| 38 38 20 62 79 0d 0a 2d | 20 2a 20 52 6f 62 6f 74 |88 by..-| * Robot|
|00001390| 69 63 73 20 50 72 69 6e | 63 69 70 6c 65 73 20 52 |ics Prin|ciples R|
|000013a0| 65 73 65 61 72 63 68 20 | 44 65 70 61 72 74 6d 65 |esearch |Departme|
|000013b0| 6e 74 2c 20 41 54 54 20 | 42 65 6c 6c 20 4c 61 62 |nt, ATT |Bell Lab|
|000013c0| 6f 72 61 74 6f 72 69 65 | 73 2e 0d 0a 2d 20 2a 20 |oratorie|s...- * |
|000013d0| 41 6c 6c 20 72 69 67 68 | 74 73 20 72 65 73 65 72 |All righ|ts reser|
|000013e0| 76 65 64 2e 0d 0a 2d 20 | 2a 20 4c 61 73 74 20 6d |ved...- |* Last m|
|000013f0| 6f 64 69 66 69 65 64 20 | 32 2f 38 2f 38 38 20 49 |odified |2/8/88 I|
|00001400| 6e 67 65 6d 61 72 20 4a | 2e 20 43 6f 78 0d 0a 2d |ngemar J|. Cox..-|
|00001410| 20 2a 20 43 20 76 65 72 | 73 69 6f 6e 20 38 2f 32 | * C ver|sion 8/2|
|00001420| 2f 38 38 20 44 65 62 6f | 72 61 68 20 41 2e 20 57 |/88 Debo|rah A. W|
|00001430| 61 6c 6c 61 63 68 0d 0a | 2d 20 2a 20 43 20 76 65 |allach..|- * C ve|
|00001440| 72 73 69 6f 6e 20 31 31 | 2f 34 2f 38 38 20 20 57 |rsion 11|/4/88  W|
|00001450| 2e 4a 20 4b 72 6f 70 66 | 6c 0d 0a 2d 20 2a 2f 0d |.J Kropf|l..- */.|
|00001460| 0a 2d 23 69 6e 63 6c 75 | 64 65 20 3c 73 74 64 69 |.-#inclu|de <stdi|
|00001470| 6f 2e 68 3e 0d 0a 2d 23 | 69 6e 63 6c 75 64 65 20 |o.h>..-#|include |
|00001480| 3c 6d 61 74 68 2e 68 3e | 0d 0a 2d 23 69 6e 63 6c |<math.h>|..-#incl|
|00001490| 75 64 65 20 22 65 64 67 | 65 5f 66 69 6e 64 65 72 |ude "edg|e_finder|
|000014a0| 2e 68 22 0d 0a 2d 0d 0a | 2d 65 78 74 65 72 6e 20 |.h"..-..|-extern |
|000014b0| 73 74 72 75 63 74 20 69 | 6d 61 67 65 20 2a 6d 79 |struct i|mage *my|
|000014c0| 5f 69 6d 61 67 65 3b 0d | 0a 2d 09 73 74 61 74 69 |_image;.|.-.stati|
|000014d0| 63 20 69 6e 74 20 69 6b | 6c 69 6d 2c 20 6e 66 2c |c int ik|lim, nf,|
|000014e0| 20 6e 66 32 3b 0d 0a 2d | 09 73 74 61 74 69 63 20 | nf2;..-|.static |
|000014f0| 69 6e 74 20 69 79 2c 20 | 6e 79 3b 0d 0a 2d 09 73 |int iy, |ny;..-.s|
|00001500| 74 61 74 69 63 20 69 6e | 74 20 6e 78 58 32 4d 31 |tatic in|t nxX2M1|
|00001510| 2c 20 69 6b 58 6e 78 2c | 20 6e 79 58 32 4d 31 58 |, ikXnx,| nyX2M1X|
|00001520| 6e 78 2c 20 6e 66 58 6e | 78 3b 0d 0a 2d 09 73 74 |nx, nfXn|x;..-.st|
|00001530| 61 74 69 63 20 69 6e 74 | 20 69 79 4d 6e 66 32 58 |atic int| iyMnf2X|
|00001540| 6e 78 2c 20 69 79 30 2c | 20 69 79 30 58 6e 78 2c |nx, iy0,| iy0Xnx,|
|00001550| 20 69 74 3b 0d 0a 2d 09 | 73 74 61 74 69 63 20 69 | it;..-.|static i|
|00001560| 6e 74 20 2a 74 65 6d 70 | 5f 69 6d 61 67 65 2c 20 |nt *temp|_image, |
|00001570| 2a 66 69 6c 74 65 72 3b | 0d 0a 2d 09 73 74 61 74 |*filter;|..-.stat|
|00001580| 69 63 20 75 6e 73 69 67 | 6e 65 64 20 63 68 61 72 |ic unsig|ned char|
|00001590| 20 2a 74 70 69 63 50 3b | 0d 0a 2d 0d 0a 2d 76 6f | *tpicP;|..-..-vo|
|000015a0| 69 64 0d 0a 2d 69 6d 61 | 67 65 5f 63 6f 6e 76 6f |id..-ima|ge_convo|
|000015b0| 6c 76 65 28 70 69 63 29 | 0d 0a 2d 20 20 20 20 20 |lve(pic)|..-     |
|000015c0| 75 6e 73 69 67 6e 65 64 | 20 63 68 61 72 20 2a 70 |unsigned| char *p|
|000015d0| 69 63 3b 0d 0a 2d 7b 0d | 0a 2d 09 72 65 67 69 73 |ic;..-{.|.-.regis|
|000015e0| 74 65 72 20 69 6e 74 20 | 72 65 73 75 6c 74 2c 20 |ter int |result, |
|000015f0| 69 6b 2c 20 69 78 2c 20 | 6e 78 3b 0d 0a 2d 09 72 |ik, ix, |nx;..-.r|
|00001600| 65 67 69 73 74 65 72 20 | 75 6e 73 69 67 6e 65 64 |egister |unsigned|
|00001610| 20 63 68 61 72 20 2a 70 | 69 63 50 3b 0d 0a 2d 09 | char *p|icP;..-.|
|00001620| 72 65 67 69 73 74 65 72 | 20 69 6e 74 20 2a 66 69 |register| int *fi|
|00001630| 6c 74 50 2c 20 2a 74 6d | 70 50 2c 20 2a 67 61 75 |ltP, *tm|pP, *gau|
|00001640| 73 73 69 61 6e 50 3b 0d | 0a 2d 0d 0a 2d 09 6e 78 |ssianP;.|.-..-.nx|
|00001650| 20 20 3d 20 6d 79 5f 69 | 6d 61 67 65 2d 3e 6e 78 |  = my_i|mage->nx|
|00001660| 3b 0d 0a 2d 09 6e 79 20 | 20 3d 20 6d 79 5f 69 6d |;..-.ny | = my_im|
|00001670| 61 67 65 2d 3e 6e 79 3b | 0d 0a 2d 09 6e 66 20 20 |age->ny;|..-.nf  |
|00001680| 3d 20 6d 79 5f 69 6d 61 | 67 65 2d 3e 6e 66 3b 09 |= my_ima|ge->nf;.|
|00001690| 6e 66 58 6e 78 20 3d 20 | 6e 66 2a 6e 78 3b 0d 0a |nfXnx = |nf*nx;..|
|000016a0| 2d 09 6e 66 32 20 3d 20 | 6d 79 5f 69 6d 61 67 65 |-.nf2 = |my_image|
|000016b0| 2d 3e 6e 66 32 3b 0d 0a | 2d 09 6e 78 58 32 4d 31 |->nf2;..|-.nxX2M1|
|000016c0| 20 3d 20 6e 78 20 2b 20 | 6e 78 20 2d 20 31 3b 0d | = nx + |nx - 1;.|
|000016d0| 0a 2d 09 6e 79 58 32 4d | 31 58 6e 78 20 3d 20 28 |.-.nyX2M|1Xnx = (|
|000016e0| 6e 79 20 2b 20 6e 79 20 | 2d 20 31 29 2a 6e 78 3b |ny + ny |- 1)*nx;|
|000016f0| 0d 0a 2d 09 66 69 6c 74 | 65 72 20 3d 20 6d 79 5f |..-.filt|er = my_|
|00001700| 69 6d 61 67 65 2d 3e 66 | 69 6c 74 65 72 3b 0d 0a |image->f|ilter;..|
|00001710| 2d 0d 0a 2d 2f 2a 0d 0a | 2d 20 2a 09 63 72 65 61 |-..-/*..|- *.crea|
|00001720| 74 65 20 74 65 6d 70 20 | 61 6e 64 20 6f 75 74 70 |te temp |and outp|
|00001730| 75 74 20 61 72 72 61 79 | 73 0d 0a 2d 20 2a 2f 0d |ut array|s..- */.|
|00001740| 0a 2d 0d 0a 2d 09 69 66 | 28 28 74 65 6d 70 5f 69 |.-..-.if|((temp_i|
|00001750| 6d 61 67 65 20 3d 20 28 | 69 6e 74 20 2a 29 6d 61 |mage = (|int *)ma|
|00001760| 6c 6c 6f 63 28 6e 79 2a | 6e 78 2a 73 69 7a 65 6f |lloc(ny*|nx*sizeo|
|00001770| 66 28 69 6e 74 29 29 29 | 3d 3d 30 29 0d 0a 2d 09 |f(int)))|==0)..-.|
|00001780| 7b 0d 0a 2d 09 09 66 70 | 72 69 6e 74 66 28 73 74 |{..-..fp|rintf(st|
|00001790| 64 65 72 72 2c 20 22 65 | 72 72 6f 72 3a 20 63 61 |derr, "e|rror: ca|
|000017a0| 6e 6e 6f 74 20 61 6c 6c | 6f 63 61 74 65 20 74 65 |nnot all|ocate te|
|000017b0| 6d 70 20 61 72 72 61 79 | 5c 6e 22 29 3b 0d 0a 2d |mp array|\n");..-|
|000017c0| 09 09 65 78 69 74 28 31 | 29 3b 0d 0a 2d 09 7d 0d |..exit(1|);..-.}.|
|000017d0| 0a 2d 0d 0a 2d 09 69 66 | 28 28 6d 79 5f 69 6d 61 |.-..-.if|((my_ima|
|000017e0| 67 65 2d 3e 67 61 75 73 | 73 69 61 6e 20 3d 20 28 |ge->gaus|sian = (|
|000017f0| 69 6e 74 20 2a 29 6d 61 | 6c 6c 6f 63 28 6e 79 2a |int *)ma|lloc(ny*|
|00001800| 6e 78 2a 73 69 7a 65 6f | 66 28 69 6e 74 29 29 29 |nx*sizeo|f(int)))|
|00001810| 3d 3d 30 29 0d 0a 2d 09 | 7b 0d 0a 2d 09 09 66 70 |==0)..-.|{..-..fp|
|00001820| 72 69 6e 74 66 28 73 74 | 64 65 72 72 2c 20 22 65 |rintf(st|derr, "e|
|00001830| 72 72 6f 72 3a 20 63 61 | 6e 6e 6f 74 20 61 6c 6c |rror: ca|nnot all|
|00001840| 6f 63 61 74 65 20 67 61 | 75 73 73 69 61 6e 20 61 |ocate ga|ussian a|
|00001850| 72 72 61 79 5c 6e 22 29 | 3b 0d 0a 2d 09 09 65 78 |rray\n")|;..-..ex|
|00001860| 69 74 28 31 29 3b 0d 0a | 2d 09 7d 0d 0a 2d 2f 2a |it(1);..|-.}..-/*|
|00001870| 0d 0a 2d 20 2a 09 66 69 | 72 73 74 20 63 6f 6e 76 |..- *.fi|rst conv|
|00001880| 6f 6c 76 65 20 69 6e 20 | 78 2d 64 69 72 6e 0d 0a |olve in |x-dirn..|
|00001890| 2d 20 2a 09 09 64 6f 20 | 6c 65 66 74 20 73 69 64 |- *..do |left sid|
|000018a0| 65 20 66 69 72 73 74 0d | 0a 2d 20 2a 09 09 74 68 |e first.|.- *..th|
|000018b0| 65 6e 20 6d 69 64 64 6c | 65 0d 0a 2d 20 2a 09 09 |en middl|e..- *..|
|000018c0| 74 68 65 6e 20 72 69 67 | 68 74 20 73 69 64 65 0d |then rig|ht side.|
|000018d0| 0a 2d 20 2a 2f 0d 0a 2d | 09 0d 0a 2d 0d 0a 2d 09 |.- */..-|...-..-.|
|000018e0| 74 6d 70 50 20 3d 20 26 | 74 65 6d 70 5f 69 6d 61 |tmpP = &|temp_ima|
|000018f0| 67 65 5b 30 5d 3b 0d 0a | 2d 09 70 69 63 50 20 3d |ge[0];..|-.picP =|
|00001900| 20 70 69 63 3b 0d 0a 2d | 09 66 6f 72 28 69 79 3d | pic;..-|.for(iy=|
|00001910| 30 3b 20 69 79 3c 6e 79 | 3b 20 69 79 2b 2b 2c 20 |0; iy<ny|; iy++, |
|00001920| 74 6d 70 50 2b 3d 6e 78 | 2c 20 70 69 63 50 2b 3d |tmpP+=nx|, picP+=|
|00001930| 6e 78 29 0d 0a 2d 09 7b | 0d 0a 2d 09 09 66 6f 72 |nx)..-.{|..-..for|
|00001940| 28 69 78 3d 30 3b 20 69 | 78 3c 6e 66 32 3b 20 69 |(ix=0; i|x<nf2; i|
|00001950| 78 2b 2b 29 0d 0a 2d 09 | 09 7b 0d 0a 2d 09 09 09 |x++)..-.|.{..-...|
|00001960| 66 69 6c 74 50 20 3d 20 | 66 69 6c 74 65 72 3b 0d |filtP = |filter;.|
|00001970| 0a 2d 09 09 09 72 65 73 | 75 6c 74 20 3d 20 30 3b |.-...res|ult = 0;|
|00001980| 0d 0a 2d 09 09 09 69 6b | 20 3d 20 69 78 20 2d 20 |..-...ik| = ix - |
|00001990| 6e 66 32 3b 0d 0a 2d 09 | 09 09 69 6b 6c 69 6d 20 |nf2;..-.|..iklim |
|000019a0| 3d 20 69 6b 20 2b 20 6e | 66 3b 0d 0a 2d 09 09 09 |= ik + n|f;..-...|
|000019b0| 66 6f 72 28 3b 20 69 6b | 3c 69 6b 6c 69 6d 3b 20 |for(; ik|<iklim; |
|000019c0| 69 6b 2b 2b 29 0d 0a 2d | 09 09 09 09 72 65 73 75 |ik++)..-|....resu|
|000019d0| 6c 74 20 2b 3d 20 28 2a | 66 69 6c 74 50 2b 2b 29 |lt += (*|filtP++)|
|000019e0| 20 2a 20 70 69 63 50 5b | 61 62 73 28 69 6b 29 5d | * picP[|abs(ik)]|
|000019f0| 3b 0d 0a 2d 0d 0a 2d 09 | 09 09 2a 28 74 6d 70 50 |;..-..-.|..*(tmpP|
|00001a00| 2b 69 78 29 20 3d 20 72 | 65 73 75 6c 74 3b 0d 0a |+ix) = r|esult;..|
|00001a10| 2d 09 09 7d 0d 0a 2d 09 | 7d 0d 0a 2d 09 0d 0a 2d |-..}..-.|}..-...-|
|00001a20| 09 74 6d 70 50 20 3d 20 | 26 74 65 6d 70 5f 69 6d |.tmpP = |&temp_im|
|00001a30| 61 67 65 5b 30 5d 3b 0d | 0a 2d 09 74 70 69 63 50 |age[0];.|.-.tpicP|
|00001a40| 20 3d 20 70 69 63 3b 0d | 0a 2d 09 66 6f 72 28 69 | = pic;.|.-.for(i|
|00001a50| 79 3d 30 3b 20 69 79 3c | 6e 79 3b 20 69 79 2b 2b |y=0; iy<|ny; iy++|
|00001a60| 2c 20 74 6d 70 50 2b 3d | 6e 78 2c 20 74 70 69 63 |, tmpP+=|nx, tpic|
|00001a70| 50 2b 3d 6e 78 29 0d 0a | 2d 09 7b 0d 0a 2d 09 09 |P+=nx)..|-.{..-..|
|00001a80| 69 6e 74 20 78 6c 69 6d | 20 3d 20 6e 78 2d 6e 66 |int xlim| = nx-nf|
|00001a90| 32 2d 31 3b 0d 0a 2d 09 | 09 66 6f 72 28 69 78 3d |2-1;..-.|.for(ix=|
|00001aa0| 6e 66 32 3b 20 69 78 3c | 78 6c 69 6d 3b 20 69 78 |nf2; ix<|xlim; ix|
|00001ab0| 2b 2b 29 0d 0a 2d 09 09 | 7b 0d 0a 2d 09 09 09 70 |++)..-..|{..-...p|
|00001ac0| 69 63 50 20 3d 20 74 70 | 69 63 50 2b 69 78 2d 6e |icP = tp|icP+ix-n|
|00001ad0| 66 32 3b 0d 0a 2d 09 09 | 09 66 69 6c 74 50 20 3d |f2;..-..|.filtP =|
|00001ae0| 20 66 69 6c 74 65 72 3b | 0d 0a 2d 09 09 09 72 65 | filter;|..-...re|
|00001af0| 73 75 6c 74 20 3d 20 30 | 3b 0d 0a 2d 09 09 09 69 |sult = 0|;..-...i|
|00001b00| 6b 3d 6e 66 3b 0d 0a 2d | 09 09 09 77 68 69 6c 65 |k=nf;..-|...while|
|00001b10| 28 69 6b 2d 2d 29 0d 0a | 2d 09 09 09 09 72 65 73 |(ik--)..|-....res|
|00001b20| 75 6c 74 20 2b 3d 20 28 | 2a 66 69 6c 74 50 2b 2b |ult += (|*filtP++|
|00001b30| 29 20 2a 20 28 2a 70 69 | 63 50 2b 2b 29 3b 0d 0a |) * (*pi|cP++);..|
|00001b40| 2d 0d 0a 2d 09 09 09 2a | 28 74 6d 70 50 2b 69 78 |-..-...*|(tmpP+ix|
|00001b50| 29 20 3d 20 72 65 73 75 | 6c 74 3b 0d 0a 2d 09 09 |) = resu|lt;..-..|
|00001b60| 7d 0d 0a 2d 09 7d 0d 0a | 2d 09 0d 0a 2d 09 74 6d |}..-.}..|-...-.tm|
|00001b70| 70 50 20 3d 20 26 74 65 | 6d 70 5f 69 6d 61 67 65 |pP = &te|mp_image|
|00001b80| 5b 30 5d 3b 0d 0a 2d 09 | 70 69 63 50 20 3d 20 70 |[0];..-.|picP = p|
|00001b90| 69 63 3b 0d 0a 2d 09 66 | 6f 72 28 69 79 3d 30 3b |ic;..-.f|or(iy=0;|
|00001ba0| 20 69 79 3c 6e 79 3b 20 | 69 79 2b 2b 2c 20 74 6d | iy<ny; |iy++, tm|
|00001bb0| 70 50 2b 3d 6e 78 2c 20 | 70 69 63 50 2b 3d 6e 78 |pP+=nx, |picP+=nx|
|00001bc0| 29 0d 0a 2d 09 7b 0d 0a | 2d 09 09 66 6f 72 28 69 |)..-.{..|-..for(i|
|00001bd0| 78 3d 6e 78 2d 6e 66 32 | 2d 31 3b 20 69 78 3c 6e |x=nx-nf2|-1; ix<n|
|00001be0| 78 3b 20 69 78 2b 2b 29 | 0d 0a 2d 09 09 7b 0d 0a |x; ix++)|..-..{..|
|00001bf0| 2d 09 09 09 66 69 6c 74 | 50 20 3d 20 66 69 6c 74 |-...filt|P = filt|
|00001c00| 65 72 3b 0d 0a 2d 09 09 | 09 72 65 73 75 6c 74 20 |er;..-..|.result |
|00001c10| 3d 20 30 3b 0d 0a 2d 09 | 09 09 69 6b 20 3d 20 69 |= 0;..-.|..ik = i|
|00001c20| 78 20 2d 20 6e 66 32 3b | 0d 0a 2d 09 09 09 69 6b |x - nf2;|..-...ik|
|00001c30| 6c 69 6d 20 3d 20 69 6b | 20 2b 20 6e 66 3b 0d 0a |lim = ik| + nf;..|
|00001c40| 2d 09 09 09 66 6f 72 28 | 3b 20 69 6b 3c 69 6b 6c |-...for(|; ik<ikl|
|00001c50| 69 6d 3b 20 69 6b 2b 2b | 29 0d 0a 2d 09 09 09 7b |im; ik++|)..-...{|
|00001c60| 0d 0a 2d 09 2f 2a 09 09 | 09 72 65 67 69 73 74 65 |..-./*..|.registe|
|00001c70| 72 20 69 6e 74 20 69 6e | 64 65 78 20 3d 20 69 78 |r int in|dex = ix|
|00001c80| 2d 6e 66 32 2b 69 6b 3b | 0d 0a 2d 09 20 2a 09 09 |-nf2+ik;|..-. *..|
|00001c90| 09 69 66 28 69 6e 64 65 | 78 3e 3d 6e 78 29 20 69 |.if(inde|x>=nx) i|
|00001ca0| 6e 64 65 78 20 3d 20 6e | 78 2b 6e 78 2d 69 6e 64 |ndex = n|x+nx-ind|
|00001cb0| 65 78 2d 31 3b 0d 0a 2d | 09 20 2a 09 09 09 72 65 |ex-1;..-|. *...re|
|00001cc0| 73 75 6c 74 20 2b 3d 20 | 28 2a 66 69 6c 74 50 2b |sult += |(*filtP+|
|00001cd0| 2b 29 20 2a 20 70 69 63 | 50 5b 69 6e 64 65 78 5d |+) * pic|P[index]|
|00001ce0| 3b 0d 0a 2d 09 20 2a 2f | 0d 0a 2d 09 09 09 09 69 |;..-. */|..-....i|
|00001cf0| 66 28 69 6b 3e 3d 6e 78 | 29 09 20 20 20 20 20 2f |f(ik>=nx|).     /|
|00001d00| 2a 20 6e 78 58 32 4d 31 | 20 3d 20 6e 78 2b 6e 78 |* nxX2M1| = nx+nx|
|00001d10| 2d 69 6b 2d 31 3b 20 2a | 2f 0d 0a 2d 09 09 09 09 |-ik-1; *|/..-....|
|00001d20| 09 72 65 73 75 6c 74 20 | 2b 3d 20 2a 66 69 6c 74 |.result |+= *filt|
|00001d30| 50 20 2a 20 70 69 63 50 | 5b 6e 78 58 32 4d 31 2d |P * picP|[nxX2M1-|
|00001d40| 69 6b 5d 3b 0d 0a 2d 09 | 09 09 09 65 6c 73 65 0d |ik];..-.|...else.|
|00001d50| 0a 2d 09 09 09 09 09 72 | 65 73 75 6c 74 20 2b 3d |.-.....r|esult +=|
|00001d60| 20 2a 66 69 6c 74 50 20 | 2a 20 70 69 63 50 5b 69 | *filtP |* picP[i|
|00001d70| 6b 5d 3b 0d 0a 2d 09 09 | 09 09 66 69 6c 74 50 2b |k];..-..|..filtP+|
|00001d80| 2b 3b 0d 0a 2d 09 09 09 | 7d 0d 0a 2d 09 09 09 2a |+;..-...|}..-...*|
|00001d90| 28 74 6d 70 50 2b 69 78 | 29 20 3d 20 72 65 73 75 |(tmpP+ix|) = resu|
|00001da0| 6c 74 3b 0d 0a 2d 09 09 | 7d 0d 0a 2d 09 7d 0d 0a |lt;..-..|}..-.}..|
|00001db0| 2d 09 0d 0a 2d 2f 2a 0d | 0a 2d 20 2a 20 74 68 65 |-...-/*.|.- * the|
|00001dc0| 6e 20 69 6e 20 79 2d 64 | 69 72 6e 0d 0a 2d 20 2a |n in y-d|irn..- *|
|00001dd0| 09 09 66 69 72 73 74 20 | 74 6f 70 0d 0a 2d 20 2a |..first |top..- *|
|00001de0| 09 09 74 68 65 6e 20 6d | 69 64 64 6c 65 0d 0a 2d |..then m|iddle..-|
|00001df0| 20 2a 09 09 74 68 65 6e | 20 62 6f 74 74 6f 6d 0d | *..then| bottom.|
|00001e00| 0a 2d 2a 2f 0d 0a 2d 09 | 69 74 20 3d 20 2d 6e 66 |.-*/..-.|it = -nf|
|00001e10| 32 2a 6e 78 3b 0d 0a 2d | 09 66 6f 72 28 69 78 3d |2*nx;..-|.for(ix=|
|00001e20| 30 3b 20 69 78 3c 6e 78 | 3b 20 69 78 2b 2b 29 0d |0; ix<nx|; ix++).|
|00001e30| 0a 2d 09 7b 0d 0a 2d 09 | 09 67 61 75 73 73 69 61 |.-.{..-.|.gaussia|
|00001e40| 6e 50 20 3d 20 26 6d 79 | 5f 69 6d 61 67 65 2d 3e |nP = &my|_image->|
|00001e50| 67 61 75 73 73 69 61 6e | 5b 69 78 5d 3b 0d 0a 2d |gaussian|[ix];..-|
|00001e60| 09 09 74 6d 70 50 20 3d | 20 26 74 65 6d 70 5f 69 |..tmpP =| &temp_i|
|00001e70| 6d 61 67 65 5b 69 78 5d | 3b 0d 0a 2d 09 09 69 79 |mage[ix]|;..-..iy|
|00001e80| 20 3d 20 30 3b 20 69 79 | 4d 6e 66 32 58 6e 78 20 | = 0; iy|Mnf2Xnx |
|00001e90| 3d 20 69 74 3b 09 09 2f | 2a 20 3d 20 28 69 79 20 |= it;../|* = (iy |
|00001ea0| 2d 20 6e 66 32 29 2a 6e | 78 3b 20 2a 2f 0d 0a 2d |- nf2)*n|x; */..-|
|00001eb0| 09 09 66 6f 72 28 3b 20 | 69 79 3c 6e 66 32 3b 20 |..for(; |iy<nf2; |
|00001ec0| 69 79 2b 2b 2c 20 69 79 | 4d 6e 66 32 58 6e 78 2b |iy++, iy|Mnf2Xnx+|
|00001ed0| 3d 6e 78 29 0d 0a 2d 09 | 09 7b 0d 0a 2d 09 09 09 |=nx)..-.|.{..-...|
|00001ee0| 66 69 6c 74 50 20 3d 20 | 66 69 6c 74 65 72 3b 0d |filtP = |filter;.|
|00001ef0| 0a 2d 09 09 09 72 65 73 | 75 6c 74 20 3d 20 30 3b |.-...res|ult = 0;|
|00001f00| 0d 0a 2d 09 2f 2a 09 09 | 69 6b 20 3d 20 69 79 20 |..-./*..|ik = iy |
|00001f10| 2d 20 6e 66 32 3b 0d 0a | 2d 09 20 2a 09 09 69 6b |- nf2;..|-. *..ik|
|00001f20| 6c 69 6d 20 3d 20 69 6b | 20 2b 20 6e 66 3b 0d 0a |lim = ik| + nf;..|
|00001f30| 2d 09 20 2a 09 09 66 6f | 72 28 3b 20 69 6b 3c 69 |-. *..fo|r(; ik<i|
|00001f40| 6b 6c 69 6d 3b 20 69 6b | 2b 2b 29 0d 0a 2d 09 20 |klim; ik|++)..-. |
|00001f50| 2a 09 09 09 72 65 73 75 | 6c 74 20 2b 3d 20 28 2a |*...resu|lt += (*|
|00001f60| 66 69 6c 74 50 2b 2b 29 | 20 2a 20 74 6d 70 50 5b |filtP++)| * tmpP[|
|00001f70| 61 62 73 28 69 6b 29 2a | 6e 78 5d 3b 0d 0a 2d 09 |abs(ik)*|nx];..-.|
|00001f80| 20 2a 2f 0d 0a 2d 09 09 | 09 69 6b 20 3d 20 69 79 | */..-..|.ik = iy|
|00001f90| 4d 6e 66 32 58 6e 78 3b | 0d 0a 2d 09 09 09 69 6b |Mnf2Xnx;|..-...ik|
|00001fa0| 6c 69 6d 20 3d 20 69 6b | 20 2b 20 6e 66 58 6e 78 |lim = ik| + nfXnx|
|00001fb0| 3b 0d 0a 2d 20 09 09 09 | 66 6f 72 28 3b 20 69 6b |;..- ...|for(; ik|
|00001fc0| 3c 69 6b 6c 69 6d 3b 20 | 69 6b 2b 3d 6e 78 29 0d |<iklim; |ik+=nx).|
|00001fd0| 0a 2d 09 09 09 09 72 65 | 73 75 6c 74 20 2b 3d 20 |.-....re|sult += |
|00001fe0| 28 2a 66 69 6c 74 50 2b | 2b 29 20 2a 20 74 6d 70 |(*filtP+|+) * tmp|
|00001ff0| 50 5b 61 62 73 28 69 6b | 29 5d 3b 0d 0a 2d 0d 0a |P[abs(ik|)];..-..|
|00002000| 2d 09 2f 2a 09 09 6d 79 | 5f 69 6d 61 67 65 2d 3e |-./*..my|_image->|
|00002010| 67 61 75 73 73 69 61 6e | 5b 69 79 2a 6e 78 2b 69 |gaussian|[iy*nx+i|
|00002020| 78 5d 20 3d 20 72 65 73 | 75 6c 74 2f 6d 79 5f 69 |x] = res|ult/my_i|
|00002030| 6d 61 67 65 2d 3e 6e 6f | 72 6d 3b 20 2a 2f 0d 0a |mage->no|rm; */..|
|00002040| 2d 09 09 09 2a 28 67 61 | 75 73 73 69 61 6e 50 29 |-...*(ga|ussianP)|
|00002050| 20 3d 20 72 65 73 75 6c | 74 2f 6d 79 5f 69 6d 61 | = resul|t/my_ima|
|00002060| 67 65 2d 3e 6e 6f 72 6d | 3b 0d 0a 2d 09 09 09 67 |ge->norm|;..-...g|
|00002070| 61 75 73 73 69 61 6e 50 | 20 2b 3d 20 6e 78 3b 0d |aussianP| += nx;.|
|00002080| 0a 2d 09 09 7d 0d 0a 2d | 09 7d 0d 0a 2d 0d 0a 2d |.-..}..-|.}..-..-|
|00002090| 09 69 79 30 20 3d 20 6e | 66 32 3b 0d 0a 2d 09 69 |.iy0 = n|f2;..-.i|
|000020a0| 79 30 58 6e 78 20 3d 20 | 69 79 30 20 2a 20 6e 78 |y0Xnx = |iy0 * nx|
|000020b0| 3b 0d 0a 2d 09 66 6f 72 | 28 69 78 3d 30 3b 20 69 |;..-.for|(ix=0; i|
|000020c0| 78 3c 6e 78 3b 20 69 78 | 2b 2b 29 0d 0a 2d 09 7b |x<nx; ix|++)..-.{|
|000020d0| 0d 0a 2d 09 09 69 6e 74 | 20 2a 63 6f 6c 5f 74 6d |..-..int| *col_tm|
|000020e0| 70 50 20 3d 20 26 74 65 | 6d 70 5f 69 6d 61 67 65 |pP = &te|mp_image|
|000020f0| 5b 69 78 5d 3b 0d 0a 2d | 09 09 69 6e 74 20 79 6c |[ix];..-|..int yl|
|00002100| 69 6d 20 3d 20 6e 79 2d | 6e 66 32 2d 31 3b 0d 0a |im = ny-|nf2-1;..|
|00002110| 2d 09 09 69 79 20 3d 20 | 69 79 30 3b 20 67 61 75 |-..iy = |iy0; gau|
|00002120| 73 73 69 61 6e 50 20 3d | 20 26 6d 79 5f 69 6d 61 |ssianP =| &my_ima|
|00002130| 67 65 2d 3e 67 61 75 73 | 73 69 61 6e 5b 69 79 30 |ge->gaus|sian[iy0|
|00002140| 58 6e 78 2b 69 78 5d 3b | 0d 0a 2d 09 09 69 79 4d |Xnx+ix];|..-..iyM|
|00002150| 6e 66 32 58 6e 78 20 3d | 20 30 3b 09 09 09 2f 2a |nf2Xnx =| 0;.../*|
|00002160| 20 3d 20 28 69 79 2d 6e | 66 32 29 2a 6e 78 3b 20 | = (iy-n|f2)*nx; |
|00002170| 2a 2f 0d 0a 2d 09 09 66 | 6f 72 28 20 3b 20 69 79 |*/..-..f|or( ; iy|
|00002180| 3c 79 6c 69 6d 3b 20 69 | 79 2b 2b 2c 20 69 79 4d |<ylim; i|y++, iyM|
|00002190| 6e 66 32 58 6e 78 2b 3d | 6e 78 29 0d 0a 2d 09 09 |nf2Xnx+=|nx)..-..|
|000021a0| 7b 0d 0a 2d 09 09 09 66 | 69 6c 74 50 20 3d 20 66 |{..-...f|iltP = f|
|000021b0| 69 6c 74 65 72 3b 0d 0a | 2d 09 09 09 74 6d 70 50 |ilter;..|-...tmpP|
|000021c0| 20 3d 20 63 6f 6c 5f 74 | 6d 70 50 20 2b 20 69 79 | = col_t|mpP + iy|
|000021d0| 4d 6e 66 32 58 6e 78 3b | 0d 0a 2d 09 09 09 72 65 |Mnf2Xnx;|..-...re|
|000021e0| 73 75 6c 74 20 3d 20 30 | 3b 0d 0a 2d 09 09 09 69 |sult = 0|;..-...i|
|000021f0| 6b 3d 6e 66 3b 0d 0a 2d | 09 09 09 77 68 69 6c 65 |k=nf;..-|...while|
|00002200| 28 69 6b 2d 2d 29 0d 0a | 2d 09 09 09 7b 0d 0a 2d |(ik--)..|-...{..-|
|00002210| 09 09 09 09 72 65 73 75 | 6c 74 20 2b 3d 20 28 2a |....resu|lt += (*|
|00002220| 66 69 6c 74 50 2b 2b 29 | 20 2a 20 28 2a 74 6d 70 |filtP++)| * (*tmp|
|00002230| 50 29 3b 0d 0a 2d 09 09 | 09 09 74 6d 70 50 20 2b |P);..-..|..tmpP +|
|00002240| 3d 20 6e 78 3b 0d 0a 2d | 09 09 09 7d 0d 0a 2d 09 |= nx;..-|...}..-.|
|00002250| 2f 2a 09 09 6d 79 5f 69 | 6d 61 67 65 2d 3e 67 61 |/*..my_i|mage->ga|
|00002260| 75 73 73 69 61 6e 5b 69 | 79 2a 6e 78 2b 69 78 5d |ussian[i|y*nx+ix]|
|00002270| 20 3d 20 72 65 73 75 6c | 74 2f 6d 79 5f 69 6d 61 | = resul|t/my_ima|
|00002280| 67 65 2d 3e 6e 6f 72 6d | 3b 20 2a 2f 0d 0a 2d 09 |ge->norm|; */..-.|
|00002290| 09 09 2a 28 67 61 75 73 | 73 69 61 6e 50 29 20 3d |..*(gaus|sianP) =|
|000022a0| 20 72 65 73 75 6c 74 2f | 6d 79 5f 69 6d 61 67 65 | result/|my_image|
|000022b0| 2d 3e 6e 6f 72 6d 3b 0d | 0a 2d 09 09 09 67 61 75 |->norm;.|.-...gau|
|000022c0| 73 73 69 61 6e 50 20 2b | 3d 20 6e 78 3b 0d 0a 2d |ssianP +|= nx;..-|
|000022d0| 09 09 7d 0d 0a 2d 09 7d | 0d 0a 2d 0d 0a 2d 09 69 |..}..-.}|..-..-.i|
|000022e0| 79 30 20 3d 20 6e 79 2d | 6e 66 32 2d 31 3b 0d 0a |y0 = ny-|nf2-1;..|
|000022f0| 2d 09 69 79 30 58 6e 78 | 20 3d 20 69 79 30 2a 6e |-.iy0Xnx| = iy0*n|
|00002300| 78 3b 0d 0a 2d 09 69 74 | 20 3d 20 28 69 79 30 2d |x;..-.it| = (iy0-|
|00002310| 6e 66 32 29 2a 6e 78 3b | 0d 0a 2d 09 66 6f 72 28 |nf2)*nx;|..-.for(|
|00002320| 69 78 3d 30 3b 20 69 78 | 3c 6e 78 3b 20 69 78 2b |ix=0; ix|<nx; ix+|
|00002330| 2b 29 0d 0a 2d 09 7b 0d | 0a 2d 09 09 74 6d 70 50 |+)..-.{.|.-..tmpP|
|00002340| 20 3d 20 26 74 65 6d 70 | 5f 69 6d 61 67 65 5b 69 | = &temp|_image[i|
|00002350| 78 5d 3b 0d 0a 2d 09 09 | 69 79 20 3d 20 69 79 30 |x];..-..|iy = iy0|
|00002360| 3b 20 67 61 75 73 73 69 | 61 6e 50 20 3d 20 26 6d |; gaussi|anP = &m|
|00002370| 79 5f 69 6d 61 67 65 2d | 3e 67 61 75 73 73 69 61 |y_image-|>gaussia|
|00002380| 6e 5b 69 79 30 58 6e 78 | 2b 69 78 5d 3b 0d 0a 2d |n[iy0Xnx|+ix];..-|
|00002390| 09 09 69 79 4d 6e 66 32 | 58 6e 78 20 3d 20 69 74 |..iyMnf2|Xnx = it|
|000023a0| 3b 09 09 09 2f 2a 20 3d | 20 28 69 79 2d 6e 66 32 |;.../* =| (iy-nf2|
|000023b0| 29 2a 6e 78 3b 20 2a 2f | 0d 0a 2d 09 09 66 6f 72 |)*nx; */|..-..for|
|000023c0| 28 20 3b 20 69 79 3c 6e | 79 3b 20 69 79 2b 2b 2c |( ; iy<n|y; iy++,|
|000023d0| 20 69 79 4d 6e 66 32 58 | 6e 78 2b 3d 6e 78 29 0d | iyMnf2X|nx+=nx).|
|000023e0| 0a 2d 09 09 7b 0d 0a 2d | 09 09 09 66 69 6c 74 50 |.-..{..-|...filtP|
|000023f0| 20 3d 20 66 69 6c 74 65 | 72 3b 0d 0a 2d 09 09 09 | = filte|r;..-...|
|00002400| 72 65 73 75 6c 74 20 3d | 20 30 3b 0d 0a 2d 09 09 |result =| 0;..-..|
|00002410| 09 69 6b 20 3d 20 69 79 | 20 2d 20 6e 66 32 3b 0d |.ik = iy| - nf2;.|
|00002420| 0a 2d 09 09 09 69 6b 6c | 69 6d 20 3d 20 69 6b 20 |.-...ikl|im = ik |
|00002430| 2b 20 6e 66 3b 0d 0a 2d | 09 09 09 69 6b 58 6e 78 |+ nf;..-|...ikXnx|
|00002440| 20 3d 20 69 79 4d 6e 66 | 32 58 6e 78 3b 0d 0a 2d | = iyMnf|2Xnx;..-|
|00002450| 09 09 09 66 6f 72 28 3b | 20 69 6b 3c 69 6b 6c 69 |...for(;| ik<ikli|
|00002460| 6d 3b 20 69 6b 2b 2b 29 | 0d 0a 2d 09 09 09 7b 0d |m; ik++)|..-...{.|
|00002470| 0a 2d 09 09 2f 2a 09 09 | 72 65 67 69 73 74 65 72 |.-../*..|register|
|00002480| 20 69 6e 74 20 69 6e 64 | 65 78 20 3d 20 69 79 4d | int ind|ex = iyM|
|00002490| 6e 66 32 2b 69 6b 3b 0d | 0a 2d 09 09 20 2a 09 09 |nf2+ik;.|.-.. *..|
|000024a0| 69 66 28 69 6e 64 65 78 | 3e 3d 6e 79 29 20 69 6e |if(index|>=ny) in|
|000024b0| 64 65 78 20 3d 20 6e 79 | 2b 6e 79 2d 69 6e 64 65 |dex = ny|+ny-inde|
|000024c0| 78 2d 31 3b 0d 0a 2d 09 | 09 20 2a 09 09 72 65 73 |x-1;..-.|. *..res|
|000024d0| 75 6c 74 20 2b 3d 20 28 | 2a 66 69 6c 74 50 2b 2b |ult += (|*filtP++|
|000024e0| 29 20 2a 20 74 6d 70 50 | 5b 69 6e 64 65 78 2a 6e |) * tmpP|[index*n|
|000024f0| 78 5d 3b 0d 0a 2d 09 09 | 2a 2f 0d 0a 2d 09 09 09 |x];..-..|*/..-...|
|00002500| 09 69 66 28 69 6b 3e 3d | 6e 79 29 7b 09 20 2f 2a |.if(ik>=|ny){. /*|
|00002510| 20 6e 79 58 32 4d 31 58 | 6e 78 20 3d 20 28 6e 79 | nyX2M1X|nx = (ny|
|00002520| 2b 6e 79 2d 31 29 2a 6e | 78 3b 20 2a 2f 0d 0a 2d |+ny-1)*n|x; */..-|
|00002530| 09 09 09 09 20 20 20 20 | 72 65 73 75 6c 74 20 2b |....    |result +|
|00002540| 3d 20 2a 66 69 6c 74 50 | 20 2a 20 74 6d 70 50 5b |= *filtP| * tmpP[|
|00002550| 6e 79 58 32 4d 31 58 6e | 78 2d 69 6b 58 6e 78 5d |nyX2M1Xn|x-ikXnx]|
|00002560| 3b 0d 0a 2d 09 09 09 09 | 7d 65 6c 73 65 0d 0a 2d |;..-....|}else..-|
|00002570| 09 09 09 09 20 20 20 20 | 72 65 73 75 6c 74 20 2b |....    |result +|
|00002580| 3d 20 2a 66 69 6c 74 50 | 20 2a 20 74 6d 70 50 5b |= *filtP| * tmpP[|
|00002590| 69 6b 58 6e 78 5d 3b 0d | 0a 2d 09 09 09 09 66 69 |ikXnx];.|.-....fi|
|000025a0| 6c 74 50 2b 2b 3b 0d 0a | 2d 09 09 09 09 69 6b 58 |ltP++;..|-....ikX|
|000025b0| 6e 78 20 2b 3d 20 6e 78 | 3b 0d 0a 2d 09 09 09 7d |nx += nx|;..-...}|
|000025c0| 0d 0a 2d 09 2f 2a 09 09 | 6d 79 5f 69 6d 61 67 65 |..-./*..|my_image|
|000025d0| 2d 3e 67 61 75 73 73 69 | 61 6e 5b 69 79 2a 6e 78 |->gaussi|an[iy*nx|
|000025e0| 2b 69 78 5d 20 3d 20 72 | 65 73 75 6c 74 2f 6d 79 |+ix] = r|esult/my|
|000025f0| 5f 69 6d 61 67 65 2d 3e | 6e 6f 72 6d 3b 20 2a 2f |_image->|norm; */|
|00002600| 0d 0a 2d 09 09 09 2a 28 | 67 61 75 73 73 69 61 6e |..-...*(|gaussian|
|00002610| 50 29 20 3d 20 72 65 73 | 75 6c 74 2f 6d 79 5f 69 |P) = res|ult/my_i|
|00002620| 6d 61 67 65 2d 3e 6e 6f | 72 6d 3b 0d 0a 2d 09 09 |mage->no|rm;..-..|
|00002630| 09 67 61 75 73 73 69 61 | 6e 50 20 2b 3d 20 6e 78 |.gaussia|nP += nx|
|00002640| 3b 0d 0a 2d 09 09 7d 0d | 0a 2d 09 7d 0d 0a 2d 09 |;..-..}.|.-.}..-.|
|00002650| 66 72 65 65 28 74 65 6d | 70 5f 69 6d 61 67 65 29 |free(tem|p_image)|
|00002660| 3b 0d 0a 2d 23 69 66 64 | 65 66 20 44 45 42 55 47 |;..-#ifd|ef DEBUG|
|00002670| 0d 0a 2d 09 69 6d 61 67 | 65 5f 57 72 69 74 65 5f |..-.imag|e_Write_|
|00002680| 69 6e 74 28 22 67 61 75 | 73 73 22 2c 20 6d 79 5f |int("gau|ss", my_|
|00002690| 69 6d 61 67 65 2d 3e 67 | 61 75 73 73 69 61 6e 29 |image->g|aussian)|
|000026a0| 3b 0d 0a 2d 23 65 6e 64 | 69 66 0d 0a 2d 7d 0d 0a |;..-#end|if..-}..|
|000026b0| 2f 2f 47 4f 2e 53 59 53 | 49 4e 20 44 44 20 63 6f |//GO.SYS|IN DD co|
|000026c0| 6e 76 2e 63 0d 0a 65 63 | 68 6f 20 63 74 6f 72 2e |nv.c..ec|ho ctor.|
|000026d0| 63 20 31 3e 26 32 0d 0a | 73 65 64 20 27 73 2f 2e |c 1>&2..|sed 's/.|
|000026e0| 2f 2f 27 20 3e 63 74 6f | 72 2e 63 20 3c 3c 27 2f |//' >cto|r.c <<'/|
|000026f0| 2f 47 4f 2e 53 59 53 49 | 4e 20 44 44 20 63 74 6f |/GO.SYSI|N DD cto|
|00002700| 72 2e 63 27 0d 0a 2d 2f | 2a 0d 0a 2d 20 2a 20 28 |r.c'..-/|*..- * (|
|00002710| 63 29 20 43 6f 70 79 72 | 69 67 68 74 20 31 39 38 |c) Copyr|ight 198|
|00002720| 38 20 62 79 0d 0a 2d 20 | 2a 20 52 6f 62 6f 74 69 |8 by..- |* Roboti|
|00002730| 63 73 20 50 72 69 6e 63 | 69 70 6c 65 73 20 52 65 |cs Princ|iples Re|
|00002740| 73 65 61 72 63 68 20 44 | 65 70 61 72 74 6d 65 6e |search D|epartmen|
|00002750| 74 2c 20 41 54 54 20 42 | 65 6c 6c 20 4c 61 62 6f |t, ATT B|ell Labo|
|00002760| 72 61 74 6f 72 69 65 73 | 2e 0d 0a 2d 20 2a 20 41 |ratories|...- * A|
|00002770| 6c 6c 20 72 69 67 68 74 | 73 20 72 65 73 65 72 76 |ll right|s reserv|
|00002780| 65 64 2e 0d 0a 2d 20 2a | 20 4c 61 73 74 20 6d 6f |ed...- *| Last mo|
|00002790| 64 69 66 69 65 64 20 32 | 2f 38 2f 38 38 20 49 6e |dified 2|/8/88 In|
|000027a0| 67 65 6d 61 72 20 4a 2e | 20 43 6f 78 0d 0a 2d 20 |gemar J.| Cox..- |
|000027b0| 2a 20 43 20 76 65 72 73 | 69 6f 6e 20 38 2f 32 2f |* C vers|ion 8/2/|
|000027c0| 38 38 20 44 65 62 6f 72 | 61 68 20 41 2e 20 57 61 |88 Debor|ah A. Wa|
|000027d0| 6c 6c 61 63 68 0d 0a 2d | 20 2a 2f 0d 0a 2d 23 69 |llach..-| */..-#i|
|000027e0| 6e 63 6c 75 64 65 20 3c | 73 74 64 69 6f 2e 68 3e |nclude <|stdio.h>|
|000027f0| 0d 0a 2d 23 69 6e 63 6c | 75 64 65 20 3c 6d 61 74 |..-#incl|ude <mat|
|00002800| 68 2e 68 3e 0d 0a 2d 23 | 69 6e 63 6c 75 64 65 20 |h.h>..-#|include |
|00002810| 22 65 64 67 65 5f 66 69 | 6e 64 65 72 2e 68 22 0d |"edge_fi|nder.h".|
|00002820| 0a 2d 0d 0a 2d 73 74 72 | 75 63 74 20 69 6d 61 67 |.-..-str|uct imag|
|00002830| 65 20 2a 6d 79 5f 69 6d | 61 67 65 3b 0d 0a 2d 0d |e *my_im|age;..-.|
|00002840| 0a 2d 69 6d 61 67 65 5f | 69 6e 69 74 28 61 72 72 |.-image_|init(arr|
|00002850| 61 79 2c 20 78 64 69 6d | 2c 20 79 64 69 6d 29 0d |ay, xdim|, ydim).|
|00002860| 0a 2d 20 20 20 20 20 75 | 6e 73 69 67 6e 65 64 20 |.-     u|nsigned |
|00002870| 63 68 61 72 20 2a 61 72 | 72 61 79 3b 0d 0a 2d 20 |char *ar|ray;..- |
|00002880| 20 20 20 20 69 6e 74 20 | 78 64 69 6d 2c 20 79 64 |    int |xdim, yd|
|00002890| 69 6d 3b 0d 0a 2d 7b 0d | 0a 2d 20 20 20 20 73 74 |im;..-{.|.-    st|
|000028a0| 61 74 69 63 20 73 74 72 | 75 63 74 20 69 6d 61 67 |atic str|uct imag|
|000028b0| 65 20 74 65 6d 70 5f 69 | 6d 61 67 65 3b 0d 0a 2d |e temp_i|mage;..-|
|000028c0| 0d 0a 2d 09 74 65 6d 70 | 5f 69 6d 61 67 65 2e 67 |..-.temp|_image.g|
|000028d0| 61 75 73 73 69 61 6e 20 | 3d 20 30 3b 0d 0a 2d 09 |aussian |= 0;..-.|
|000028e0| 74 65 6d 70 5f 69 6d 61 | 67 65 2e 66 69 6c 74 65 |temp_ima|ge.filte|
|000028f0| 72 20 3d 20 30 3b 0d 0a | 2d 09 74 65 6d 70 5f 69 |r = 0;..|-.temp_i|
|00002900| 6d 61 67 65 2e 69 64 78 | 3d 30 3b 0d 0a 2d 09 74 |mage.idx|=0;..-.t|
|00002910| 65 6d 70 5f 69 6d 61 67 | 65 2e 69 64 79 3d 30 3b |emp_imag|e.idy=0;|
|00002920| 0d 0a 2d 09 74 65 6d 70 | 5f 69 6d 61 67 65 2e 69 |..-.temp|_image.i|
|00002930| 64 34 35 3d 30 3b 0d 0a | 2d 09 74 65 6d 70 5f 69 |d45=0;..|-.temp_i|
|00002940| 6d 61 67 65 2e 69 64 31 | 33 35 3d 30 3b 0d 0a 2d |mage.id1|35=0;..-|
|00002950| 09 74 65 6d 70 5f 69 6d | 61 67 65 2e 65 64 67 65 |.temp_im|age.edge|
|00002960| 5f 6d 61 70 3d 30 3b 0d | 0a 2d 09 74 65 6d 70 5f |_map=0;.|.-.temp_|
|00002970| 69 6d 61 67 65 2e 61 6e | 5f 69 6d 61 67 65 20 3d |image.an|_image =|
|00002980| 20 61 72 72 61 79 3b 0d | 0a 2d 09 74 65 6d 70 5f | array;.|.-.temp_|
|00002990| 69 6d 61 67 65 2e 6e 78 | 3d 78 64 69 6d 3b 0d 0a |image.nx|=xdim;..|
|000029a0| 2d 09 74 65 6d 70 5f 69 | 6d 61 67 65 2e 6e 79 3d |-.temp_i|mage.ny=|
|000029b0| 79 64 69 6d 3b 0d 0a 2d | 20 20 20 20 20 20 20 20 |ydim;..-|        |
|000029c0| 6d 79 5f 69 6d 61 67 65 | 20 3d 20 26 74 65 6d 70 |my_image| = &temp|
|000029d0| 5f 69 6d 61 67 65 3b 0d | 0a 2d 20 20 0d 0a 2d 7d |_image;.|.-  ..-}|
|000029e0| 0d 0a 2d 0d 0a 2f 2f 47 | 4f 2e 53 59 53 49 4e 20 |..-..//G|O.SYSIN |
|000029f0| 44 44 20 63 74 6f 72 2e | 63 0d 0a 65 63 68 6f 20 |DD ctor.|c..echo |
|00002a00| 64 65 74 2e 63 20 31 3e | 26 32 0d 0a 73 65 64 20 |det.c 1>|&2..sed |
|00002a10| 27 73 2f 2e 2f 2f 27 20 | 3e 64 65 74 2e 63 20 3c |'s/.//' |>det.c <|
|00002a20| 3c 27 2f 2f 47 4f 2e 53 | 59 53 49 4e 20 44 44 20 |<'//GO.S|YSIN DD |
|00002a30| 64 65 74 2e 63 27 0d 0a | 2d 2f 2a 0d 0a 2d 20 2a |det.c'..|-/*..- *|
|00002a40| 20 28 63 29 20 43 6f 70 | 79 72 69 67 68 74 20 31 | (c) Cop|yright 1|
|00002a50| 39 38 38 20 62 79 0d 0a | 2d 20 2a 20 52 6f 62 6f |988 by..|- * Robo|
|00002a60| 74 69 63 73 20 50 72 69 | 6e 63 69 70 6c 65 73 20 |tics Pri|nciples |
|00002a70| 52 65 73 65 61 72 63 68 | 20 44 65 70 61 72 74 6d |Research| Departm|
|00002a80| 65 6e 74 2c 20 41 54 54 | 20 42 65 6c 6c 20 4c 61 |ent, ATT| Bell La|
|00002a90| 62 6f 72 61 74 6f 72 69 | 65 73 2e 0d 0a 2d 20 2a |boratori|es...- *|
|00002aa0| 20 41 6c 6c 20 72 69 67 | 68 74 73 20 72 65 73 65 | All rig|hts rese|
|00002ab0| 72 76 65 64 2e 0d 0a 2d | 20 2a 20 4c 61 73 74 20 |rved...-| * Last |
|00002ac0| 6d 6f 64 69 66 69 65 64 | 20 32 2f 38 2f 38 38 20 |modified| 2/8/88 |
|00002ad0| 49 6e 67 65 6d 61 72 20 | 4a 2e 20 43 6f 78 0d 0a |Ingemar |J. Cox..|
|00002ae0| 2d 20 2a 20 43 20 76 65 | 72 73 69 6f 6e 20 38 2f |- * C ve|rsion 8/|
|00002af0| 32 2f 38 38 20 44 65 62 | 6f 72 61 68 20 41 2e 20 |2/88 Deb|orah A. |
|00002b00| 57 61 6c 6c 61 63 68 0d | 0a 2d 20 2a 20 43 20 76 |Wallach.|.- * C v|
|00002b10| 65 72 73 69 6f 6e 20 31 | 31 2f 34 2f 38 38 20 57 |ersion 1|1/4/88 W|
|00002b20| 2e 20 4a 2e 20 4b 72 6f | 70 66 6c 0d 0a 2d 20 2a |. J. Kro|pfl..- *|
|00002b30| 2f 0d 0a 2d 23 69 6e 63 | 6c 75 64 65 20 3c 73 74 |/..-#inc|lude <st|
|00002b40| 64 69 6f 2e 68 3e 0d 0a | 2d 23 69 6e 63 6c 75 64 |dio.h>..|-#includ|
|00002b50| 65 20 3c 6d 61 74 68 2e | 68 3e 0d 0a 2d 23 69 6e |e <math.|h>..-#in|
|00002b60| 63 6c 75 64 65 20 22 65 | 64 67 65 5f 66 69 6e 64 |clude "e|dge_find|
|00002b70| 65 72 2e 68 22 0d 0a 2d | 0d 0a 2d 65 78 74 65 72 |er.h"..-|..-exter|
|00002b80| 6e 20 73 74 72 75 63 74 | 20 69 6d 61 67 65 20 2a |n struct| image *|
|00002b90| 6d 79 5f 69 6d 61 67 65 | 3b 0d 0a 2d 0d 0a 2d 76 |my_image|;..-..-v|
|00002ba0| 6f 69 64 0d 0a 2d 69 6d | 61 67 65 5f 64 65 74 78 |oid..-im|age_detx|
|00002bb0| 28 29 0d 0a 2d 7b 0d 0a | 2d 09 72 65 67 69 73 74 |()..-{..|-.regist|
|00002bc0| 65 72 20 69 6e 74 20 2a | 64 50 2c 20 2a 67 50 3b |er int *|dP, *gP;|
|00002bd0| 0d 0a 2d 09 72 65 67 69 | 73 74 65 72 20 69 6e 74 |..-.regi|ster int|
|00002be0| 20 69 78 2c 20 69 79 3b | 0d 0a 2d 09 72 65 67 69 | ix, iy;|..-.regi|
|00002bf0| 73 74 65 72 20 69 6e 74 | 20 6e 78 2c 20 6e 79 2c |ster int| nx, ny,|
|00002c00| 20 6e 78 4d 31 3b 0d 0a | 2d 0d 0a 2d 09 6e 78 20 | nxM1;..|-..-.nx |
|00002c10| 3d 20 6d 79 5f 69 6d 61 | 67 65 2d 3e 6e 78 3b 09 |= my_ima|ge->nx;.|
|00002c20| 6e 78 4d 31 20 3d 20 6e | 78 2d 31 3b 0d 0a 2d 09 |nxM1 = n|x-1;..-.|
|00002c30| 6e 79 20 3d 20 6d 79 5f | 69 6d 61 67 65 2d 3e 6e |ny = my_|image->n|
|00002c40| 79 3b 0d 0a 2d 09 0d 0a | 2d 09 69 66 28 28 6d 79 |y;..-...|-.if((my|
|00002c50| 5f 69 6d 61 67 65 2d 3e | 69 64 78 20 3d 20 28 69 |_image->|idx = (i|
|00002c60| 6e 74 20 2a 29 6d 61 6c | 6c 6f 63 28 6e 79 2a 6e |nt *)mal|loc(ny*n|
|00002c70| 78 2a 73 69 7a 65 6f 66 | 28 69 6e 74 29 29 29 3d |x*sizeof|(int)))=|
|00002c80| 3d 30 29 0d 0a 2d 09 7b | 0d 0a 2d 09 09 66 70 72 |=0)..-.{|..-..fpr|
|00002c90| 69 6e 74 66 28 73 74 64 | 65 72 72 2c 20 22 65 72 |intf(std|err, "er|
|00002ca0| 72 6f 72 3a 20 63 61 6e | 6e 6f 74 20 61 6c 6c 6f |ror: can|not allo|
|00002cb0| 63 61 74 65 20 69 64 78 | 20 61 72 72 61 79 5c 6e |cate idx| array\n|
|00002cc0| 22 29 3b 0d 0a 2d 09 09 | 65 78 69 74 28 31 29 3b |");..-..|exit(1);|
|00002cd0| 0d 0a 2d 09 7d 0d 0a 2d | 09 0d 0a 2d 09 64 50 20 |..-.}..-|...-.dP |
|00002ce0| 3d 20 6d 79 5f 69 6d 61 | 67 65 2d 3e 69 64 78 3b |= my_ima|ge->idx;|
|00002cf0| 0d 0a 2d 09 67 50 20 3d | 20 6d 79 5f 69 6d 61 67 |..-.gP =| my_imag|
|00002d00| 65 2d 3e 67 61 75 73 73 | 69 61 6e 3b 0d 0a 2d 09 |e->gauss|ian;..-.|
|00002d10| 66 6f 72 28 69 79 3d 30 | 3b 20 69 79 3c 6e 79 3b |for(iy=0|; iy<ny;|
|00002d20| 20 69 79 2b 2b 29 0d 0a | 2d 09 7b 0d 0a 2d 2f 2a | iy++)..|-.{..-/*|
|00002d30| 09 09 64 5b 69 79 2a 6e | 78 5d 20 3d 20 67 5b 69 |..d[iy*n|x] = g[i|
|00002d40| 79 2a 6e 78 5d 20 2d 20 | 67 5b 69 79 2a 6e 78 2b |y*nx] - |g[iy*nx+|
|00002d50| 31 5d 3b 20 2a 2f 0d 0a | 2d 09 09 2a 64 50 20 3d |1]; */..|-..*dP =|
|00002d60| 20 2a 67 50 20 2d 20 2a | 28 67 50 2b 31 29 3b 0d | *gP - *|(gP+1);.|
|00002d70| 0a 2d 09 09 64 50 2b 3d | 6e 78 3b 20 67 50 2b 3d |.-..dP+=|nx; gP+=|
|00002d80| 6e 78 3b 0d 0a 2d 09 7d | 0d 0a 2d 09 66 6f 72 28 |nx;..-.}|..-.for(|
|00002d90| 69 79 3d 30 3b 20 69 79 | 3c 6e 79 3b 20 69 79 2b |iy=0; iy|<ny; iy+|
|00002da0| 2b 29 0d 0a 2d 09 7b 0d | 0a 2d 09 09 64 50 20 3d |+)..-.{.|.-..dP =|
|00002db0| 20 6d 79 5f 69 6d 61 67 | 65 2d 3e 69 64 78 2b 69 | my_imag|e->idx+i|
|00002dc0| 79 2a 6e 78 2b 31 3b 0d | 0a 2d 09 09 67 50 20 3d |y*nx+1;.|.-..gP =|
|00002dd0| 20 6d 79 5f 69 6d 61 67 | 65 2d 3e 67 61 75 73 73 | my_imag|e->gauss|
|00002de0| 69 61 6e 2b 69 79 2a 6e | 78 2b 32 3b 0d 0a 2d 09 |ian+iy*n|x+2;..-.|
|00002df0| 09 66 6f 72 28 69 78 3d | 31 3b 20 69 78 3c 6e 78 |.for(ix=|1; ix<nx|
|00002e00| 4d 31 3b 20 69 78 2b 2b | 29 0d 0a 2d 09 09 7b 0d |M1; ix++|)..-..{.|
|00002e10| 0a 2d 2f 2a 09 09 09 64 | 5b 69 78 5d 20 3d 20 67 |.-/*...d|[ix] = g|
|00002e20| 5b 69 78 2b 31 5d 20 2d | 20 67 5b 69 78 2d 31 5d |[ix+1] -| g[ix-1]|
|00002e30| 3b 20 2a 2f 0d 0a 2d 09 | 09 09 2a 64 50 20 3d 20 |; */..-.|..*dP = |
|00002e40| 2a 67 50 20 2d 20 2a 28 | 67 50 2d 32 29 3b 0d 0a |*gP - *(|gP-2);..|
|00002e50| 2d 09 09 09 64 50 2b 2b | 3b 20 67 50 2b 2b 3b 0d |-...dP++|; gP++;.|
|00002e60| 0a 2d 09 09 7d 0d 0a 2d | 09 7d 0d 0a 2d 09 64 50 |.-..}..-|.}..-.dP|
|00002e70| 20 3d 20 6d 79 5f 69 6d | 61 67 65 2d 3e 69 64 78 | = my_im|age->idx|
|00002e80| 2b 6e 78 2d 31 3b 0d 0a | 2d 09 67 50 20 3d 20 6d |+nx-1;..|-.gP = m|
|00002e90| 79 5f 69 6d 61 67 65 2d | 3e 67 61 75 73 73 69 61 |y_image-|>gaussia|
|00002ea0| 6e 2b 6e 78 2d 31 3b 0d | 0a 2d 09 66 6f 72 28 69 |n+nx-1;.|.-.for(i|
|00002eb0| 79 3d 30 3b 20 69 79 3c | 6e 79 3b 20 69 79 2b 2b |y=0; iy<|ny; iy++|
|00002ec0| 29 0d 0a 2d 09 7b 0d 0a | 2d 2f 2a 09 09 64 5b 69 |)..-.{..|-/*..d[i|
|00002ed0| 79 2a 6e 78 2b 6e 78 2d | 31 5d 20 3d 20 67 5b 69 |y*nx+nx-|1] = g[i|
|00002ee0| 79 2a 6e 78 2b 6e 78 2d | 31 5d 20 2d 20 67 5b 69 |y*nx+nx-|1] - g[i|
|00002ef0| 79 2a 6e 78 2b 6e 78 2d | 32 5d 3b 20 2a 2f 0d 0a |y*nx+nx-|2]; */..|
|00002f00| 2d 09 09 2a 64 50 20 3d | 20 2a 67 50 20 2d 20 2a |-..*dP =| *gP - *|
|00002f10| 28 67 50 2d 31 29 3b 0d | 0a 2d 09 09 64 50 2b 3d |(gP-1);.|.-..dP+=|
|00002f20| 6e 78 3b 20 67 50 2b 3d | 6e 78 3b 0d 0a 2d 09 7d |nx; gP+=|nx;..-.}|
|00002f30| 0d 0a 2d 23 69 66 64 65 | 66 20 44 45 42 55 47 0d |..-#ifde|f DEBUG.|
|00002f40| 0a 2d 09 69 6d 61 67 65 | 5f 57 72 69 74 65 5f 69 |.-.image|_Write_i|
|00002f50| 6e 74 28 22 64 65 74 78 | 22 2c 20 6d 79 5f 69 6d |nt("detx|", my_im|
|00002f60| 61 67 65 2d 3e 69 64 78 | 29 3b 0d 0a 2d 23 65 6e |age->idx|);..-#en|
|00002f70| 64 69 66 0d 0a 2d 7d 0d | 0a 2d 0d 0a 2d 76 6f 69 |dif..-}.|.-..-voi|
|00002f80| 64 0d 0a 2d 69 6d 61 67 | 65 5f 64 65 74 79 28 29 |d..-imag|e_dety()|
|00002f90| 0d 0a 2d 7b 0d 0a 2d 09 | 72 65 67 69 73 74 65 72 |..-{..-.|register|
|00002fa0| 20 69 6e 74 20 2a 64 50 | 2c 20 2a 67 50 3b 0d 0a | int *dP|, *gP;..|
|00002fb0| 2d 09 72 65 67 69 73 74 | 65 72 20 69 6e 74 20 69 |-.regist|er int i|
|00002fc0| 78 2c 20 69 79 2c 20 6e | 78 2c 20 6e 79 2c 20 6e |x, iy, n|x, ny, n|
|00002fd0| 79 4d 31 2c 20 6e 78 50 | 6e 78 3b 0d 0a 2d 0d 0a |yM1, nxP|nx;..-..|
|00002fe0| 2d 09 6e 78 20 3d 20 6d | 79 5f 69 6d 61 67 65 2d |-.nx = m|y_image-|
|00002ff0| 3e 6e 78 3b 09 6e 78 50 | 6e 78 20 3d 20 6e 78 2b |>nx;.nxP|nx = nx+|
|00003000| 6e 78 3b 0d 0a 2d 09 6e | 79 20 3d 20 6d 79 5f 69 |nx;..-.n|y = my_i|
|00003010| 6d 61 67 65 2d 3e 6e 79 | 3b 09 6e 79 4d 31 20 3d |mage->ny|;.nyM1 =|
|00003020| 20 6e 79 2d 31 3b 0d 0a | 2d 09 0d 0a 2d 09 69 66 | ny-1;..|-...-.if|
|00003030| 28 28 6d 79 5f 69 6d 61 | 67 65 2d 3e 69 64 79 20 |((my_ima|ge->idy |
|00003040| 3d 20 28 69 6e 74 20 2a | 29 6d 61 6c 6c 6f 63 28 |= (int *|)malloc(|
|00003050| 6e 79 2a 6e 78 2a 73 69 | 7a 65 6f 66 28 69 6e 74 |ny*nx*si|zeof(int|
|00003060| 29 29 29 3d 3d 30 29 0d | 0a 2d 09 7b 0d 0a 2d 09 |)))==0).|.-.{..-.|
|00003070| 09 66 70 72 69 6e 74 66 | 28 73 74 64 65 72 72 2c |.fprintf|(stderr,|
|00003080| 20 22 65 72 72 6f 72 3a | 20 63 61 6e 6e 6f 74 20 | "error:| cannot |
|00003090| 61 6c 6c 6f 63 61 74 65 | 20 69 64 79 20 61 72 72 |allocate| idy arr|
|000030a0| 61 79 5c 6e 22 29 3b 0d | 0a 2d 09 09 65 78 69 74 |ay\n");.|.-..exit|
|000030b0| 28 31 29 3b 0d 0a 2d 09 | 7d 0d 0a 2d 09 0d 0a 2d |(1);..-.|}..-...-|
|000030c0| 09 64 50 20 3d 20 6d 79 | 5f 69 6d 61 67 65 2d 3e |.dP = my|_image->|
|000030d0| 69 64 79 3b 0d 0a 2d 09 | 67 50 20 3d 20 6d 79 5f |idy;..-.|gP = my_|
|000030e0| 69 6d 61 67 65 2d 3e 67 | 61 75 73 73 69 61 6e 3b |image->g|aussian;|
|000030f0| 0d 0a 2d 09 66 6f 72 28 | 69 78 3d 30 3b 20 69 78 |..-.for(|ix=0; ix|
|00003100| 3c 6e 78 3b 20 69 78 2b | 2b 29 0d 0a 2d 09 7b 0d |<nx; ix+|+)..-.{.|
|00003110| 0a 2d 09 09 64 50 5b 69 | 78 5d 20 3d 20 67 50 5b |.-..dP[i|x] = gP[|
|00003120| 69 78 5d 20 2d 20 67 50 | 5b 69 78 2b 6e 78 5d 3b |ix] - gP|[ix+nx];|
|00003130| 20 0d 0a 2d 2f 2a 09 09 | 2a 64 50 20 3d 20 2a 67 | ..-/*..|*dP = *g|
|00003140| 50 20 2d 20 2a 28 67 50 | 2b 6e 78 29 3b 09 2a 2f |P - *(gP|+nx);.*/|
|00003150| 0d 0a 2d 2f 2a 09 09 64 | 50 2b 2b 3b 20 67 50 2b |..-/*..d|P++; gP+|
|00003160| 2b 3b 09 2a 2f 0d 0a 2d | 09 7d 0d 0a 2d 09 66 6f |+;.*/..-|.}..-.fo|
|00003170| 72 28 69 78 3d 30 3b 20 | 69 78 3c 6e 78 3b 20 69 |r(ix=0; |ix<nx; i|
|00003180| 78 2b 2b 29 0d 0a 2d 09 | 7b 0d 0a 2d 09 09 64 50 |x++)..-.|{..-..dP|
|00003190| 20 3d 20 6d 79 5f 69 6d | 61 67 65 2d 3e 69 64 79 | = my_im|age->idy|
|000031a0| 2b 69 78 3b 0d 0a 2d 09 | 09 67 50 20 3d 20 6d 79 |+ix;..-.|.gP = my|
|000031b0| 5f 69 6d 61 67 65 2d 3e | 67 61 75 73 73 69 61 6e |_image->|gaussian|
|000031c0| 2b 69 78 3b 0d 0a 2d 09 | 09 66 6f 72 28 69 79 3d |+ix;..-.|.for(iy=|
|000031d0| 31 3b 20 69 79 3c 6e 79 | 4d 31 3b 20 69 79 2b 2b |1; iy<ny|M1; iy++|
|000031e0| 29 0d 0a 2d 09 09 7b 0d | 0a 2d 09 09 09 64 50 20 |)..-..{.|.-...dP |
|000031f0| 2b 3d 20 6e 78 3b 0d 0a | 2d 09 09 09 2a 64 50 20 |+= nx;..|-...*dP |
|00003200| 3d 20 2a 28 67 50 2b 6e | 78 50 6e 78 29 2d 20 2a |= *(gP+n|xPnx)- *|
|00003210| 67 50 3b 0d 0a 2d 09 09 | 09 67 50 20 2b 3d 20 6e |gP;..-..|.gP += n|
|00003220| 78 3b 0d 0a 2d 09 09 7d | 0d 0a 2d 09 7d 0d 0a 2d |x;..-..}|..-.}..-|
|00003230| 09 64 50 20 3d 20 6d 79 | 5f 69 6d 61 67 65 2d 3e |.dP = my|_image->|
|00003240| 69 64 79 2b 28 6e 79 2d | 31 29 2a 6e 78 3b 0d 0a |idy+(ny-|1)*nx;..|
|00003250| 2d 09 67 50 20 3d 20 6d | 79 5f 69 6d 61 67 65 2d |-.gP = m|y_image-|
|00003260| 3e 67 61 75 73 73 69 61 | 6e 2b 28 6e 79 2d 31 29 |>gaussia|n+(ny-1)|
|00003270| 2a 6e 78 3b 0d 0a 2d 09 | 66 6f 72 28 69 78 3d 30 |*nx;..-.|for(ix=0|
|00003280| 3b 20 69 78 3c 6e 78 3b | 20 69 78 2b 2b 29 0d 0a |; ix<nx;| ix++)..|
|00003290| 2d 09 7b 0d 0a 2d 09 09 | 64 50 5b 69 78 5d 20 3d |-.{..-..|dP[ix] =|
|000032a0| 20 67 50 5b 69 78 5d 20 | 2d 20 67 50 5b 69 78 2d | gP[ix] |- gP[ix-|
|000032b0| 6e 78 5d 3b 0d 0a 2d 2f | 2a 09 09 2a 64 50 20 3d |nx];..-/|*..*dP =|
|000032c0| 20 2a 67 50 20 2d 20 2a | 28 67 50 2d 6e 78 29 3b | *gP - *|(gP-nx);|
|000032d0| 09 2a 2f 0d 0a 2d 2f 2a | 09 09 64 50 2b 2b 3b 20 |.*/..-/*|..dP++; |
|000032e0| 67 50 2b 2b 3b 20 2a 2f | 0d 0a 2d 09 7d 0d 0a 2d |gP++; */|..-.}..-|
|000032f0| 23 69 66 64 65 66 20 44 | 45 42 55 47 0d 0a 2d 09 |#ifdef D|EBUG..-.|
|00003300| 69 6d 61 67 65 5f 57 72 | 69 74 65 5f 69 6e 74 28 |image_Wr|ite_int(|
|00003310| 22 64 65 74 79 22 2c 20 | 6d 79 5f 69 6d 61 67 65 |"dety", |my_image|
|00003320| 2d 3e 69 64 79 29 3b 0d | 0a 2d 23 65 6e 64 69 66 |->idy);.|.-#endif|
|00003330| 0d 0a 2d 7d 0d 0a 2d 76 | 6f 69 64 0d 0a 2d 69 6d |..-}..-v|oid..-im|
|00003340| 61 67 65 5f 64 65 74 34 | 35 28 29 0d 0a 2d 7b 0d |age_det4|5()..-{.|
|00003350| 0a 2d 09 72 65 67 69 73 | 74 65 72 20 69 6e 74 20 |.-.regis|ter int |
|00003360| 2a 64 50 2c 20 2a 67 50 | 3b 0d 0a 2d 09 72 65 67 |*dP, *gP|;..-.reg|
|00003370| 69 73 74 65 72 20 69 6e | 74 20 69 2c 20 69 78 2c |ister in|t i, ix,|
|00003380| 20 69 79 2c 20 6e 78 2c | 69 6e 64 65 78 2c 20 6e | iy, nx,|index, n|
|00003390| 79 3b 0d 0a 2d 09 72 65 | 67 69 73 74 65 72 20 6e |y;..-.re|gister n|
|000033a0| 78 50 31 2c 20 6e 78 4d | 31 2c 20 6e 79 4d 31 3b |xP1, nxM|1, nyM1;|
|000033b0| 0d 0a 2d 0d 0a 2d 09 6e | 78 20 3d 20 6d 79 5f 69 |..-..-.n|x = my_i|
|000033c0| 6d 61 67 65 2d 3e 6e 78 | 3b 09 6e 78 4d 31 20 3d |mage->nx|;.nxM1 =|
|000033d0| 20 6e 78 2d 31 3b 09 6e | 78 50 31 20 3d 20 6e 78 | nx-1;.n|xP1 = nx|
|000033e0| 2b 31 3b 0d 0a 2d 09 6e | 79 20 3d 20 6d 79 5f 69 |+1;..-.n|y = my_i|
|000033f0| 6d 61 67 65 2d 3e 6e 79 | 3b 09 6e 79 4d 31 20 3d |mage->ny|;.nyM1 =|
|00003400| 20 6e 79 2d 31 3b 0d 0a | 2d 09 0d 0a 2d 09 69 66 | ny-1;..|-...-.if|
|00003410| 28 28 6d 79 5f 69 6d 61 | 67 65 2d 3e 69 64 34 35 |((my_ima|ge->id45|
|00003420| 20 3d 20 28 69 6e 74 20 | 2a 29 6d 61 6c 6c 6f 63 | = (int |*)malloc|
|00003430| 28 6e 79 2a 6e 78 2a 73 | 69 7a 65 6f 66 28 69 6e |(ny*nx*s|izeof(in|
|00003440| 74 29 29 29 3d 3d 30 29 | 0d 0a 2d 09 7b 0d 0a 2d |t)))==0)|..-.{..-|
|00003450| 09 09 66 70 72 69 6e 74 | 66 28 73 74 64 65 72 72 |..fprint|f(stderr|
|00003460| 2c 20 22 65 72 72 6f 72 | 3a 20 63 61 6e 6e 6f 74 |, "error|: cannot|
|00003470| 20 61 6c 6c 6f 63 61 74 | 65 20 69 64 34 35 20 61 | allocat|e id45 a|
|00003480| 72 72 61 79 5c 6e 22 29 | 3b 0d 0a 2d 09 09 65 78 |rray\n")|;..-..ex|
|00003490| 69 74 28 31 29 3b 0d 0a | 2d 09 7d 0d 0a 2d 2f 2a |it(1);..|-.}..-/*|
|000034a0| 0d 0a 2d 20 2a 09 64 6f | 20 74 6f 70 0d 0a 2d 20 |..- *.do| top..- |
|000034b0| 2a 09 74 68 65 6e 20 6d | 69 64 64 6c 65 0d 0a 2d |*.then m|iddle..-|
|000034c0| 20 2a 09 74 68 65 6e 20 | 72 69 67 68 74 20 73 69 | *.then |right si|
|000034d0| 64 65 0d 0a 2d 2a 2f 0d | 0a 2d 09 64 50 20 3d 20 |de..-*/.|.-.dP = |
|000034e0| 6d 79 5f 69 6d 61 67 65 | 2d 3e 69 64 34 35 2b 6e |my_image|->id45+n|
|000034f0| 78 2d 32 3b 0d 0a 2d 09 | 67 50 20 3d 20 6d 79 5f |x-2;..-.|gP = my_|
|00003500| 69 6d 61 67 65 2d 3e 67 | 61 75 73 73 69 61 6e 2b |image->g|aussian+|
|00003510| 6e 78 2d 32 3b 0d 0a 2d | 09 66 6f 72 28 69 78 3d |nx-2;..-|.for(ix=|
|00003520| 6e 78 2d 32 3b 20 69 78 | 3e 2d 31 3b 20 69 78 2d |nx-2; ix|>-1; ix-|
|00003530| 2d 29 0d 0a 2d 09 7b 0d | 0a 2d 2f 2a 09 09 64 5b |-)..-.{.|.-/*..d[|
|00003540| 69 78 5d 20 3d 20 67 5b | 69 78 5d 20 2d 20 67 5b |ix] = g[|ix] - g[|
|00003550| 69 78 2b 6e 78 2b 31 5d | 3b 20 2a 2f 0d 0a 2d 09 |ix+nx+1]|; */..-.|
|00003560| 09 2a 64 50 20 3d 20 2a | 67 50 20 2d 20 2a 28 67 |.*dP = *|gP - *(g|
|00003570| 50 2b 6e 78 50 31 29 3b | 0d 0a 2d 09 09 64 50 2d |P+nxP1);|..-..dP-|
|00003580| 2d 3b 20 67 50 2d 2d 3b | 0d 0a 2d 09 7d 0d 0a 2d |-; gP--;|..-.}..-|
|00003590| 09 64 50 20 3d 20 6d 79 | 5f 69 6d 61 67 65 2d 3e |.dP = my|_image->|
|000035a0| 69 64 34 35 3b 0d 0a 2d | 09 67 50 20 3d 20 6d 79 |id45;..-|.gP = my|
|000035b0| 5f 69 6d 61 67 65 2d 3e | 67 61 75 73 73 69 61 6e |_image->|gaussian|
|000035c0| 3b 0d 0a 2d 09 66 6f 72 | 28 69 3d 6e 78 2d 32 3b |;..-.for|(i=nx-2;|
|000035d0| 20 69 3e 2d 31 3b 20 69 | 2d 2d 29 0d 0a 2d 09 7b | i>-1; i|--)..-.{|
|000035e0| 0d 0a 2d 09 09 69 78 3d | 69 3b 0d 0a 2d 09 09 69 |..-..ix=|i;..-..i|
|000035f0| 79 3d 31 3b 0d 0a 2d 09 | 09 69 6e 64 65 78 20 3d |y=1;..-.|.index =|
|00003600| 20 6e 78 2b 69 78 3b 0d | 0a 2d 09 09 77 68 69 6c | nx+ix;.|.-..whil|
|00003610| 65 28 69 78 3e 2d 31 20 | 26 26 20 69 79 3c 6e 79 |e(ix>-1 |&& iy<ny|
|00003620| 4d 31 29 0d 0a 2d 09 09 | 7b 0d 0a 2d 2f 2a 09 09 |M1)..-..|{..-/*..|
|00003630| 09 64 50 5b 69 79 2a 6e | 78 2b 69 78 5d 20 3d 20 |.dP[iy*n|x+ix] = |
|00003640| 67 50 5b 28 69 79 2b 31 | 29 2a 6e 78 2b 69 78 2d |gP[(iy+1|)*nx+ix-|
|00003650| 31 5d 2d 67 50 5b 28 69 | 79 2d 31 29 2a 6e 78 2b |1]-gP[(i|y-1)*nx+|
|00003660| 69 78 2b 31 5d 3b 20 2a | 2f 0d 0a 2d 2f 2a 09 09 |ix+1]; *|/..-/*..|
|00003670| 09 64 50 5b 69 79 2a 6e | 78 2b 69 78 5d 20 3d 20 |.dP[iy*n|x+ix] = |
|00003680| 67 50 5b 69 79 2a 6e 78 | 2b 69 78 2b 6e 78 2d 31 |gP[iy*nx|+ix+nx-1|
|00003690| 5d 20 2d 67 50 5b 69 79 | 2a 6e 78 2b 69 78 2d 6e |] -gP[iy|*nx+ix-n|
|000036a0| 78 2b 31 5d 3b 20 20 2a | 2f 0d 0a 2d 09 09 09 64 |x+1];  *|/..-...d|
|000036b0| 50 5b 69 6e 64 65 78 5d | 20 3d 20 67 50 5b 69 6e |P[index]| = gP[in|
|000036c0| 64 65 78 2b 6e 78 4d 31 | 5d 20 2d 20 67 50 5b 69 |dex+nxM1|] - gP[i|
|000036d0| 6e 64 65 78 2d 6e 78 4d | 31 5d 3b 0d 0a 2d 09 09 |ndex-nxM|1];..-..|
|000036e0| 09 69 6e 64 65 78 20 2b | 3d 20 6e 78 4d 31 3b 0d |.index +|= nxM1;.|
|000036f0| 0a 2d 09 09 09 69 78 2d | 2d 3b 0d 0a 2d 09 09 09 |.-...ix-|-;..-...|
|00003700| 69 79 2b 2b 3b 0d 0a 2d | 09 09 7d 0d 0a 2d 09 7d |iy++;..-|..}..-.}|
|00003710| 0d 0a 2d 09 66 6f 72 28 | 69 3d 31 3b 20 69 3c 6e |..-.for(|i=1; i<n|
|00003720| 79 3b 20 69 2b 2b 29 0d | 0a 2d 09 7b 0d 0a 2d 09 |y; i++).|.-.{..-.|
|00003730| 09 69 78 3d 6e 78 2d 32 | 3b 0d 0a 2d 09 09 69 79 |.ix=nx-2|;..-..iy|
|00003740| 3d 69 3b 0d 0a 2d 09 09 | 69 6e 64 65 78 20 3d 20 |=i;..-..|index = |
|00003750| 69 79 2a 6e 78 2b 69 78 | 3b 0d 0a 2d 09 09 77 68 |iy*nx+ix|;..-..wh|
|00003760| 69 6c 65 28 69 78 3e 2d | 31 20 26 26 20 69 79 3c |ile(ix>-|1 && iy<|
|00003770| 6e 79 4d 31 29 0d 0a 2d | 09 09 7b 0d 0a 2d 2f 2a |nyM1)..-|..{..-/*|
|00003780| 09 09 09 64 50 5b 69 79 | 2a 6e 78 2b 69 78 5d 20 |...dP[iy|*nx+ix] |
|00003790| 3d 20 67 50 5b 28 69 79 | 2b 31 29 2a 6e 78 2b 69 |= gP[(iy|+1)*nx+i|
|000037a0| 78 2d 31 5d 2d 67 50 5b | 28 69 79 2d 31 29 2a 6e |x-1]-gP[|(iy-1)*n|
|000037b0| 78 2b 69 78 2b 31 5d 3b | 20 2a 2f 0d 0a 2d 2f 2a |x+ix+1];| */..-/*|
|000037c0| 09 09 09 64 50 5b 69 79 | 2a 6e 78 2b 69 78 5d 20 |...dP[iy|*nx+ix] |
|000037d0| 3d 20 67 50 5b 69 79 2a | 6e 78 2b 69 78 2b 6e 78 |= gP[iy*|nx+ix+nx|
|000037e0| 2d 31 5d 20 2d 67 50 5b | 69 79 2a 6e 78 2b 69 78 |-1] -gP[|iy*nx+ix|
|000037f0| 2d 6e 78 2b 31 5d 3b 20 | 20 2a 2f 0d 0a 2d 09 09 |-nx+1]; | */..-..|
|00003800| 09 64 50 5b 69 6e 64 65 | 78 5d 20 3d 20 67 50 5b |.dP[inde|x] = gP[|
|00003810| 69 6e 64 65 78 2b 6e 78 | 4d 31 5d 20 2d 20 67 50 |index+nx|M1] - gP|
|00003820| 5b 69 6e 64 65 78 2d 6e | 78 4d 31 5d 3b 0d 0a 2d |[index-n|xM1];..-|
|00003830| 09 09 09 69 6e 64 65 78 | 20 2b 3d 20 6e 78 4d 31 |...index| += nxM1|
|00003840| 3b 0d 0a 2d 09 09 09 69 | 78 2d 2d 3b 0d 0a 2d 09 |;..-...i|x--;..-.|
|00003850| 09 09 69 79 2b 2b 3b 0d | 0a 2d 09 09 7d 0d 0a 2d |..iy++;.|.-..}..-|
|00003860| 09 7d 0d 0a 2d 09 64 50 | 20 3d 20 6d 79 5f 69 6d |.}..-.dP| = my_im|
|00003870| 61 67 65 2d 3e 69 64 34 | 35 2b 6e 78 2b 6e 78 2d |age->id4|5+nx+nx-|
|00003880| 31 3b 0d 0a 2d 09 67 50 | 20 3d 20 6d 79 5f 69 6d |1;..-.gP| = my_im|
|00003890| 61 67 65 2d 3e 67 61 75 | 73 73 69 61 6e 2b 6e 78 |age->gau|ssian+nx|
|000038a0| 2b 6e 78 2d 31 3b 0d 0a | 2d 09 66 6f 72 28 69 79 |+nx-1;..|-.for(iy|
|000038b0| 3d 31 3b 20 69 79 3c 6e | 79 3b 20 69 79 2b 2b 29 |=1; iy<n|y; iy++)|
|000038c0| 0d 0a 2d 09 7b 0d 0a 2d | 2f 2a 09 09 64 5b 69 79 |..-.{..-|/*..d[iy|
|000038d0| 2a 6e 78 2b 6e 78 2d 31 | 5d 20 3d 20 67 5b 69 79 |*nx+nx-1|] = g[iy|
|000038e0| 2a 6e 78 2b 6e 78 2d 31 | 5d 20 2d 20 67 5b 28 69 |*nx+nx-1|] - g[(i|
|000038f0| 79 2d 31 29 2a 6e 78 2b | 6e 78 2d 32 5d 3b 20 2a |y-1)*nx+|nx-2]; *|
|00003900| 2f 0d 0a 2d 2f 2a 09 09 | 64 5b 69 79 2a 6e 78 2b |/..-/*..|d[iy*nx+|
|00003910| 6e 78 2d 31 5d 20 3d 20 | 67 5b 69 79 2a 6e 78 2b |nx-1] = |g[iy*nx+|
|00003920| 6e 78 2d 31 5d 20 2d 20 | 67 5b 69 79 2a 6e 78 2b |nx-1] - |g[iy*nx+|
|00003930| 6e 78 2d 31 2d 6e 78 2d | 31 5d 3b 20 2a 2f 0d 0a |nx-1-nx-|1]; */..|
|00003940| 2d 09 09 2a 64 50 20 3d | 20 2a 67 50 20 2d 20 2a |-..*dP =| *gP - *|
|00003950| 28 67 50 2d 6e 78 50 31 | 29 3b 0d 0a 2d 09 09 64 |(gP-nxP1|);..-..d|
|00003960| 50 2b 3d 6e 78 3b 20 67 | 50 2b 3d 6e 78 3b 0d 0a |P+=nx; g|P+=nx;..|
|00003970| 2d 09 7d 0d 0a 2d 09 64 | 50 20 3d 20 6d 79 5f 69 |-.}..-.d|P = my_i|
|00003980| 6d 61 67 65 2d 3e 69 64 | 34 35 2b 28 6e 79 2d 31 |mage->id|45+(ny-1|
|00003990| 29 2a 6e 78 3b 0d 0a 2d | 09 67 50 20 3d 20 6d 79 |)*nx;..-|.gP = my|
|000039a0| 5f 69 6d 61 67 65 2d 3e | 67 61 75 73 73 69 61 6e |_image->|gaussian|
|000039b0| 2b 28 6e 79 2d 31 29 2a | 6e 78 3b 0d 0a 2d 09 66 |+(ny-1)*|nx;..-.f|
|000039c0| 6f 72 28 69 78 3d 30 3b | 20 69 78 3c 6e 78 4d 31 |or(ix=0;| ix<nxM1|
|000039d0| 3b 20 69 78 2b 2b 29 0d | 0a 2d 09 7b 0d 0a 2d 2f |; ix++).|.-.{..-/|
|000039e0| 2a 09 09 64 50 5b 28 6e | 79 2d 31 29 2a 6e 78 2b |*..dP[(n|y-1)*nx+|
|000039f0| 69 78 5d 20 3d 20 67 50 | 5b 28 6e 79 2d 31 29 2a |ix] = gP|[(ny-1)*|
|00003a00| 6e 78 2b 69 78 5d 20 2d | 20 67 50 5b 28 6e 79 2d |nx+ix] -| gP[(ny-|
|00003a10| 32 29 2a 6e 78 2b 69 78 | 2b 31 5d 3b 20 20 20 20 |2)*nx+ix|+1];    |
|00003a20| 2a 2f 0d 0a 2d 2f 2a 09 | 09 64 50 5b 28 6e 79 2d |*/..-/*.|.dP[(ny-|
|00003a30| 31 29 2a 6e 78 2b 69 78 | 5d 20 3d 20 67 50 5b 28 |1)*nx+ix|] = gP[(|
|00003a40| 6e 79 2d 31 29 2a 6e 78 | 2b 69 78 5d 20 2d 20 67 |ny-1)*nx|+ix] - g|
|00003a50| 50 5b 28 6e 79 2d 31 29 | 2a 6e 78 2b 69 78 2d 6e |P[(ny-1)|*nx+ix-n|
|00003a60| 78 2b 31 5d 3b 20 2a 2f | 0d 0a 2d 09 09 2a 64 50 |x+1]; */|..-..*dP|
|00003a70| 20 3d 20 2a 67 50 20 2d | 20 2a 28 67 50 2d 6e 78 | = *gP -| *(gP-nx|
|00003a80| 4d 31 29 3b 0d 0a 2d 09 | 09 64 50 2b 2b 3b 20 67 |M1);..-.|.dP++; g|
|00003a90| 50 2b 2b 3b 0d 0a 2d 09 | 7d 0d 0a 2d 09 64 50 20 |P++;..-.|}..-.dP |
|00003aa0| 3d 20 6d 79 5f 69 6d 61 | 67 65 2d 3e 69 64 34 35 |= my_ima|ge->id45|
|00003ab0| 2b 28 6e 79 2d 31 29 2a | 6e 78 3b 0d 0a 2d 09 67 |+(ny-1)*|nx;..-.g|
|00003ac0| 50 20 3d 20 6d 79 5f 69 | 6d 61 67 65 2d 3e 67 61 |P = my_i|mage->ga|
|00003ad0| 75 73 73 69 61 6e 2b 28 | 6e 79 2d 31 29 2a 6e 78 |ussian+(|ny-1)*nx|
|00003ae0| 3b 0d 0a 2d 2f 2a 09 64 | 50 5b 28 6e 79 2d 31 29 |;..-/*.d|P[(ny-1)|
|00003af0| 2a 6e 78 2b 6e 78 2d 31 | 5d 20 3d 20 67 50 5b 28 |*nx+nx-1|] = gP[(|
|00003b00| 6e 79 2d 31 29 2a 6e 78 | 2b 6e 78 2d 31 5d 20 2d |ny-1)*nx|+nx-1] -|
|00003b10| 20 67 50 5b 28 6e 79 2d | 32 29 2a 6e 78 2b 6e 78 | gP[(ny-|2)*nx+nx|
|00003b20| 2d 32 5d 3b 20 20 20 20 | 20 20 2a 2f 0d 0a 2d 2f |-2];    |  */..-/|
|00003b30| 2a 09 64 50 5b 28 6e 79 | 2d 31 29 2a 6e 78 2b 6e |*.dP[(ny|-1)*nx+n|
|00003b40| 78 2d 31 5d 20 3d 20 67 | 50 5b 28 6e 79 2d 31 29 |x-1] = g|P[(ny-1)|
|00003b50| 2a 6e 78 2b 6e 78 2d 31 | 5d 20 2d 20 67 50 5b 28 |*nx+nx-1|] - gP[(|
|00003b60| 6e 79 2d 31 29 2a 6e 78 | 2b 6e 78 2d 31 2d 6e 78 |ny-1)*nx|+nx-1-nx|
|00003b70| 2d 31 5d 3b 20 2a 2f 0d | 0a 2d 09 64 50 5b 6e 78 |-1]; */.|.-.dP[nx|
|00003b80| 4d 31 5d 20 3d 20 67 50 | 5b 6e 78 4d 31 5d 20 2d |M1] = gP|[nxM1] -|
|00003b90| 20 67 50 5b 2d 32 5d 3b | 0d 0a 2d 09 64 50 20 3d | gP[-2];|..-.dP =|
|00003ba0| 20 6d 79 5f 69 6d 61 67 | 65 2d 3e 69 64 34 35 3b | my_imag|e->id45;|
|00003bb0| 0d 0a 2d 09 67 50 20 3d | 20 6d 79 5f 69 6d 61 67 |..-.gP =| my_imag|
|00003bc0| 65 2d 3e 67 61 75 73 73 | 69 61 6e 3b 0d 0a 2d 09 |e->gauss|ian;..-.|
|00003bd0| 64 50 5b 6e 78 4d 31 5d | 20 3d 20 67 50 5b 6e 78 |dP[nxM1]| = gP[nx|
|00003be0| 4d 31 5d 2d 67 50 5b 6e | 78 2b 6e 78 2d 32 5d 3b |M1]-gP[n|x+nx-2];|
|00003bf0| 0d 0a 2d 23 69 66 64 65 | 66 20 44 45 42 55 47 0d |..-#ifde|f DEBUG.|
|00003c00| 0a 2d 09 69 6d 61 67 65 | 5f 57 72 69 74 65 5f 69 |.-.image|_Write_i|
|00003c10| 6e 74 28 22 64 65 74 34 | 35 22 2c 20 6d 79 5f 69 |nt("det4|5", my_i|
|00003c20| 6d 61 67 65 2d 3e 69 64 | 34 35 29 3b 0d 0a 2d 23 |mage->id|45);..-#|
|00003c30| 65 6e 64 69 66 0d 0a 2d | 7d 0d 0a 2d 0d 0a 2d 76 |endif..-|}..-..-v|
|00003c40| 6f 69 64 0d 0a 2d 69 6d | 61 67 65 5f 64 65 74 31 |oid..-im|age_det1|
|00003c50| 33 35 28 29 0d 0a 2d 7b | 0d 0a 2d 09 72 65 67 69 |35()..-{|..-.regi|
|00003c60| 73 74 65 72 20 69 6e 74 | 20 2a 64 50 2c 20 2a 67 |ster int| *dP, *g|
|00003c70| 50 3b 0d 0a 2d 09 72 65 | 67 69 73 74 65 72 20 69 |P;..-.re|gister i|
|00003c80| 6e 74 20 69 2c 20 69 78 | 2c 20 69 79 2c 20 6e 78 |nt i, ix|, iy, nx|
|00003c90| 2c 69 6e 64 65 78 2c 20 | 6e 79 3b 0d 0a 2d 09 72 |,index, |ny;..-.r|
|00003ca0| 65 67 69 73 74 65 72 20 | 6e 78 50 31 2c 20 6e 78 |egister |nxP1, nx|
|00003cb0| 4d 31 2c 20 6e 79 4d 31 | 2c 20 6e 78 4d 32 3b 0d |M1, nyM1|, nxM2;.|
|00003cc0| 0a 2d 0d 0a 2d 09 6e 78 | 20 3d 20 6d 79 5f 69 6d |.-..-.nx| = my_im|
|00003cd0| 61 67 65 2d 3e 6e 78 3b | 09 6e 78 4d 31 20 3d 20 |age->nx;|.nxM1 = |
|00003ce0| 6e 78 2d 31 3b 09 6e 78 | 50 31 20 3d 20 6e 78 2b |nx-1;.nx|P1 = nx+|
|00003cf0| 31 3b 09 6e 78 4d 32 20 | 3d 20 6e 78 2d 32 3b 0d |1;.nxM2 |= nx-2;.|
|00003d00| 0a 2d 09 6e 79 20 3d 20 | 6d 79 5f 69 6d 61 67 65 |.-.ny = |my_image|
|00003d10| 2d 3e 6e 79 3b 09 6e 79 | 4d 31 20 3d 20 6e 79 2d |->ny;.ny|M1 = ny-|
|00003d20| 31 3b 0d 0a 2d 09 0d 0a | 2d 09 69 66 28 28 6d 79 |1;..-...|-.if((my|
|00003d30| 5f 69 6d 61 67 65 2d 3e | 69 64 31 33 35 20 3d 20 |_image->|id135 = |
|00003d40| 28 69 6e 74 20 2a 29 6d | 61 6c 6c 6f 63 28 6e 79 |(int *)m|alloc(ny|
|00003d50| 2a 6e 78 2a 73 69 7a 65 | 6f 66 28 69 6e 74 29 29 |*nx*size|of(int))|
|00003d60| 29 3d 3d 30 29 0d 0a 2d | 09 7b 0d 0a 2d 09 09 66 |)==0)..-|.{..-..f|
|00003d70| 70 72 69 6e 74 66 28 73 | 74 64 65 72 72 2c 20 22 |printf(s|tderr, "|
|00003d80| 65 72 72 6f 72 3a 20 63 | 61 6e 6e 6f 74 20 61 6c |error: c|annot al|
|00003d90| 6c 6f 63 61 74 65 20 69 | 64 31 33 35 20 61 72 72 |locate i|d135 arr|
|00003da0| 61 79 5c 6e 22 29 3b 0d | 0a 2d 09 09 65 78 69 74 |ay\n");.|.-..exit|
|00003db0| 28 31 29 3b 0d 0a 2d 09 | 7d 0d 0a 2d 2f 2a 0d 0a |(1);..-.|}..-/*..|
|00003dc0| 2d 20 2a 09 64 6f 20 74 | 6f 70 0d 0a 2d 20 2a 09 |- *.do t|op..- *.|
|00003dd0| 74 68 65 6e 20 6d 69 64 | 64 6c 65 0d 0a 2d 20 2a |then mid|dle..- *|
|00003de0| 09 74 68 65 6e 20 72 69 | 67 68 74 20 73 69 64 65 |.then ri|ght side|
|00003df0| 0d 0a 2d 2a 2f 0d 0a 2d | 09 64 50 20 3d 20 6d 79 |..-*/..-|.dP = my|
|00003e00| 5f 69 6d 61 67 65 2d 3e | 69 64 31 33 35 2b 31 3b |_image->|id135+1;|
|00003e10| 0d 0a 2d 09 67 50 20 3d | 20 6d 79 5f 69 6d 61 67 |..-.gP =| my_imag|
|00003e20| 65 2d 3e 67 61 75 73 73 | 69 61 6e 2b 31 3b 0d 0a |e->gauss|ian+1;..|
|00003e30| 2d 09 66 6f 72 28 69 78 | 3d 31 3b 20 69 78 3c 6e |-.for(ix|=1; ix<n|
|00003e40| 78 3b 20 69 78 2b 2b 29 | 0d 0a 2d 09 7b 0d 0a 2d |x; ix++)|..-.{..-|
|00003e50| 2f 2a 09 09 64 5b 69 78 | 5d 20 3d 20 67 5b 69 78 |/*..d[ix|] = g[ix|
|00003e60| 5d 20 2d 20 67 5b 69 78 | 2b 6e 78 2d 31 5d 3b 20 |] - g[ix|+nx-1]; |
|00003e70| 2a 2f 0d 0a 2d 09 09 2a | 64 50 20 3d 20 2a 67 50 |*/..-..*|dP = *gP|
|00003e80| 20 2d 20 2a 28 67 50 2b | 6e 78 4d 31 29 3b 0d 0a | - *(gP+|nxM1);..|
|00003e90| 2d 09 09 64 50 2b 2b 3b | 20 67 50 2b 2b 3b 0d 0a |-..dP++;| gP++;..|
|00003ea0| 2d 09 7d 0d 0a 2d 09 64 | 50 20 3d 20 6d 79 5f 69 |-.}..-.d|P = my_i|
|00003eb0| 6d 61 67 65 2d 3e 69 64 | 31 33 35 3b 0d 0a 2d 09 |mage->id|135;..-.|
|00003ec0| 67 50 20 3d 20 6d 79 5f | 69 6d 61 67 65 2d 3e 67 |gP = my_|image->g|
|00003ed0| 61 75 73 73 69 61 6e 3b | 0d 0a 2d 09 66 6f 72 28 |aussian;|..-.for(|
|00003ee0| 69 3d 31 3b 20 69 3c 6e | 78 4d 31 3b 20 69 2b 2b |i=1; i<n|xM1; i++|
|00003ef0| 29 0d 0a 2d 09 7b 0d 0a | 2d 09 09 69 78 3d 69 3b |)..-.{..|-..ix=i;|
|00003f00| 0d 0a 2d 09 09 69 79 3d | 31 3b 0d 0a 2d 09 09 69 |..-..iy=|1;..-..i|
|00003f10| 6e 64 65 78 20 3d 20 6e | 78 2b 69 78 3b 0d 0a 2d |ndex = n|x+ix;..-|
|00003f20| 09 09 77 68 69 6c 65 28 | 69 78 3c 6e 78 4d 31 20 |..while(|ix<nxM1 |
|00003f30| 26 26 20 69 79 3c 6e 79 | 4d 31 29 0d 0a 2d 09 09 |&& iy<ny|M1)..-..|
|00003f40| 7b 0d 0a 2d 2f 2a 09 09 | 09 64 50 5b 69 79 2a 6e |{..-/*..|.dP[iy*n|
|00003f50| 78 2b 69 78 5d 20 3d 20 | 67 50 5b 28 69 79 2b 31 |x+ix] = |gP[(iy+1|
|00003f60| 29 2a 6e 78 2b 69 78 2b | 31 5d 2d 67 50 5b 28 69 |)*nx+ix+|1]-gP[(i|
|00003f70| 79 2d 31 29 2a 6e 78 2b | 69 78 2d 31 5d 3b 2a 2f |y-1)*nx+|ix-1];*/|
|00003f80| 0d 0a 2d 2f 2a 09 09 09 | 64 50 5b 69 79 2a 6e 78 |..-/*...|dP[iy*nx|
|00003f90| 2b 69 78 5d 20 3d 20 67 | 50 5b 69 79 2a 6e 78 2b |+ix] = g|P[iy*nx+|
|00003fa0| 69 78 2b 6e 78 2b 31 5d | 2d 67 50 5b 69 79 2a 6e |ix+nx+1]|-gP[iy*n|
|00003fb0| 78 2b 69 78 2d 6e 78 2d | 31 5d 3b 20 20 2a 2f 0d |x+ix-nx-|1];  */.|
|00003fc0| 0a 2d 09 09 09 64 50 5b | 69 6e 64 65 78 5d 20 3d |.-...dP[|index] =|
|00003fd0| 20 67 50 5b 69 6e 64 65 | 78 2b 6e 78 50 31 5d 20 | gP[inde|x+nxP1] |
|00003fe0| 2d 20 67 50 5b 69 6e 64 | 65 78 2d 6e 78 50 31 5d |- gP[ind|ex-nxP1]|
|00003ff0| 3b 0d 0a 2d 09 09 09 69 | 6e 64 65 78 20 2b 3d 20 |;..-...i|ndex += |
|00004000| 6e 78 50 31 3b 0d 0a 2d | 09 09 09 69 78 2b 2b 3b |nxP1;..-|...ix++;|
|00004010| 0d 0a 2d 09 09 09 69 79 | 2b 2b 3b 0d 0a 2d 09 09 |..-...iy|++;..-..|
|00004020| 7d 0d 0a 2d 09 7d 0d 0a | 2d 09 66 6f 72 28 69 3d |}..-.}..|-.for(i=|
|00004030| 31 3b 20 69 3c 6e 79 3b | 20 69 2b 2b 29 0d 0a 2d |1; i<ny;| i++)..-|
|00004040| 09 7b 0d 0a 2d 09 09 69 | 78 3d 31 3b 0d 0a 2d 09 |.{..-..i|x=1;..-.|
|00004050| 09 69 79 3d 69 3b 0d 0a | 2d 09 09 69 6e 64 65 78 |.iy=i;..|-..index|
|00004060| 20 3d 20 69 79 2a 6e 78 | 2b 69 78 3b 0d 0a 2d 09 | = iy*nx|+ix;..-.|
|00004070| 09 77 68 69 6c 65 28 69 | 78 3c 6e 78 4d 31 20 26 |.while(i|x<nxM1 &|
|00004080| 26 20 69 79 3c 6e 79 4d | 31 29 0d 0a 2d 09 09 7b |& iy<nyM|1)..-..{|
|00004090| 0d 0a 2d 2f 2a 09 09 09 | 64 50 5b 69 79 2a 6e 78 |..-/*...|dP[iy*nx|
|000040a0| 2b 69 78 5d 20 3d 20 67 | 50 5b 28 69 79 2b 31 29 |+ix] = g|P[(iy+1)|
|000040b0| 2a 6e 78 2b 69 78 2b 31 | 5d 2d 20 67 50 5b 28 69 |*nx+ix+1|]- gP[(i|
|000040c0| 79 2d 31 29 2a 6e 78 2b | 69 78 2d 31 5d 3b 2a 2f |y-1)*nx+|ix-1];*/|
|000040d0| 0d 0a 2d 2f 2a 09 09 09 | 64 50 5b 69 79 2a 6e 78 |..-/*...|dP[iy*nx|
|000040e0| 2b 69 78 5d 20 3d 20 67 | 50 5b 69 79 2a 6e 78 2b |+ix] = g|P[iy*nx+|
|000040f0| 69 78 2b 6e 78 2b 31 5d | 20 2d 20 67 50 5b 69 79 |ix+nx+1]| - gP[iy|
|00004100| 2a 6e 78 2b 69 78 2d 6e | 78 2d 31 5d 3b 20 2a 2f |*nx+ix-n|x-1]; */|
|00004110| 0d 0a 2d 09 09 09 64 50 | 5b 69 6e 64 65 78 5d 20 |..-...dP|[index] |
|00004120| 3d 20 67 50 5b 69 6e 64 | 65 78 2b 6e 78 50 31 5d |= gP[ind|ex+nxP1]|
|00004130| 20 2d 20 67 50 5b 69 6e | 64 65 78 2d 6e 78 50 31 | - gP[in|dex-nxP1|
|00004140| 5d 3b 0d 0a 2d 09 09 09 | 69 6e 64 65 78 20 2b 3d |];..-...|index +=|
|00004150| 20 6e 78 50 31 3b 0d 0a | 2d 09 09 09 69 78 2b 2b | nxP1;..|-...ix++|
|00004160| 3b 0d 0a 2d 09 09 09 69 | 79 2b 2b 3b 0d 0a 2d 09 |;..-...i|y++;..-.|
|00004170| 09 7d 0d 0a 2d 09 7d 0d | 0a 2d 09 64 50 20 3d 20 |.}..-.}.|.-.dP = |
|00004180| 6d 79 5f 69 6d 61 67 65 | 2d 3e 69 64 31 33 35 2b |my_image|->id135+|
|00004190| 6e 78 3b 0d 0a 2d 09 67 | 50 20 3d 20 6d 79 5f 69 |nx;..-.g|P = my_i|
|000041a0| 6d 61 67 65 2d 3e 67 61 | 75 73 73 69 61 6e 2b 6e |mage->ga|ussian+n|
|000041b0| 78 3b 0d 0a 2d 09 66 6f | 72 28 69 79 3d 31 3b 20 |x;..-.fo|r(iy=1; |
|000041c0| 69 79 3c 6e 79 3b 20 69 | 79 2b 2b 29 0d 0a 2d 09 |iy<ny; i|y++)..-.|
|000041d0| 7b 0d 0a 2d 2f 2a 09 09 | 64 5b 69 79 2a 6e 78 5d |{..-/*..|d[iy*nx]|
|000041e0| 20 3d 20 67 5b 69 79 2a | 6e 78 5d 20 2d 20 67 5b | = g[iy*|nx] - g[|
|000041f0| 28 69 79 2d 31 29 2a 6e | 78 2b 32 5d 3b 20 2a 2f |(iy-1)*n|x+2]; */|
|00004200| 0d 0a 2d 2f 2a 09 09 64 | 5b 69 79 2a 6e 78 5d 20 |..-/*..d|[iy*nx] |
|00004210| 3d 20 67 5b 69 79 2a 6e | 78 5d 20 2d 20 67 5b 69 |= g[iy*n|x] - g[i|
|00004220| 79 2a 6e 78 2d 6e 78 2b | 32 5d 3b 20 2a 2f 0d 0a |y*nx-nx+|2]; */..|
|00004230| 2d 09 09 2a 64 50 20 3d | 20 2a 67 50 20 2d 20 2a |-..*dP =| *gP - *|
|00004240| 28 67 50 2d 6e 78 4d 32 | 29 3b 0d 0a 2d 09 09 64 |(gP-nxM2|);..-..d|
|00004250| 50 2b 3d 6e 78 3b 20 67 | 50 2b 3d 6e 78 3b 0d 0a |P+=nx; g|P+=nx;..|
|00004260| 2d 09 7d 0d 0a 2d 09 64 | 50 20 3d 20 6d 79 5f 69 |-.}..-.d|P = my_i|
|00004270| 6d 61 67 65 2d 3e 69 64 | 31 33 35 2b 6e 78 2b 6e |mage->id|135+nx+n|
|00004280| 78 2d 31 3b 0d 0a 2d 09 | 67 50 20 3d 20 6d 79 5f |x-1;..-.|gP = my_|
|00004290| 69 6d 61 67 65 2d 3e 67 | 61 75 73 73 69 61 6e 2b |image->g|aussian+|
|000042a0| 6e 78 2b 6e 78 2d 31 3b | 0d 0a 2d 09 66 6f 72 28 |nx+nx-1;|..-.for(|
|000042b0| 69 79 3d 31 3b 20 69 79 | 3c 6e 79 3b 20 69 79 2b |iy=1; iy|<ny; iy+|
|000042c0| 2b 29 0d 0a 2d 09 7b 0d | 0a 2d 2f 2a 09 09 64 5b |+)..-.{.|.-/*..d[|
|000042d0| 69 79 2a 6e 78 2b 6e 78 | 2d 31 5d 20 3d 20 67 5b |iy*nx+nx|-1] = g[|
|000042e0| 69 79 2a 6e 78 2b 6e 78 | 2d 31 5d 20 2d 20 67 5b |iy*nx+nx|-1] - g[|
|000042f0| 28 69 79 2d 31 29 2a 6e | 78 2b 6e 78 2d 32 5d 3b |(iy-1)*n|x+nx-2];|
|00004300| 20 20 2a 2f 0d 0a 2d 2f | 2a 09 09 64 5b 69 79 2a |  */..-/|*..d[iy*|
|00004310| 6e 78 2b 6e 78 2d 31 5d | 20 3d 20 67 5b 69 79 2a |nx+nx-1]| = g[iy*|
|00004320| 6e 78 2b 6e 78 2d 31 5d | 20 2d 20 67 5b 69 79 2a |nx+nx-1]| - g[iy*|
|00004330| 6e 78 2b 6e 78 2d 31 2d | 6e 78 2d 31 5d 3b 20 2a |nx+nx-1-|nx-1]; *|
|00004340| 2f 0d 0a 2d 09 09 2a 64 | 50 20 3d 20 2a 67 50 20 |/..-..*d|P = *gP |
|00004350| 2d 20 2a 28 67 50 2d 6e | 78 50 31 29 3b 0d 0a 2d |- *(gP-n|xP1);..-|
|00004360| 09 09 64 50 2b 3d 6e 78 | 3b 20 67 50 2b 3d 6e 78 |..dP+=nx|; gP+=nx|
|00004370| 3b 0d 0a 2d 09 7d 0d 0a | 2d 0d 0a 2d 09 64 50 20 |;..-.}..|-..-.dP |
|00004380| 3d 20 6d 79 5f 69 6d 61 | 67 65 2d 3e 69 64 31 33 |= my_ima|ge->id13|
|00004390| 35 2b 28 6e 79 2d 31 29 | 2a 6e 78 2b 31 3b 0d 0a |5+(ny-1)|*nx+1;..|
|000043a0| 2d 09 67 50 20 3d 20 6d | 79 5f 69 6d 61 67 65 2d |-.gP = m|y_image-|
|000043b0| 3e 67 61 75 73 73 69 61 | 6e 2b 28 6e 79 2d 31 29 |>gaussia|n+(ny-1)|
|000043c0| 2a 6e 78 2b 31 3b 0d 0a | 2d 09 66 6f 72 28 69 78 |*nx+1;..|-.for(ix|
|000043d0| 3d 31 3b 20 69 78 3c 6e | 78 4d 31 3b 20 69 78 2b |=1; ix<n|xM1; ix+|
|000043e0| 2b 29 0d 0a 2d 09 7b 0d | 0a 2d 2f 2a 09 09 64 5b |+)..-.{.|.-/*..d[|
|000043f0| 28 6e 79 2d 31 29 2a 6e | 78 2b 69 78 5d 20 3d 20 |(ny-1)*n|x+ix] = |
|00004400| 67 5b 28 6e 79 2d 31 29 | 2a 6e 78 2b 69 78 5d 20 |g[(ny-1)|*nx+ix] |
|00004410| 2d 20 67 5b 28 6e 79 2d | 32 29 2a 6e 78 2b 69 78 |- g[(ny-|2)*nx+ix|
|00004420| 2d 31 5d 3b 20 20 20 20 | 2a 2f 0d 0a 2d 2f 2a 09 |-1];    |*/..-/*.|
|00004430| 09 64 5b 28 6e 79 2d 31 | 29 2a 6e 78 2b 69 78 5d |.d[(ny-1|)*nx+ix]|
|00004440| 20 3d 20 67 5b 28 6e 79 | 2d 31 29 2a 6e 78 2b 69 | = g[(ny|-1)*nx+i|
|00004450| 78 5d 20 2d 20 67 5b 28 | 6e 79 2d 31 29 2a 6e 78 |x] - g[(|ny-1)*nx|
|00004460| 2b 69 78 2d 6e 78 2d 31 | 5d 3b 20 2a 2f 0d 0a 2d |+ix-nx-1|]; */..-|
|00004470| 09 09 2a 64 50 20 3d 20 | 2a 67 50 20 2d 20 2a 28 |..*dP = |*gP - *(|
|00004480| 67 50 2d 6e 78 50 31 29 | 3b 0d 0a 2d 09 09 64 50 |gP-nxP1)|;..-..dP|
|00004490| 2b 2b 3b 20 67 50 2b 2b | 3b 0d 0a 2d 09 7d 0d 0a |++; gP++|;..-.}..|
|000044a0| 2d 09 64 50 20 3d 20 6d | 79 5f 69 6d 61 67 65 2d |-.dP = m|y_image-|
|000044b0| 3e 69 64 31 33 35 3b 0d | 0a 2d 09 67 50 20 3d 20 |>id135;.|.-.gP = |
|000044c0| 6d 79 5f 69 6d 61 67 65 | 2d 3e 67 61 75 73 73 69 |my_image|->gaussi|
|000044d0| 61 6e 3b 0d 0a 2d 09 64 | 50 5b 6e 78 4d 31 5d 20 |an;..-.d|P[nxM1] |
|000044e0| 3d 20 67 50 5b 6e 78 4d | 31 5d 20 2d 20 67 50 5b |= gP[nxM|1] - gP[|
|000044f0| 6e 78 2b 6e 78 2d 32 5d | 3b 0d 0a 2d 09 64 50 5b |nx+nx-2]|;..-.dP[|
|00004500| 30 5d 20 3d 20 67 50 5b | 30 5d 20 2d 20 67 50 5b |0] = gP[|0] - gP[|
|00004510| 6e 78 2b 31 5d 3b 0d 0a | 2d 23 69 66 64 65 66 20 |nx+1];..|-#ifdef |
|00004520| 44 45 42 55 47 0d 0a 2d | 09 69 6d 61 67 65 5f 57 |DEBUG..-|.image_W|
|00004530| 72 69 74 65 5f 69 6e 74 | 28 22 64 65 74 31 33 35 |rite_int|("det135|
|00004540| 22 2c 20 6d 79 5f 69 6d | 61 67 65 2d 3e 69 64 31 |", my_im|age->id1|
|00004550| 33 35 29 3b 0d 0a 2d 23 | 65 6e 64 69 66 0d 0a 2d |35);..-#|endif..-|
|00004560| 7d 0d 0a 2f 2f 47 4f 2e | 53 59 53 49 4e 20 44 44 |}..//GO.|SYSIN DD|
|00004570| 20 64 65 74 2e 63 0d 0a | 65 63 68 6f 20 64 74 6f | det.c..|echo dto|
|00004580| 72 2e 63 20 31 3e 26 32 | 0d 0a 73 65 64 20 27 73 |r.c 1>&2|..sed 's|
|00004590| 2f 2e 2f 2f 27 20 3e 64 | 74 6f 72 2e 63 20 3c 3c |/.//' >d|tor.c <<|
|000045a0| 27 2f 2f 47 4f 2e 53 59 | 53 49 4e 20 44 44 20 64 |'//GO.SY|SIN DD d|
|000045b0| 74 6f 72 2e 63 27 0d 0a | 2d 2f 2a 0d 0a 2d 20 2a |tor.c'..|-/*..- *|
|000045c0| 20 28 63 29 20 43 6f 70 | 79 72 69 67 68 74 20 31 | (c) Cop|yright 1|
|000045d0| 39 38 38 20 62 79 0d 0a | 2d 20 2a 20 52 6f 62 6f |988 by..|- * Robo|
|000045e0| 74 69 63 73 20 50 72 69 | 6e 63 69 70 6c 65 73 20 |tics Pri|nciples |
|000045f0| 52 65 73 65 61 72 63 68 | 20 44 65 70 61 72 74 6d |Research| Departm|
|00004600| 65 6e 74 2c 20 41 54 54 | 20 42 65 6c 6c 20 4c 61 |ent, ATT| Bell La|
|00004610| 62 6f 72 61 74 6f 72 69 | 65 73 2e 0d 0a 2d 20 2a |boratori|es...- *|
|00004620| 20 41 6c 6c 20 72 69 67 | 68 74 73 20 72 65 73 65 | All rig|hts rese|
|00004630| 72 76 65 64 2e 0d 0a 2d | 20 2a 20 4c 61 73 74 20 |rved...-| * Last |
|00004640| 6d 6f 64 69 66 69 65 64 | 20 32 2f 38 2f 38 38 20 |modified| 2/8/88 |
|00004650| 49 6e 67 65 6d 61 72 20 | 4a 2e 20 43 6f 78 0d 0a |Ingemar |J. Cox..|
|00004660| 2d 20 2a 20 43 20 76 65 | 72 73 69 6f 6e 20 38 2f |- * C ve|rsion 8/|
|00004670| 32 2f 38 38 20 44 65 62 | 6f 72 61 68 20 41 2e 20 |2/88 Deb|orah A. |
|00004680| 57 61 6c 6c 61 63 68 0d | 0a 2d 20 2a 2f 0d 0a 2d |Wallach.|.- */..-|
|00004690| 23 69 6e 63 6c 75 64 65 | 20 3c 73 74 64 69 6f 2e |#include| <stdio.|
|000046a0| 68 3e 0d 0a 2d 23 69 6e | 63 6c 75 64 65 20 3c 6d |h>..-#in|clude <m|
|000046b0| 61 74 68 2e 68 3e 0d 0a | 2d 23 69 6e 63 6c 75 64 |ath.h>..|-#includ|
|000046c0| 65 20 22 65 64 67 65 5f | 66 69 6e 64 65 72 2e 68 |e "edge_|finder.h|
|000046d0| 22 0d 0a 2d 0d 0a 2d 65 | 78 74 65 72 6e 20 73 74 |"..-..-e|xtern st|
|000046e0| 72 75 63 74 20 69 6d 61 | 67 65 20 2a 6d 79 5f 69 |ruct ima|ge *my_i|
|000046f0| 6d 61 67 65 3b 0d 0a 2d | 0d 0a 2d 69 6d 61 67 65 |mage;..-|..-image|
|00004700| 5f 63 6c 65 61 6e 75 70 | 28 29 0d 0a 2d 7b 0d 0a |_cleanup|()..-{..|
|00004710| 2d 20 20 20 20 66 72 65 | 65 28 6d 79 5f 69 6d 61 |-    fre|e(my_ima|
|00004720| 67 65 2d 3e 67 61 75 73 | 73 69 61 6e 29 3b 0d 0a |ge->gaus|sian);..|
|00004730| 2d 20 20 20 20 66 72 65 | 65 28 6d 79 5f 69 6d 61 |-    fre|e(my_ima|
|00004740| 67 65 2d 3e 66 69 6c 74 | 65 72 29 3b 0d 0a 2d 20 |ge->filt|er);..- |
|00004750| 20 20 20 66 72 65 65 28 | 6d 79 5f 69 6d 61 67 65 |   free(|my_image|
|00004760| 2d 3e 69 64 78 29 3b 0d | 0a 2d 20 20 20 20 66 72 |->idx);.|.-    fr|
|00004770| 65 65 28 6d 79 5f 69 6d | 61 67 65 2d 3e 69 64 79 |ee(my_im|age->idy|
|00004780| 29 3b 0d 0a 2d 20 20 20 | 20 66 72 65 65 28 6d 79 |);..-   | free(my|
|00004790| 5f 69 6d 61 67 65 2d 3e | 69 64 34 35 29 3b 0d 0a |_image->|id45);..|
|000047a0| 2d 20 20 20 20 66 72 65 | 65 28 6d 79 5f 69 6d 61 |-    fre|e(my_ima|
|000047b0| 67 65 2d 3e 69 64 31 33 | 35 29 3b 0d 0a 2d 20 20 |ge->id13|5);..-  |
|000047c0| 20 20 66 72 65 65 28 6d | 79 5f 69 6d 61 67 65 2d |  free(m|y_image-|
|000047d0| 3e 65 64 67 65 5f 6d 61 | 70 5f 68 69 29 3b 0d 0a |>edge_ma|p_hi);..|
|000047e0| 2d 7d 0d 0a 2d 0d 0a 2f | 2f 47 4f 2e 53 59 53 49 |-}..-../|/GO.SYSI|
|000047f0| 4e 20 44 44 20 64 74 6f | 72 2e 63 0d 0a 65 63 68 |N DD dto|r.c..ech|
|00004800| 6f 20 65 64 67 65 5f 66 | 69 6e 64 65 72 2e 68 20 |o edge_f|inder.h |
|00004810| 31 3e 26 32 0d 0a 73 65 | 64 20 27 73 2f 2e 2f 2f |1>&2..se|d 's/.//|
|00004820| 27 20 3e 65 64 67 65 5f | 66 69 6e 64 65 72 2e 68 |' >edge_|finder.h|
|00004830| 20 3c 3c 27 2f 2f 47 4f | 2e 53 59 53 49 4e 20 44 | <<'//GO|.SYSIN D|
|00004840| 44 20 65 64 67 65 5f 66 | 69 6e 64 65 72 2e 68 27 |D edge_f|inder.h'|
|00004850| 0d 0a 2d 2f 2a 0d 0a 2d | 20 2a 20 28 63 29 20 43 |..-/*..-| * (c) C|
|00004860| 6f 70 79 72 69 67 68 74 | 20 31 39 38 38 20 62 79 |opyright| 1988 by|
|00004870| 0d 0a 2d 20 2a 20 52 6f | 62 6f 74 69 63 73 20 50 |..- * Ro|botics P|
|00004880| 72 69 6e 63 69 70 6c 65 | 73 20 52 65 73 65 61 72 |rinciple|s Resear|
|00004890| 63 68 20 44 65 70 61 72 | 74 6d 65 6e 74 2c 20 41 |ch Depar|tment, A|
|000048a0| 54 54 20 42 65 6c 6c 20 | 4c 61 62 6f 72 61 74 6f |TT Bell |Laborato|
|000048b0| 72 69 65 73 2e 0d 0a 2d | 20 2a 20 41 6c 6c 20 72 |ries...-| * All r|
|000048c0| 69 67 68 74 73 20 72 65 | 73 65 72 76 65 64 2e 0d |ights re|served..|
|000048d0| 0a 2d 20 2a 20 4c 61 73 | 74 20 6d 6f 64 69 66 69 |.- * Las|t modifi|
|000048e0| 65 64 20 32 2f 38 2f 38 | 38 20 49 6e 67 65 6d 61 |ed 2/8/8|8 Ingema|
|000048f0| 72 20 4a 2e 20 43 6f 78 | 0d 0a 2d 20 2a 20 43 20 |r J. Cox|..- * C |
|00004900| 76 65 72 73 69 6f 6e 20 | 38 2f 32 2f 38 38 20 44 |version |8/2/88 D|
|00004910| 65 62 6f 72 61 68 20 41 | 2e 20 57 61 6c 6c 61 63 |eborah A|. Wallac|
|00004920| 68 0d 0a 2d 20 2a 2f 0d | 0a 2d 0d 0a 2d 23 64 65 |h..- */.|.-..-#de|
|00004930| 66 69 6e 65 20 61 62 73 | 28 78 29 20 28 28 78 29 |fine abs|(x) ((x)|
|00004940| 3c 30 20 3f 20 2d 28 78 | 29 20 3a 20 28 78 29 29 |<0 ? -(x|) : (x))|
|00004950| 0d 0a 2d 23 64 65 66 69 | 6e 65 20 6d 61 78 28 78 |..-#defi|ne max(x|
|00004960| 2c 20 79 29 20 28 28 78 | 29 3e 28 79 29 20 3f 20 |, y) ((x|)>(y) ? |
|00004970| 28 78 29 20 3a 20 28 79 | 29 29 0d 0a 2d 23 64 65 |(x) : (y|))..-#de|
|00004980| 66 69 6e 65 20 6d 69 6e | 28 78 2c 20 79 29 20 28 |fine min|(x, y) (|
|00004990| 28 78 29 3c 28 79 29 20 | 3f 20 28 78 29 20 3a 20 |(x)<(y) |? (x) : |
|000049a0| 28 79 29 29 0d 0a 2d 0d | 0a 2d 23 64 65 66 69 6e |(y))..-.|.-#defin|
|000049b0| 65 20 45 44 47 45 5f 58 | 20 31 0d 0a 2d 23 64 65 |e EDGE_X| 1..-#de|
|000049c0| 66 69 6e 65 20 45 44 47 | 45 5f 59 20 32 0d 0a 2d |fine EDG|E_Y 2..-|
|000049d0| 23 64 65 66 69 6e 65 20 | 45 44 47 45 5f 34 35 20 |#define |EDGE_45 |
|000049e0| 33 0d 0a 2d 23 64 65 66 | 69 6e 65 20 45 44 47 45 |3..-#def|ine EDGE|
|000049f0| 5f 31 33 35 20 34 0d 0a | 2d 0d 0a 2d 23 64 65 66 |_135 4..|-..-#def|
|00004a00| 69 6e 65 20 46 57 44 20 | 30 0d 0a 2d 23 64 65 66 |ine FWD |0..-#def|
|00004a10| 69 6e 65 20 42 57 44 20 | 31 0d 0a 2d 0d 0a 2d 23 |ine BWD |1..-..-#|
|00004a20| 64 65 66 69 6e 65 20 64 | 69 61 67 5f 73 63 61 6c |define d|iag_scal|
|00004a30| 65 20 31 2e 34 31 34 20 | 2f 2a 20 73 63 61 6c 65 |e 1.414 |/* scale|
|00004a40| 20 66 61 63 74 6f 72 20 | 66 6f 72 20 64 69 61 67 | factor |for diag|
|00004a50| 6f 6e 61 6c 73 20 2a 2f | 0d 0a 2d 0d 0a 2d 73 74 |onals */|..-..-st|
|00004a60| 72 75 63 74 20 69 6d 61 | 67 65 0d 0a 2d 7b 0d 0a |ruct ima|ge..-{..|
|00004a70| 2d 09 75 6e 73 69 67 6e | 65 64 20 63 68 61 72 20 |-.unsign|ed char |
|00004a80| 2a 61 6e 5f 69 6d 61 67 | 65 3b 20 2f 2a 20 70 6f |*an_imag|e; /* po|
|00004a90| 69 6e 74 65 72 20 74 6f | 20 69 6d 61 67 65 20 61 |inter to| image a|
|00004aa0| 72 72 61 79 20 2a 2f 0d | 0a 2d 09 69 6e 74 20 6e |rray */.|.-.int n|
|00004ab0| 78 2c 20 6e 79 3b 20 2f | 2a 20 64 69 6d 65 6e 73 |x, ny; /|* dimens|
|00004ac0| 69 6f 6e 73 20 6f 66 20 | 69 6d 61 67 65 20 61 72 |ions of |image ar|
|00004ad0| 72 61 79 20 2a 2f 0d 0a | 2d 0d 0a 2d 09 69 6e 74 |ray */..|-..-.int|
|00004ae0| 2a 09 67 61 75 73 73 69 | 61 6e 3b 20 2f 2a 20 70 |*.gaussi|an; /* p|
|00004af0| 6f 69 6e 74 65 72 20 74 | 6f 20 6f 75 74 70 75 74 |ointer t|o output|
|00004b00| 20 69 6d 61 67 65 20 61 | 66 74 65 72 20 63 6f 6e | image a|fter con|
|00004b10| 76 6f 6c 75 74 69 6f 6e | 20 2a 2f 0d 0a 2d 09 69 |volution| */..-.i|
|00004b20| 6e 74 20 2a 69 64 78 2c | 20 2a 69 64 79 2c 20 2a |nt *idx,| *idy, *|
|00004b30| 69 64 34 35 2c 20 2a 69 | 64 31 33 35 3b 20 2f 2a |id45, *i|d135; /*|
|00004b40| 20 70 6f 69 6e 74 65 72 | 73 20 74 6f 20 6f 2f 70 | pointer|s to o/p|
|00004b50| 20 66 72 6f 6d 20 64 65 | 74 6e 20 66 69 6c 74 65 | from de|tn filte|
|00004b60| 72 73 20 2a 2f 0d 0a 2d | 0d 0a 2d 09 69 6e 74 20 |rs */..-|..-.int |
|00004b70| 2a 66 69 6c 74 65 72 3b | 20 2f 2a 20 70 6f 69 6e |*filter;| /* poin|
|00004b80| 74 65 72 20 74 6f 20 31 | 44 20 61 72 72 61 79 20 |ter to 1|D array |
|00004b90| 6f 66 20 66 69 6c 74 65 | 72 20 63 6f 65 66 66 69 |of filte|r coeffi|
|00004ba0| 63 69 65 6e 74 73 20 2a | 2f 0d 0a 2d 09 69 6e 74 |cients *|/..-.int|
|00004bb0| 20 6e 6f 72 6d 3b 20 2f | 2a 20 6e 6f 72 6d 61 6c | norm; /|* normal|
|00004bc0| 69 7a 61 74 69 6f 6e 20 | 66 61 63 74 6f 72 20 61 |ization |factor a|
|00004bd0| 66 74 65 72 20 63 6f 6e | 76 6f 6c 75 74 69 6f 6e |fter con|volution|
|00004be0| 20 2a 2f 0d 0a 2d 09 69 | 6e 74 20 6e 66 2c 20 6e | */..-.i|nt nf, n|
|00004bf0| 66 32 3b 20 2f 2a 20 6e | 66 20 3d 20 73 69 7a 65 |f2; /* n|f = size|
|00004c00| 20 6f 66 20 66 69 6c 74 | 65 72 20 28 69 6e 20 70 | of filt|er (in p|
|00004c10| 69 78 65 6c 73 29 2c 20 | 6e 66 32 3d 6e 66 2f 32 |ixels), |nf2=nf/2|
|00004c20| 20 2a 2f 0d 0a 2d 09 69 | 6e 74 20 6f 72 69 65 6e | */..-.i|nt orien|
|00004c30| 74 5f 66 6c 61 67 3b 0d | 0a 2d 09 75 6e 73 69 67 |t_flag;.|.-.unsig|
|00004c40| 6e 65 64 20 63 68 61 72 | 20 2a 65 64 67 65 5f 6d |ned char| *edge_m|
|00004c50| 61 70 2c 20 2a 65 64 67 | 65 5f 6d 61 70 5f 68 69 |ap, *edg|e_map_hi|
|00004c60| 3b 0d 0a 2d 09 69 6e 74 | 20 6c 6f 5f 74 68 72 65 |;..-.int| lo_thre|
|00004c70| 73 68 6f 6c 64 3b 0d 0a | 2d 7d 3b 0d 0a 2d 0d 0a |shold;..|-};..-..|
|00004c80| 2d 0d 0a 2d 0d 0a 2f 2f | 47 4f 2e 53 59 53 49 4e |-..-..//|GO.SYSIN|
|00004c90| 20 44 44 20 65 64 67 65 | 5f 66 69 6e 64 65 72 2e | DD edge|_finder.|
|00004ca0| 68 0d 0a 65 63 68 6f 20 | 65 64 67 65 5f 6d 61 70 |h..echo |edge_map|
|00004cb0| 73 2e 63 20 31 3e 26 32 | 0d 0a 73 65 64 20 27 73 |s.c 1>&2|..sed 's|
|00004cc0| 2f 2e 2f 2f 27 20 3e 65 | 64 67 65 5f 6d 61 70 73 |/.//' >e|dge_maps|
|00004cd0| 2e 63 20 3c 3c 27 2f 2f | 47 4f 2e 53 59 53 49 4e |.c <<'//|GO.SYSIN|
|00004ce0| 20 44 44 20 65 64 67 65 | 5f 6d 61 70 73 2e 63 27 | DD edge|_maps.c'|
|00004cf0| 0d 0a 2d 2f 2a 0d 0a 2d | 20 2a 20 28 63 29 20 43 |..-/*..-| * (c) C|
|00004d00| 6f 70 79 72 69 67 68 74 | 20 31 39 38 38 20 62 79 |opyright| 1988 by|
|00004d10| 0d 0a 2d 20 2a 20 52 6f | 62 6f 74 69 63 73 20 50 |..- * Ro|botics P|
|00004d20| 72 69 6e 63 69 70 6c 65 | 73 20 52 65 73 65 61 72 |rinciple|s Resear|
|00004d30| 63 68 20 44 65 70 61 72 | 74 6d 65 6e 74 2c 20 41 |ch Depar|tment, A|
|00004d40| 54 54 20 42 65 6c 6c 20 | 4c 61 62 6f 72 61 74 6f |TT Bell |Laborato|
|00004d50| 72 69 65 73 2e 0d 0a 2d | 20 2a 20 41 6c 6c 20 72 |ries...-| * All r|
|00004d60| 69 67 68 74 73 20 72 65 | 73 65 72 76 65 64 2e 0d |ights re|served..|
|00004d70| 0a 2d 20 2a 20 4c 61 73 | 74 20 6d 6f 64 69 66 69 |.- * Las|t modifi|
|00004d80| 65 64 20 32 2f 38 2f 38 | 38 20 49 6e 67 65 6d 61 |ed 2/8/8|8 Ingema|
|00004d90| 72 20 4a 2e 20 43 6f 78 | 0d 0a 2d 20 2a 20 43 20 |r J. Cox|..- * C |
|00004da0| 76 65 72 73 69 6f 6e 20 | 38 2f 32 2f 38 38 20 44 |version |8/2/88 D|
|00004db0| 65 62 6f 72 61 68 20 41 | 2e 20 57 61 6c 6c 61 63 |eborah A|. Wallac|
|00004dc0| 68 0d 0a 2d 20 2a 20 43 | 20 76 65 72 73 69 6f 6e |h..- * C| version|
|00004dd0| 20 31 31 2f 32 2f 38 38 | 20 57 2e 20 4a 2e 20 4b | 11/2/88| W. J. K|
|00004de0| 72 6f 70 66 6c 0d 0a 2d | 20 2a 2f 0d 0a 2d 23 69 |ropfl..-| */..-#i|
|00004df0| 6e 63 6c 75 64 65 20 3c | 73 74 64 69 6f 2e 68 3e |nclude <|stdio.h>|
|00004e00| 0d 0a 2d 23 69 6e 63 6c | 75 64 65 20 3c 6d 61 74 |..-#incl|ude <mat|
|00004e10| 68 2e 68 3e 0d 0a 2d 23 | 69 6e 63 6c 75 64 65 20 |h.h>..-#|include |
|00004e20| 22 65 64 67 65 5f 66 69 | 6e 64 65 72 2e 68 22 0d |"edge_fi|nder.h".|
|00004e30| 0a 2d 0d 0a 2d 65 78 74 | 65 72 6e 20 73 74 72 75 |.-..-ext|ern stru|
|00004e40| 63 74 20 69 6d 61 67 65 | 20 2a 6d 79 5f 69 6d 61 |ct image| *my_ima|
|00004e50| 67 65 3b 0d 0a 2d 65 78 | 74 65 72 6e 20 69 6e 74 |ge;..-ex|tern int|
|00004e60| 20 6e 78 58 6e 79 3b 09 | 09 2f 2a 3f 3f 3f 2a 2f | nxXny;.|./*???*/|
|00004e70| 0d 0a 2d 0d 0a 2d 23 64 | 65 66 69 6e 65 20 5f 6e |..-..-#d|efine _n|
|00004e80| 78 20 6d 79 5f 69 6d 61 | 67 65 2d 3e 6e 78 0d 0a |x my_ima|ge->nx..|
|00004e90| 2d 23 64 65 66 69 6e 65 | 20 69 6d 61 67 65 5f 7a |-#define| image_z|
|00004ea0| 63 5f 78 28 69 78 2c 69 | 79 2c 69 29 20 5c 0d 0a |c_x(ix,i|y,i) \..|
|00004eb0| 2d 28 28 28 69 6d 61 67 | 65 5f 6c 6f 63 78 28 69 |-(((imag|e_locx(i|
|00004ec0| 79 2c 69 78 2c 69 29 3e | 30 29 20 21 3d 20 28 69 |y,ix,i)>|0) != (i|
|00004ed0| 6d 61 67 65 5f 6c 6f 63 | 78 28 69 79 2c 69 78 2b |mage_loc|x(iy,ix+|
|00004ee0| 31 2c 69 2b 31 29 3e 30 | 29 29 20 3f 45 44 47 45 |1,i+1)>0|)) ?EDGE|
|00004ef0| 5f 58 3a 30 29 0d 0a 2d | 0d 0a 2d 23 64 65 66 69 |_X:0)..-|..-#defi|
|00004f00| 6e 65 20 69 6d 61 67 65 | 5f 7a 63 5f 79 28 69 78 |ne image|_zc_y(ix|
|00004f10| 2c 69 79 2c 69 29 20 5c | 0d 0a 2d 28 28 28 69 6d |,iy,i) \|..-(((im|
|00004f20| 61 67 65 5f 6c 6f 63 79 | 28 69 79 2b 31 2c 69 78 |age_locy|(iy+1,ix|
|00004f30| 2c 69 2b 5f 6e 78 29 3e | 30 29 20 21 3d 20 28 69 |,i+_nx)>|0) != (i|
|00004f40| 6d 61 67 65 5f 6c 6f 63 | 79 28 69 79 2c 69 78 2c |mage_loc|y(iy,ix,|
|00004f50| 69 29 3e 30 29 29 20 3f | 45 44 47 45 5f 59 3a 30 |i)>0)) ?|EDGE_Y:0|
|00004f60| 29 0d 0a 2d 0d 0a 2d 23 | 64 65 66 69 6e 65 20 69 |)..-..-#|define i|
|00004f70| 6d 61 67 65 5f 7a 63 5f | 31 33 35 28 69 78 2c 69 |mage_zc_|135(ix,i|
|00004f80| 79 2c 69 29 20 5c 0d 0a | 2d 28 28 28 69 6d 61 67 |y,i) \..|-(((imag|
|00004f90| 65 5f 6c 6f 63 31 33 35 | 28 69 79 2b 31 2c 69 78 |e_loc135|(iy+1,ix|
|00004fa0| 2b 31 2c 69 2b 5f 6e 78 | 2b 31 29 3e 30 29 20 21 |+1,i+_nx|+1)>0) !|
|00004fb0| 3d 20 28 69 6d 61 67 65 | 5f 6c 6f 63 31 33 35 28 |= (image|_loc135(|
|00004fc0| 69 79 2c 69 78 2c 69 29 | 3e 30 29 29 20 3f 45 44 |iy,ix,i)|>0)) ?ED|
|00004fd0| 47 45 5f 31 33 35 3a 30 | 29 0d 0a 2d 0d 0a 2d 23 |GE_135:0|)..-..-#|
|00004fe0| 64 65 66 69 6e 65 20 69 | 6d 61 67 65 5f 7a 63 5f |define i|mage_zc_|
|00004ff0| 34 35 28 69 78 2c 69 79 | 2c 69 29 20 5c 0d 0a 2d |45(ix,iy|,i) \..-|
|00005000| 28 28 28 69 6d 61 67 65 | 5f 6c 6f 63 34 35 28 69 |(((image|_loc45(i|
|00005010| 79 2b 31 2c 69 78 2d 31 | 2c 69 2b 5f 6e 78 2d 31 |y+1,ix-1|,i+_nx-1|
|00005020| 29 3e 30 29 20 21 3d 20 | 28 69 6d 61 67 65 5f 6c |)>0) != |(image_l|
|00005030| 6f 63 34 35 28 69 79 2c | 69 78 2c 69 29 3e 30 29 |oc45(iy,|ix,i)>0)|
|00005040| 29 20 3f 45 44 47 45 5f | 34 35 3a 30 29 0d 0a 2d |) ?EDGE_|45:0)..-|
|00005050| 0d 0a 2d 75 6e 73 69 67 | 6e 65 64 20 63 68 61 72 |..-unsig|ned char|
|00005060| 2a 20 69 6d 61 67 65 5f | 65 64 67 65 73 28 74 68 |* image_|edges(th|
|00005070| 72 65 73 68 6f 6c 64 29 | 0d 0a 2d 20 20 20 20 20 |reshold)|..-     |
|00005080| 69 6e 74 20 74 68 72 65 | 73 68 6f 6c 64 3b 0d 0a |int thre|shold;..|
|00005090| 2d 7b 0d 0a 2d 09 75 6e | 73 69 67 6e 65 64 20 63 |-{..-.un|signed c|
|000050a0| 68 61 72 20 2a 6d 61 70 | 3b 0d 0a 2d 09 72 65 67 |har *map|;..-.reg|
|000050b0| 69 73 74 65 72 20 69 6e | 74 20 69 78 2c 20 69 79 |ister in|t ix, iy|
|000050c0| 2c 20 6e 78 2c 20 6f 72 | 69 65 6e 74 5f 66 6c 61 |, nx, or|ient_fla|
|000050d0| 67 2c 20 69 6e 64 65 78 | 2c 20 6e 79 3b 0d 0a 2d |g, index|, ny;..-|
|000050e0| 0d 0a 2d 09 72 65 67 69 | 73 74 65 72 20 69 6e 74 |..-.regi|ster int|
|000050f0| 20 2a 70 78 20 3d 20 6d | 79 5f 69 6d 61 67 65 2d | *px = m|y_image-|
|00005100| 3e 69 64 78 3b 0d 0a 2d | 09 72 65 67 69 73 74 65 |>idx;..-|.registe|
|00005110| 72 20 69 6e 74 20 2a 70 | 79 20 3d 20 6d 79 5f 69 |r int *p|y = my_i|
|00005120| 6d 61 67 65 2d 3e 69 64 | 79 3b 0d 0a 2d 09 72 65 |mage->id|y;..-.re|
|00005130| 67 69 73 74 65 72 20 69 | 6e 74 20 2a 70 34 35 20 |gister i|nt *p45 |
|00005140| 20 3d 20 6d 79 5f 69 6d | 61 67 65 2d 3e 69 64 34 | = my_im|age->id4|
|00005150| 35 3b 0d 0a 2d 09 72 65 | 67 69 73 74 65 72 20 69 |5;..-.re|gister i|
|00005160| 6e 74 20 2a 70 31 33 35 | 20 3d 20 6d 79 5f 69 6d |nt *p135| = my_im|
|00005170| 61 67 65 2d 3e 69 64 31 | 33 35 3b 0d 0a 2d 09 72 |age->id1|35;..-.r|
|00005180| 65 67 69 73 74 65 72 20 | 75 6e 73 69 67 6e 65 64 |egister |unsigned|
|00005190| 20 63 68 61 72 20 2a 6d | 61 70 70 3b 0d 0a 2d 0d | char *m|app;..-.|
|000051a0| 0a 2d 09 6e 78 20 3d 20 | 6d 79 5f 69 6d 61 67 65 |.-.nx = |my_image|
|000051b0| 2d 3e 6e 78 3b 0d 0a 2d | 09 6e 79 20 3d 20 6d 79 |->nx;..-|.ny = my|
|000051c0| 5f 69 6d 61 67 65 2d 3e | 6e 79 3b 0d 0a 2d 09 6e |_image->|ny;..-.n|
|000051d0| 78 58 6e 79 20 3d 20 6e | 78 2a 6e 79 3b 0d 0a 2d |xXny = n|x*ny;..-|
|000051e0| 0d 0a 2d 09 69 66 28 28 | 6d 61 70 20 3d 20 28 75 |..-.if((|map = (u|
|000051f0| 6e 73 69 67 6e 65 64 20 | 63 68 61 72 20 2a 29 6d |nsigned |char *)m|
|00005200| 61 6c 6c 6f 63 28 6e 79 | 2a 6e 78 2a 73 69 7a 65 |alloc(ny|*nx*size|
|00005210| 6f 66 28 75 6e 73 69 67 | 6e 65 64 20 63 68 61 72 |of(unsig|ned char|
|00005220| 29 29 29 3d 3d 30 29 0d | 0a 2d 09 7b 0d 0a 2d 09 |)))==0).|.-.{..-.|
|00005230| 09 66 70 72 69 6e 74 66 | 28 73 74 64 65 72 72 2c |.fprintf|(stderr,|
|00005240| 20 22 65 72 72 6f 72 3a | 20 63 61 6e 6e 6f 74 20 | "error:| cannot |
|00005250| 61 6c 6c 6f 63 61 74 65 | 20 6d 61 70 20 61 72 72 |allocate| map arr|
|00005260| 61 79 5c 6e 22 29 3b 0d | 0a 2d 09 09 65 78 69 74 |ay\n");.|.-..exit|
|00005270| 28 31 29 3b 0d 0a 2d 09 | 7d 0d 0a 2d 09 6d 61 70 |(1);..-.|}..-.map|
|00005280| 70 20 3d 20 6d 61 70 3b | 09 0d 0a 2d 0d 0a 2d 09 |p = map;|...-..-.|
|00005290| 66 6f 72 28 69 79 3d 30 | 3b 20 69 79 3c 6e 79 3b |for(iy=0|; iy<ny;|
|000052a0| 20 69 79 2b 2b 29 0d 0a | 2d 09 7b 0d 0a 2d 09 09 | iy++)..|-.{..-..|
|000052b0| 66 6f 72 28 69 78 3d 30 | 3b 20 69 78 3c 6e 78 3b |for(ix=0|; ix<nx;|
|000052c0| 20 69 78 2b 2b 29 0d 0a | 2d 09 09 7b 0d 0a 2d 09 | ix++)..|-..{..-.|
|000052d0| 09 09 72 65 67 69 73 74 | 65 72 20 69 6e 74 20 73 |..regist|er int s|
|000052e0| 78 2c 20 73 79 2c 20 73 | 34 35 2c 20 73 31 33 35 |x, sy, s|45, s135|
|000052f0| 3b 0d 0a 2d 09 09 09 72 | 65 67 69 73 74 65 72 20 |;..-...r|egister |
|00005300| 69 6e 74 20 62 69 67 3b | 0d 0a 2d 09 09 09 73 78 |int big;|..-...sx|
|00005310| 20 3d 20 61 62 73 28 2a | 70 78 29 3b 20 20 70 78 | = abs(*|px);  px|
|00005320| 2b 2b 3b 0d 0a 2d 09 09 | 09 73 79 20 3d 20 61 62 |++;..-..|.sy = ab|
|00005330| 73 28 2a 70 79 29 3b 20 | 09 70 79 2b 2b 3b 0d 0a |s(*py); |.py++;..|
|00005340| 2d 09 09 09 73 34 35 20 | 3d 20 28 69 6e 74 29 28 |-...s45 |= (int)(|
|00005350| 28 66 6c 6f 61 74 29 28 | 61 62 73 28 2a 70 34 35 |(float)(|abs(*p45|
|00005360| 29 29 2f 64 69 61 67 5f | 73 63 61 6c 65 29 3b 20 |))/diag_|scale); |
|00005370| 70 34 35 2b 2b 3b 0d 0a | 2d 09 09 09 73 31 33 35 |p45++;..|-...s135|
|00005380| 20 3d 20 28 69 6e 74 29 | 28 28 66 6c 6f 61 74 29 | = (int)|((float)|
|00005390| 28 61 62 73 28 2a 70 31 | 33 35 29 29 2f 64 69 61 |(abs(*p1|35))/dia|
|000053a0| 67 5f 73 63 61 6c 65 29 | 3b 20 70 31 33 35 2b 2b |g_scale)|; p135++|
|000053b0| 3b 0d 0a 2d 0d 0a 2d 09 | 09 09 62 69 67 20 3d 20 |;..-..-.|..big = |
|000053c0| 73 78 3b 0d 0a 2d 0d 0a | 2d 09 09 09 6f 72 69 65 |sx;..-..|-...orie|
|000053d0| 6e 74 5f 66 6c 61 67 20 | 3d 20 45 44 47 45 5f 58 |nt_flag |= EDGE_X|
|000053e0| 3b 0d 0a 2d 0d 0a 2d 09 | 09 09 69 66 28 73 79 3e |;..-..-.|..if(sy>|
|000053f0| 62 69 67 29 0d 0a 2d 09 | 09 09 7b 0d 0a 2d 09 09 |big)..-.|..{..-..|
|00005400| 09 09 62 69 67 3d 73 79 | 3b 0d 0a 2d 09 09 09 09 |..big=sy|;..-....|
|00005410| 6f 72 69 65 6e 74 5f 66 | 6c 61 67 20 3d 20 45 44 |orient_f|lag = ED|
|00005420| 47 45 5f 59 3b 0d 0a 2d | 09 09 09 7d 0d 0a 2d 09 |GE_Y;..-|...}..-.|
|00005430| 09 09 69 66 28 73 34 35 | 3e 62 69 67 29 0d 0a 2d |..if(s45|>big)..-|
|00005440| 09 09 09 7b 0d 0a 2d 09 | 09 09 09 62 69 67 3d 73 |...{..-.|...big=s|
|00005450| 34 35 3b 0d 0a 2d 09 09 | 09 09 6f 72 69 65 6e 74 |45;..-..|..orient|
|00005460| 5f 66 6c 61 67 20 3d 20 | 45 44 47 45 5f 34 35 3b |_flag = |EDGE_45;|
|00005470| 0d 0a 2d 09 09 09 7d 0d | 0a 2d 09 09 09 69 66 28 |..-...}.|.-...if(|
|00005480| 73 31 33 35 3e 62 69 67 | 29 0d 0a 2d 09 09 09 7b |s135>big|)..-...{|
|00005490| 0d 0a 2d 09 09 09 09 62 | 69 67 3d 73 31 33 35 3b |..-....b|ig=s135;|
|000054a0| 0d 0a 2d 09 09 09 09 6f | 72 69 65 6e 74 5f 66 6c |..-....o|rient_fl|
|000054b0| 61 67 20 3d 20 45 44 47 | 45 5f 31 33 35 3b 0d 0a |ag = EDG|E_135;..|
|000054c0| 2d 09 09 09 7d 0d 0a 2d | 0d 0a 2d 09 09 09 6d 79 |-...}..-|..-...my|
|000054d0| 5f 69 6d 61 67 65 2d 3e | 6f 72 69 65 6e 74 5f 66 |_image->|orient_f|
|000054e0| 6c 61 67 20 3d 20 6f 72 | 69 65 6e 74 5f 66 6c 61 |lag = or|ient_fla|
|000054f0| 67 3b 20 0d 0a 2d 0d 0a | 2d 09 09 09 69 66 28 62 |g; ..-..|-...if(b|
|00005500| 69 67 20 3e 20 74 68 72 | 65 73 68 6f 6c 64 29 0d |ig > thr|eshold).|
|00005510| 0a 2d 09 09 09 7b 0d 0a | 2d 09 09 09 09 69 6e 64 |.-...{..|-....ind|
|00005520| 65 78 20 3d 20 6d 61 70 | 70 20 2d 20 6d 61 70 3b |ex = map|p - map;|
|00005530| 0d 0a 2d 09 09 09 09 73 | 77 69 74 63 68 20 28 6f |..-....s|witch (o|
|00005540| 72 69 65 6e 74 5f 66 6c | 61 67 29 0d 0a 2d 09 09 |rient_fl|ag)..-..|
|00005550| 09 09 7b 0d 0a 2d 09 09 | 09 09 20 20 20 63 61 73 |..{..-..|..   cas|
|00005560| 65 20 45 44 47 45 5f 58 | 3a 0d 0a 2d 09 09 09 09 |e EDGE_X|:..-....|
|00005570| 09 2a 6d 61 70 70 2b 2b | 20 3d 20 69 6d 61 67 65 |.*mapp++| = image|
|00005580| 5f 7a 63 5f 78 28 69 78 | 2c 20 69 79 2c 20 69 6e |_zc_x(ix|, iy, in|
|00005590| 64 65 78 29 3b 0d 0a 2d | 09 09 09 09 09 62 72 65 |dex);..-|.....bre|
|000055a0| 61 6b 3b 0d 0a 2d 09 09 | 09 09 20 20 20 63 61 73 |ak;..-..|..   cas|
|000055b0| 65 20 45 44 47 45 5f 59 | 3a 0d 0a 2d 09 09 09 09 |e EDGE_Y|:..-....|
|000055c0| 09 2a 6d 61 70 70 2b 2b | 20 3d 20 69 6d 61 67 65 |.*mapp++| = image|
|000055d0| 5f 7a 63 5f 79 28 69 78 | 2c 20 69 79 2c 20 69 6e |_zc_y(ix|, iy, in|
|000055e0| 64 65 78 29 3b 0d 0a 2d | 09 09 09 09 09 62 72 65 |dex);..-|.....bre|
|000055f0| 61 6b 3b 0d 0a 2d 09 09 | 09 09 20 20 20 63 61 73 |ak;..-..|..   cas|
|00005600| 65 20 45 44 47 45 5f 34 | 35 3a 0d 0a 2d 09 09 09 |e EDGE_4|5:..-...|
|00005610| 09 09 2a 6d 61 70 70 2b | 2b 20 3d 20 69 6d 61 67 |..*mapp+|+ = imag|
|00005620| 65 5f 7a 63 5f 34 35 28 | 69 78 2c 20 69 79 2c 20 |e_zc_45(|ix, iy, |
|00005630| 69 6e 64 65 78 29 3b 0d | 0a 2d 09 09 09 09 09 62 |index);.|.-.....b|
|00005640| 72 65 61 6b 3b 0d 0a 2d | 09 09 09 09 20 20 20 63 |reak;..-|....   c|
|00005650| 61 73 65 20 45 44 47 45 | 5f 31 33 35 3a 0d 0a 2d |ase EDGE|_135:..-|
|00005660| 09 09 09 09 09 2a 6d 61 | 70 70 2b 2b 20 3d 20 69 |.....*ma|pp++ = i|
|00005670| 6d 61 67 65 5f 7a 63 5f | 31 33 35 28 69 78 2c 20 |mage_zc_|135(ix, |
|00005680| 69 79 2c 20 69 6e 64 65 | 78 29 3b 0d 0a 2d 09 09 |iy, inde|x);..-..|
|00005690| 09 09 09 62 72 65 61 6b | 3b 0d 0a 2d 09 09 09 09 |...break|;..-....|
|000056a0| 20 20 20 64 65 66 61 75 | 6c 74 3a 0d 0a 2d 09 09 |   defau|lt:..-..|
|000056b0| 09 09 09 66 70 72 69 6e | 74 66 28 73 74 64 65 72 |...fprin|tf(stder|
|000056c0| 72 2c 20 22 65 72 72 6f | 72 20 69 6e 20 63 61 73 |r, "erro|r in cas|
|000056d0| 65 20 73 74 61 74 65 6d | 65 6e 74 73 5c 6e 22 29 |e statem|ents\n")|
|000056e0| 3b 0d 0a 2d 09 09 09 09 | 09 65 78 69 74 28 31 29 |;..-....|.exit(1)|
|000056f0| 3b 0d 0a 2d 09 09 09 09 | 09 62 72 65 61 6b 3b 0d |;..-....|.break;.|
|00005700| 0a 2d 09 09 09 09 7d 0d | 0a 2d 09 09 09 7d 0d 0a |.-....}.|.-...}..|
|00005710| 2d 09 09 09 65 6c 73 65 | 0d 0a 2d 09 09 09 7b 0d |-...else|..-...{.|
|00005720| 0a 2d 09 09 09 09 2a 6d | 61 70 70 2b 2b 20 3d 20 |.-....*m|app++ = |
|00005730| 30 3b 0d 0a 2d 09 09 09 | 7d 0d 0a 2d 09 09 20 20 |0;..-...|}..-..  |
|00005740| 20 20 7d 0d 0a 2d 09 7d | 0d 0a 2d 23 69 66 64 65 |  }..-.}|..-#ifde|
|00005750| 66 20 44 45 42 55 47 0d | 0a 2d 09 69 6d 61 67 65 |f DEBUG.|.-.image|
|00005760| 5f 57 72 69 74 65 5f 63 | 68 61 72 28 22 68 69 5f |_Write_c|har("hi_|
|00005770| 6d 61 70 22 2c 20 6d 61 | 70 29 3b 0d 0a 2d 23 65 |map", ma|p);..-#e|
|00005780| 6e 64 69 66 0d 0a 2d 09 | 72 65 74 75 72 6e 28 6d |ndif..-.|return(m|
|00005790| 61 70 29 3b 0d 0a 2d 7d | 0d 0a 2d 0d 0a 2d 69 6d |ap);..-}|..-..-im|
|000057a0| 61 67 65 5f 65 64 67 65 | 5f 6d 61 70 5f 6c 6f 28 |age_edge|_map_lo(|
|000057b0| 69 78 2c 20 69 79 2c 20 | 69 6e 64 65 78 29 0d 0a |ix, iy, |index)..|
|000057c0| 2d 72 65 67 69 73 74 65 | 72 20 69 6e 74 20 69 78 |-registe|r int ix|
|000057d0| 2c 20 69 79 2c 20 69 6e | 64 65 78 3b 09 09 2f 2a |, iy, in|dex;../*|
|000057e0| 20 69 6e 64 65 78 20 3d | 20 69 79 2a 6e 78 2b 69 | index =| iy*nx+i|
|000057f0| 78 3b 20 2a 2f 0d 0a 2d | 7b 0d 0a 2d 09 72 65 67 |x; */..-|{..-.reg|
|00005800| 69 73 74 65 72 20 69 6e | 74 20 62 69 67 2c 20 73 |ister in|t big, s|
|00005810| 78 2c 20 73 79 2c 20 73 | 34 35 2c 20 73 31 33 35 |x, sy, s|45, s135|
|00005820| 3b 0d 0a 2d 0d 0a 2d 09 | 73 78 20 3d 20 61 62 73 |;..-..-.|sx = abs|
|00005830| 28 6d 79 5f 69 6d 61 67 | 65 2d 3e 69 64 78 5b 69 |(my_imag|e->idx[i|
|00005840| 6e 64 65 78 5d 29 3b 0d | 0a 2d 09 73 79 20 3d 20 |ndex]);.|.-.sy = |
|00005850| 61 62 73 28 6d 79 5f 69 | 6d 61 67 65 2d 3e 69 64 |abs(my_i|mage->id|
|00005860| 79 5b 69 6e 64 65 78 5d | 29 3b 0d 0a 2d 09 73 34 |y[index]|);..-.s4|
|00005870| 35 20 20 3d 20 28 69 6e | 74 29 28 28 66 6c 6f 61 |5  = (in|t)((floa|
|00005880| 74 29 28 61 62 73 28 6d | 79 5f 69 6d 61 67 65 2d |t)(abs(m|y_image-|
|00005890| 3e 69 64 34 35 5b 69 6e | 64 65 78 5d 29 29 2f 64 |>id45[in|dex]))/d|
|000058a0| 69 61 67 5f 73 63 61 6c | 65 29 3b 0d 0a 2d 09 73 |iag_scal|e);..-.s|
|000058b0| 31 33 35 20 3d 20 28 69 | 6e 74 29 28 28 66 6c 6f |135 = (i|nt)((flo|
|000058c0| 61 74 29 28 61 62 73 28 | 6d 79 5f 69 6d 61 67 65 |at)(abs(|my_image|
|000058d0| 2d 3e 69 64 31 33 35 5b | 69 6e 64 65 78 5d 29 29 |->id135[|index]))|
|000058e0| 2f 64 69 61 67 5f 73 63 | 61 6c 65 29 3b 0d 0a 2d |/diag_sc|ale);..-|
|000058f0| 09 62 69 67 3d 73 78 3b | 0d 0a 2d 0d 0a 2d 09 6d |.big=sx;|..-..-.m|
|00005900| 79 5f 69 6d 61 67 65 2d | 3e 6f 72 69 65 6e 74 5f |y_image-|>orient_|
|00005910| 66 6c 61 67 20 3d 20 45 | 44 47 45 5f 58 3b 0d 0a |flag = E|DGE_X;..|
|00005920| 2d 09 0d 0a 2d 09 69 66 | 28 73 79 3e 62 69 67 29 |-...-.if|(sy>big)|
|00005930| 0d 0a 2d 09 7b 0d 0a 2d | 09 09 62 69 67 3d 73 79 |..-.{..-|..big=sy|
|00005940| 3b 0d 0a 2d 09 09 6d 79 | 5f 69 6d 61 67 65 2d 3e |;..-..my|_image->|
|00005950| 6f 72 69 65 6e 74 5f 66 | 6c 61 67 20 3d 20 45 44 |orient_f|lag = ED|
|00005960| 47 45 5f 59 3b 0d 0a 2d | 09 7d 0d 0a 2d 09 69 66 |GE_Y;..-|.}..-.if|
|00005970| 28 73 34 35 3e 62 69 67 | 29 0d 0a 2d 09 7b 0d 0a |(s45>big|)..-.{..|
|00005980| 2d 09 09 62 69 67 3d 73 | 34 35 3b 0d 0a 2d 09 09 |-..big=s|45;..-..|
|00005990| 6d 79 5f 69 6d 61 67 65 | 2d 3e 6f 72 69 65 6e 74 |my_image|->orient|
|000059a0| 5f 66 6c 61 67 20 3d 20 | 45 44 47 45 5f 34 35 3b |_flag = |EDGE_45;|
|000059b0| 0d 0a 2d 09 7d 0d 0a 2d | 09 69 66 28 73 31 33 35 |..-.}..-|.if(s135|
|000059c0| 3e 62 69 67 29 0d 0a 2d | 09 7b 0d 0a 2d 09 09 62 |>big)..-|.{..-..b|
|000059d0| 69 67 3d 73 31 33 35 3b | 0d 0a 2d 09 09 6d 79 5f |ig=s135;|..-..my_|
|000059e0| 69 6d 61 67 65 2d 3e 6f | 72 69 65 6e 74 5f 66 6c |image->o|rient_fl|
|000059f0| 61 67 20 3d 20 45 44 47 | 45 5f 31 33 35 3b 0d 0a |ag = EDG|E_135;..|
|00005a00| 2d 09 7d 0d 0a 2d 0d 0a | 2d 09 69 66 28 62 69 67 |-.}..-..|-.if(big|
|00005a10| 20 3e 20 6d 79 5f 69 6d | 61 67 65 2d 3e 6c 6f 5f | > my_im|age->lo_|
|00005a20| 74 68 72 65 73 68 6f 6c | 64 29 0d 0a 2d 09 7b 0d |threshol|d)..-.{.|
|00005a30| 0a 2d 09 09 73 77 69 74 | 63 68 20 28 6d 79 5f 69 |.-..swit|ch (my_i|
|00005a40| 6d 61 67 65 2d 3e 6f 72 | 69 65 6e 74 5f 66 6c 61 |mage->or|ient_fla|
|00005a50| 67 29 0d 0a 2d 09 09 7b | 0d 0a 2d 09 09 09 63 61 |g)..-..{|..-...ca|
|00005a60| 73 65 20 45 44 47 45 5f | 58 3a 0d 0a 2d 09 09 09 |se EDGE_|X:..-...|
|00005a70| 09 72 65 74 75 72 6e 28 | 69 6d 61 67 65 5f 7a 63 |.return(|image_zc|
|00005a80| 5f 78 28 69 78 2c 20 69 | 79 2c 20 69 6e 64 65 78 |_x(ix, i|y, index|
|00005a90| 29 29 3b 0d 0a 2d 09 09 | 09 63 61 73 65 20 45 44 |));..-..|.case ED|
|00005aa0| 47 45 5f 59 3a 0d 0a 2d | 09 09 09 09 72 65 74 75 |GE_Y:..-|....retu|
|00005ab0| 72 6e 28 69 6d 61 67 65 | 5f 7a 63 5f 79 28 69 78 |rn(image|_zc_y(ix|
|00005ac0| 2c 20 69 79 2c 20 69 6e | 64 65 78 29 29 3b 0d 0a |, iy, in|dex));..|
|00005ad0| 2d 09 09 09 63 61 73 65 | 20 45 44 47 45 5f 34 35 |-...case| EDGE_45|
|00005ae0| 3a 0d 0a 2d 09 09 09 09 | 72 65 74 75 72 6e 28 69 |:..-....|return(i|
|00005af0| 6d 61 67 65 5f 7a 63 5f | 34 35 28 69 78 2c 20 69 |mage_zc_|45(ix, i|
|00005b00| 79 2c 20 69 6e 64 65 78 | 29 29 3b 0d 0a 2d 09 09 |y, index|));..-..|
|00005b10| 09 63 61 73 65 20 45 44 | 47 45 5f 31 33 35 3a 0d |.case ED|GE_135:.|
|00005b20| 0a 2d 09 09 09 09 72 65 | 74 75 72 6e 28 69 6d 61 |.-....re|turn(ima|
|00005b30| 67 65 5f 7a 63 5f 31 33 | 35 28 69 78 2c 20 69 79 |ge_zc_13|5(ix, iy|
|00005b40| 2c 20 69 6e 64 65 78 29 | 29 3b 0d 0a 2d 09 09 09 |, index)|);..-...|
|00005b50| 64 65 66 61 75 6c 74 3a | 0d 0a 2d 09 09 09 09 66 |default:|..-....f|
|00005b60| 70 72 69 6e 74 66 28 73 | 74 64 65 72 72 2c 20 22 |printf(s|tderr, "|
|00005b70| 65 72 72 6f 72 20 69 6e | 20 63 61 73 65 20 73 74 |error in| case st|
|00005b80| 61 74 65 6d 65 6e 74 73 | 5c 6e 22 29 3b 0d 0a 2d |atements|\n");..-|
|00005b90| 09 09 09 09 65 78 69 74 | 28 31 29 3b 0d 0a 2d 09 |....exit|(1);..-.|
|00005ba0| 09 09 09 62 72 65 61 6b | 3b 0d 0a 2d 09 09 7d 0d |...break|;..-..}.|
|00005bb0| 0a 2d 09 7d 0d 0a 2d 09 | 65 6c 73 65 0d 0a 2d 09 |.-.}..-.|else..-.|
|00005bc0| 09 72 65 74 75 72 6e 28 | 30 29 3b 0d 0a 2d 7d 0d |.return(|0);..-}.|
|00005bd0| 0a 2f 2f 47 4f 2e 53 59 | 53 49 4e 20 44 44 20 65 |.//GO.SY|SIN DD e|
|00005be0| 64 67 65 5f 6d 61 70 73 | 2e 63 0d 0a 65 63 68 6f |dge_maps|.c..echo|
|00005bf0| 20 66 69 6c 74 65 72 2e | 63 20 31 3e 26 32 0d 0a | filter.|c 1>&2..|
|00005c00| 73 65 64 20 27 73 2f 2e | 2f 2f 27 20 3e 66 69 6c |sed 's/.|//' >fil|
|00005c10| 74 65 72 2e 63 20 3c 3c | 27 2f 2f 47 4f 2e 53 59 |ter.c <<|'//GO.SY|
|00005c20| 53 49 4e 20 44 44 20 66 | 69 6c 74 65 72 2e 63 27 |SIN DD f|ilter.c'|
|00005c30| 0d 0a 2d 2f 2a 0d 0a 2d | 20 2a 20 28 63 29 20 43 |..-/*..-| * (c) C|
|00005c40| 6f 70 79 72 69 67 68 74 | 20 31 39 38 38 20 62 79 |opyright| 1988 by|
|00005c50| 0d 0a 2d 20 2a 20 52 6f | 62 6f 74 69 63 73 20 50 |..- * Ro|botics P|
|00005c60| 72 69 6e 63 69 70 6c 65 | 73 20 52 65 73 65 61 72 |rinciple|s Resear|
|00005c70| 63 68 20 44 65 70 61 72 | 74 6d 65 6e 74 2c 20 41 |ch Depar|tment, A|
|00005c80| 54 54 20 42 65 6c 6c 20 | 4c 61 62 6f 72 61 74 6f |TT Bell |Laborato|
|00005c90| 72 69 65 73 2e 0d 0a 2d | 20 2a 20 41 6c 6c 20 72 |ries...-| * All r|
|00005ca0| 69 67 68 74 73 20 72 65 | 73 65 72 76 65 64 2e 0d |ights re|served..|
|00005cb0| 0a 2d 20 2a 20 4c 61 73 | 74 20 6d 6f 64 69 66 69 |.- * Las|t modifi|
|00005cc0| 65 64 20 32 2f 38 2f 38 | 38 20 49 6e 67 65 6d 61 |ed 2/8/8|8 Ingema|
|00005cd0| 72 20 4a 2e 20 43 6f 78 | 0d 0a 2d 20 2a 20 43 20 |r J. Cox|..- * C |
|00005ce0| 76 65 72 73 69 6f 6e 20 | 38 2f 32 2f 38 38 20 44 |version |8/2/88 D|
|00005cf0| 65 62 6f 72 61 68 20 41 | 2e 20 57 61 6c 6c 61 63 |eborah A|. Wallac|
|00005d00| 68 0d 0a 2d 20 2a 2f 0d | 0a 2d 23 69 6e 63 6c 75 |h..- */.|.-#inclu|
|00005d10| 64 65 20 3c 73 74 64 69 | 6f 2e 68 3e 0d 0a 2d 23 |de <stdi|o.h>..-#|
|00005d20| 69 6e 63 6c 75 64 65 20 | 3c 6d 61 74 68 2e 68 3e |include |<math.h>|
|00005d30| 0d 0a 2d 23 69 6e 63 6c | 75 64 65 20 22 65 64 67 |..-#incl|ude "edg|
|00005d40| 65 5f 66 69 6e 64 65 72 | 2e 68 22 0d 0a 2d 0d 0a |e_finder|.h"..-..|
|00005d50| 2d 65 78 74 65 72 6e 20 | 73 74 72 75 63 74 20 69 |-extern |struct i|
|00005d60| 6d 61 67 65 20 2a 6d 79 | 5f 69 6d 61 67 65 3b 0d |mage *my|_image;.|
|00005d70| 0a 2d 0d 0a 2d 76 6f 69 | 64 20 67 65 6e 65 72 61 |.-..-voi|d genera|
|00005d80| 74 65 5f 66 69 6c 74 65 | 72 28 73 69 67 6d 61 29 |te_filte|r(sigma)|
|00005d90| 0d 0a 2d 20 20 20 20 20 | 66 6c 6f 61 74 20 73 69 |..-     |float si|
|00005da0| 67 6d 61 3b 0d 0a 2d 7b | 0d 0a 2d 09 69 6e 74 20 |gma;..-{|..-.int |
|00005db0| 69 3b 0d 0a 2d 09 66 6c | 6f 61 74 20 74 77 6f 5f |i;..-.fl|oat two_|
|00005dc0| 73 69 67 6d 61 5f 73 71 | 3b 0d 0a 2d 0d 0a 2d 09 |sigma_sq|;..-..-.|
|00005dd0| 6d 79 5f 69 6d 61 67 65 | 2d 3e 6e 66 20 3d 20 28 |my_image|->nf = (|
|00005de0| 69 6e 74 29 28 33 2a 32 | 2a 73 69 67 6d 61 29 3b |int)(3*2|*sigma);|
|00005df0| 0d 0a 2d 09 6d 79 5f 69 | 6d 61 67 65 2d 3e 6e 66 |..-.my_i|mage->nf|
|00005e00| 20 3d 20 28 6d 79 5f 69 | 6d 61 67 65 2d 3e 6e 66 | = (my_i|mage->nf|
|00005e10| 25 32 3e 30 3f 6d 79 5f | 69 6d 61 67 65 2d 3e 6e |%2>0?my_|image->n|
|00005e20| 66 3a 6d 79 5f 69 6d 61 | 67 65 2d 3e 6e 66 2b 31 |f:my_ima|ge->nf+1|
|00005e30| 29 3b 0d 0a 2d 09 6d 79 | 5f 69 6d 61 67 65 2d 3e |);..-.my|_image->|
|00005e40| 6e 66 32 20 3d 09 6d 79 | 5f 69 6d 61 67 65 2d 3e |nf2 =.my|_image->|
|00005e50| 6e 66 2f 32 3b 0d 0a 2d | 09 0d 0a 2d 09 69 66 28 |nf/2;..-|...-.if(|
|00005e60| 28 6d 79 5f 69 6d 61 67 | 65 2d 3e 66 69 6c 74 65 |(my_imag|e->filte|
|00005e70| 72 20 3d 20 28 69 6e 74 | 20 2a 29 6d 61 6c 6c 6f |r = (int| *)mallo|
|00005e80| 63 28 6d 79 5f 69 6d 61 | 67 65 2d 3e 6e 66 2a 73 |c(my_ima|ge->nf*s|
|00005e90| 69 7a 65 6f 66 28 69 6e | 74 29 29 29 3d 3d 30 29 |izeof(in|t)))==0)|
|00005ea0| 0d 0a 2d 09 7b 0d 0a 2d | 09 09 66 70 72 69 6e 74 |..-.{..-|..fprint|
|00005eb0| 66 28 73 74 64 65 72 72 | 2c 20 22 65 72 72 6f 72 |f(stderr|, "error|
|00005ec0| 3a 20 63 61 6e 6e 6f 74 | 20 61 6c 6c 6f 63 61 74 |: cannot| allocat|
|00005ed0| 65 20 66 69 6c 74 65 72 | 20 61 72 72 61 79 5c 6e |e filter| array\n|
|00005ee0| 22 29 3b 0d 0a 2d 09 09 | 65 78 69 74 28 31 29 3b |");..-..|exit(1);|
|00005ef0| 0d 0a 2d 09 7d 0d 0a 2d | 09 74 77 6f 5f 73 69 67 |..-.}..-|.two_sig|
|00005f00| 6d 61 5f 73 71 20 3d 20 | 32 2e 30 2a 73 69 67 6d |ma_sq = |2.0*sigm|
|00005f10| 61 2a 73 69 67 6d 61 3b | 0d 0a 2d 09 6d 79 5f 69 |a*sigma;|..-.my_i|
|00005f20| 6d 61 67 65 2d 3e 6e 6f | 72 6d 3d 30 3b 0d 0a 2d |mage->no|rm=0;..-|
|00005f30| 09 66 6f 72 28 69 3d 30 | 3b 69 3c 6d 79 5f 69 6d |.for(i=0|;i<my_im|
|00005f40| 61 67 65 2d 3e 6e 66 3b | 69 2b 2b 29 0d 0a 2d 09 |age->nf;|i++)..-.|
|00005f50| 7b 0d 0a 2d 09 09 66 6c | 6f 61 74 20 78 20 3d 20 |{..-..fl|oat x = |
|00005f60| 69 2d 6d 79 5f 69 6d 61 | 67 65 2d 3e 6e 66 32 3b |i-my_ima|ge->nf2;|
|00005f70| 0d 0a 2d 09 09 6d 79 5f | 69 6d 61 67 65 2d 3e 66 |..-..my_|image->f|
|00005f80| 69 6c 74 65 72 5b 69 5d | 20 3d 20 28 69 6e 74 29 |ilter[i]| = (int)|
|00005f90| 28 32 35 35 2e 30 2a 65 | 78 70 28 2d 78 2a 78 2f |(255.0*e|xp(-x*x/|
|00005fa0| 74 77 6f 5f 73 69 67 6d | 61 5f 73 71 29 29 3b 0d |two_sigm|a_sq));.|
|00005fb0| 0a 2d 09 09 6d 79 5f 69 | 6d 61 67 65 2d 3e 6e 6f |.-..my_i|mage->no|
|00005fc0| 72 6d 20 2b 3d 20 6d 79 | 5f 69 6d 61 67 65 2d 3e |rm += my|_image->|
|00005fd0| 66 69 6c 74 65 72 5b 69 | 5d 3b 0d 0a 2d 09 7d 0d |filter[i|];..-.}.|
|00005fe0| 0a 2d 09 6d 79 5f 69 6d | 61 67 65 2d 3e 6e 6f 72 |.-.my_im|age->nor|
|00005ff0| 6d 3d 61 62 73 28 6d 79 | 5f 69 6d 61 67 65 2d 3e |m=abs(my|_image->|
|00006000| 6e 6f 72 6d 29 3b 0d 0a | 2d 23 69 66 64 65 66 20 |norm);..|-#ifdef |
|00006010| 44 45 42 55 47 0d 0a 2d | 09 66 70 72 69 6e 74 66 |DEBUG..-|.fprintf|
|00006020| 28 73 74 64 65 72 72 2c | 20 22 66 69 6c 74 65 72 |(stderr,| "filter|
|00006030| 20 6c 65 6e 67 74 68 20 | 3d 20 25 64 5c 6e 22 2c | length |= %d\n",|
|00006040| 20 6d 79 5f 69 6d 61 67 | 65 2d 3e 6e 66 29 3b 0d | my_imag|e->nf);.|
|00006050| 0a 2d 09 66 6f 72 28 69 | 3d 30 3b 69 3c 6d 79 5f |.-.for(i|=0;i<my_|
|00006060| 69 6d 61 67 65 2d 3e 6e | 66 3b 69 2b 2b 29 0d 0a |image->n|f;i++)..|
|00006070| 2d 09 09 66 70 72 69 6e | 74 66 28 73 74 64 65 72 |-..fprin|tf(stder|
|00006080| 72 2c 22 66 5b 25 64 5d | 20 3d 20 25 64 5c 6e 22 |r,"f[%d]| = %d\n"|
|00006090| 2c 20 69 2c 20 6d 79 5f | 69 6d 61 67 65 2d 3e 66 |, i, my_|image->f|
|000060a0| 69 6c 74 65 72 5b 69 5d | 29 3b 0d 0a 2d 23 65 6e |ilter[i]|);..-#en|
|000060b0| 64 69 66 0d 0a 2d 7d 0d | 0a 2f 2f 47 4f 2e 53 59 |dif..-}.|.//GO.SY|
|000060c0| 53 49 4e 20 44 44 20 66 | 69 6c 74 65 72 2e 63 0d |SIN DD f|ilter.c.|
|000060d0| 0a 65 63 68 6f 20 66 69 | 6e 64 5f 65 64 67 65 73 |.echo fi|nd_edges|
|000060e0| 2e 63 20 31 3e 26 32 0d | 0a 73 65 64 20 27 73 2f |.c 1>&2.|.sed 's/|
|000060f0| 2e 2f 2f 27 20 3e 66 69 | 6e 64 5f 65 64 67 65 73 |.//' >fi|nd_edges|
|00006100| 2e 63 20 3c 3c 27 2f 2f | 47 4f 2e 53 59 53 49 4e |.c <<'//|GO.SYSIN|
|00006110| 20 44 44 20 66 69 6e 64 | 5f 65 64 67 65 73 2e 63 | DD find|_edges.c|
|00006120| 27 0d 0a 2d 2f 2a 0d 0a | 2d 20 2a 20 28 63 29 20 |'..-/*..|- * (c) |
|00006130| 43 6f 70 79 72 69 67 68 | 74 20 31 39 38 38 20 62 |Copyrigh|t 1988 b|
|00006140| 79 0d 0a 2d 20 2a 20 52 | 6f 62 6f 74 69 63 73 20 |y..- * R|obotics |
|00006150| 50 72 69 6e 63 69 70 6c | 65 73 20 52 65 73 65 61 |Principl|es Resea|
|00006160| 72 63 68 20 44 65 70 61 | 72 74 6d 65 6e 74 2c 20 |rch Depa|rtment, |
|00006170| 41 54 54 20 42 65 6c 6c | 20 4c 61 62 6f 72 61 74 |ATT Bell| Laborat|
|00006180| 6f 72 69 65 73 2e 0d 0a | 2d 20 2a 20 41 6c 6c 20 |ories...|- * All |
|00006190| 72 69 67 68 74 73 20 72 | 65 73 65 72 76 65 64 2e |rights r|eserved.|
|000061a0| 0d 0a 2d 20 2a 20 4c 61 | 73 74 20 6d 6f 64 69 66 |..- * La|st modif|
|000061b0| 69 65 64 20 32 2f 38 2f | 38 38 20 49 6e 67 65 6d |ied 2/8/|88 Ingem|
|000061c0| 61 72 20 4a 2e 20 43 6f | 78 0d 0a 2d 20 2a 20 43 |ar J. Co|x..- * C|
|000061d0| 20 76 65 72 73 69 6f 6e | 20 38 2f 32 2f 38 38 20 | version| 8/2/88 |
|000061e0| 44 65 62 6f 72 61 68 20 | 41 2e 20 57 61 6c 6c 61 |Deborah |A. Walla|
|000061f0| 63 68 0d 0a 2d 20 2a 20 | 41 75 74 6f 6d 61 74 69 |ch..- * |Automati|
|00006200| 63 20 74 68 72 65 73 68 | 6f 6c 64 69 6e 67 20 39 |c thresh|olding 9|
|00006210| 2f 33 30 2f 38 38 20 49 | 6e 67 65 6d 61 72 20 4a |/30/88 I|ngemar J|
|00006220| 2e 20 43 6f 78 0d 0a 2d | 20 2a 2f 0d 0a 2d 23 69 |. Cox..-| */..-#i|
|00006230| 6e 63 6c 75 64 65 20 3c | 73 74 64 69 6f 2e 68 3e |nclude <|stdio.h>|
|00006240| 0d 0a 2d 23 69 6e 63 6c | 75 64 65 20 3c 6d 61 74 |..-#incl|ude <mat|
|00006250| 68 2e 68 3e 0d 0a 2d 23 | 69 6e 63 6c 75 64 65 20 |h.h>..-#|include |
|00006260| 22 65 64 67 65 5f 66 69 | 6e 64 65 72 2e 68 22 0d |"edge_fi|nder.h".|
|00006270| 0a 2d 0d 0a 2d 65 78 74 | 65 72 6e 20 73 74 72 75 |.-..-ext|ern stru|
|00006280| 63 74 20 69 6d 61 67 65 | 20 2a 6d 79 5f 69 6d 61 |ct image| *my_ima|
|00006290| 67 65 3b 0d 0a 2d 0d 0a | 2d 75 6e 73 69 67 6e 65 |ge;..-..|-unsigne|
|000062a0| 64 20 63 68 61 72 2a 0d | 0a 2d 69 6d 61 67 65 5f |d char*.|.-image_|
|000062b0| 66 69 6e 64 5f 65 64 67 | 65 73 28 73 69 67 6d 61 |find_edg|es(sigma|
|000062c0| 2c 20 74 68 72 65 73 68 | 6f 6c 64 29 0d 0a 2d 20 |, thresh|old)..- |
|000062d0| 20 20 20 20 66 6c 6f 61 | 74 20 73 69 67 6d 61 3b |    floa|t sigma;|
|000062e0| 0d 0a 2d 20 20 20 20 20 | 69 6e 74 20 2a 74 68 72 |..-     |int *thr|
|000062f0| 65 73 68 6f 6c 64 3b 0d | 0a 2d 7b 0d 0a 2d 20 20 |eshold;.|.-{..-  |
|00006300| 20 20 75 6e 73 69 67 6e | 65 64 20 63 68 61 72 20 |  unsign|ed char |
|00006310| 2a 69 6d 61 67 65 5f 65 | 64 67 65 73 28 29 3b 0d |*image_e|dges();.|
|00006320| 0a 2d 0d 0a 2d 09 66 70 | 72 69 6e 74 66 28 73 74 |.-..-.fp|rintf(st|
|00006330| 64 65 72 72 2c 20 22 5c | 6e 67 65 6e 65 72 61 74 |derr, "\|ngenerat|
|00006340| 69 6e 67 20 66 69 6c 74 | 65 72 20 73 69 67 6d 61 |ing filt|er sigma|
|00006350| 20 3d 20 25 2e 31 66 22 | 2c 20 73 69 67 6d 61 29 | = %.1f"|, sigma)|
|00006360| 3b 0d 0a 2d 09 67 65 6e | 65 72 61 74 65 5f 66 69 |;..-.gen|erate_fi|
|00006370| 6c 74 65 72 28 73 69 67 | 6d 61 29 3b 0d 0a 2d 09 |lter(sig|ma);..-.|
|00006380| 66 70 72 69 6e 74 66 28 | 73 74 64 65 72 72 2c 20 |fprintf(|stderr, |
|00006390| 22 20 28 20 6c 65 6e 67 | 74 68 20 3d 20 25 64 29 |" ( leng|th = %d)|
|000063a0| 5c 6e 22 2c 20 6d 79 5f | 69 6d 61 67 65 2d 3e 6e |\n", my_|image->n|
|000063b0| 66 29 3b 0d 0a 2d 0d 0a | 2d 09 66 70 72 69 6e 74 |f);..-..|-.fprint|
|000063c0| 66 28 73 74 64 65 72 72 | 2c 20 22 5c 6e 63 6f 6e |f(stderr|, "\ncon|
|000063d0| 76 6f 6c 76 69 6e 67 20 | 69 6d 61 67 65 5c 6e 22 |volving |image\n"|
|000063e0| 29 3b 0d 0a 2d 09 69 6d | 61 67 65 5f 63 6f 6e 76 |);..-.im|age_conv|
|000063f0| 6f 6c 76 65 28 6d 79 5f | 69 6d 61 67 65 2d 3e 61 |olve(my_|image->a|
+--------+-------------------------+-------------------------+--------+--------+
Only 25.0 KB of data is shown above.