[libav.git] / vhook / fish.c

/*
 * Fish Detector Hook
 * Copyright (c) 2002 Philip Gladstone
 *
 * This file implements a fish detector. It is used to see when a 
 * goldfish passes in front of the camera. It does this by counting
 * the number of input pixels that fall within a particular HSV
 * range.
 *
 * It takes a multitude of arguments:
 *
 * -h <num>-<num>    the range of H values that are fish
 * -s <num>-<num>    the range of S values that are fish
 * -v <num>-<num>    the range of V values that are fish
 * -z                zap all non-fish values to black
 * -l <num>          limit the number of saved files to <num>
 * -i <num>          only check frames every <num> seconds
 * -t <num>          the threshold for the amount of fish pixels (range 0-1)
 * -d                turn debugging on
 * -D <directory>    where to put the fish images
 *
 * This library is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2 of the License, or (at your option) any later version.
 *
 * This library is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with this library; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 */
#include <stdlib.h>
#include <fcntl.h>
#include <unistd.h>
#include <stdarg.h>
#include <string.h>
#include <time.h>
#include <stdio.h>
#include <dirent.h>

#include "framehook.h"
#include "dsputil.h"
#include "avformat.h"

#define SCALEBITS 10
#define ONE_HALF  (1 << (SCALEBITS - 1))
#define FIX(x)    ((int) ((x) * (1<<SCALEBITS) + 0.5))

#define YUV_TO_RGB1_CCIR(cb1, cr1)\
{\
    cb = (cb1) - 128;\
    cr = (cr1) - 128;\
    r_add = FIX(1.40200*255.0/224.0) * cr + ONE_HALF;\
    g_add = - FIX(0.34414*255.0/224.0) * cb - FIX(0.71414*255.0/224.0) * cr + \
                    ONE_HALF;\
    b_add = FIX(1.77200*255.0/224.0) * cb + ONE_HALF;\
}

#define YUV_TO_RGB2_CCIR(r, g, b, y1)\
{\
    yt = ((y1) - 16) * FIX(255.0/219.0);\
    r = cm[(yt + r_add) >> SCALEBITS];\
    g = cm[(yt + g_add) >> SCALEBITS];\
    b = cm[(yt + b_add) >> SCALEBITS];\
}


typedef struct {
    int h;  /* 0 .. 360 */
    int s;  /* 0 .. 255 */
    int v;  /* 0 .. 255 */
} HSV;
              
typedef struct {
    int zapping;
    int threshold;
    HSV dark, bright;
    char *dir;
    int file_limit;
    int debug;
    int min_interval;
    int64_t next_pts;
    int inset;
    int min_width;
} ContextInfo;

static void dorange(const char *s, int *first, int *second, int maxval)
{
    sscanf(s, "%d-%d", first, second);
    if (*first > maxval)
        *first = maxval;
    if (*second > maxval)
        *second = maxval;
}

void Release(void *ctx)
{
    if (ctx)
        av_free(ctx);
}

int Configure(void **ctxp, int argc, char *argv[])
{
    ContextInfo *ci;
    int c;

    *ctxp = av_mallocz(sizeof(ContextInfo));
    ci = (ContextInfo *) *ctxp;

    optind = 0;

    ci->dir = "/tmp";
    ci->threshold = 100;
    ci->file_limit = 100;
    ci->min_interval = 1000000;
    ci->inset = 10;     /* Percent */

    while ((c = getopt(argc, argv, "w:i:dh:s:v:zl:t:D:")) > 0) {
        switch (c) {
            case 'h':
                dorange(optarg, &ci->dark.h, &ci->bright.h, 360);
                break;
            case 's':
                dorange(optarg, &ci->dark.s, &ci->bright.s, 255);
                break;
            case 'v':
                dorange(optarg, &ci->dark.v, &ci->bright.v, 255);
                break;
            case 'z':
                ci->zapping = 1;
                break;
            case 'l':
                ci->file_limit = atoi(optarg);
                break;
            case 'i':
                ci->min_interval = 1000000 * atof(optarg);
                break;
            case 't':
                ci->threshold = atof(optarg) * 1000;
                if (ci->threshold > 1000 || ci->threshold < 0) {
                    fprintf(stderr, "Invalid threshold value '%s' (range is 0-1)\n", optarg);
                    return -1;
                }
                break;
            case 'w':
                ci->min_width = atoi(optarg);
                break;
            case 'd':
                ci->debug++;
                break;
            case 'D':
                ci->dir = av_strdup(optarg);
                break;
            default:
                fprintf(stderr, "Unrecognized argument '%s'\n", argv[optind]);
                return -1;
        }
    }

    fprintf(stderr, "Fish detector configured:\n");
    fprintf(stderr, "    HSV range: %d,%d,%d - %d,%d,%d\n",
                        ci->dark.h,
                        ci->dark.s,
                        ci->dark.v,
                        ci->bright.h,
                        ci->bright.s,
                        ci->bright.v);
    fprintf(stderr, "    Threshold is %d%% pixels\n", ci->threshold / 10);


    return 0;
}

