1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
|
/*
* OpenSlide, a library for reading whole slide image files
*
* Copyright (c) 2007-2015 Carnegie Mellon University
* All rights reserved.
*
* OpenSlide 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, version 2.1.
*
* OpenSlide 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 OpenSlide. If not, see
* <http://www.gnu.org/licenses/>.
*
*/
#include <stdio.h>
#include <stdlib.h>
#include <stdarg.h>
#include <inttypes.h>
#include <glib.h>
#include <openslide.h>
#include "config.h"
#ifdef HAVE_VALGRIND
#include <callgrind.h>
#endif
#define BUFWIDTH 1000
#define BUFHEIGHT 1000
#define MAXWIDTH 10000
#define MAXHEIGHT 10000
static void fail(const char *str, ...) {
va_list ap;
va_start(ap, str);
vfprintf(stderr, str, ap);
fprintf(stderr, "\n");
va_end(ap);
exit(1);
}
int main(int argc, char **argv) {
if (argc != 3) {
fail("Usage: %s <slide> <level>", argv[0]);
}
const char *path = argv[1];
int level = atoi(argv[2]);
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);
}
if (level >= openslide_get_level_count(osr)) {
fail("No such level: %d", level);
}
// Get dimensions
int64_t x = 0;
int64_t y = 0;
int64_t w, h;
openslide_get_level_dimensions(osr, level, &w, &h);
// Read from active region, if denoted
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);
const char *bounds_w = openslide_get_property_value(osr, OPENSLIDE_PROPERTY_NAME_BOUNDS_WIDTH);
const char *bounds_h = openslide_get_property_value(osr, OPENSLIDE_PROPERTY_NAME_BOUNDS_HEIGHT);
if (bounds_x && bounds_y) {
x = g_ascii_strtoll(bounds_x, NULL, 10);
y = g_ascii_strtoll(bounds_y, NULL, 10);
}
if (bounds_w && bounds_h) {
double downsample = openslide_get_level_downsample(osr, level);
w = g_ascii_strtoll(bounds_w, NULL, 10) / downsample;
h = g_ascii_strtoll(bounds_h, NULL, 10) / downsample;
}
w = MIN(w, MAXWIDTH);
h = MIN(h, MAXHEIGHT);
uint32_t *buf = g_new(uint32_t, BUFWIDTH * BUFHEIGHT);
printf("Reading (%"PRId64", %"PRId64") in level %d for "
"%"PRId64" x %"PRId64"\n\n", x, y, level, w, h);
#ifdef HAVE_VALGRIND
CALLGRIND_START_INSTRUMENTATION;
#endif
for (int yy = 0; yy < h; yy += BUFHEIGHT) {
for (int xx = 0; xx < w; xx += BUFWIDTH) {
openslide_read_region(osr, buf, x + xx, y + yy, level,
MIN(BUFWIDTH, w - xx), MIN(BUFHEIGHT, h - yy));
}
}
#ifdef HAVE_VALGRIND
CALLGRIND_STOP_INSTRUMENTATION;
#endif
err = openslide_get_error(osr);
if (err) {
fail("Read failed: %s", err);
}
g_free(buf);
openslide_close(osr);
return 0;
}
|
