tac2nii.c File Reference

Create NIfTI format PET image with contents from specified TAC file. More...

Go to the source code of this file.

Functions
int	simRCDims (const long long N, const int mxdim, const int mydim, const int mzdim, int *rnx, int *rny, int *rnz, int *rcxdim, int *rcydim, int *rczdim)

Detailed Description

Create NIfTI format PET image with contents from specified TAC file.

Copyright

(c) Turku PET Centre

Author

Vesa Oikonen

Definition in file tac2nii.c.

Function Documentation

simRCDims()

int simRCDims	(	const long long	N,
		const int	mxdim,
		const int	mydim,
		const int	mzdim,
		int *	rnx,
		int *	rny,
		int *	rnz,
		int *	rcxdim,
		int *	rcydim,
		int *	rczdim )

Determine dimensions of N rectangular cuboids fitting specified 3D matrix.

Returns

Returns 0 when successful.

Parameters

N	Nr of rectangular cuboids.
mxdim	3D matrix xdim.
mydim	3D matrix ydim.
mzdim	3D matrix zdim.
rnx	Nr of cuboids in x dim.
rny	Nr of cuboids in y dim.
rnz	Nr of cuboids in z dim.
rcxdim	Rectangular cuboid xdim.
rcydim	Rectangular cuboid ydim.
rczdim	Rectangular cuboid zdim.

Definition at line 476 of file tac2nii.c.

  {
  if(rnx!=NULL) *rnx=0;
  if(rny!=NULL) *rny=0;
  if(rnz!=NULL) *rnz=0;
  if(rcxdim!=NULL) *rcxdim=0;
  if(rcydim!=NULL) *rcydim=0;
  if(rczdim!=NULL) *rczdim=0;
  if(N<1 || mxdim<1 || mydim<1 || mzdim<1) return(1);
 
  int xd, yd, zd;
  xd=mxdim; yd=mydim; zd=mzdim;
 
  long long nx=1, ny=1, nz=1;
  long long maxn=nx*ny*nz;
  long long n=N-maxn;
  while(n>0 && (nz<mzdim || ny<mydim || nx<mxdim)) {
//printf("  n=%lld nx=%lld ny=%lld nz=%lld\n", n, nx, ny, nz);
    if(n>0 && nx<mxdim) {
      nx++; maxn=nx*ny*nz; n=N-maxn;
    }
    if(n>0 && ny<mydim) {
      ny++; maxn=nx*ny*nz; n=N-maxn;
    }
    if(n>0 && nz<mzdim) {
      nz++; maxn=nx*ny*nz; n=N-maxn;
    }
  }
  if(n>0) return(2);
 
  xd/=nx; yd/=ny; zd/=nz;
 
  if(rnx!=NULL) *rnx=nx;
  if(rny!=NULL) *rny=ny;
  if(rnz!=NULL) *rnz=nz;
  if(rcxdim!=NULL) *rcxdim=xd;
  if(rcydim!=NULL) *rcydim=yd;
  if(rczdim!=NULL) *rczdim=zd;
  return(0);
}