static void get_hsv(HSV *hsv, int r, int g, int b)
{
    int i, v, x, f;
         
    x = (r < g) ? r : g;
    if (b < x)
        x = b;
    v = (r > g) ? r : g;
    if (b > v)
        v = b;
          
    if (v == x) {
        hsv->h = 0;
        hsv->s = 0;
        hsv->v = v;
        return;
    }
       
    if (r == v) {
        f = g - b;
        i = 0;
    } else if (g == v) {
        f = b - r;
        i = 2 * 60;
    } else {
        f = r - g;
        i = 4 * 60;
    }
        
    hsv->h = i + (60 * f) / (v - x);
    if (hsv->h < 0)
        hsv->h += 360;

    hsv->s = (255 * (v - x)) / v;
    hsv->v = v;
         
    return;
}                                                                               

void Process(void *ctx, AVPicture *picture, enum PixelFormat pix_fmt, int width, int height, int64_t pts)
{
    ContextInfo *ci = (ContextInfo *) ctx;
    uint8_t *cm = cropTbl + MAX_NEG_CROP;                                         
    int rowsize = picture->linesize[0];

#if 0
    printf("pix_fmt = %d, width = %d, pts = %lld, ci->next_pts = %lld\n",
        pix_fmt, width, pts, ci->next_pts);
#endif

    if (pts < ci->next_pts)
        return;

    if (width < ci->min_width)
        return;

    ci->next_pts = pts + 1000000;    

    if (pix_fmt == PIX_FMT_YUV420P) {
        uint8_t *y, *u, *v;
        int width2 = width >> 1;
        int inrange = 0;
        int pixcnt;
        int h;
        int h_start, h_end;
        int w_start, w_end;

        h_end = 2 * ((ci->inset * height) / 200);
        h_start = height - h_end;

        w_end = (ci->inset * width2) / 100;
        w_start = width2 - w_end;

        pixcnt = ((h_start - h_end) >> 1) * (w_start - w_end);

        y = picture->data[0] + h_end * picture->linesize[0] + w_end * 2;
        u = picture->data[1] + h_end * picture->linesize[1] / 2 + w_end;
        v = picture->data[2] + h_end * picture->linesize[2] / 2 + w_end;

        for (h = h_start; h > h_end; h -= 2) {
            int w;

            for (w = w_start; w > w_end; w--) {
                unsigned int r,g,b;
                HSV hsv;
                int cb, cr, yt, r_add, g_add, b_add;

                YUV_TO_RGB1_CCIR(u[0], v[0]);
                YUV_TO_RGB2_CCIR(r, g, b, y[0]);

                get_hsv(&hsv, r, g, b);

                if (ci->debug > 1) 
                    fprintf(stderr, "(%d,%d,%d) -> (%d,%d,%d)\n",
                        r,g,b,hsv.h,hsv.s,hsv.v);


                if (hsv.h >= ci->dark.h && hsv.h <= ci->bright.h &&
                    hsv.s >= ci->dark.s && hsv.s <= ci->bright.s &&
                    hsv.v >= ci->dark.v && hsv.v <= ci->bright.v) {            
                    inrange++;
                } else if (ci->zapping) {
                    y[0] = y[1] = y[rowsize] = y[rowsize + 1] = 16;
                    u[0] = 128;
                    v[0] = 128;
                }

                y+= 2;
                u++;
                v++;
            }

            y += picture->linesize[0] * 2 - (w_start - w_end) * 2;
            u += picture->linesize[1] - (w_start - w_end);
            v += picture->linesize[2] - (w_start - w_end);
        }

        if (ci->debug) 
            fprintf(stderr, "Fish: Inrange=%d of %d = %d threshold\n", inrange, pixcnt, 1000 * inrange / pixcnt);

        if (inrange * 1000 / pixcnt >= ci->threshold) {
            /* Save to file */
            int size;
            char *buf;
            AVPicture picture1;
            static int frame_counter;
            static int foundfile;

            if ((frame_counter++ % 20) == 0) {
                /* Check how many files we have */
                DIR *d;

                foundfile = 0;

                d = opendir(ci->dir);
                if (d) {
                    struct dirent *dent;

                    while ((dent = readdir(d))) {
                        if (strncmp("fishimg", dent->d_name, 7) == 0) {
                            if (strcmp(".ppm", dent->d_name + strlen(dent->d_name) - 4) == 0) {
                                foundfile++;
                            }
                        }
                    }
                    closedir(d);
                }
            }

            if (foundfile < ci->file_limit) {
                size = avpicture_get_size(PIX_FMT_RGB24, width, height);
                buf = av_malloc(size);

                avpicture_fill(&picture1, buf, PIX_FMT_RGB24, width, height);
                if (img_convert(&picture1, PIX_FMT_RGB24, 
                                picture, pix_fmt, width, height) >= 0) {
                    /* Write out the PPM file */

                    FILE *f;
                    char fname[256];

                    snprintf(fname, sizeof(fname), "%s/fishimg%ld_%"PRId64".ppm", ci->dir, (long)(av_gettime() / 1000000), pts);
                    f = fopen(fname, "w");
                    if (f) {
                        fprintf(f, "P6 %d %d 255\n", width, height);
                        fwrite(buf, width * height * 3, 1, f);
                        fclose(f);
                    }
                }

                av_free(buf);
                ci->next_pts = pts + ci->min_interval;    
            }
        }
    }
}
Commit	Line	Data
26b4bb70 PG	1	/*
	2	* Fish Detector Hook
	3	* Copyright (c) 2002 Philip Gladstone
	4	*
	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
	8	* range.
	9	*
	10	* It takes a multitude of arguments:
	11	*
	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
	21	*
	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.
	26	*
	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.
	31	*
	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
	35	*/
	36	#include <stdlib.h>
	37	#include <fcntl.h>
