/***************************************************************************** Major portions of this software are copyrighted by the Medical College of Wisconsin, 1994-2000, and are released under the Gnu General Public License, Version 2. See the file README.Copyright for details. ******************************************************************************/ #include "mrilib.h" #include "thd.h" static int use_3D_format = 0 ; void THD_use_3D_format( int uu ){ use_3D_format = uu; } /*---------------------------------------------------------------- this routine writes all the data from the dataset into the datablock attributes, then writes the datablock to disk 29 April 1998: erase attributes that are unused, so that they won't be left over from a previous life ------------------------------------------------------------------*/ #define IFILL_DIM 100 #define FFILL_DIM (MAPPING_LINEAR_FSIZE*20) #define IFILL -999 #define FFILL -999999.0 #define ITFILL(p,q) for(ii=p;iidblk) || ! ISVALID_DISKPTR(dset->dblk->diskptr) ) RETURN(False) ; /* Can only write AFNI formatted datasets */ if( DSET_IS_MINC(dset) ) RETURN(False) ; /* 29 Oct 2001 */ if( DSET_IS_MASTERED(dset) ) RETURN(False) ; /* 11 Jan 1999 */ if( DSET_IS_ANALYZE(dset) ) RETURN(False) ; /* 27 Aug 2002 */ if( DSET_IS_NIFTI(dset) ) RETURN(False) ; /* 28 Aug 2003 */ if( DSET_IS_CTFMRI(dset) ) RETURN(False) ; /* 05 Dec 2002 */ if( DSET_IS_CTFSAM(dset) ) RETURN(False) ; /* 05 Dec 2002 */ if( DSET_IS_TCAT(dset) ) RETURN(False) ; /* 05 Aug 2004 */ if( DSET_IS_VOLUMES(dset) && write_brick ) RETURN(False) ; /* 20 Jun 2002 */ if( DSET_IS_1D(dset) || ( DSET_NY(dset)==1 && DSET_NZ(dset)==1 ) ){ /* 04 Mar 2003 */ THD_write_1D( new_sessname , new_prefixname , dset ) ; RETURN(True) ; } blk = dset->dblk ; daxes = dset->daxes ; /* always used fixed daxes */ /*------------------------------*/ /*----- change filenames? -----*/ THD_init_diskptr_names( blk->diskptr , new_sessname , NULL , new_prefixname , dset->view_type , True ) ; /*----- write TYPESTRING attribute -----*/ THD_set_string_atr( blk , ATRNAME_TYPESTRING , DATASET_typestr[dset->type] ) ; /*----- write IDCODE attributes -----*/ THD_set_string_atr( blk , ATRNAME_IDSTRING , dset->idcode.str ) ; THD_set_string_atr( blk , ATRNAME_IDDATE , dset->idcode.date ) ; if( ! ISZERO_IDCODE(dset->anat_parent_idcode) ) THD_set_string_atr( blk , ATRNAME_IDANATPAR , dset->anat_parent_idcode.str ) ; else THD_erase_one_atr( blk , ATRNAME_IDANATPAR ) ; if( ! ISZERO_IDCODE(dset->warp_parent_idcode) ) THD_set_string_atr( blk , ATRNAME_IDWARPPAR , dset->warp_parent_idcode.str ) ; else THD_erase_one_atr( blk , ATRNAME_IDWARPPAR ) ; /*----- write SCENE_TYPE attribute -----*/ itemp[0] = dset->view_type ; itemp[1] = dset->func_type ; itemp[2] = dset->type ; ITFILL(3,ATRSIZE_SCENE_TYPE) ; THD_set_int_atr( blk , ATRNAME_SCENE_TYPE , ATRSIZE_SCENE_TYPE , itemp ) ; /*----- write data labels -----*/ if( strlen(dset->self_name) == 0 ) DSET_FIX_NAMES(dset) ; THD_set_string_atr( blk , ATRNAME_LABEL1 , dset->label1 ) ; THD_set_string_atr( blk , ATRNAME_LABEL2 , dset->label2 ) ; THD_set_string_atr( blk , ATRNAME_DATANAME , dset->self_name ) ; if( dset->keywords != NULL ) THD_set_string_atr( blk , ATRNAME_KEYWORDS , dset->keywords ) ; else THD_erase_one_atr( blk , ATRNAME_KEYWORDS ) ; /*----- write parent names, if they exist -----*/ if( strlen(dset->warp_parent_name) > 0 ) THD_set_string_atr( blk , ATRNAME_WARP_PARENT , dset->warp_parent_name ) ; else THD_erase_one_atr( blk , ATRNAME_WARP_PARENT ) ; if( strlen(dset->anat_parent_name) > 0 ) THD_set_string_atr( blk , ATRNAME_ANATOMY_PARENT , dset->anat_parent_name ) ; else THD_erase_one_atr( blk , ATRNAME_ANATOMY_PARENT ) ; /*----- write axes orientation -----*/ itemp[0] = daxes->xxorient ; itemp[1] = daxes->yyorient ; itemp[2] = daxes->zzorient ; ITFILL(3,ATRSIZE_ORIENT_SPECIFIC) ; THD_set_int_atr( blk , ATRNAME_ORIENT_SPECIFIC , ATRSIZE_ORIENT_SPECIFIC , itemp ) ; /*----- write axes origin -----*/ ftemp[0] = daxes->xxorg ; ftemp[1] = daxes->yyorg ; ftemp[2] = daxes->zzorg ; FTFILL(3,ATRSIZE_ORIGIN) ; THD_set_float_atr( blk , ATRNAME_ORIGIN , ATRSIZE_ORIGIN , ftemp ) ; /*----- write axes spacings -----*/ ftemp[0] = daxes->xxdel ; ftemp[1] = daxes->yydel ; ftemp[2] = daxes->zzdel ; FTFILL(3,ATRSIZE_DELTA) ; THD_set_float_atr( blk , ATRNAME_DELTA , ATRSIZE_DELTA , ftemp ) ; /*----- write markers, if present -----*/ if( dset->markers != NULL ){ for( ii=0 ; ii < MARKS_MAXNUM ; ii++ ){ /* put bad data in */ if( !dset->markers->valid[ii] ) /* invalid markers */ dset->markers->xyz[ii][0] = dset->markers->xyz[ii][1] = dset->markers->xyz[ii][2] = FFILL ; } THD_set_float_atr( blk , ATRNAME_MARKSXYZ , ATRSIZE_MARKSXYZ , &(dset->markers->xyz[0][0]) ) ; THD_set_char_atr( blk , ATRNAME_MARKSLAB , ATRSIZE_MARKSLAB , &(dset->markers->label[0][0]) ) ; THD_set_char_atr( blk , ATRNAME_MARKSHELP , ATRSIZE_MARKSHELP , &(dset->markers->help[0][0]) ) ; THD_set_int_atr( blk , ATRNAME_MARKSFLAG , ATRSIZE_MARKSFLAG , &(dset->markers->aflags[0]) ) ; } else { THD_erase_one_atr( blk , ATRNAME_MARKSXYZ ) ; THD_erase_one_atr( blk , ATRNAME_MARKSLAB ) ; THD_erase_one_atr( blk , ATRNAME_MARKSHELP ) ; THD_erase_one_atr( blk , ATRNAME_MARKSFLAG ) ; } /*----- write warp, if present -----*/ if( dset->warp != NULL ){ int wdata_size = 0 ; switch( dset->warp->type ){ case WARP_AFFINE_TYPE:{ THD_affine_warp * ww = (THD_affine_warp *) dset->warp ; itemp[0] = WARP_AFFINE_TYPE ; itemp[1] = ww->resam_type ; wdata_size = MAPPING_LINEAR_FSIZE ; COPY_FROM_STRUCT( ww->warp , MAPPING_LINEAR_FSTART , float , ftemp , MAPPING_LINEAR_FSIZE ) ; } break ; case WARP_TALAIRACH_12_TYPE:{ THD_talairach_12_warp * ww = (THD_talairach_12_warp *) dset->warp ; int iw , ioff ; itemp[0] = WARP_TALAIRACH_12_TYPE ; itemp[1] = ww->resam_type ; wdata_size = WARP_TALAIRACH_12_SIZE ; for( iw=0 ; iw < 12 ; iw++ ){ ioff = iw * MAPPING_LINEAR_FSIZE ; COPY_FROM_STRUCT( ww->warp[iw] , MAPPING_LINEAR_FSTART , float , &(ftemp[ioff]) , MAPPING_LINEAR_FSIZE ) ; } } break ; } /* end of switch on warp type */ ITFILL(2,ATRSIZE_WARP_TYPE) ; THD_set_int_atr( blk , ATRNAME_WARP_TYPE , ATRSIZE_WARP_TYPE , itemp ) ; THD_set_float_atr( blk , ATRNAME_WARP_DATA , wdata_size , ftemp ) ; } else { /* no warp exists */ THD_erase_one_atr( blk , ATRNAME_WARP_TYPE ) ; THD_erase_one_atr( blk , ATRNAME_WARP_DATA ) ; } /* end of if warp exists */ /*----- if statistics exist, write them out (modern style of floats) -----*/ if( ISVALID_STATISTIC(dset->stats) ){ float * tfil ; int qq ; tfil = (float *) malloc( sizeof(float) * 2*blk->nvals ) ; if( tfil != NULL ){ for( qq=0 ; qq < dset->stats->nbstat ; qq++ ){ tfil[2*qq] = dset->stats->bstat[qq].min ; tfil[2*qq+1] = dset->stats->bstat[qq].max ; } THD_set_float_atr( blk, ATRNAME_BRICK_STATS, 2*dset->stats->nbstat, tfil ) ; free( tfil ) ; } } else { THD_erase_one_atr( blk , ATRNAME_BRICK_STATS ) ; } /*----- if auxiliary statistics data exists, write them out too -----*/ for( ii=MAX_STAT_AUX-1 ; ii>=0 && dset->stat_aux[ii]==0.0 ; ii-- ) ; /* nada */ if( ii >= 0 ){ /* ii = largest index where stat_aux != 0 */ ii++ ; /* number of stat_aux values to save */ THD_set_float_atr( blk , ATRNAME_STAT_AUX , ii , dset->stat_aux ) ; } else { THD_erase_one_atr( blk , ATRNAME_STAT_AUX ) ; } /*----- if time-dependent data, write that stuff out too -----*/ if( dset->taxis != NULL ){ itemp[0] = dset->taxis->ntt ; itemp[1] = dset->taxis->nsl ; itemp[2] = dset->taxis->units_type ; /* 21 Oct 1996 */ ITFILL(3,ATRSIZE_TAXIS_NUMS) ; THD_set_int_atr( blk , ATRNAME_TAXIS_NUMS , ATRSIZE_TAXIS_NUMS , itemp ) ; ftemp[0] = dset->taxis->ttorg ; ftemp[1] = dset->taxis->ttdel ; ftemp[2] = dset->taxis->ttdur ; ftemp[3] = dset->taxis->zorg_sl ; ftemp[4] = dset->taxis->dz_sl ; FTFILL(5,ATRSIZE_TAXIS_FLOATS) ; THD_set_float_atr( blk , ATRNAME_TAXIS_FLOATS , ATRSIZE_TAXIS_FLOATS , ftemp ) ; if( dset->taxis->toff_sl != NULL ) THD_set_float_atr( blk , ATRNAME_TAXIS_OFFSETS , dset->taxis->nsl , dset->taxis->toff_sl ) ; } else { THD_erase_one_atr( blk , ATRNAME_TAXIS_NUMS ) ; THD_erase_one_atr( blk , ATRNAME_TAXIS_FLOATS ) ; THD_erase_one_atr( blk , ATRNAME_TAXIS_OFFSETS ) ; } /*----- 23 Oct 1998: write out tagset, if present -----*/ #undef NFPER #define NFPER 5 #define TF(i,j) vtag[(i)*NFPER+(j)] if( dset->tagset != NULL && dset->tagset->num > 0 ){ int ii , ntag=dset->tagset->num , tlen,ilen,jj ; float * vtag = (float *) malloc( sizeof(float) * (NFPER*ntag) ) ; char * ctag ; /* set the counts */ itemp[0] = ntag ; itemp[1] = NFPER ; THD_set_int_atr( blk , ATRNAME_TAGSET_NUM , 2 , itemp ) ; /* set the values */ for( ii=0 ; ii < ntag ; ii++ ){ TF(ii,0) = dset->tagset->tag[ii].x ; TF(ii,1) = dset->tagset->tag[ii].y ; TF(ii,2) = dset->tagset->tag[ii].z ; TF(ii,3) = dset->tagset->tag[ii].val ; if( dset->tagset->tag[ii].set ) TF(ii,4) = dset->tagset->tag[ii].ti ; else TF(ii,4) = -1.0 ; } THD_set_float_atr( blk , ATRNAME_TAGSET_FLOATS , NFPER*ntag , vtag ) ; free(vtag) ; /* set the labels */ tlen = 4 ; /* a little slop space */ for( ii=0 ; ii < ntag ; ii++ ) tlen += strlen( dset->tagset->tag[ii].label ) + 1 ; /* +1 for the '\0' */ ctag = (char *) malloc( sizeof(char) * tlen ) ; /* to hold all labels */ jj = 0 ; for( ii=0 ; ii < ntag ; ii++ ){ ilen = strlen( dset->tagset->tag[ii].label ) + 1 ; memcpy( ctag+jj , dset->tagset->tag[ii].label , ilen ) ; jj += ilen ; } THD_set_char_atr( blk , ATRNAME_TAGSET_LABELS , tlen , ctag ) ; free(ctag) ; } else { THD_erase_one_atr( blk , ATRNAME_TAGSET_NUM ) ; THD_erase_one_atr( blk , ATRNAME_TAGSET_LABELS ) ; THD_erase_one_atr( blk , ATRNAME_TAGSET_FLOATS ) ; } #undef NFPER #undef TF if( DSET_IS_3D(dset) || use_3D_format ){ /* 21 Mar 2003 */ THD_write_3D( NULL, NULL, dset ) ; RETURN(True) ; } /*----- write datablock to disk -----*/ RETURN( THD_write_datablock(blk,write_brick) ) ; }