Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036 #include <stdio.h>
00037 #include <stdlib.h>
00038 #include <math.h>
00039 #include <time.h>
00040 #include <string.h>
00041 #include <ctype.h>
00042
00043 #include "halflife.h"
00044
00045 #include "include/img.h"
00046 #include "include/imgdecay.h"
00047
00048
00049
00059 int imgDecayCorrection(IMG *image, int mode) {
00060 int pi, fi, i, j;
00061 float lambda;
00062 float cf, dur;
00063
00064
00065 if(image->status!=IMG_STATUS_OCCUPIED) return(1);
00066 if(image->isotopeHalflife<=0.0) return(1);
00067
00068 if(mode==1 && image->decayCorrected!=0) return(2);
00069 if(mode==0 && image->decayCorrected==0) return(2);
00070
00071
00072 for(fi=0; fi<image->dimt; fi++) {
00073 dur=image->end[fi]-image->start[fi];
00074 if(image->end[fi]>0.0) {
00075 if(mode==0 && image->decayCorrFactor[fi]>1.000001) {
00076
00077 cf=1.0/image->decayCorrFactor[fi];
00078 } else {
00079 lambda=hl2lambda(image->isotopeHalflife); if(lambda<0.0) return(1);
00080 if(mode==0) lambda=-lambda;
00081 if(fi==image->dimt-1 && image->end[fi]<=0.0) return(3);
00082 cf=hlLambda2factor_float(lambda, image->start[fi], dur);
00083 }
00084 if(IMG_TEST) printf("applied_dc_factor[%d] := %g\n", fi+1, cf);
00085
00086 if(mode==0) {
00087 image->decayCorrFactor[fi]=1.0;
00088 } else {
00089 image->decayCorrFactor[fi]=cf;
00090 }
00091
00092 for(pi=0; pi<image->dimz; pi++)
00093 for(i=0; i<image->dimy; i++)
00094 for(j=0; j<image->dimx; j++)
00095 image->m[pi][i][j][fi]*=cf;
00096 image->decayCorrected=mode;
00097 }
00098 }
00099 return(0);
00100 }
00101
00102
00103
00110 char *imgIsotope(IMG *img) {
00111 return(hlIsotopeCode(hlIsotopeFromHalflife(img->isotopeHalflife/60.0)));
00112 }
00113
00114
00115
00126 int imgSetDecayCorrFactors(IMG *image, int mode) {
00127 int fi;
00128 float lambda, cf, dur;
00129
00130
00131 if(image->status!=IMG_STATUS_OCCUPIED) return(1);
00132 if(image->isotopeHalflife<=0.0) return(1);
00133
00134
00135 if(mode!=0 && image->end[image->dimt-1]<=0.0) return(3);
00136
00137
00138 for(fi=0; fi<image->dimt; fi++) {
00139 if(mode==0) {
00140 image->decayCorrFactor[fi]=1.0;
00141 } else {
00142 dur=image->end[fi]-image->start[fi];
00143 if(image->end[fi]>0.0) {
00144 lambda=hl2lambda(image->isotopeHalflife); if(lambda<0.0) return(2);
00145 cf=hlLambda2factor_float(lambda, image->start[fi], dur);
00146 image->decayCorrFactor[fi]=cf;
00147 }
00148 }
00149 }
00150 image->decayCorrected=mode;
00151 return(0);
00152 }
00153
00154
00155
00156