c_nifti 0.0.1

data/ext/nifticlib/nifti_image_converters.c

@@ -0,0 +1,65 @@
+#include <ruby.h> // Ruby C extensions
+#include <nifti1_io.h> // standard NIfTI C library
+#include <stdbool.h> // Boolean type handling
+#include <nifti_image_converters.h>
+
+nifti_image * to_nifti_image(VALUE rb_nifti_image){
+  nifti_image *img = NULL;
+
+  Data_Get_Struct(rb_nifti_image, nifti_image, img);
+
+  return img;
+}
+
+VALUE voxel_data_to_value(nifti_image *img, int index){
+  // Types nifti1.h
+  switch(img->datatype){
+    case DT_BINARY: // DT_UINT8
+      if (((bool *) img->data)[index]){
+        return Qtrue;
+      }else{
+        return Qfalse;
+      }
+    case DT_UNSIGNED_CHAR:
+      return rb_sprintf("%u", ((unsigned char *) img->data)[index]);
+    case DT_SIGNED_SHORT: // DT_INT16
+      return INT2NUM((int) ((short *) img->data)[index]);
+    case DT_SIGNED_INT: // DT_INT32
+      return INT2NUM(((int *) img->data)[index]);
+    case DT_FLOAT: // DT_FLOAT32
+      return rb_float_new((double) ((float *) img->data)[index]);
+    case DT_DOUBLE: // DT_FLOAT64
+      return rb_float_new(((double *) img->data)[index]);
+    case DT_INT8:
+      return rb_sprintf("%c", ((char *) img->data)[index]);
+    case DT_UINT16:
+      return INT2NUM((int) ((unsigned short *) img->data)[index]);
+    case DT_UINT32:
+      return INT2NUM((int) ((unsigned int *) img->data)[index]);
+    case DT_INT64:
+      return INT2NUM((int) ((long long *) img->data)[index]);
+    case DT_UINT64:
+      return INT2NUM((int) ((unsigned long long *) img->data)[index]);
+    case DT_FLOAT128:
+      return rb_float_new(((double) ((long double *) img->data)[index]));
+    default: // Unsupported types: DT_COMPLEX, DT_RGB, DT_ALL, DT_COMPLEX128, DT_COMPLEX256, DT_RGBA32
+      return Qnil;
+  }
+}
+
+VALUE mat44_to_rb_array(mat44 matrix){
+  VALUE rb_data = rb_ary_new2(4);
+  int x, y;
+
+  for(x = 0; x < 4; x++){
+    VALUE line = rb_ary_new2(4);
+
+    for(y = 0; y < 4; y++){
+      rb_ary_store(line, y, rb_float_new((double) matrix.m[x][y]));
+    }
+
+    rb_ary_store(rb_data, x, line);
+  }
+
+  return rb_data;
+}

data/ext/nifticlib/nifti_image_dimensions.c

@@ -0,0 +1,113 @@
+#include <ruby.h> // Ruby C extensions
+#include <nifti1_io.h> // standard NIfTI C library
+
+#include <nifti_image_converters.h>
+#include <nifti_image_dimensions.h>
+
+VALUE nifti_image_ndim(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->ndim);
+}
+
+VALUE nifti_image_nvox(VALUE self){
+  return INT2NUM(to_nifti_image(self)->nvox);
+}
+
+VALUE nifti_image_nx(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->nx);
+}
+
+VALUE nifti_image_ny(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->ny);
+}
+
+VALUE nifti_image_nz(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->nz);
+}
+
+VALUE nifti_image_nt(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->nt);
+}
+
+VALUE nifti_image_nu(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->nu);
+}
+
+VALUE nifti_image_nv(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->nv);
+}
+
+VALUE nifti_image_nw(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->nw);
+}
+
+VALUE nifti_image_dim(VALUE self){
+  int i = 0;
+  nifti_image *img = to_nifti_image(self);
+
+  VALUE rb_data = rb_ary_new2(img->nvox);
+
+  for(i = 0; i < 8; i++){
+    rb_ary_store(rb_data, i, INT2NUM(img->dim[i]));
+  }
+
+  return rb_data;
+}
+
+VALUE nifti_image_nbyper(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->nbyper);
+}
+
+VALUE nifti_image_datatype(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->datatype);
+}
+
+VALUE nifti_image_xyz_units(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->xyz_units);
+}
+
+VALUE nifti_image_time_units(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->time_units);
+}
+
+VALUE nifti_image_dimensions_init(VALUE klass){
+  rb_define_method(klass, "ndim", nifti_image_ndim, 0);
+  rb_define_method(klass, "nx", nifti_image_nx, 0);
+  rb_define_method(klass, "ny", nifti_image_ny, 0);
+  rb_define_method(klass, "nz", nifti_image_nz, 0);
+  rb_define_method(klass, "nt", nifti_image_nt, 0);
+  rb_define_method(klass, "nu", nifti_image_nu, 0);
+  rb_define_method(klass, "nv", nifti_image_nv, 0);
+  rb_define_method(klass, "nw", nifti_image_nw, 0);
+  rb_define_method(klass, "dim", nifti_image_dim, 0);
+  rb_define_method(klass, "nvox", nifti_image_nvox, 0);
+  rb_define_method(klass, "nbyper", nifti_image_nbyper, 0);
+  rb_define_method(klass, "datatype", nifti_image_datatype, 0);
+  rb_define_method(klass, "xyz_units", nifti_image_xyz_units, 0);
+  rb_define_method(klass, "time_units", nifti_image_time_units, 0);
+
+  return klass;
+}

