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
|
#if HAVE_CONFIG_H
#include "config.h"
#endif
/* Create a grid transform from an existing MINC volume.
*/
#include <stdio.h>
#include <volume_io.h>
#include "time_stamp.h"
int main( int ac, char* av[] )
{
VIO_Volume v;
minc_input_options mio;
VIO_General_transform t;
VIO_Status st;
if ( ac != 3 ) {
fprintf( stderr, "usage: %s in_grid.mnc out.xfm\n", av[0] );
return 1;
}
/* The displacement volume must retain the vector dimension,
* so we turn off "vector -> scalar" conversion.
*/
set_default_minc_input_options( &mio );
set_minc_input_vector_to_scalar_flag( &mio, 0 );
st = input_volume( av[1],
0, NULL,
MI_ORIGINAL_TYPE, FALSE, 0.0, 0.0,
TRUE, &v,
&mio );
if ( st != VIO_OK ) {
fprintf( stderr, "failed to read grid volume \"%s\"\n", av[1] );
return 1;
}
create_grid_transform( &t, v, NULL );
st = output_transform_file( av[2], time_stamp(ac,av), &t );
if ( st != VIO_OK ) {
fprintf( stderr, "error writing to xfm file \"%s\"\n", av[2] );
return 1;
}
return 0;
}
|
