|
|
compute the overall minimum and maximum voxel values for a dataset
Definition at line 31 of file 3dinfo_OLD.c.
References THD_3dim_dataset::anat_parent_idcode, THD_3dim_dataset::anat_parent_name, argc, THD_datablock::brick_lab, THD_statistics::bstat, MCW_idcode::date, THD_3dim_dataset::daxes, DAXES_NUM, THD_3dim_dataset::dblk, DSET_BRICK_FACTOR, DSET_BRICK_KEYWORDS, DSET_BRICK_STATCODE, DSET_BRICK_STATPAR, DSET_BRICK_TYPE, DSET_KEYWORDS, DSET_NUM_TIMES, DSET_NVALS_PER_TIME, THD_timeaxis::dz_sl, FUNC_IS_STAT, THD_3dim_dataset::func_type, THD_3dim_dataset::idcode, ISANAT, ISFUNC, ISVALID_STATISTIC, ISZERO_IDCODE, THD_3dim_dataset::label1, LOAD_FVEC3, THD_brick_stats::max, THD_brick_stats::min, n1, n2, THD_timeaxis::nsl, THD_datablock::nvals, THD_dataxes::nxx, THD_dataxes::nyy, THD_dataxes::nzz, THD_3dim_dataset::self_name, THD_3dim_dataset::stat_aux, THD_3dim_dataset::stats, MCW_idcode::str, Syntax(), THD_3dim_dataset::taxis, THD_3dmm_to_dicomm(), THD_delete_3dim_dataset(), THD_open_one_dataset(), THD_timeaxis::toff_sl, THD_timeaxis::ttdel, THD_timeaxis::units_type, UNITS_TYPE_LABEL, THD_3dim_dataset::warp_parent_idcode, THD_3dim_dataset::warp_parent_name, THD_dataxes::xxdel, THD_dataxes::xxorg, THD_dataxes::xxorient, THD_fvec3::xyz, THD_dataxes::yydel, THD_dataxes::yyorg, THD_dataxes::yyorient, THD_dataxes::zzdel, THD_dataxes::zzorg, and THD_dataxes::zzorient.
00032 {
00033 THD_3dim_dataset * dset ;
00034 THD_dataxes * daxes ;
00035 THD_fvec3 fv1 , fv2 , fv3 ;
00036 THD_ivec3 iv ;
00037 int iarg , ival , ntimes , nval_per , n1,n2,n3 , verbose = 0 , kv,npar ;
00038 float tf ;
00039
00040 static char * RR="[R]" , * LL="[L]" ,
00041 * PP="[P]" , * AA="[A]" ,
00042 * SS="[S]" , * II="[I]" , * ZZ=" " ;
00043 char * xlbot , * xltop , * ylbot , * yltop , * zlbot , * zltop , * cpt ;
00044 char str[256] ;
00045 int nstr ;
00046
00047 #define XLAB(xl,xv) ((xl) = ((xv)==0.0) ? (ZZ) : ( ((xv)<0.0) ? (RR) : (LL) ))
00048 #define YLAB(yl,yv) ((yl) = ((yv)==0.0) ? (ZZ) : ( ((yv)<0.0) ? (AA) : (PP) ))
00049 #define ZLAB(zl,zv) ((zl) = ((zv)==0.0) ? (ZZ) : ( ((zv)<0.0) ? (II) : (SS) ))
00050
00051 #define FSWAP(x,y) (tf=(x),(x)=(y),(y)=tf)
00052
00053 if( argc < 2 || strncmp(argv[1],"-help",4) == 0 ) Syntax() ;
00054
00055 iarg = 1 ;
00056
00057 if( strcmp(argv[iarg],"-v") == 0 ){ verbose = 1 ; iarg++ ; }
00058
00059 for( ; iarg < argc ; iarg++ ){
00060 dset = THD_open_one_dataset( argv[iarg] ) ;
00061 if( dset == NULL ) continue ;
00062 daxes = dset->daxes ;
00063
00064 printf("\n"
00065 "Dataset File: %s\n"
00066 #if 0
00067 "Dataset Name: %s\n"
00068 "Dataset Label: %s\n"
00069 #endif
00070 , argv[iarg]
00071 #if 0
00072 , dset->self_name , dset->label1
00073 #endif
00074 ) ;
00075
00076 #ifndef OMIT_DATASET_IDCODES
00077 printf("Identifier Code: %s Creation Date: %s\n" ,
00078 dset->idcode.str , dset->idcode.date ) ;
00079 #endif
00080
00081 if( ISANAT(dset) ){
00082 printf("Dataset Type: %s (-%s)\n",
00083 ANAT_typestr[dset->func_type] , ANAT_prefixstr[dset->func_type] ) ;
00084 } else {
00085 printf("Dataset Type: %s (-%s)\n",
00086 FUNC_typestr[dset->func_type] , FUNC_prefixstr[dset->func_type] ) ;
00087 }
00088
00089 cpt = DSET_KEYWORDS(dset) ;
00090 if( cpt != NULL && cpt[0] != '\0' )
00091 printf("Keywords: %s\n" , cpt ) ;
00092
00093 #ifdef OMIT_DATASET_IDCODES
00094 if( strlen(dset->anat_parent_name) > 0 )
00095 printf("Anatomy Parent: %s\n" , dset->anat_parent_name ) ;
00096
00097 if( strlen(dset->warp_parent_name) > 0 )
00098 printf("Warp Parent: %s\n" , dset->warp_parent_name ) ;
00099 #else
00100 if( ! ISZERO_IDCODE(dset->anat_parent_idcode) )
00101 printf("Anatomy Parent: %s [%s]\n" ,
00102 dset->anat_parent_name , dset->anat_parent_idcode.str ) ;
00103 else if( strlen(dset->anat_parent_name) > 0 )
00104 printf("Anatomy Parent: %s\n" , dset->anat_parent_name ) ;
00105
00106 if( ! ISZERO_IDCODE(dset->warp_parent_idcode) )
00107 printf("Warp Parent: %s [%s]\n" ,
00108 dset->warp_parent_name , dset->warp_parent_idcode.str) ;
00109 else if( strlen(dset->warp_parent_name) > 0 )
00110 printf("Warp Parent: %s\n" , dset->warp_parent_name ) ;
00111 #endif
00112
00113 printf(
00114 "Data Axes Orientation:\n"
00115 " first (x) = %s\n"
00116 " second (y) = %s\n"
00117 " third (z) = %s [-orient %c%c%c]\n" ,
00118 ORIENT_typestr[daxes->xxorient] ,
00119 ORIENT_typestr[daxes->yyorient] ,
00120 ORIENT_typestr[daxes->zzorient] ,
00121 ORIENT_typestr[daxes->xxorient][0] ,
00122 ORIENT_typestr[daxes->yyorient][0] ,
00123 ORIENT_typestr[daxes->zzorient][0] ) ;
00124
00125 LOAD_FVEC3(fv1 , daxes->xxorg , daxes->yyorg , daxes->zzorg) ;
00126 fv1 = THD_3dmm_to_dicomm( dset , fv1 ) ;
00127
00128 LOAD_FVEC3(fv2 , daxes->xxorg + (daxes->nxx-1)*daxes->xxdel ,
00129 daxes->yyorg + (daxes->nyy-1)*daxes->yydel ,
00130 daxes->zzorg + (daxes->nzz-1)*daxes->zzdel ) ;
00131 fv2 = THD_3dmm_to_dicomm( dset , fv2 ) ;
00132
00133 if( fv1.xyz[0] > fv2.xyz[0] ) FSWAP( fv1.xyz[0] , fv2.xyz[0] ) ;
00134 if( fv1.xyz[1] > fv2.xyz[1] ) FSWAP( fv1.xyz[1] , fv2.xyz[1] ) ;
00135 if( fv1.xyz[2] > fv2.xyz[2] ) FSWAP( fv1.xyz[2] , fv2.xyz[2] ) ;
00136
00137 LOAD_FVEC3(fv3 , daxes->xxdel , daxes->yydel , daxes->zzdel) ;
00138 fv3 = THD_3dmm_to_dicomm( dset , fv3 ) ;
00139
00140 XLAB(xlbot,fv1.xyz[0]) ; YLAB(ylbot,fv1.xyz[1]) ; ZLAB(zlbot,fv1.xyz[2]) ;
00141 XLAB(xltop,fv2.xyz[0]) ; YLAB(yltop,fv2.xyz[1]) ; ZLAB(zltop,fv2.xyz[2]) ;
00142
00143 n1 = DAXES_NUM(daxes,daxes->xxorient) ;
00144 n2 = DAXES_NUM(daxes,daxes->yyorient) ;
00145 n3 = DAXES_NUM(daxes,daxes->zzorient) ;
00146
00147 printf(
00148 "R-to-L extent: %9.3f %s -to- %9.3f %s -step- %9.3f mm [%3d voxels]\n"
00149 "A-to-P extent: %9.3f %s -to- %9.3f %s -step- %9.3f mm [%3d voxels]\n"
00150 "I-to-S extent: %9.3f %s -to- %9.3f %s -step- %9.3f mm [%3d voxels]\n" ,
00151 fv1.xyz[0],xlbot , fv2.xyz[0],xltop , fabs(fv3.xyz[0]) , n1 ,
00152 fv1.xyz[1],ylbot , fv2.xyz[1],yltop , fabs(fv3.xyz[1]) , n2 ,
00153 fv1.xyz[2],zlbot , fv2.xyz[2],zltop , fabs(fv3.xyz[2]) , n3 ) ;
00154
00155 ntimes = DSET_NUM_TIMES(dset) ;
00156 nval_per = DSET_NVALS_PER_TIME(dset) ;
00157 if( ntimes > 1 ){
00158 printf(
00159 "Number of time steps = %d Number of values at each pixel = %d\n",
00160 ntimes , nval_per ) ;
00161
00162 printf( "Time step = %.3f (%s)" ,
00163 dset->taxis->ttdel ,
00164 UNITS_TYPE_LABEL(dset->taxis->units_type) ) ;
00165 if( dset->taxis->nsl > 0 )
00166 printf(" Number time-offset slices = %d Thickness = %.3f",
00167 dset->taxis->nsl , fabs(dset->taxis->dz_sl) ) ;
00168 printf("\n") ;
00169
00170 if( verbose && dset->taxis->nsl > 0 ){
00171 printf("Time-offsets per slice:") ;
00172 for( ival=0 ; ival < dset->taxis->nsl ; ival++ )
00173 printf( " %.3f" , dset->taxis->toff_sl[ival] ) ;
00174 printf("\n") ;
00175 }
00176 } else {
00177 printf(
00178 "Number of values stored at each pixel = %d\n" , nval_per ) ;
00179 }
00180
00181 if( verbose && ntimes > 1 ) nval_per = dset->dblk->nvals ;
00182
00183 for( ival=0 ; ival < nval_per ; ival++ ){
00184 sprintf( str ,
00185 " -- At sub-brick #%d [%s] datum type is %s" ,
00186 ival , dset->dblk->brick_lab[ival] ,
00187 MRI_TYPE_name[DSET_BRICK_TYPE(dset,ival)] ) ;
00188 nstr = strlen(str) ;
00189
00190 tf = DSET_BRICK_FACTOR(dset,ival) ;
00191
00192 if( ISVALID_STATISTIC(dset->stats) ){
00193 if( tf != 0.0 )
00194 sprintf( str+nstr ,
00195 ":%13.6g to %13.6g [internal]\n"
00196 "%*s[*%13.6g] %13.6g to %13.6g [scaled]\n" ,
00197 dset->stats->bstat[ival].min/tf ,
00198 dset->stats->bstat[ival].max/tf ,
00199 nstr-16," " , tf ,
00200 dset->stats->bstat[ival].min , dset->stats->bstat[ival].max ) ;
00201 else
00202 sprintf( str+nstr , ":%13.6g to %13.6g\n" ,
00203 dset->stats->bstat[ival].min , dset->stats->bstat[ival].max ) ;
00204 } else if( tf != 0.0 ){
00205 sprintf( str+nstr , " [*%g]\n",tf) ;
00206 } else {
00207 sprintf( str+nstr , "\n") ;
00208 }
00209 printf(str) ;
00210
00211
00212
00213 kv = DSET_BRICK_STATCODE(dset,ival) ;
00214 if( FUNC_IS_STAT(kv) ){
00215 printf(" statcode = %s",FUNC_prefixstr[kv] ) ;
00216 npar = FUNC_need_stat_aux[kv] ;
00217 if( npar > 0 ){
00218 printf("; statpar =") ;
00219 for( kv=0 ; kv < npar ; kv++ )
00220 printf(" %g",DSET_BRICK_STATPAR(dset,ival,kv)) ;
00221 }
00222 printf("\n") ;
00223 }
00224
00225 cpt = DSET_BRICK_KEYWORDS(dset,ival) ;
00226 if( cpt != NULL && cpt[0] != '\0' )
00227 printf(" keywords = %s\n",cpt) ;
00228 }
00229
00230
00231
00232 if( ISFUNC(dset) && FUNC_need_stat_aux[dset->func_type] > 0 ){
00233 printf("Auxiliary functional statistical parameters:\n %s\n",
00234 FUNC_label_stat_aux[dset->func_type] ) ;
00235 for( ival=0 ; ival < FUNC_need_stat_aux[dset->func_type] ; ival++ )
00236 printf(" %g",dset->stat_aux[ival]) ;
00237 printf("\n") ;
00238 }
00239
00240 THD_delete_3dim_dataset( dset , False ) ;
00241 }
00242 exit(0) ;
00243 }
|
