sergiotarxz
/
libvips
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity
libvips/libvips/colour/icc_transform.c

1514 lines
38 KiB
C
Raw Permalink Blame History

/* Transform images with little cms
*
* 26/4/02 JC
* 26/8/05
* - attach profiles and intents to output images
* - added im_icc_import_embedded() to import with an embedded profile
* 12/5/06
* - lock around cmsDoTransform
* 23/1/07
* - set RGB16 on 16-bit RGB export
* 6/4/09
* - catch lcms error messages
* 2/11/09
* - gtkdoc
* - small cleanups
* - call attach_profile() before im_wrapone() so the profile will get
* written if we are wrinting to a file
* 2/8/10
* - add lcms2
* 12/7/11
* - import and export cast @in to an appropriate format for you
* 25/9/12
* - redo as a class
* 14/5/13
* - import and export would segv on very wide images
* 12/11/13
* - support XYZ as an alternative PCS
* 10/9/14
* - support GRAY as an input and output space
* 29/9/14
* - check input profiles for compatibility with the input image, thanks
* James
* 26/6/15
* - better profile sanity checking for icc import
* 2/8/17
* - remove lcms1 support, it was untested
* 10/10/17
* - more input profile sanity tests
* 8/3/18
* - attach fallback profile on import if we used it
* 28/12/18
* - remove warning messages from vips_icc_is_compatible_profile() since
* they can be triggered under normal circumstances
* 17/4/19 kleisauke
* - better rejection of broken embedded profiles
* 29/3/21 [hanssonrickard]
* - add black_point_compensation
*/
/*
This file is part of VIPS.
VIPS 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 program 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 program; 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>
#ifdef HAVE_LCMS2
#include <stdio.h>
#include <math.h>
/* Has to be before VIPS to avoid nameclashes.
*/
#include <lcms2.h>
#include <vips/vips.h>
#include "pcolour.h"
/* Call lcms with up to this many pixels at once.
*/
#define PIXEL_BUFFER_SIZE (10000)
/**
* VipsIntent:
* @VIPS_INTENT_PERCEPTUAL: perceptual rendering intent
* @VIPS_INTENT_RELATIVE: relative colorimetric rendering intent
* @VIPS_INTENT_SATURATION: saturation rendering intent
* @VIPS_INTENT_ABSOLUTE: absolute colorimetric rendering intent
*
* The rendering intent. #VIPS_INTENT_ABSOLUTE is best for
* scientific work, #VIPS_INTENT_RELATIVE is usually best for
* accurate communication with other imaging libraries.
*/
/**
* VipsPCS:
* @VIPS_PCS_LAB: use CIELAB D65 as the Profile Connection Space
* @VIPS_PCS_XYZ: use XYZ as the Profile Connection Space
*
* Pick a Profile Connection Space for vips_icc_import() and
* vips_icc_export(). LAB is usually best, XYZ can be more convenient in some
* cases.
*/
/**
* vips_icc_present:
*
* VIPS can optionally be built without the ICC library. Use this function to
* test for its availability.
*
* Returns: non-zero if the ICC library is present.
*/
int
vips_icc_present( void )
{
return( 1 );
}
#define VIPS_TYPE_ICC (vips_icc_get_type())
#define VIPS_ICC( obj ) \
(G_TYPE_CHECK_INSTANCE_CAST( (obj), \
VIPS_TYPE_ICC, VipsIcc ))
#define VIPS_ICC_CLASS( klass ) \
(G_TYPE_CHECK_CLASS_CAST( (klass), \
VIPS_TYPE_ICC, VipsIccClass))
#define VIPS_IS_ICC( obj ) \
(G_TYPE_CHECK_INSTANCE_TYPE( (obj), VIPS_TYPE_ICC ))
#define VIPS_IS_ICC_CLASS( klass ) \
(G_TYPE_CHECK_CLASS_TYPE( (klass), VIPS_TYPE_ICC ))
#define VIPS_ICC_GET_CLASS( obj ) \
(G_TYPE_INSTANCE_GET_CLASS( (obj), \
VIPS_TYPE_ICC, VipsIccClass ))
typedef struct _VipsIcc {
VipsColourCode parent_instance;
VipsIntent intent;
VipsPCS pcs;
int depth;
gboolean black_point_compensation;
VipsBlob *in_blob;
cmsHPROFILE in_profile;
VipsBlob *out_blob;
cmsHPROFILE out_profile;
cmsUInt32Number in_icc_format;
cmsUInt32Number out_icc_format;
cmsHTRANSFORM trans;
gboolean non_standard_input_profile;
} VipsIcc;
typedef VipsColourCodeClass VipsIccClass;
G_DEFINE_ABSTRACT_TYPE( VipsIcc, vips_icc, VIPS_TYPE_COLOUR_CODE );
/* Error from lcms.
*/
static void
icc_error( cmsContext context, cmsUInt32Number code, const char *text )
{
vips_error( "VipsIcc", "%s", text );
}
static void
vips_icc_dispose( GObject *gobject )
{
VipsIcc *icc = (VipsIcc *) gobject;
VIPS_FREEF( cmsDeleteTransform, icc->trans );
VIPS_FREEF( cmsCloseProfile, icc->in_profile );
VIPS_FREEF( cmsCloseProfile, icc->out_profile );
if( icc->in_blob ) {
vips_area_unref( (VipsArea *) icc->in_blob );
icc->in_blob = NULL;
}
if( icc->out_blob ) {
vips_area_unref( (VipsArea *) icc->out_blob );
icc->out_blob = NULL;
}
G_OBJECT_CLASS( vips_icc_parent_class )->dispose( gobject );
}
/* Is a profile just a pcs stub.
*/
static gboolean
is_pcs( cmsHPROFILE profile )
{
return( cmsGetColorSpace( profile ) == cmsSigLabData ||
cmsGetColorSpace( profile ) == cmsSigXYZData );
}
/* Info for all supported lcms profile signatures.
*/
typedef struct _VipsIccInfo {
int signature;
int bands;
guint lcms_type8;
guint lcms_type16;
} VipsIccInfo;
static VipsIccInfo vips_icc_info_table[] = {
{ cmsSigGrayData, 1, TYPE_GRAY_8, TYPE_GRAY_16 },
{ cmsSigRgbData, 3, TYPE_RGB_8, TYPE_RGB_16 },
{ cmsSigLabData, 3, TYPE_Lab_8, TYPE_Lab_16 },
{ cmsSigXYZData, 3, -1, TYPE_XYZ_16 },
{ cmsSigCmykData, 4, TYPE_CMYK_8, TYPE_CMYK_16 },
{ cmsSig4colorData, 4, TYPE_CMYK_8, TYPE_CMYK_16 },
{ cmsSig5colorData, 5, TYPE_CMYK5_8, TYPE_CMYK5_16 },
{ cmsSig6colorData, 6, TYPE_CMYK6_8, TYPE_CMYK6_16 },
{ cmsSig7colorData, 7, TYPE_CMYK7_8, TYPE_CMYK7_16 },
{ cmsSig8colorData, 8, TYPE_CMYK8_8, TYPE_CMYK8_16 },
{ cmsSig9colorData, 9, TYPE_CMYK9_8, TYPE_CMYK9_16 },
{ cmsSig10colorData, 10, TYPE_CMYK10_8, TYPE_CMYK10_16 },
{ cmsSig11colorData, 11, TYPE_CMYK11_8, TYPE_CMYK11_16 },
{ cmsSig12colorData, 12, TYPE_CMYK12_8, TYPE_CMYK12_16 },
};
static VipsIccInfo *
vips_icc_info( int signature )
{
int i;
for( i = 0; i < VIPS_NUMBER( vips_icc_info_table ); i++ )
if( vips_icc_info_table[i].signature == signature )
return( &vips_icc_info_table[i] );
return( NULL );
}
static int
vips_icc_build( VipsObject *object )
{
VipsObjectClass *class = VIPS_OBJECT_GET_CLASS( object );
VipsColour *colour = (VipsColour *) object;
VipsColourCode *code = (VipsColourCode *) object;
VipsIcc *icc = (VipsIcc *) object;
cmsUInt32Number flags;
if( icc->depth != 8 &&
icc->depth != 16 ) {
vips_error( class->nickname,
"%s", _( "depth must be 8 or 16" ) );
return( -1 );
}
if( icc->in_profile &&
code->in ) {
int signature;
VipsIccInfo *info;
signature = cmsGetColorSpace( icc->in_profile );
if( !(info = vips_icc_info( signature )) ) {
vips_error( class->nickname,
_( "unimplemented input color space 0x%x" ),
signature );
return( -1 );
}
colour->input_bands = info->bands;
switch( signature ) {
case cmsSigGrayData:
code->input_format =
code->in->BandFmt == VIPS_FORMAT_USHORT ?
VIPS_FORMAT_USHORT : VIPS_FORMAT_UCHAR;
icc->in_icc_format =
code->in->BandFmt == VIPS_FORMAT_USHORT ?
info->lcms_type16 : info->lcms_type8;
break;
case cmsSigRgbData:
code->input_format =
code->in->BandFmt == VIPS_FORMAT_USHORT ?
VIPS_FORMAT_USHORT : VIPS_FORMAT_UCHAR;
icc->in_icc_format =
code->in->BandFmt == VIPS_FORMAT_USHORT ?
info->lcms_type16 : info->lcms_type8;
break;
case cmsSigLabData:
code->input_format = VIPS_FORMAT_FLOAT;
code->input_interpretation =
VIPS_INTERPRETATION_LAB;
icc->in_icc_format = info->lcms_type16;
break;
case cmsSigXYZData:
code->input_format = VIPS_FORMAT_FLOAT;
code->input_interpretation =
VIPS_INTERPRETATION_XYZ;
icc->in_icc_format = info->lcms_type16;
break;
case cmsSigCmykData:
case cmsSig5colorData:
case cmsSig6colorData:
case cmsSig7colorData:
case cmsSig8colorData:
case cmsSig9colorData:
case cmsSig10colorData:
case cmsSig11colorData:
case cmsSig12colorData:
/* Treat as forms of CMYK.
*/
info = vips_icc_info(
cmsGetColorSpace( icc->in_profile ) );
code->input_format =
code->in->BandFmt == VIPS_FORMAT_USHORT ?
VIPS_FORMAT_USHORT : VIPS_FORMAT_UCHAR;
icc->in_icc_format =
code->in->BandFmt == VIPS_FORMAT_USHORT ?
info->lcms_type16 : info->lcms_type8;
break;
default:
g_assert_not_reached();
return( -1 );
}
}
if( icc->out_profile ) {
int signature;
VipsIccInfo *info;
signature = cmsGetColorSpace( icc->out_profile );
if( !(info = vips_icc_info( signature )) ) {
vips_error( class->nickname,
_( "unimplemented output color space 0x%x" ),
signature );
return( -1 );
}
colour->bands = info->bands;
switch( signature ) {
case cmsSigGrayData:
colour->interpretation =
icc->depth == 8 ?
VIPS_INTERPRETATION_B_W :
VIPS_INTERPRETATION_GREY16;
colour->format =
icc->depth == 8 ?
VIPS_FORMAT_UCHAR : VIPS_FORMAT_USHORT;
icc->out_icc_format =
icc->depth == 16 ?
info->lcms_type16 : info->lcms_type8;
break;
case cmsSigRgbData:
colour->interpretation =
icc->depth == 8 ?
VIPS_INTERPRETATION_sRGB :
VIPS_INTERPRETATION_RGB16;
colour->format =
icc->depth == 8 ?
VIPS_FORMAT_UCHAR : VIPS_FORMAT_USHORT;
icc->out_icc_format =
icc->depth == 16 ?
info->lcms_type16 : info->lcms_type8;
break;
case cmsSigLabData:
colour->interpretation = VIPS_INTERPRETATION_LAB;
colour->format = VIPS_FORMAT_FLOAT;
icc->out_icc_format = info->lcms_type16;
break;
case cmsSigXYZData:
colour->interpretation = VIPS_INTERPRETATION_XYZ;
colour->format = VIPS_FORMAT_FLOAT;
icc->out_icc_format = info->lcms_type16;
break;
case cmsSigCmykData:
case cmsSig5colorData:
case cmsSig6colorData:
case cmsSig7colorData:
case cmsSig8colorData:
case cmsSig9colorData:
case cmsSig10colorData:
case cmsSig11colorData:
case cmsSig12colorData:
/* Treat as forms of CMYK.
*/
colour->interpretation = VIPS_INTERPRETATION_CMYK;
colour->format =
icc->depth == 8 ?
VIPS_FORMAT_UCHAR : VIPS_FORMAT_USHORT;
icc->out_icc_format =
icc->depth == 16 ?
info->lcms_type16 : info->lcms_type8;
break;
default:
g_assert_not_reached();
return( -1 );
}
}
/* At least one must be a device profile.
*/
if( icc->in_profile &&
icc->out_profile &&
is_pcs( icc->in_profile ) &&
is_pcs( icc->out_profile ) ) {
vips_error( class->nickname,
"%s", _( "no device profile" ) );
return( -1 );
}
/* Use cmsFLAGS_NOCACHE to disable the 1-pixel cache and make
* calling cmsDoTransform() from multiple threads safe.
*/
flags = cmsFLAGS_NOCACHE;
if( icc->black_point_compensation )
flags |= cmsFLAGS_BLACKPOINTCOMPENSATION;
if( !(icc->trans = cmsCreateTransform(
icc->in_profile, icc->in_icc_format,
icc->out_profile, icc->out_icc_format,
icc->intent, flags )) )
return( -1 );
if( VIPS_OBJECT_CLASS( vips_icc_parent_class )->
build( object ) )
return( -1 );
return( 0 );
}
/* Get from an image.
*/
static VipsBlob *
vips_icc_get_profile_image( VipsImage *image )
{
const void *data;
size_t size;
if( !vips_image_get_typeof( image, VIPS_META_ICC_NAME ) )
return( NULL );
if( vips_image_get_blob( image, VIPS_META_ICC_NAME, &data, &size ) )
return( NULL );
return( vips_blob_new( NULL, data, size ) );
}
#ifdef DEBUG
static void
vips_icc_print_profile( const char *name, cmsHPROFILE profile )
{
static const cmsInfoType info_types[] = {
cmsInfoDescription,
cmsInfoManufacturer,
cmsInfoModel,
cmsInfoCopyright
};
static const char *info_names[] = {
"description",
"manufacturer",
"model",
"copyright"
};
int i;
cmsUInt32Number n_bytes;
cmsUInt32Number n_intents;
cmsUInt32Number *intent_codes;
char **intent_descriptions;
printf( "icc profile %s: %p\n", name, profile );
for( i = 0; i < VIPS_NUMBER( info_types ); i++ ) {
if( (n_bytes = cmsGetProfileInfoASCII( profile,
info_types[i], "en", "US",
NULL, 0 )) ) {
char *buffer;
buffer = VIPS_ARRAY( NULL, n_bytes, char );
(void) cmsGetProfileInfoASCII( profile,
info_types[i], "en", "US",
buffer, n_bytes );
printf( "%s: %s\n", info_names[i], buffer );
g_free( buffer );
}
}
printf( "profile class: %#x\n", cmsGetDeviceClass( profile ) );
{
cmsColorSpaceSignature pcs = cmsGetPCS( profile );
guint pcs_as_be = GUINT_TO_BE( pcs );
char name[5];
vips_strncpy( name, (const char *) &pcs_as_be, 5 );
printf( "PCS: %#x (%s)\n", pcs, name );
}
printf( "matrix shaper: %d\n", cmsIsMatrixShaper( profile ) );
printf( "version: %g\n", cmsGetProfileVersion( profile ) );
n_intents = cmsGetSupportedIntents( 0, NULL, NULL );
printf( "n_intents = %u\n", n_intents );
intent_codes = VIPS_ARRAY( NULL, n_intents, cmsUInt32Number );
intent_descriptions = VIPS_ARRAY( NULL, n_intents, char * );
(void) cmsGetSupportedIntents( n_intents,
intent_codes, intent_descriptions );
for( i = 0; i < n_intents; i++ ) {
printf( " %#x: %s, in CLUT = %d, out CLUT = %d\n",
intent_codes[i], intent_descriptions[i],
cmsIsCLUT( profile,
intent_codes[i], LCMS_USED_AS_INPUT ),
cmsIsCLUT( profile,
intent_codes[i], LCMS_USED_AS_OUTPUT ) );
}
g_free( intent_codes );
g_free( intent_descriptions );
}
#endif /*DEBUG*/
/* Load a profile from a blob and check compatibility with image, intent and
* direction.
*
* Don't set any errors since this is used to test compatibility.
*/
static cmsHPROFILE
vips_icc_load_profile_blob( VipsBlob *blob,
VipsImage *image, VipsIntent intent, int direction )
{
const void *data;
size_t size;
cmsHPROFILE profile;
VipsIccInfo *info;
#ifdef DEBUG
printf( "loading %s profile, intent %s, from blob %p\n",
direction == LCMS_USED_AS_INPUT ? _( "input" ) : _( "output" ),
vips_enum_nick( VIPS_TYPE_INTENT, intent ),
blob );
#endif /*DEBUG*/
data = vips_blob_get( blob, &size );
if( !(profile = cmsOpenProfileFromMem( data, size )) ) {
g_warning( "%s", _( "corrupt profile" ) );
return( NULL );
}
#ifdef DEBUG
vips_icc_print_profile( "loaded from blob to make", profile );
#endif /*DEBUG*/
if( !(info = vips_icc_info( cmsGetColorSpace( profile ) )) ) {
VIPS_FREEF( cmsCloseProfile, profile );
g_warning( "%s", _( "unsupported profile" ) );
return( NULL );
}
if( image &&
image->Bands < info->bands ) {
VIPS_FREEF( cmsCloseProfile, profile );
g_warning( "%s", _( "profile incompatible with image" ) );
return( NULL );
}
/* Spot the common error of an RGB profile and a CMYK image.
*/
if( image &&
image->Type == VIPS_INTERPRETATION_CMYK &&
info->bands == 3 ) {
VIPS_FREEF( cmsCloseProfile, profile );
g_warning( "%s", _( "profile incompatible with image" ) );
return( NULL );
}
if( !cmsIsIntentSupported( profile, intent, direction ) ) {
VIPS_FREEF( cmsCloseProfile, profile );
g_warning( _( "profile does not support %s %s intent" ),
vips_enum_nick( VIPS_TYPE_INTENT, intent ),
direction == LCMS_USED_AS_INPUT ?
_( "input" ) : _( "output" ) );
return( NULL );
}
return( profile );
}
/* Verify that a blob is not corrupt and is compatible with this image.
*
* unref the blob if it's useless.
*/
static cmsHPROFILE
vips_icc_verify_blob( VipsBlob **blob,
VipsImage *image, VipsIntent intent, int direction )
{
if( *blob ) {
cmsHPROFILE profile;
if( !(profile = vips_icc_load_profile_blob( *blob,
image, intent, direction )) ) {
vips_area_unref( (VipsArea *) *blob );
*blob = NULL;
}
return( profile );
}
return( NULL );
}
/* Try to set the inport profile. We read the input profile like this:
*
* embedded filename action
* 0 0 image
* 1 0 image
* 0 1 file
* 1 1 image, then fall back to file
*
* If averything fails, we fall back to our built-in profiles, either
* srgb or cmyk, depending on the input image.
*
* We set attach_input_profile if we used a non-emdedded profile. The profile
* in in_blob will need to be attached to the output image in some way.
*/
static int
vips_icc_set_import( VipsIcc *icc,
gboolean embedded, const char *input_profile_filename )
{
VipsObjectClass *class = VIPS_OBJECT_GET_CLASS( icc );
VipsColourCode *code = (VipsColourCode *) icc;
icc->non_standard_input_profile = FALSE;
/* Try embedded profile.
*/
if( code->in &&
(embedded || !input_profile_filename) ) {
icc->in_blob = vips_icc_get_profile_image( code->in );
icc->in_profile = vips_icc_verify_blob( &icc->in_blob,
code->in, icc->intent, LCMS_USED_AS_INPUT );
}
/* Try profile from filename.
*/
if( code->in &&
!icc->in_blob &&
input_profile_filename ) {
if(
!vips_profile_load( input_profile_filename,
&icc->in_blob, NULL ) &&
(icc->in_profile = vips_icc_verify_blob( &icc->in_blob,
code->in, icc->intent, LCMS_USED_AS_INPUT )) )
icc->non_standard_input_profile = TRUE;
}
/* Try built-in profile.
*/
if( code->in &&
!icc->in_profile ) {
const char *name = code->in->Type == VIPS_INTERPRETATION_CMYK ?
"cmyk" : "srgb";
if(
!vips_profile_load( name, &icc->in_blob, NULL ) &&
(icc->in_profile = vips_icc_verify_blob( &icc->in_blob,
code->in, icc->intent, LCMS_USED_AS_INPUT )) )
icc->non_standard_input_profile = TRUE;
}
if( !icc->in_profile ) {
vips_error( class->nickname, "%s", _( "unable to load or "
"find any compatible input profile" ) );
return( -1 );
}
return( 0 );
}
static void
vips_icc_class_init( VipsIccClass *class )
{
GObjectClass *gobject_class = G_OBJECT_CLASS( class );
VipsObjectClass *object_class = (VipsObjectClass *) class;
gobject_class->dispose = vips_icc_dispose;
gobject_class->set_property = vips_object_set_property;
gobject_class->get_property = vips_object_get_property;
object_class->nickname = "icc";
object_class->description = _( "transform using ICC profiles" );
object_class->build = vips_icc_build;
VIPS_ARG_ENUM( class, "intent", 6,
_( "Intent" ),
_( "Rendering intent" ),
VIPS_ARGUMENT_OPTIONAL_INPUT,
G_STRUCT_OFFSET( VipsIcc, intent ),
VIPS_TYPE_INTENT, VIPS_INTENT_RELATIVE );
VIPS_ARG_ENUM( class, "pcs", 6,
_( "PCS" ),
_( "Set Profile Connection Space" ),
VIPS_ARGUMENT_OPTIONAL_INPUT,
G_STRUCT_OFFSET( VipsIcc, pcs ),
VIPS_TYPE_PCS, VIPS_PCS_LAB );
VIPS_ARG_BOOL( class, "black_point_compensation", 7,
_( "Black point compensation" ),
_( "Enable black point compensation" ),
VIPS_ARGUMENT_OPTIONAL_INPUT,
G_STRUCT_OFFSET( VipsIcc, black_point_compensation ),
FALSE );
cmsSetLogErrorHandler( icc_error );
}
static void
vips_icc_init( VipsIcc *icc )
{
icc->intent = VIPS_INTENT_RELATIVE;
icc->pcs = VIPS_PCS_LAB;
icc->depth = 8;
}
typedef struct _VipsIccImport {
VipsIcc parent_instance;
gboolean embedded;
char *input_profile_filename;
} VipsIccImport;
typedef VipsIccClass VipsIccImportClass;
G_DEFINE_TYPE( VipsIccImport, vips_icc_import, VIPS_TYPE_ICC );
static int
vips_icc_import_build( VipsObject *object )
{
VipsColour *colour = (VipsColour *) object;
VipsIcc *icc = (VipsIcc *) object;
VipsIccImport *import = (VipsIccImport *) object;
if( vips_icc_set_import( icc,
import->embedded, import->input_profile_filename ) )
return( -1 );
if( icc->pcs == VIPS_PCS_LAB ) {
cmsCIExyY white;
cmsWhitePointFromTemp( &white, 6500 );
icc->out_profile = cmsCreateLab4Profile( &white );
}
else
icc->out_profile = cmsCreateXYZProfile();
if( VIPS_OBJECT_CLASS( vips_icc_import_parent_class )->build( object ) )
return( -1 );
/* If we used the fallback profile, we need to attach it to the PCS
* image, since the PCS image needs a route back to device space.
*
* In the same way, we don't remove the embedded input profile on
* import.
*/
if( icc->non_standard_input_profile &&
icc->in_blob ) {
const void *data;
size_t size;
data = vips_blob_get( icc->in_blob, &size );
vips_image_set_blob( colour->out, VIPS_META_ICC_NAME,
NULL, data, size );
}
return( 0 );
}
static void
decode_lab( guint16 *fixed, float *lab, int n )
{
int i;
for( i = 0; i < n; i++ ) {
lab[0] = (double) fixed[0] / 652.800;
lab[1] = ((double) fixed[1] / 256.0) - 128.0;
lab[2] = ((double) fixed[2] / 256.0) - 128.0;
lab += 3;
fixed += 3;
}
}
#define X_FAC (VIPS_D50_X0 * 32768 / (VIPS_D65_X0 * 100))
#define Y_FAC (VIPS_D50_Y0 * 32768 / (VIPS_D65_Y0 * 100))
#define Z_FAC (VIPS_D50_Z0 * 32768 / (VIPS_D65_Z0 * 100))
static void
decode_xyz( guint16 *fixed, float *xyz, int n )
{
int i;
for( i = 0; i < n; i++ ) {
xyz[0] = (double) fixed[0] / X_FAC;
xyz[1] = (double) fixed[1] / Y_FAC;
xyz[2] = (double) fixed[2] / Z_FAC;
xyz += 3;
fixed += 3;
}
}
/* Process a buffer of data.
*/
static void
vips_icc_import_line( VipsColour *colour,
VipsPel *out, VipsPel **in, int width )
{
VipsIcc *icc = (VipsIcc *) colour;
VipsPel *p;
float *q;
int i;
/* Buffer of encoded 16-bit pixels we transform.
*/
guint16 encoded[3 * PIXEL_BUFFER_SIZE];
p = (VipsPel *) in[0];
q = (float *) out;
for( i = 0; i < width; i += PIXEL_BUFFER_SIZE ) {
const int chunk = VIPS_MIN( width - i, PIXEL_BUFFER_SIZE );
cmsDoTransform( icc->trans, p, encoded, chunk );
if( icc->pcs == VIPS_PCS_LAB )
decode_lab( encoded, q, chunk );
else
decode_xyz( encoded, q, chunk );
p += PIXEL_BUFFER_SIZE * VIPS_IMAGE_SIZEOF_PEL( colour->in[0] );
q += PIXEL_BUFFER_SIZE * 3;
}
}
static void
vips_icc_import_class_init( VipsIccImportClass *class )
{
GObjectClass *gobject_class = G_OBJECT_CLASS( class );
VipsObjectClass *object_class = (VipsObjectClass *) class;
VipsColourClass *colour_class = VIPS_COLOUR_CLASS( class );
gobject_class->set_property = vips_object_set_property;
gobject_class->get_property = vips_object_get_property;
object_class->nickname = "icc_import";
object_class->description = _( "import from device with ICC profile" );
object_class->build = vips_icc_import_build;
colour_class->process_line = vips_icc_import_line;
VIPS_ARG_BOOL( class, "embedded", 110,
_( "Embedded" ),
_( "Use embedded input profile, if available" ),
VIPS_ARGUMENT_OPTIONAL_INPUT,
G_STRUCT_OFFSET( VipsIccImport, embedded ),
FALSE );
VIPS_ARG_STRING( class, "input_profile", 120,
_( "Input profile" ),
_( "Filename to load input profile from" ),
VIPS_ARGUMENT_OPTIONAL_INPUT,
G_STRUCT_OFFSET( VipsIccImport, input_profile_filename ),
NULL );
}
static void
vips_icc_import_init( VipsIccImport *import )
{
}
typedef struct _VipsIccExport {
VipsIcc parent_instance;
char *output_profile_filename;
} VipsIccExport;
typedef VipsIccClass VipsIccExportClass;
G_DEFINE_TYPE( VipsIccExport, vips_icc_export, VIPS_TYPE_ICC );
static int
vips_icc_export_build( VipsObject *object )
{
VipsObjectClass *class = VIPS_OBJECT_GET_CLASS( object );
VipsColour *colour = (VipsColour *) object;
VipsColourCode *code = (VipsColourCode *) object;
VipsIcc *icc = (VipsIcc *) object;
VipsIccExport *export = (VipsIccExport *) object;
/* If icc->pcs hasn't been set and this image is tagged as XYZ, swap
* to XYZ pcs. This will save a XYZ->LAB conversion when we chain up.
*/
if( !vips_object_argument_isset( object, "pcs" ) &&
code->in &&
code->in->Type == VIPS_INTERPRETATION_XYZ )
icc->pcs = VIPS_PCS_XYZ;
if( icc->pcs == VIPS_PCS_LAB ) {
cmsCIExyY white;
cmsWhitePointFromTemp( &white, 6500 );
icc->in_profile = cmsCreateLab4Profile( &white );
}
else
icc->in_profile = cmsCreateXYZProfile();
if( code->in &&
!export->output_profile_filename )
icc->out_blob = vips_icc_get_profile_image( code->in );
if( !icc->out_blob &&
export->output_profile_filename ) {
if( vips_profile_load( export->output_profile_filename,
&icc->out_blob, NULL ) )
return( -1 );
colour->profile_filename = export->output_profile_filename;
}
if( icc->out_blob &&
!(icc->out_profile = vips_icc_load_profile_blob( icc->out_blob,
NULL, icc->intent, LCMS_USED_AS_OUTPUT )) ) {
vips_error( class->nickname, "%s", _( "no output profile" ) );
return( -1 );
}
if( VIPS_OBJECT_CLASS( vips_icc_export_parent_class )->build( object ) )
return( -1 );
return( 0 );
}
/* Pack a buffer of floats into lcms's fixed-point formats. Cut from
* lcms-1.0.8.
*/
static void
encode_lab( float *lab, guint16 *fixed, int n )
{
int i;
for( i = 0; i < n; i++ ) {
float L = lab[0];
float a = lab[1];
float b = lab[2];
if( L < 0 )
L = 0;
if( L > 100. )
L = 100.;
if( a < -128. )
a = -128;
if( a > 127.9961 )
a = 127.9961;
if( b < -128. )
b = -128;
if( b > 127.9961 )
b = 127.9961;
fixed[0] = L * 652.800 + 0.5;
fixed[1] = (a + 128.0) * 256.0 + 0.5;
fixed[2] = (b + 128.0) * 256.0 + 0.5;
lab += 3;
fixed += 3;
}
}
#define MAX_ENCODEABLE_XYZ (100 * (1.0 + 32767.0 / 32768.0))
// 1.15 fixed point for XYZ
static void
encode_xyz( float *xyz, guint16 *fixed, int n )
{
int i;
for( i = 0; i < n; i++ ) {
float X = xyz[0];
float Y = xyz[1];
float Z = xyz[2];
if( X < 0 )
X = 0;
if( X > MAX_ENCODEABLE_XYZ )
X = MAX_ENCODEABLE_XYZ;
if( Y < 0 )
Y = 0;
if( Y > MAX_ENCODEABLE_XYZ )
Y = MAX_ENCODEABLE_XYZ;
if( Z < 0 )
Z = 0;
if( Z > MAX_ENCODEABLE_XYZ )
Z = MAX_ENCODEABLE_XYZ;
fixed[0] = X * X_FAC + 0.5;
fixed[1] = Y * Y_FAC + 0.5;
fixed[2] = Z * Z_FAC + 0.5;
xyz += 3;
fixed += 3;
}
}
/* Process a buffer of data.
*/
static void
vips_icc_export_line( VipsColour *colour,
VipsPel *out, VipsPel **in, int width )
{
VipsIcc *icc = (VipsIcc *) colour;
float *p;
VipsPel *q;
int x;
/* Buffer of encoded 16-bit pixels we transform.
*/
guint16 encoded[3 * PIXEL_BUFFER_SIZE];
p = (float *) in[0];
q = (VipsPel *) out;
for( x = 0; x < width; x += PIXEL_BUFFER_SIZE ) {
const int chunk = VIPS_MIN( width - x, PIXEL_BUFFER_SIZE );
if( icc->pcs == VIPS_PCS_LAB )
encode_lab( p, encoded, chunk );
else
encode_xyz( p, encoded, chunk );
cmsDoTransform( icc->trans, encoded, q, chunk );
p += PIXEL_BUFFER_SIZE * 3;
q += PIXEL_BUFFER_SIZE * VIPS_IMAGE_SIZEOF_PEL( colour->out );
}
}
static void
vips_icc_export_class_init( VipsIccExportClass *class )
{
GObjectClass *gobject_class = G_OBJECT_CLASS( class );
VipsObjectClass *object_class = (VipsObjectClass *) class;
VipsColourClass *colour_class = VIPS_COLOUR_CLASS( class );
gobject_class->set_property = vips_object_set_property;
gobject_class->get_property = vips_object_get_property;
object_class->nickname = "icc_export";
object_class->description = _( "output to device with ICC profile" );
object_class->build = vips_icc_export_build;
colour_class->process_line = vips_icc_export_line;
VIPS_ARG_STRING( class, "output_profile", 110,
_( "Output profile" ),
_( "Filename to load output profile from" ),
VIPS_ARGUMENT_OPTIONAL_INPUT,
G_STRUCT_OFFSET( VipsIccExport, output_profile_filename ),
NULL );
VIPS_ARG_INT( class, "depth", 130,
_( "Depth" ),
_( "Output device space depth in bits" ),
VIPS_ARGUMENT_OPTIONAL_INPUT,
G_STRUCT_OFFSET( VipsIcc, depth ),
8, 16, 8 );
}
static void
vips_icc_export_init( VipsIccExport *export )
{
}
typedef struct _VipsIccTransform {
VipsIcc parent_instance;
gboolean embedded;
char *input_profile_filename;
char *output_profile_filename;
} VipsIccTransform;
typedef VipsIccClass VipsIccTransformClass;
G_DEFINE_TYPE( VipsIccTransform, vips_icc_transform, VIPS_TYPE_ICC );
static int
vips_icc_transform_build( VipsObject *object )
{
VipsObjectClass *class = VIPS_OBJECT_GET_CLASS( object );
VipsColour *colour = (VipsColour *) object;
VipsIcc *icc = (VipsIcc *) object;
VipsIccTransform *transform = (VipsIccTransform *) object;
if( vips_icc_set_import( icc,
transform->embedded, transform->input_profile_filename ) )
return( -1 );
if( transform->output_profile_filename ) {
if( vips_profile_load( transform->output_profile_filename,
&icc->out_blob, NULL ) )
return( -1 );
colour->profile_filename = transform->output_profile_filename;
}
if( icc->out_blob )
icc->out_profile = vips_icc_load_profile_blob( icc->out_blob,
NULL, icc->intent, LCMS_USED_AS_OUTPUT );
if( !icc->out_profile ) {
vips_error( class->nickname, "%s", _( "no output profile" ) );
return( -1 );
}
if( VIPS_OBJECT_CLASS( vips_icc_transform_parent_class )->
build( object ) )
return( -1 );
return( 0 );
}
/* Process a buffer of data.
*/
static void
vips_icc_transform_line( VipsColour *colour,
VipsPel *out, VipsPel **in, int width )
{
VipsIcc *icc = (VipsIcc *) colour;
cmsDoTransform( icc->trans, in[0], out, width );
}
static void
vips_icc_transform_class_init( VipsIccImportClass *class )
{
GObjectClass *gobject_class = G_OBJECT_CLASS( class );
VipsObjectClass *object_class = (VipsObjectClass *) class;
VipsColourClass *colour_class = VIPS_COLOUR_CLASS( class );
gobject_class->set_property = vips_object_set_property;
gobject_class->get_property = vips_object_get_property;
object_class->nickname = "icc_transform";
object_class->description =
_( "transform between devices with ICC profiles" );
object_class->build = vips_icc_transform_build;
colour_class->process_line = vips_icc_transform_line;
VIPS_ARG_STRING( class, "output_profile", 110,
_( "Output profile" ),
_( "Filename to load output profile from" ),
VIPS_ARGUMENT_REQUIRED_INPUT,
G_STRUCT_OFFSET( VipsIccTransform, output_profile_filename ),
NULL );
VIPS_ARG_BOOL( class, "embedded", 120,
_( "Embedded" ),
_( "Use embedded input profile, if available" ),
VIPS_ARGUMENT_OPTIONAL_INPUT,
G_STRUCT_OFFSET( VipsIccTransform, embedded ),
FALSE );
VIPS_ARG_STRING( class, "input_profile", 130,
_( "Input profile" ),
_( "Filename to load input profile from" ),
VIPS_ARGUMENT_OPTIONAL_INPUT,
G_STRUCT_OFFSET( VipsIccTransform, input_profile_filename ),
NULL );
VIPS_ARG_INT( class, "depth", 140,
_( "Depth" ),
_( "Output device space depth in bits" ),
VIPS_ARGUMENT_OPTIONAL_INPUT,
G_STRUCT_OFFSET( VipsIcc, depth ),
8, 16, 8 );
}
static void
vips_icc_transform_init( VipsIccTransform *transform )
{
}
/**
* vips_icc_ac2rc: (method)
* @in: input image
* @out: (out): output image
* @profile_filename: use this profile
*
* Transform an image from absolute to relative colorimetry using the
* MediaWhitePoint stored in the ICC profile.
*
* See also: vips_icc_transform(), vips_icc_import().
*
* Returns: 0 on success, -1 on error.
*/
int
vips_icc_ac2rc( VipsImage *in, VipsImage **out, const char *profile_filename )
{
VipsImage *t;
cmsHPROFILE profile;
cmsCIEXYZ *media;
double X, Y, Z;
double *add;
double *mul;
int i;
#ifdef G_OS_WIN32
if( !(profile = cmsOpenProfileFromFile( profile_filename, "r" )) )
#else /*!G_OS_WIN32*/
if( !(profile = cmsOpenProfileFromFile( profile_filename, "re" )) )
#endif /*G_OS_WIN32*/
return( -1 );
#ifdef DEBUG
vips_icc_print_profile( profile_filename, profile );
#endif /*DEBUG*/
if( !(media = cmsReadTag( profile, cmsSigMediaWhitePointTag )) ) {
vips_error( "vips_icc_ac2rc",
"%s", _( "unable to get media white point" ) );
return( -1 );
}
X = media->X;
Y = media->Y;
Z = media->Z;
cmsCloseProfile( profile );
/* We need XYZ so we can adjust the white balance.
*/
if( vips_colourspace( in, &t, VIPS_INTERPRETATION_XYZ, NULL ) )
return( -1 );
in = t;
if( !(add = VIPS_ARRAY( in, in->Bands, double )) ||
!(mul = VIPS_ARRAY( in, in->Bands, double )) )
return( -1 );
/* There might be extra bands off to the right somewhere.
*/
for( i = 0; i < in->Bands; i++ )
add[i] = 0.0;
mul[0] = VIPS_D50_X0 / (X * 100.0);
mul[1] = VIPS_D50_Y0 / (Y * 100.0);
mul[2] = VIPS_D50_Z0 / (Z * 100.0);
for( i = 3; i < in->Bands; i++ )
mul[i] = 1.0;
if( vips_linear( in, &t, add, mul, in->Bands, NULL ) ) {
g_object_unref( in );
return( -1 );
}
g_object_unref( in );
in = t;
*out = in;
return( 0 );
}
/* TRUE if a profile is sane and is compatible with an image.
*/
gboolean
vips_icc_is_compatible_profile( VipsImage *image,
const void *data, size_t data_length )
{
cmsHPROFILE profile;
VipsIccInfo *info;
if( !(profile = cmsOpenProfileFromMem( data, data_length )) )
/* Corrupt profile.
*/
return( FALSE );
#ifdef DEBUG
vips_icc_print_profile( "from memory", profile );
#endif /*DEBUG*/
if( !(info = vips_icc_info( cmsGetColorSpace( profile ) )) ) {
/* Unsupported profile.
*/
VIPS_FREEF( cmsCloseProfile, profile );
return( FALSE );
}
if( image->Bands < info->bands ) {
/* Too few bands,
*/
VIPS_FREEF( cmsCloseProfile, profile );
return( FALSE );
}
VIPS_FREEF( cmsCloseProfile, profile );
return( TRUE );
}
#else /*!HAVE_LCMS2*/
#include <vips/vips.h>
int
vips_icc_present( void )
{
return( 0 );
}
int
vips_icc_ac2rc( VipsImage *in, VipsImage **out, const char *profile_filename )
{
vips_error( "VipsIcc", "%s",
_( "libvips configured without lcms support" ) );
return( -1 );
}
gboolean
vips_icc_is_compatible_profile( VipsImage *image,
const void *data, size_t data_length )
{
return( TRUE );
}
#endif /*HAVE_LCMS2*/
/**
* vips_icc_import: (method)
* @in: input image
* @out: (out): output image
* @...: %NULL-terminated list of optional named arguments
*
* Optional arguments:
*
* * @pcs: #VipsPCS, use XYZ or LAB PCS
* * @intent: #VipsIntent, transform with this intent
* * @black_point_compensation: %gboolean, enable black point compensation
* * @embedded: %gboolean, use profile embedded in input image
* * @input_profile: %gchararray, get the input profile from here
*
* Import an image from device space to D65 LAB with an ICC profile. If @pcs is
* set to #VIPS_PCS_XYZ, use CIE XYZ PCS instead.
*
* If @embedded is set, the input profile is taken from the input image
* metadata. If there is no embedded profile,
* @input_profile_filename is used as a fall-back.
* You can test for the
* presence of an embedded profile with
* vips_image_get_typeof() with #VIPS_META_ICC_NAME as an argument. This will
* return %GType 0 if there is no profile.
*
* If @embedded is not set, the input profile is taken from
* @input_profile. If @input_profile is not supplied, the
* metadata profile, if any, is used as a fall-back.
*
* If @black_point_compensation is set, LCMS black point compensation is
* enabled.
*
* Returns: 0 on success, -1 on error.
*/
int
vips_icc_import( VipsImage *in, VipsImage **out, ... )
{
va_list ap;
int result;
va_start( ap, out );
result = vips_call_split( "icc_import", ap, in, out );
va_end( ap );
return( result );
}
/**
* vips_icc_export: (method)
* @in: input image
* @out: (out): output image
* @...: %NULL-terminated list of optional named arguments
*
* Optional arguments:
*
* * @pcs: #VipsPCS, use XYZ or LAB PCS
* * @intent: #VipsIntent, transform with this intent
* * @black_point_compensation: %gboolean, enable black point compensation
* * @output_profile: %gchararray, get the output profile from here
* * @depth: %gint, depth of output image in bits
*
* Export an image from D65 LAB to device space with an ICC profile.
* If @pcs is
* set to #VIPS_PCS_XYZ, use CIE XYZ PCS instead.
* If @output_profile is not set, use the embedded profile, if any.
* If @output_profile is set, export with that and attach it to the output
* image.
*
* If @black_point_compensation is set, LCMS black point compensation is
* enabled.
*
* Returns: 0 on success, -1 on error.
*/
int
vips_icc_export( VipsImage *in, VipsImage **out, ... )
{
va_list ap;
int result;
va_start( ap, out );
result = vips_call_split( "icc_export", ap, in, out );
va_end( ap );
return( result );
}
/**
* vips_icc_transform: (method)
* @in: input image
* @out: (out): output image
* @output_profile: get the output profile from here
* @...: %NULL-terminated list of optional named arguments
*
* Optional arguments:
*
* * @pcs: #VipsPCS, use XYZ or LAB PCS
* * @intent: #VipsIntent, transform with this intent
* * @black_point_compensation: %gboolean, enable black point compensation
* * @embedded: %gboolean, use profile embedded in input image
* * @input_profile: %gchararray, get the input profile from here
* * @depth: %gint, depth of output image in bits
*
* Transform an image with a pair of ICC profiles. The input image is moved to
* profile-connection space with the input profile and then to the output
* space with the output profile.
*
* If @embedded is set, the input profile is taken from the input image
* metadata, if present. If there is no embedded profile,
* @input_profile is used as a fall-back.
* You can test for the
* presence of an embedded profile with
* vips_image_get_typeof() with #VIPS_META_ICC_NAME as an argument. This will
* return %GType 0 if there is no profile.
*
* If @embedded is not set, the input profile is taken from
* @input_profile. If @input_profile is not supplied, the
* metadata profile, if any, is used as a fall-back.
*
* If @black_point_compensation is set, LCMS black point compensation is
* enabled.
*
* The output image has the output profile attached to the #VIPS_META_ICC_NAME
* field.
*
* Use vips_icc_import() and vips_icc_export() to do either the first or
* second half of this operation in isolation.
*
* Returns: 0 on success, -1 on error.
*/
int
vips_icc_transform( VipsImage *in, VipsImage **out,
const char *output_profile, ... )
{
va_list ap;
int result;
va_start( ap, output_profile );
result = vips_call_split( "icc_transform", ap,
in, out, output_profile );
va_end( ap );
return( result );
}
Reference in New Issue View Git Blame Copy Permalink