static void test_image_fetch(openslide_t *osr,
int64_t x, int64_t y,
int64_t w, int64_t h) {
uint32_t *buf = g_new(uint32_t, w * h);
for (int32_t level = 0; level < openslide_get_level_count(osr); level++) {
openslide_read_region(osr, buf, x, y, level, w, h);
}
g_free(buf);
const char *err = openslide_get_error(osr);
if (err) {
fail("Read failed: %"G_GINT64_FORMAT" %"G_GINT64_FORMAT
" %"G_GINT64_FORMAT" %"G_GINT64_FORMAT": %s", x, y, w, h, err);
}
}
static ReadSlide *
readslide_new( const char *filename, VipsImage *out,
int level, gboolean autocrop, const char *associated )
{
ReadSlide *rslide;
int64_t w, h;
const char *error;
const char *background;
const char * const *properties;
char *associated_names;
if( level &&
associated ) {
vips_error( "openslide2vips",
"%s", _( "specify only one of level or associated "
"image" ) );
return( NULL );
}
rslide = VIPS_NEW( out, ReadSlide );
memset( rslide, 0, sizeof( *rslide ) );
g_signal_connect( out, "close", G_CALLBACK( readslide_destroy_cb ),
rslide );
rslide->level = level;
rslide->autocrop = autocrop;
rslide->associated = g_strdup( associated );
/* Non-crazy defaults, override below if we can.
*/
rslide->tile_width = 256;
rslide->tile_height = 256;
rslide->osr = openslide_open( filename );
if( rslide->osr == NULL ) {
vips_error( "openslide2vips",
"%s", _( "unsupported slide format" ) );
return( NULL );
}
error = openslide_get_error( rslide->osr );
if( error ) {
vips_error( "openslide2vips",
_( "opening slide: %s" ), error );
return( NULL );
}
if( level < 0 ||
level >= openslide_get_level_count( rslide->osr ) ) {
vips_error( "openslide2vips",
"%s", _( "invalid slide level" ) );
return( NULL );
}
if( associated &&
check_associated_image( rslide->osr, associated ) )
return( NULL );
if( associated ) {
openslide_get_associated_image_dimensions( rslide->osr,
associated, &w, &h );
vips_image_set_string( out, "slide-associated-image",
associated );
vips_image_pipelinev( out, VIPS_DEMAND_STYLE_THINSTRIP, NULL );
}
else {
char buf[256];
const char *value;
openslide_get_level_dimensions( rslide->osr,
level, &w, &h );
rslide->downsample = openslide_get_level_downsample(
rslide->osr, level );
vips_image_set_int( out, "slide-level", level );
vips_image_pipelinev( out, VIPS_DEMAND_STYLE_SMALLTILE, NULL );
/* Try to get tile width/height. An undocumented, experimental
* feature.
*/
vips_snprintf( buf, 256,
"openslide.level[%d].tile-width", level );
if( (value = openslide_get_property_value( rslide->osr, buf )) )
rslide->tile_width = atoi( value );
vips_snprintf( buf, 256,
"openslide.level[%d].tile-height", level );
if( (value = openslide_get_property_value( rslide->osr, buf )) )
rslide->tile_height = atoi( value );
if( value )
VIPS_DEBUG_MSG( "readslide_new: found tile-size\n" );
/* Some images have a bounds in the header. Crop to
* that if autocrop is set.
*/
if( rslide->autocrop )
if( !get_bounds( rslide->osr, &rslide->bounds ) )
rslide->autocrop = FALSE;
if( rslide->autocrop ) {
VipsRect image;
rslide->bounds.left /= rslide->downsample;
rslide->bounds.top /= rslide->downsample;
rslide->bounds.width /= rslide->downsample;
rslide->bounds.height /= rslide->downsample;
/* Clip against image size.
*/
image.left = 0;
image.top = 0;
image.width = w;
image.height = h;
vips_rect_intersectrect( &rslide->bounds, &image,
&rslide->bounds );
/* If we've clipped to nothing, ignore bounds.
*/
if( vips_rect_isempty( &rslide->bounds ) )
rslide->autocrop = FALSE;
}
if( rslide->autocrop ) {
w = rslide->bounds.width;
h = rslide->bounds.height;
}
}
rslide->bg = 0xffffff;
if( (background = openslide_get_property_value( rslide->osr,
OPENSLIDE_PROPERTY_NAME_BACKGROUND_COLOR )) )
rslide->bg = strtoul( background, NULL, 16 );
if( w <= 0 ||
h <= 0 ||
rslide->downsample < 0 ) {
vips_error( "openslide2vips", _( "getting dimensions: %s" ),
openslide_get_error( rslide->osr ) );
return( NULL );
}
if( w > INT_MAX ||
h > INT_MAX ) {
vips_error( "openslide2vips",
"%s", _( "image dimensions overflow int" ) );
return( NULL );
}
if( !rslide->autocrop ) {
rslide->bounds.left = 0;
rslide->bounds.top = 0;
rslide->bounds.width = w;
rslide->bounds.height = h;
}
vips_image_init_fields( out, w, h, 4, VIPS_FORMAT_UCHAR,
VIPS_CODING_NONE, VIPS_INTERPRETATION_RGB, 1.0, 1.0 );
for( properties = openslide_get_property_names( rslide->osr );
*properties != NULL; properties++ )
vips_image_set_string( out, *properties,
openslide_get_property_value( rslide->osr,
*properties ) );
associated_names = g_strjoinv( ", ", (char **)
openslide_get_associated_image_names( rslide->osr ) );
vips_image_set_string( out,
"slide-associated-images", associated_names );
VIPS_FREE( associated_names );
return( rslide );
}
int main(int argc, char **argv) {
if (!g_thread_supported()) {
g_thread_init(NULL);
}
if (argc != 2) {
fail("No file specified");
}
const char *path = argv[1];
if (g_str_equal(path, "--leak-check--")) {
child_check_open_fds();
return 0;
}
openslide_get_version();
if (!openslide_detect_vendor(path)) {
fail("No vendor for %s", path);
}
openslide_t *osr = openslide_open(path);
if (!osr) {
fail("Couldn't open %s", path);
}
const char *err = openslide_get_error(osr);
if (err) {
fail("Open failed: %s", err);
}
openslide_close(osr);
osr = openslide_open(path);
if (!osr || openslide_get_error(osr)) {
fail("Reopen failed");
}
int64_t w, h;
openslide_get_level0_dimensions(osr, &w, &h);
int32_t levels = openslide_get_level_count(osr);
for (int32_t i = -1; i < levels + 1; i++) {
int64_t ww, hh;
openslide_get_level_dimensions(osr, i, &ww, &hh);
openslide_get_level_downsample(osr, i);
}
openslide_get_best_level_for_downsample(osr, 0.8);
openslide_get_best_level_for_downsample(osr, 1.0);
openslide_get_best_level_for_downsample(osr, 1.5);
openslide_get_best_level_for_downsample(osr, 2.0);
openslide_get_best_level_for_downsample(osr, 3.0);
openslide_get_best_level_for_downsample(osr, 3.1);
openslide_get_best_level_for_downsample(osr, 10);
openslide_get_best_level_for_downsample(osr, 20);
openslide_get_best_level_for_downsample(osr, 25);
openslide_get_best_level_for_downsample(osr, 100);
openslide_get_best_level_for_downsample(osr, 1000);
openslide_get_best_level_for_downsample(osr, 10000);
// NULL buffer
openslide_read_region(osr, NULL, 0, 0, 0, 1000, 1000);
// empty region
openslide_read_region(osr, NULL, 0, 0, 0, 0, 0);
// read properties
const char * const *property_names = openslide_get_property_names(osr);
while (*property_names) {
const char *name = *property_names;
openslide_get_property_value(osr, name);
property_names++;
}
// read associated images
const char * const *associated_image_names =
openslide_get_associated_image_names(osr);
while (*associated_image_names) {
int64_t w, h;
const char *name = *associated_image_names;
openslide_get_associated_image_dimensions(osr, name, &w, &h);
uint32_t *buf = g_new(uint32_t, w * h);
openslide_read_associated_image(osr, name, buf);
g_free(buf);
associated_image_names++;
}
test_image_fetch(osr, -10, -10, 200, 200);
test_image_fetch(osr, w/2, h/2, 500, 500);
test_image_fetch(osr, w - 200, h - 100, 500, 400);
test_image_fetch(osr, w*2, h*2, 400, 400);
test_image_fetch(osr, w - 20, 0, 40, 100);
test_image_fetch(osr, 0, h - 20, 100, 40);
// active region
const char *bounds_x = openslide_get_property_value(osr, OPENSLIDE_PROPERTY_NAME_BOUNDS_X);
const char *bounds_y = openslide_get_property_value(osr, OPENSLIDE_PROPERTY_NAME_BOUNDS_Y);
int64_t bounds_xx = 0;
int64_t bounds_yy = 0;
if (bounds_x && bounds_y) {
bounds_xx = g_ascii_strtoll(bounds_x, NULL, 10);
bounds_yy = g_ascii_strtoll(bounds_y, NULL, 10);
test_image_fetch(osr, bounds_xx, bounds_yy, 200, 200);
}
openslide_close(osr);
check_cloexec_leaks(path, argv[0], bounds_xx, bounds_yy);
return 0;
}