libvips/libvips/conversion/bandbool.c

354 lines
8.3 KiB
C

/* bandbool.c --- bool op across image bands
*
* 7/12/12
* - from boolean.c
*/
/*
Copyright (C) 1991-2005 The National Gallery
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.1 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., 51 Franklin Street, Fifth Floor, Boston, MA
02110-1301 USA
*/
/*
These files are distributed with VIPS - http://www.vips.ecs.soton.ac.uk
*/
/*
#define DEBUG
*/
#ifdef HAVE_CONFIG_H
#include <config.h>
#endif /*HAVE_CONFIG_H*/
#include <glib/gi18n-lib.h>
#include <stdio.h>
#include <stdlib.h>
#include <vips/vips.h>
#include "bandary.h"
typedef struct _VipsBandbool {
VipsBandary parent_instance;
VipsImage *in;
VipsOperationBoolean operation;
} VipsBandbool;
typedef VipsBandaryClass VipsBandboolClass;
G_DEFINE_TYPE( VipsBandbool, vips_bandbool, VIPS_TYPE_BANDARY );
static int
vips_bandbool_build( VipsObject *object )
{
VipsObjectClass *class = VIPS_OBJECT_GET_CLASS( object );
VipsBandary *bandary = (VipsBandary *) object;
VipsBandbool *bandbool = (VipsBandbool *) object;
/* << and >> don't work over bands.
*/
if( bandbool->operation == VIPS_OPERATION_BOOLEAN_LSHIFT ||
bandbool->operation == VIPS_OPERATION_BOOLEAN_RSHIFT ) {
vips_error( class->nickname,
_( "operator %s not supported across image bands" ),
vips_enum_nick( VIPS_TYPE_OPERATION_BOOLEAN,
bandbool->operation ) );
return( -1 );
}
if( bandbool->in ) {
if( vips_check_noncomplex( class->nickname, bandbool->in ) )
return( -1 );
bandary->n = 1;
bandary->in = &bandbool->in;
if( bandbool->in->Bands == 1 )
return( vips_bandary_copy( bandary ) );
}
bandary->out_bands = 1;
if( VIPS_OBJECT_CLASS( vips_bandbool_parent_class )->
build( object ) )
return( -1 );
return( 0 );
}
#define SWITCH( I, F, OP ) \
switch( vips_image_get_format( im ) ) { \
case VIPS_FORMAT_UCHAR: I( unsigned char, OP ); break; \
case VIPS_FORMAT_CHAR: I( signed char, OP ); break; \
case VIPS_FORMAT_USHORT: I( unsigned short, OP ); break; \
case VIPS_FORMAT_SHORT: I( signed short, OP ); break; \
case VIPS_FORMAT_UINT: I( unsigned int, OP ); break; \
case VIPS_FORMAT_INT: I( signed int, OP ); break; \
case VIPS_FORMAT_FLOAT: F( float, OP ); break; \
case VIPS_FORMAT_DOUBLE: F( double, OP ); break;\
\
default: \
g_assert_not_reached(); \
}
#define LOOPB( TYPE, OP ) { \
TYPE *p = (TYPE *) in[0]; \
TYPE *q = (TYPE *) out; \
\
for( x = 0; x < width; x++ ) { \
TYPE acc; \
\
acc = p[0]; \
for( b = 1; b < bands; b++ ) \
acc = acc OP p[b]; \
\
q[x] = acc; \
p += bands; \
} \
}
#define FLOOPB( TYPE, OP ) { \
TYPE *p = (TYPE *) in[0]; \
int *q = (int *) out; \
\
for( x = 0; x < width; x++ ) { \
int acc; \
\
acc = (int) p[0]; \
for( b = 1; b < bands; b++ ) \
acc = acc OP ((int) p[b]); \
\
q[x] = acc; \
p += bands; \
} \
}
static void
vips_bandbool_buffer( VipsBandarySequence *seq,
VipsPel *out, VipsPel **in, int width )
{
VipsBandary *bandary = seq->bandary;
VipsBandbool *bandbool = (VipsBandbool *) bandary;
VipsImage *im = bandary->ready[0];
int bands = im->Bands;
int x, b;
switch( bandbool->operation ) {
case VIPS_OPERATION_BOOLEAN_AND:
SWITCH( LOOPB, FLOOPB, & );
break;
case VIPS_OPERATION_BOOLEAN_OR:
SWITCH( LOOPB, FLOOPB, | );
break;
case VIPS_OPERATION_BOOLEAN_EOR:
SWITCH( LOOPB, FLOOPB, ^ );
break;
default:
g_assert_not_reached();
}
}
/* Save a bit of typing.
*/
#define UC VIPS_FORMAT_UCHAR
#define C VIPS_FORMAT_CHAR
#define US VIPS_FORMAT_USHORT
#define S VIPS_FORMAT_SHORT
#define UI VIPS_FORMAT_UINT
#define I VIPS_FORMAT_INT
#define F VIPS_FORMAT_FLOAT
#define X VIPS_FORMAT_COMPLEX
#define D VIPS_FORMAT_DOUBLE
#define DX VIPS_FORMAT_DPCOMPLEX
/* Format conversions for boolean.
*/
static const VipsBandFormat vips_bandbool_format_table[10] = {
/* Band format: UC C US S UI I F X D DX */
/* Promotion: */ UC, C, US, S, UI, I, I, I, I, I
};
static void
vips_bandbool_class_init( VipsBandboolClass *class )
{
GObjectClass *gobject_class = G_OBJECT_CLASS( class );
VipsObjectClass *object_class = (VipsObjectClass *) class;
VipsBandaryClass *bandary_class = VIPS_BANDARY_CLASS( class );
gobject_class->set_property = vips_object_set_property;
gobject_class->get_property = vips_object_get_property;
object_class->nickname = "bandbool";
object_class->description = _( "boolean operation across image bands" );
object_class->build = vips_bandbool_build;
bandary_class->process_line = vips_bandbool_buffer;
bandary_class->format_table = vips_bandbool_format_table;
VIPS_ARG_IMAGE( class, "in", 0,
_( "Input" ),
_( "Input image argument" ),
VIPS_ARGUMENT_REQUIRED_INPUT,
G_STRUCT_OFFSET( VipsBandbool, in ) );
VIPS_ARG_ENUM( class, "boolean", 200,
_( "Operation" ),
_( "Boolean to perform" ),
VIPS_ARGUMENT_REQUIRED_INPUT,
G_STRUCT_OFFSET( VipsBandbool, operation ),
VIPS_TYPE_OPERATION_BOOLEAN,
VIPS_OPERATION_BOOLEAN_AND );
}
static void
vips_bandbool_init( VipsBandbool *bandbool )
{
bandbool->operation = VIPS_OPERATION_BOOLEAN_AND;
}
static int
vips_bandboolv( VipsImage *in, VipsImage **out,
VipsOperationBoolean operation, va_list ap )
{
return( vips_call_split( "bandbool", ap, in, out, operation ) );
}
/**
* vips_bandbool: (method)
* @in: left-hand input #VipsImage
* @out: (out): output #VipsImage
* @boolean: boolean operation to perform
* @...: %NULL-terminated list of optional named arguments
*
* Perform various boolean operations across the bands of an image. For
* example, a three-band uchar image operated on with
* #VIPS_OPERATION_BOOLEAN_AND will produce a one-band uchar image where each
* pixel is the bitwise and of the band elements of the corresponding pixel in
* the input image.
*
* The output image is the same format as the input image for integer
* types. Float types are cast to int before processing. Complex types are not
* supported.
*
* The output image always has one band.
*
* This operation is useful in conjuction with vips_relational(). You can use
* it to see if all image bands match exactly.
*
* See also: vips_boolean_const().
*
* Returns: 0 on success, -1 on error
*/
int
vips_bandbool( VipsImage *in, VipsImage **out,
VipsOperationBoolean boolean, ... )
{
va_list ap;
int result;
va_start( ap, boolean );
result = vips_bandboolv( in, out, boolean, ap );
va_end( ap );
return( result );
}
/**
* vips_bandand: (method)
* @in: left-hand input #VipsImage
* @out: (out): output #VipsImage
* @...: %NULL-terminated list of optional named arguments
*
* Perform #VIPS_OPERATION_BOOLEAN_AND on an image. See
* vips_bandbool().
*
* Returns: 0 on success, -1 on error
*/
int
vips_bandand( VipsImage *in, VipsImage **out, ... )
{
va_list ap;
int result;
va_start( ap, out );
result = vips_bandboolv( in, out, VIPS_OPERATION_BOOLEAN_AND, ap );
va_end( ap );
return( result );
}
/**
* vips_bandor: (method)
* @in: left-hand input #VipsImage
* @out: (out): output #VipsImage
* @...: %NULL-terminated list of optional named arguments
*
* Perform #VIPS_OPERATION_BOOLEAN_OR on an image. See
* vips_bandbool().
*
* Returns: 0 on success, -1 on error
*/
int
vips_bandor( VipsImage *in, VipsImage **out, ... )
{
va_list ap;
int result;
va_start( ap, out );
result = vips_bandboolv( in, out, VIPS_OPERATION_BOOLEAN_OR, ap );
va_end( ap );
return( result );
}
/**
* vips_bandeor: (method)
* @in: left-hand input #VipsImage
* @out: (out): output #VipsImage
* @...: %NULL-terminated list of optional named arguments
*
* Perform #VIPS_OPERATION_BOOLEAN_EOR on an image. See
* vips_bandbool().
*
* Returns: 0 on success, -1 on error
*/
int
vips_bandeor( VipsImage *in, VipsImage **out, ... )
{
va_list ap;
int result;
va_start( ap, out );
result = vips_bandboolv( in, out, VIPS_OPERATION_BOOLEAN_EOR, ap );
va_end( ap );
return( result );
}