3 * Copyright (c) 2002 Philip Gladstone
5 * This file implements a fish detector. It is used to see when a
6 * goldfish passes in front of the camera. It does this by counting
7 * the number of input pixels that fall within a particular HSV
10 * It takes a multitude of arguments:
12 * -h <num>-<num> the range of H values that are fish
13 * -s <num>-<num> the range of S values that are fish
14 * -v <num>-<num> the range of V values that are fish
15 * -z zap all non-fish values to black
16 * -l <num> limit the number of saved files to <num>
17 * -i <num> only check frames every <num> seconds
18 * -t <num> the threshold for the amount of fish pixels (range 0-1)
19 * -d turn debugging on
20 * -D <directory> where to put the fish images
22 * This library is free software; you can redistribute it and/or
23 * modify it under the terms of the GNU Lesser General Public
24 * License as published by the Free Software Foundation; either
25 * version 2 of the License, or (at your option) any later version.
27 * This library is distributed in the hope that it will be useful,
28 * but WITHOUT ANY WARRANTY; without even the implied warranty of
29 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
30 * Lesser General Public License for more details.
32 * You should have received a copy of the GNU Lesser General Public
33 * License along with this library; if not, write to the Free Software
34 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
45 #include "framehook.h"
49 #define ONE_HALF (1 << (SCALEBITS - 1))
50 #define FIX(x) ((int) ((x) * (1<<SCALEBITS) + 0.5))
52 #define YUV_TO_RGB1_CCIR(cb1, cr1)\
56 r_add = FIX(1.40200*255.0/224.0) * cr + ONE_HALF;\
57 g_add = - FIX(0.34414*255.0/224.0) * cb - FIX(0.71414*255.0/224.0) * cr + \
59 b_add = FIX(1.77200*255.0/224.0) * cb + ONE_HALF;\
62 #define YUV_TO_RGB2_CCIR(r, g, b, y1)\
64 yt = ((y1) - 16) * FIX(255.0/219.0);\
65 r = cm[(yt + r_add) >> SCALEBITS];\
66 g = cm[(yt + g_add) >> SCALEBITS];\
67 b = cm[(yt + b_add) >> SCALEBITS];\
92 static void dorange(const char *s
, int *first
, int *second
, int maxval
)
94 sscanf(s
, "%d-%d", first
, second
);
101 void Release(void *ctx
)
107 int Configure(void **ctxp
, int argc
, char *argv
[])
112 *ctxp
= av_mallocz(sizeof(ContextInfo
));
113 ci
= (ContextInfo
*) *ctxp
;
119 ci
->file_limit
= 100;
120 ci
->min_interval
= 1000000;
121 ci
->inset
= 10; /* Percent */
123 while ((c
= getopt(argc
, argv
, "w:i:dh:s:v:zl:t:D:")) > 0) {
126 dorange(optarg
, &ci
->dark
.h
, &ci
->bright
.h
, 360);
129 dorange(optarg
, &ci
->dark
.s
, &ci
->bright
.s
, 255);
132 dorange(optarg
, &ci
->dark
.v
, &ci
->bright
.v
, 255);
138 ci
->file_limit
= atoi(optarg
);
141 ci
->min_interval
= 1000000 * atof(optarg
);
144 ci
->threshold
= atof(optarg
) * 1000;
145 if (ci
->threshold
> 1000 || ci
->threshold
< 0) {
146 fprintf(stderr
, "Invalid threshold value '%s' (range is 0-1)\n", optarg
);
151 ci
->min_width
= atoi(optarg
);
157 ci
->dir
= av_strdup(optarg
);
160 fprintf(stderr
, "Unrecognized argument '%s'\n", argv
[optind
]);
165 fprintf(stderr
, "Fish detector configured:\n");
166 fprintf(stderr
, " HSV range: %d,%d,%d - %d,%d,%d\n",
173 fprintf(stderr
, " Threshold is %d%% pixels\n", ci
->threshold
/ 10);
179 static void get_hsv(HSV
*hsv
, int r
, int g
, int b
)
208 hsv
->h
= i
+ (60 * f
) / (v
- x
);
212 hsv
->s
= (255 * (v
- x
)) / v
;
218 void Process(void *ctx
, AVPicture
*picture
, enum PixelFormat pix_fmt
, int width
, int height
, int64_t pts
)
220 ContextInfo
*ci
= (ContextInfo
*) ctx
;
221 uint8_t *cm
= cropTbl
+ MAX_NEG_CROP
;
222 int rowsize
= picture
->linesize
[0];
225 printf("pix_fmt = %d, width = %d, pts = %lld, ci->next_pts = %lld\n",
226 pix_fmt
, width
, pts
, ci
->next_pts
);
229 if (pts
< ci
->next_pts
)
232 if (width
< ci
->min_width
)
235 ci
->next_pts
= pts
+ 1000000;
237 if (pix_fmt
== PIX_FMT_YUV420P
) {
239 int width2
= width
>> 1;
246 h_end
= 2 * ((ci
->inset
* height
) / 200);
247 h_start
= height
- h_end
;
249 w_end
= (ci
->inset
* width2
) / 100;
250 w_start
= width2
- w_end
;
252 pixcnt
= ((h_start
- h_end
) >> 1) * (w_start
- w_end
);
254 y
= picture
->data
[0] + h_end
* picture
->linesize
[0] + w_end
* 2;
255 u
= picture
->data
[1] + h_end
* picture
->linesize
[1] / 2 + w_end
;
256 v
= picture
->data
[2] + h_end
* picture
->linesize
[2] / 2 + w_end
;
258 for (h
= h_start
; h
> h_end
; h
-= 2) {
261 for (w
= w_start
; w
> w_end
; w
--) {
264 int cb
, cr
, yt
, r_add
, g_add
, b_add
;
266 YUV_TO_RGB1_CCIR(u
[0], v
[0]);
267 YUV_TO_RGB2_CCIR(r
, g
, b
, y
[0]);
269 get_hsv(&hsv
, r
, g
, b
);
272 fprintf(stderr
, "(%d,%d,%d) -> (%d,%d,%d)\n",
273 r
,g
,b
,hsv
.h
,hsv
.s
,hsv
.v
);
276 if (hsv
.h
>= ci
->dark
.h
&& hsv
.h
<= ci
->bright
.h
&&
277 hsv
.s
>= ci
->dark
.s
&& hsv
.s
<= ci
->bright
.s
&&
278 hsv
.v
>= ci
->dark
.v
&& hsv
.v
<= ci
->bright
.v
) {
280 } else if (ci
->zapping
) {
281 y
[0] = y
[1] = y
[rowsize
] = y
[rowsize
+ 1] = 16;
291 y
+= picture
->linesize
[0] * 2 - (w_start
- w_end
) * 2;
292 u
+= picture
->linesize
[1] - (w_start
- w_end
);
293 v
+= picture
->linesize
[2] - (w_start
- w_end
);
297 fprintf(stderr
, "Fish: Inrange=%d of %d = %d threshold\n", inrange
, pixcnt
, 1000 * inrange
/ pixcnt
);
299 if (inrange
* 1000 / pixcnt
>= ci
->threshold
) {
304 static int frame_counter
;
305 static int foundfile
;
307 if ((frame_counter
++ % 20) == 0) {
308 /* Check how many files we have */
313 d
= opendir(ci
->dir
);
317 while ((dent
= readdir(d
))) {
318 if (strncmp("fishimg", dent
->d_name
, 7) == 0) {
319 if (strcmp(".ppm", dent
->d_name
+ strlen(dent
->d_name
) - 4) == 0) {
328 if (foundfile
< ci
->file_limit
) {
329 size
= avpicture_get_size(PIX_FMT_RGB24
, width
, height
);
330 buf
= av_malloc(size
);
332 avpicture_fill(&picture1
, buf
, PIX_FMT_RGB24
, width
, height
);
333 if (img_convert(&picture1
, PIX_FMT_RGB24
,
334 picture
, pix_fmt
, width
, height
) >= 0) {
335 /* Write out the PPM file */
340 sprintf(fname
, "%s/fishimg%ld_%lld.ppm", ci
->dir
, time(0), pts
);
341 f
= fopen(fname
, "w");
343 fprintf(f
, "P6 %d %d 255\n", width
, height
);
344 fwrite(buf
, width
* height
* 3, 1, f
);
350 ci
->next_pts
= pts
+ ci
->min_interval
;