data/ext/nifticlib/nifti_image_intents.c

@@ -0,0 +1,45 @@
+#include <ruby.h> // Ruby C extensions
+#include <nifti1_io.h> // standard NIfTI C library
+
+#include <nifti_image_converters.h>
+#include <nifti_image_intents.h>
+
+VALUE nifti_image_intent_code(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->intent_code);
+}
+
+VALUE nifti_image_intent_p1(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_float_new((double) img->intent_p1);
+}
+
+VALUE nifti_image_intent_p2(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_float_new((double) img->intent_p2);
+}
+
+VALUE nifti_image_intent_p3(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_float_new((double) img->intent_p3);
+}
+
+VALUE nifti_image_intent_name(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_str_new2(img->intent_name);
+}
+
+VALUE nifti_image_intents_init(VALUE klass){
+  rb_define_method(klass, "intent_code", nifti_image_intent_code, 0);
+  rb_define_method(klass, "intent_p1", nifti_image_intent_p1, 0);
+  rb_define_method(klass, "intent_p2", nifti_image_intent_p2, 0);
+  rb_define_method(klass, "intent_p3", nifti_image_intent_p3, 0);
+  rb_define_method(klass, "intent_name", nifti_image_intent_name, 0);
+
+  return klass;
+}

data/ext/nifticlib/nifti_image_metadata.c

@@ -0,0 +1,133 @@
+#include <ruby.h> // Ruby C extensions
+#include <nifti1_io.h> // standard NIfTI C library
+
+#include <nifti_image_converters.h>
+#include <nifti_image_intents.h>
+
+#define MAX(a,b) (((a)>(b))?(a):(b))
+
+/***********/
+/* Getters */
+/***********/
+
+VALUE nifti_image_descrip(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_str_new2(img->descrip);
+}
+
+VALUE nifti_image_aux_file(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_str_new2(img->aux_file);
+}
+
+VALUE nifti_image_fname(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  if(img->fname){
+    return rb_str_new2(img->fname);
+  }else{
+    return Qnil;
+  }
+}
+
+VALUE nifti_image_iname(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_str_new2(img->iname);
+}
+
+VALUE nifti_image_iname_offset(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->iname_offset);
+}
+
+VALUE nifti_image_swapsize(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->swapsize);
+}
+
+VALUE nifti_image_byteorder(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->byteorder);
+}
+
+VALUE nifti_image_analyze75_orient(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM((int) img->analyze75_orient);
+}
+
+/***********/
+/* Setters */
+/***********/
+
+VALUE nifti_image_set_descrip(VALUE self, VALUE r_value){
+  nifti_image *img = to_nifti_image(self);
+  char *value;
+  int i = 0;
+
+  StringValue(r_value);
+  value = StringValuePtr(r_value);
+
+  for(i = 0; i < MAX(80, RSTRING_LEN(r_value)); i++){
+    img->descrip[i] = value[i];
+  }
+
+  return nifti_image_descrip(self);
+}
+
+VALUE nifti_image_set_aux_file(VALUE self, VALUE r_value){
+  nifti_image *img = to_nifti_image(self);
+  char *value;
+  int i = 0;
+
+  StringValue(r_value);
+  value = StringValuePtr(r_value);
+
+  for(i = 0; i < MAX(24, RSTRING_LEN(r_value)); i++){
+    img->aux_file[i] = value[i];
+  }
+
+  return nifti_image_aux_file(self);
+}
+
+VALUE nifti_image_set_filenames(VALUE self, VALUE r_value){
+  nifti_image *img = to_nifti_image(self);
+  char *value;
+
+  StringValue(r_value);
+  value = StringValuePtr(r_value);
+
+  nifti_set_filenames(img, value, 1, 1);
+
+  return r_value;
+}
+
+/********/
+/* Init */
+/********/
+
+VALUE nifti_image_metadata_init(VALUE klass){
+  // Getters
+  rb_define_method(klass, "descrip", nifti_image_descrip, 0);
+  rb_define_method(klass, "aux_file", nifti_image_aux_file, 0);
+  rb_define_method(klass, "fname", nifti_image_fname, 0);
+  rb_define_method(klass, "iname", nifti_image_iname, 0);
+  rb_define_method(klass, "iname_offset", nifti_image_iname_offset, 0);
+  rb_define_method(klass, "swapsize", nifti_image_swapsize, 0);
+  rb_define_method(klass, "byteorder", nifti_image_byteorder, 0);
+  rb_define_method(klass, "analyze75_orient", nifti_image_analyze75_orient, 0);
+
+  // Setters
+  // byteorder and swapsize cannot be setted
+  rb_define_method(klass, "descrip=", nifti_image_set_descrip, 1);
+  rb_define_method(klass, "aux_file=", nifti_image_set_aux_file, 1);
+  rb_define_method(klass, "filenames=", nifti_image_set_filenames, 1);
+
+  return klass;
+}