47930f09	38	#include <unistd.h>
26b4bb70 PG	39	#include <stdarg.h>
26b4bb70 PG	40	#include <string.h>
b55e4ef4	41	#include <time.h>
26b4bb70 PG	42	#include <stdio.h>
	43	#include <dirent.h>
	44
	45	#include "framehook.h"
	46	#include "dsputil.h"
cf344ec7	47	#include "avformat.h"
26b4bb70	48
25a2ee7d PG	49	#define SCALEBITS 10
	50	#define ONE_HALF (1 << (SCALEBITS - 1))
	51	#define FIX(x) ((int) ((x) * (1<<SCALEBITS) + 0.5))
	52
	53	#define YUV_TO_RGB1_CCIR(cb1, cr1)\
	54	{\
	55	cb = (cb1) - 128;\
	56	cr = (cr1) - 128;\
	57	r_add = FIX(1.40200255.0/224.0) cr + ONE_HALF;\
	58	g_add = - FIX(0.34414255.0/224.0) cb - FIX(0.71414255.0/224.0) cr + \
	59	ONE_HALF;\
	60	b_add = FIX(1.77200255.0/224.0) cb + ONE_HALF;\
	61	}
	62
	63	#define YUV_TO_RGB2_CCIR(r, g, b, y1)\
	64	{\
	65	yt = ((y1) - 16) * FIX(255.0/219.0);\
	66	r = cm[(yt + r_add) >> SCALEBITS];\
	67	g = cm[(yt + g_add) >> SCALEBITS];\
	68	b = cm[(yt + b_add) >> SCALEBITS];\
	69	}
	70
26b4bb70	71
25a2ee7d	72
26b4bb70 PG	73
	74	typedef struct {
	75	int h; /* 0 .. 360 */
	76	int s; /* 0 .. 255 */
	77	int v; /* 0 .. 255 */
	78	} HSV;
	79
	80	typedef struct {
	81	int zapping;
	82	int threshold;
	83	HSV dark, bright;
	84	char *dir;
	85	int file_limit;
	86	int debug;
	87	int min_interval;
0c1a9eda	88	int64_t next_pts;
26b4bb70 PG	89	int inset;
	90	int min_width;
	91	} ContextInfo;
	92
	93	static void dorange(const char s, int first, int *second, int maxval)
	94	{
	95	sscanf(s, "%d-%d", first, second);
	96	if (*first > maxval)
	97	*first = maxval;
	98	if (*second > maxval)
	99	*second = maxval;
	100	}
	101
6c11d48c PG	102	void Release(void *ctx)
	103	{
	104	if (ctx)
	105	av_free(ctx);
	106	}
