File: minc_rw_test2.cpp

package info (click to toggle)
libminc 2.4.07-2
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 4,288 kB
  • sloc: ansic: 57,268; cpp: 3,654; sh: 100; makefile: 23; ruby: 18
file content (42 lines) | stat: -rw-r--r-- 1,058 bytes parent folder | download | duplicates (6)
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
#include "minc_1_rw.h"
#include <iostream>
#include "minc_1_simple.h"
#include "minc_io_simple_volume.h"

using namespace minc;

int main(int argc,char **argv)
{
  try
  {
    if(argc<3) {
      std::cerr<<"Usage: "<<argv[0]<<" <input.mnc> <output.mnc>"<<std::endl;
      return 1;
    }
    minc_1_reader rdr;
    rdr.open(argv[1],true);
    rdr.setup_read_float();
    simple_volume<float> vol;
    vol.resize(rdr.ndim(1),rdr.ndim(2),rdr.ndim(3));
    load_non_standard_volume<float>(rdr,vol.c_buf());
    //rdr.close();
    for(size_t z=0;z<vol.dim(2);z++)
      for(size_t y=0;y<vol.dim(1);y++)
        for(size_t x=0;x<vol.dim(0);x++)
        {
          vol(x,y,z)=x+y+z;
        }
    minc_1_writer wrt;
    wrt.open(argv[2],rdr.info(),3,NC_FLOAT,false);
    wrt.setup_write_float();
    save_non_standard_volume<float>(wrt,vol.c_buf());
    
  } catch (const minc::generic_error & err) {
    std::cerr << "Got an error at:" << err.file () << ":" << err.line () << std::endl;
    std::cerr << err.msg()<<std::endl;
    return 1;
  }
  
  return 0;
}