Doxygen Source Code Documentation
arentest.c File Reference
#include "aren.h"
#include <stdlib.h>
Go to the source code of this file.
Functions | |
int | main (int argc, char *argv[]) |
Function Documentation
|
convert DTIStudio fiber format data to SUMA segment data Definition at line 4 of file arentest.c. References AREN_be_verbose(), AREN_render(), AREN_set_graybrick(), AREN_set_opabrick(), AREN_set_rgbbricks(), AREN_set_size(), AREN_set_viewpoint(), argc, DSET_BRICK, DSET_delete, DSET_DX, DSET_DY, DSET_DZ, DSET_load, DSET_LOADED, DSET_NVOX, DSET_NX, DSET_NY, DSET_NZ, MRI_BYTE_PTR, mri_free(), mri_write_pnm(), new_AREN_renderer(), strtod(), and THD_open_dataset().
00005 { 00006 THD_3dim_dataset * dset , * oset ; 00007 MRI_IMAGE * grim , * opim , * renim ; 00008 float theta , phi , sx,sy,sz ; 00009 void * arh ; 00010 char fname[128] ; 00011 int nopt ; 00012 00013 dset = THD_open_dataset( argv[1] ) ; 00014 if( dset == NULL ) exit(1) ; 00015 DSET_load(dset) ; if( !DSET_LOADED(dset) ) exit(1) ; 00016 grim = DSET_BRICK(dset,0) ; 00017 sx = DSET_NX(dset) * DSET_DX(dset) ; 00018 sy = DSET_NY(dset) * DSET_DY(dset) ; 00019 sz = DSET_NZ(dset) * DSET_DZ(dset) ; 00020 00021 oset = THD_open_dataset( argv[2] ) ; 00022 if( oset == NULL ) exit(1) ; 00023 DSET_load(oset) ; if( !DSET_LOADED(oset) ) exit(1) ; 00024 opim = DSET_BRICK(oset,0) ; 00025 00026 { int nvox = DSET_NVOX(oset) , ii ; 00027 byte * opar = MRI_BYTE_PTR(opim) ; 00028 for( ii=0 ; ii < nvox/4 ; ii++ ) opar[ii] = 0 ; 00029 for( ii=3*nvox/4 ; ii < nvox ; ii++ ) opar[ii] = 0 ; 00030 } 00031 00032 arh = new_AREN_renderer() ; AREN_be_verbose(arh) ; 00033 00034 #if 1 00035 AREN_set_graybrick( arh , grim ) ; 00036 #else 00037 AREN_set_rgbbricks( arh , grim,grim,grim ) ; 00038 #endif 00039 00040 AREN_set_opabrick ( arh , opim ) ; 00041 AREN_set_size ( arh , sx,sy,sz ) ; 00042 00043 #if 0 00044 DSET_delete(dset) ; DSET_delete(oset) ; 00045 #endif 00046 00047 for( nopt=3 ; nopt < argc-1 ; nopt+=2 ){ 00048 theta = strtod( argv[nopt] , NULL ) ; 00049 phi = strtod( argv[nopt+1] , NULL ) ; 00050 AREN_set_viewpoint( arh , theta , phi ) ; 00051 renim = AREN_render( arh , 256 ) ; 00052 sprintf(fname,"aren.%03d.pnm",(nopt-1)/2) ; 00053 mri_write_pnm( fname , renim ) ; 00054 mri_free(renim) ; 00055 } 00056 00057 exit(0) ; 00058 } |