26b4bb70 PG	107
	108	int Configure(void *ctxp, int argc, char argv[])
	109	{
	110	ContextInfo *ci;
	111	int c;
	112
	113	*ctxp = av_mallocz(sizeof(ContextInfo));
	114	ci = (ContextInfo ) ctxp;
	115
	116	optind = 0;
	117
	118	ci->dir = "/tmp";
d9974a48	119	ci->threshold = 100;
26b4bb70 PG	120	ci->file_limit = 100;
	121	ci->min_interval = 1000000;
	122	ci->inset = 10; /* Percent */
	123
	124	while ((c = getopt(argc, argv, "w:i:dh:s:v:zl:t:D:")) > 0) {
	125	switch (c) {
	126	case 'h':
	127	dorange(optarg, &ci->dark.h, &ci->bright.h, 360);
	128	break;
	129	case 's':
	130	dorange(optarg, &ci->dark.s, &ci->bright.s, 255);
	131	break;
	132	case 'v':
	133	dorange(optarg, &ci->dark.v, &ci->bright.v, 255);
	134	break;
	135	case 'z':
	136	ci->zapping = 1;
	137	break;
	138	case 'l':
	139	ci->file_limit = atoi(optarg);
	140	break;
	141	case 'i':
	142	ci->min_interval = 1000000 * atof(optarg);
	143	break;
	144	case 't':
	145	ci->threshold = atof(optarg) * 1000;
d9974a48 PG	146	if (ci->threshold > 1000 \|\| ci->threshold < 0) {
	147	fprintf(stderr, "Invalid threshold value '%s' (range is 0-1)\n", optarg);
	148	return -1;
	149	}
26b4bb70 PG	150	break;
	151	case 'w':
	152	ci->min_width = atoi(optarg);
	153	break;
	154	case 'd':
	155	ci->debug++;
	156	break;
	157	case 'D':
e9a9e0c2	158	ci->dir = av_strdup(optarg);
26b4bb70 PG	159	break;
	160	default:
	161	fprintf(stderr, "Unrecognized argument '%s'\n", argv[optind]);
	162	return -1;
	163	}
	164	}
	165
	166	fprintf(stderr, "Fish detector configured:\n");
	167	fprintf(stderr, " HSV range: %d,%d,%d - %d,%d,%d\n",
	168	ci->dark.h,
	169	ci->dark.s,
	170	ci->dark.v,
	171	ci->bright.h,
	172	ci->bright.s,
	173	ci->bright.v);
25a2ee7d PG	174	fprintf(stderr, " Threshold is %d%% pixels\n", ci->threshold / 10);
25a2ee7d PG	175
26b4bb70 PG	176
	177	return 0;
	178	}
	179
	180	static void get_hsv(HSV *hsv, int r, int g, int b)
	181	{
	182	int i, v, x, f;
	183
	184	x = (r < g) ? r : g;
	185	if (b < x)
	186	x = b;
	187	v = (r > g) ? r : g;
	188	if (b > v)
	189	v = b;
	190
	191	if (v == x) {
	192	hsv->h = 0;
	193	hsv->s = 0;
	194	hsv->v = v;
	195	return;
	196	}
	197
	198	if (r == v) {
	199	f = g - b;
	200	i = 0;
	201	} else if (g == v) {
	202	f = b - r;
	203	i = 2 * 60;
	204	} else {
	205	f = r - g;
	206	i = 4 * 60;
	207	}
	208
	209	hsv->h = i + (60 * f) / (v - x);
	210	if (hsv->h < 0)
	211	hsv->h += 360;
	212
	213	hsv->s = (255 * (v - x)) / v;
	214	hsv->v = v;
	215
	216	return;
	217	}
	218
