static int
vips_foreign_load_fits_header( VipsForeignLoad *load )
{
VipsForeignLoadFits *fits = (VipsForeignLoadFits *) load;
if( vips__fits_read_header( fits->filename, load->out ) )
return( -1 );
VIPS_SETSTR( load->out->filename, fits->filename );
return( 0 );
}
int
vips__fits_read( const char *filename, VipsImage *out )
{
VipsImage *t;
int n_bands;
VIPS_DEBUG_MSG( "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.
*/
t = vips_image_new();
if( vips__fits_read_header( filename, t ) ) {
g_object_unref( t );
return( -1 );
}
n_bands = t->Bands;
g_object_unref( t );
if( n_bands == 1 ) {
if( fits2vips( filename, out, 0 ) )
return( -1 );
}
else {
VipsImage **x;
int i;
t = vips_image_new();
x = (VipsImage **) vips_object_local_array( VIPS_OBJECT( t ),
n_bands + 1 );
for( i = 0; i < n_bands; i++ ) {
x[i] = vips_image_new();
if( fits2vips( filename, x[i], i ) ) {
g_object_unref( t );
return( -1 );
}
}
if( vips_bandjoin( x, &x[n_bands], n_bands, NULL ) ||
vips_image_write( x[n_bands], out ) ) {
g_object_unref( t );
return( -1 );
}
g_object_unref( t );
}
return( 0 );
}