data/ext/nifticlib/nifti_image_quaternions.c

@@ -0,0 +1,59 @@
+#include <ruby.h> // Ruby C extensions
+#include <nifti1_io.h> // standard NIfTI C library
+
+#include <nifti_image_converters.h>
+#include <nifti_image_quaternions.h>
+
+VALUE nifti_image_quatern_b(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_float_new((double) img->quatern_b);
+}
+
+VALUE nifti_image_quatern_c(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_float_new((double) img->quatern_c);
+}
+
+VALUE nifti_image_quatern_d(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_float_new((double) img->quatern_d);
+}
+
+VALUE nifti_image_qoffset_x(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_float_new((double) img->qoffset_x);
+}
+
+VALUE nifti_image_qoffset_y(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_float_new((double) img->qoffset_y);
+}
+
+VALUE nifti_image_qoffset_z(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_float_new((double) img->qoffset_z);
+}
+
+VALUE nifti_image_qfac(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_float_new((double) img->qfac);
+}
+
+VALUE nifti_image_quaternions_init(VALUE klass){
+  rb_define_method(klass, "quatern_b", nifti_image_quatern_b, 0);
+  rb_define_method(klass, "quatern_c", nifti_image_quatern_c, 0);
+  rb_define_method(klass, "quatern_d", nifti_image_quatern_d, 0);
+  rb_define_method(klass, "qoffset_x", nifti_image_qoffset_x, 0);
+  rb_define_method(klass, "qoffset_y", nifti_image_qoffset_y, 0);
+  rb_define_method(klass, "qoffset_z", nifti_image_qoffset_z, 0);
+  rb_define_method(klass, "qfac", nifti_image_qfac, 0);
+
+  return klass;
+}

data/ext/nifticlib/nifti_image_spacings.c

@@ -0,0 +1,87 @@
+#include <ruby.h> // Ruby C extensions
+#include <nifti1_io.h> // standard NIfTI C library
+
+#include <nifti_image_converters.h>
+#include <nifti_image_spacings.h>
+
+VALUE nifti_image_dx(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_float_new((double) img->dx);
+}
+
+VALUE nifti_image_dy(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_float_new((double) img->dy);
+}
+
+VALUE nifti_image_dz(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_float_new((double) img->dz);
+}
+
+VALUE nifti_image_dt(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_float_new((double) img->dt);
+}
+
+VALUE nifti_image_du(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_float_new((double) img->du);
+}
+
+VALUE nifti_image_dv(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_float_new((double) img->dv);
+}
+
+VALUE nifti_image_dw(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return rb_float_new((double) img->dw);
+}
+
+VALUE nifti_image_pixdim(VALUE self){
+  int i = 0;
+  nifti_image *img = to_nifti_image(self);
+
+  VALUE rb_data = rb_ary_new2(img->nvox);
+
+  for(i = 0; i < 8; i++){
+    rb_ary_store(rb_data, i, rb_float_new((double) img->pixdim[i]));
+  }
+
+  return rb_data;
+}
+
+VALUE nifti_image_qform_code(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->qform_code);
+}
+
+VALUE nifti_image_sform_code(VALUE self){
+  nifti_image *img = to_nifti_image(self);
+
+  return INT2NUM(img->sform_code);
+}
+
+VALUE nifti_image_spacings_init(VALUE klass){
+  rb_define_method(klass, "dx", nifti_image_dx, 0);
+  rb_define_method(klass, "dy", nifti_image_dy, 0);
+  rb_define_method(klass, "dz", nifti_image_dz, 0);
+  rb_define_method(klass, "dt", nifti_image_dt, 0);
+  rb_define_method(klass, "du", nifti_image_du, 0);
+  rb_define_method(klass, "dv", nifti_image_dv, 0);
+  rb_define_method(klass, "dw", nifti_image_dw, 0);
+  rb_define_method(klass, "pixdim", nifti_image_pixdim, 0);
+  rb_define_method(klass, "qform_code", nifti_image_qform_code, 0);
+  rb_define_method(klass, "sform_code", nifti_image_sform_code, 0);
+
+  return klass;
+}