0c1a9eda	219	void Process(void ctx, AVPicture picture, enum PixelFormat pix_fmt, int width, int height, int64_t pts)
26b4bb70 PG	220	{
26b4bb70 PG	221	ContextInfo ci = (ContextInfo ) ctx;
0c1a9eda	222	uint8_t *cm = cropTbl + MAX_NEG_CROP;
26b4bb70 PG	223	int rowsize = picture->linesize[0];
26b4bb70 PG	224
d9974a48 PG	225	#if 0
	226	printf("pix_fmt = %d, width = %d, pts = %lld, ci->next_pts = %lld\n",
	227	pix_fmt, width, pts, ci->next_pts);
	228	#endif
	229
26b4bb70 PG	230	if (pts < ci->next_pts)
	231	return;
	232
	233	if (width < ci->min_width)
	234	return;
	235
	236	ci->next_pts = pts + 1000000;
	237
	238	if (pix_fmt == PIX_FMT_YUV420P) {
0c1a9eda	239	uint8_t y, u, *v;
26b4bb70 PG	240	int width2 = width >> 1;
	241	int inrange = 0;
	242	int pixcnt;
	243	int h;
	244	int h_start, h_end;
	245	int w_start, w_end;
	246
	247	h_end = 2 * ((ci->inset * height) / 200);
	248	h_start = height - h_end;
	249
	250	w_end = (ci->inset * width2) / 100;
	251	w_start = width2 - w_end;
	252
	253	pixcnt = ((h_start - h_end) >> 1) * (w_start - w_end);
	254
d9974a48	255	y = picture->data[0] + h_end * picture->linesize[0] + w_end * 2;
25a2ee7d PG	256	u = picture->data[1] + h_end * picture->linesize[1] / 2 + w_end;
25a2ee7d PG	257	v = picture->data[2] + h_end * picture->linesize[2] / 2 + w_end;
26b4bb70 PG	258
	259	for (h = h_start; h > h_end; h -= 2) {
	260	int w;
	261
	262	for (w = w_start; w > w_end; w--) {
25a2ee7d	263	unsigned int r,g,b;
26b4bb70	264	HSV hsv;
25a2ee7d	265	int cb, cr, yt, r_add, g_add, b_add;
26b4bb70	266
25a2ee7d PG	267	YUV_TO_RGB1_CCIR(u[0], v[0]);
25a2ee7d PG	268	YUV_TO_RGB2_CCIR(r, g, b, y[0]);
26b4bb70 PG	269
	270	get_hsv(&hsv, r, g, b);
	271
	272	if (ci->debug > 1)
	273	fprintf(stderr, "(%d,%d,%d) -> (%d,%d,%d)\n",
	274	r,g,b,hsv.h,hsv.s,hsv.v);
	275
	276
	277	if (hsv.h >= ci->dark.h && hsv.h <= ci->bright.h &&
	278	hsv.s >= ci->dark.s && hsv.s <= ci->bright.s &&
	279	hsv.v >= ci->dark.v && hsv.v <= ci->bright.v) {
	280	inrange++;
	281	} else if (ci->zapping) {
25a2ee7d PG	282	y[0] = y[1] = y[rowsize] = y[rowsize + 1] = 16;
	283	u[0] = 128;
	284	v[0] = 128;
26b4bb70 PG	285	}
	286
	287	y+= 2;
	288	u++;
	289	v++;
	290	}
	291
d9974a48 PG	292	y += picture->linesize[0] * 2 - (w_start - w_end) * 2;
	293	u += picture->linesize[1] - (w_start - w_end);
	294	v += picture->linesize[2] - (w_start - w_end);
26b4bb70 PG	295	}
26b4bb70 PG	296
d9974a48 PG	297	if (ci->debug)
	298	fprintf(stderr, "Fish: Inrange=%d of %d = %d threshold\n", inrange, pixcnt, 1000 * inrange / pixcnt);
	299
26b4bb70 PG	300	if (inrange * 1000 / pixcnt >= ci->threshold) {
	301	/* Save to file */
	302	int size;
	303	char *buf;
	304	AVPicture picture1;
	305	static int frame_counter;
	306	static int foundfile;
	307
26b4bb70 PG	308	if ((frame_counter++ % 20) == 0) {
	309	/* Check how many files we have */
	310	DIR *d;
	311
	312	foundfile = 0;
	313
	314	d = opendir(ci->dir);
	315	if (d) {
	316	struct dirent *dent;
	317
	318	while ((dent = readdir(d))) {
	319	if (strncmp("fishimg", dent->d_name, 7) == 0) {
	320	if (strcmp(".ppm", dent->d_name + strlen(dent->d_name) - 4) == 0) {
	321	foundfile++;
	322	}
	323	}
	324	}
	325	closedir(d);
	326	}
	327	}
	328
	329	if (foundfile < ci->file_limit) {
	330	size = avpicture_get_size(PIX_FMT_RGB24, width, height);
	331	buf = av_malloc(size);
	332
	333	avpicture_fill(&picture1, buf, PIX_FMT_RGB24, width, height);
	334	if (img_convert(&picture1, PIX_FMT_RGB24,
	335	picture, pix_fmt, width, height) >= 0) {
	336	/* Write out the PPM file */
	337
	338	FILE *f;
	339	char fname[256];
	340
4733abcb	341	snprintf(fname, sizeof(fname), "%s/fishimg%ld_%"PRId64".ppm", ci->dir, (long)(av_gettime() / 1000000), pts);
26b4bb70 PG	342	f = fopen(fname, "w");
	343	if (f) {
	344	fprintf(f, "P6 %d %d 255\n", width, height);
	345	fwrite(buf, width * height * 3, 1, f);
	346	fclose(f);
	347	}
	348	}
	349
	350	av_free(buf);
	351	ci->next_pts = pts + ci->min_interval;
	352	}
	353	}
	354	}
	355	}
	356