/**
* im_fits2vips:
* @filename: file to load
* @out: image to write to
*
* Read a FITS image file into a VIPS image.
*
* See also: im_vips2fits(), #VipsFormat.
*
* Returns: 0 on success, -1 on error.
*/
int
im_fits2vips( const char *filename, VipsImage *out )
{
VipsImage *t;
int n_bands;
VIPS_DEBUG_MSG( "im_fits2vips: reading \"%s\"\n", filename );
/* fits is naturally a band-separated format. For single-band images,
* we can just read out. For many bands, we read each band out
* separately then join them.
*/
if( !(t = vips_image_new()) ||
vips_object_local( out, t ) ||
fits2vips_header( filename, t ) )
return( -1 );
n_bands = t->Bands;
if( n_bands == 1 ) {
if( !(t = vips_image_new()) ||
vips_object_local( out, t ) ||
fits2vips( filename, t, 0 ) )
return( -1 );
}
else {
VipsImage *acc;
int i;
acc = NULL;
for( i = 0; i < n_bands; i++ ) {
if( !(t = vips_image_new()) ||
vips_object_local( out, t ) ||
fits2vips( filename, t, i ) )
return( -1 );
if( !acc )
acc = t;
else {
VipsImage *t2;
if( !(t2 = vips_image_new()) ||
vips_object_local( out, t2 ) ||
im_bandjoin( acc, t, t2 ) )
return( -1 );
acc = t2;
}
}
t = acc;
}
/* fits has inverted y.
*/
if( im_flipver( t, out ) )
return( -1 );
return( 0 );
}
/**
* im_tone_map:
* @in: input image
* @out: output image
* @lut: look-up table
*
* Map the first channel of @in through @lut. If @in is IM_CODING_LABQ, unpack
* to LABS, map L and then repack.
*
* @in should be a LABS or LABQ image for this to work
* sensibly.
*
* See also: im_maplut().
*
* Returns: 0 on success, -1 on error
*/
int
im_tone_map( IMAGE *in, IMAGE *out, IMAGE *lut )
{
IMAGE *t[8];
if( im_check_hist( "im_tone_map", lut ) ||
im_open_local_array( out, t, 8, "im_tone_map", "p" ) )
return( -1 );
/* If in is IM_CODING_LABQ, unpack.
*/
if( in->Coding == IM_CODING_LABQ ) {
if( im_LabQ2LabS( in, t[0] ) )
return( -1 );
}
else
t[0] = in;
/* Split into bands.
*/
if( im_extract_band( t[0], t[1], 0 ) )
return( -1 );
if( t[0]->Bands > 1 ) {
if( im_extract_bands( t[0], t[2], 1, t[0]->Bands - 1 ) )
return( -1 );
}
/* Map L.
*/
if( im_maplut( t[1], t[3], lut ) )
return( -1 );
/* Recombine bands.
*/
if( t[0]->Bands > 1 ) {
if( im_bandjoin( t[3], t[2], t[4] ) )
return( -1 );
}
else
t[4] = t[3];
/* If input was LabQ, repack.
*/
if( in->Coding == IM_CODING_LABQ ) {
if( im_LabS2LabQ( t[4], t[5] ) )
return( -1 );
}
else
t[5] = t[4];
return( im_copy( t[4], out ) );
}
/* Transform an n-band image with a 1-band processing function.
*/
int
im__fftproc( IMAGE *dummy, IMAGE *in, IMAGE *out, im__fftproc_fn fn )
{
if( im_pincheck( in ) || im_outcheck( out ) )
return( -1 );
if( in->Bands == 1 ) {
if( fn( dummy, in, out ) )
return( -1 );
}
else {
IMAGE *acc;
int b;
for( acc = NULL, b = 0; b < in->Bands; b++ ) {
IMAGE *t1 = im_open_local( dummy,
"fwfftn:1", "p" );
IMAGE *t2 = im_open_local( dummy,
"fwfftn:2", "p" );
if( !t1 || !t2 ||
im_extract_band( in, t1, b ) ||
fn( dummy, t1, t2 ) )
return( -1 );
if( !acc )
acc = t2;
else {
IMAGE *t3 = im_open_local( dummy,
"fwfftn:3", "p" );
if( !t3 || im_bandjoin( acc, t2, t3 ) )
return( -1 );
acc = t3;
}
}
if( im_copy( acc, out ) )
return( -1 );
}
return( 0 );
}
/* Call im_bandjoin via arg vector.
*/
static int
bandjoin_vec( im_object *argv )
{
return( im_bandjoin( argv[0], argv[1], argv[2] ) );
}
