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
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
|
#include <stdio.h>
#include <stdlib.h>
#include "minc2.h"
/* This test will attempt to create a few different test images
* which can be tested with minc2.0
*/
#define TESTRPT(msg, val) (error_cnt++, fprintf(stderr, \
"Error reported on line #%d, %s: %d\n", \
__LINE__, msg, val))
#define CZ 142
#define CY 245
#define CX 120
#define NDIMS 3
int create_real_as_int_image(void)
{
int r;
double start_values[NDIMS]={-6.96, -12.453, -9.48};
double separations[NDIMS]={0.09,0.09,0.09};
midimhandle_t hdim[NDIMS];
mihandle_t hvol;
int *buf = ( int *) malloc(CX * CY * CZ * sizeof(int));
int i;
long count[NDIMS];
long start[NDIMS];
miboolean_t flag=1;
double min = -1.0;
double max = 1.0;
r = micreate_dimension("yspace", MI_DIMCLASS_SPATIAL,
MI_DIMATTR_REGULARLY_SAMPLED, CY, &hdim[0]);
r = micreate_dimension("xspace", MI_DIMCLASS_SPATIAL,
MI_DIMATTR_REGULARLY_SAMPLED, CX, &hdim[1]);
r = micreate_dimension("zspace", MI_DIMCLASS_SPATIAL,
MI_DIMATTR_REGULARLY_SAMPLED, CZ, &hdim[2]);
r = miset_dimension_starts(hdim, NDIMS, start_values);
r = miset_dimension_separations(hdim, NDIMS, separations);
r = micreate_volume("real_as_int_image.mnc", NDIMS, hdim, MI_TYPE_INT,
MI_CLASS_REAL, NULL, &hvol);
/* set slice scaling flag to true */
r = miset_slice_scaling_flag(hvol, flag);
r = micreate_volume_image(hvol);
for (i = 0; i < CY*CX*CZ; i++) {
buf[i] = (int) i * 0.001;
}
start[0] = start[1] = start[2] = 0;
count[0] = CY; count[1] = CX; count[2] = CZ;
r = miset_voxel_value_hyperslab(hvol, MI_TYPE_INT, start, count, buf);
/* Set random values to slice min and max for slice scaling*/
start[0] =start[1]=start[2]=0;
for (i=0; i < CY; i++) {
start[0] = i;
min += -0.1;
max += 0.1;
r = miset_slice_range(hvol,start,NDIMS , max, min);
}
r = miclose_volume(hvol);
}
int create_real_as_float_image(void)
{
int r;
double start_values[NDIMS]={-6.96, -12.453, -9.48};
double separations[NDIMS]={0.09,0.09,0.09};
midimhandle_t hdim[NDIMS];
mihandle_t hvol;
float *buf = (float *) malloc(CX * CY * CZ * sizeof(float));
int i;
long count[NDIMS];
long start[NDIMS];
miboolean_t flag=1;
double min = -1.0;
double max = 1.0;
r = micreate_dimension("yspace", MI_DIMCLASS_SPATIAL,
MI_DIMATTR_REGULARLY_SAMPLED, CY, &hdim[0]);
r = micreate_dimension("xspace", MI_DIMCLASS_SPATIAL,
MI_DIMATTR_REGULARLY_SAMPLED, CX, &hdim[1]);
r = micreate_dimension("zspace", MI_DIMCLASS_SPATIAL,
MI_DIMATTR_REGULARLY_SAMPLED, CZ, &hdim[2]);
r = miset_dimension_starts(hdim, NDIMS, start_values);
r = miset_dimension_separations(hdim, NDIMS, separations);
r = micreate_volume("real_as_float_image.mnc", NDIMS, hdim, MI_TYPE_FLOAT,
MI_CLASS_REAL, NULL, &hvol);
/* set slice scaling flag to true */
r = miset_slice_scaling_flag(hvol, flag);
r = micreate_volume_image(hvol);
for (i = 0; i < CY*CX*CZ; i++) {
buf[i] = i * 0.001;
}
start[0] = start[1] = start[2] = 0;
count[0] = CY; count[1] = CX; count[2] = CZ;
r = miset_voxel_value_hyperslab(hvol, MI_TYPE_FLOAT, start, count, buf);
/* Set random values to slice min and max for slice scaling*/
start[0] =start[1]=start[2]=0;
for (i=0; i < CY; i++) {
start[0] = i;
min += -0.1;
max += 0.1;
r = miset_slice_range(hvol,start,NDIMS , max, min);
}
r = miclose_volume(hvol);
}
int main(int argc, char **argv)
{
int r = 0;
printf("Creating 3D image REAL stored as INT w/ slice scaling!! (real_as_int_image.mnc)\n");
r +=create_real_as_int_image();
printf("Creating 3D image REAL stored as FLOAT w/ slice scaling!! (real_as_float_image.mnc)\n"); r +=create_real_as_float_image();
if (r != 0) {
fprintf(stderr, "%d error%s reported\n",
r, (r == 1) ? "" : "s");
}
else {
fprintf(stderr, "\n No errors\n");
}
return (r);
}
|
