/* ----------------------------- MNI Header -----------------------------------
@NAME : setup_image_variables
@INPUT : inmincid - id of input minc file (MI_ERROR if no file)
mincid - id of output minc file
ndims - number of dimensions
dim - list of dimension ids
@OUTPUT : (nothing)
@RETURNS : (nothing)
@DESCRIPTION: Routine to set up image, image-max and image-min variables
in the output minc file
@METHOD :
@GLOBALS :
@CALLS :
@CREATED : August 26, 1993 (Peter Neelin)
@MODIFIED :
---------------------------------------------------------------------------- */
static void setup_image_variables(int inmincid, int mincid,
int ndims, int dim[])
{
int imgid, maxid, minid;
static double valid_range[2] = {0.0, 255.0};
/* Create the image max and min variables (varying over slices) */
maxid = micreate_std_variable(mincid, MIimagemax, NC_DOUBLE, 1, dim);
minid = micreate_std_variable(mincid, MIimagemin, NC_DOUBLE, 1, dim);
if (inmincid != MI_ERROR) {
(void) micopy_all_atts(inmincid, ncvarid(inmincid, MIimagemax),
mincid, maxid);
(void) micopy_all_atts(inmincid, ncvarid(inmincid, MIimagemin),
mincid, minid);
}
/* Create the image variable, copy attributes, set the signtype attribute,
set the valid range attribute and delete valid max/min attributes */
imgid = micreate_std_variable(mincid, MIimage, NC_BYTE, ndims, dim);
if (inmincid != MI_ERROR) {
(void) micopy_all_atts(inmincid, ncvarid(inmincid, MIimage),
mincid, imgid);
ncopts = 0;
(void) ncattdel(mincid, imgid, MIvalid_max);
(void) ncattdel(mincid, imgid, MIvalid_min);
ncopts = NC_OPTS_VAL;
}
(void) miattputstr(mincid, imgid, MIsigntype, MI_UNSIGNED);
(void) miset_valid_range(mincid, imgid, valid_range);
}
void
copy_init(struct conversion_info *ci_ptr)
{
ci_ptr->frame_buffer = malloc(ci_ptr->frame_nbytes);
if (ci_ptr->frame_buffer == NULL) {
message(MSG_FATAL, "Out of memory\n");
exit(-1);
}
/* Create the image, imagemax, and imagemin variables.
*/
micreate_std_variable(ci_ptr->mnc_fd, MIimagemax, NC_DOUBLE,
1, ci_ptr->dim_ids);
micreate_std_variable(ci_ptr->mnc_fd, MIimagemin, NC_DOUBLE,
1, ci_ptr->dim_ids);
micreate_std_variable(ci_ptr->mnc_fd, MIimage, ci_ptr->minc_type,
ci_ptr->dim_count + 1, ci_ptr->dim_ids);
/* Set up the dimension step and start values. Because of the microPET
* data orientation, we set Z and Y to be the inverse of the norm, to
* put the animal's nose at the top of the display.
* TODO: allow this behavior to be controlled on the command line.
*/
miattputdbl(ci_ptr->mnc_fd, ncvarid(ci_ptr->mnc_fd, _dimnames[DIM_Z]),
MIstep, -ci_ptr->dim_steps[DIM_Z]);
miattputdbl(ci_ptr->mnc_fd, ncvarid(ci_ptr->mnc_fd, _dimnames[DIM_X]),
MIstep, ci_ptr->dim_steps[DIM_X]);
miattputdbl(ci_ptr->mnc_fd, ncvarid(ci_ptr->mnc_fd, _dimnames[DIM_Y]),
MIstep, -ci_ptr->dim_steps[DIM_Y]);
miattputdbl(ci_ptr->mnc_fd, ncvarid(ci_ptr->mnc_fd, _dimnames[DIM_Z]),
MIstart, ci_ptr->dim_steps[DIM_Z] * ci_ptr->dim_lengths[DIM_Z]);
miattputdbl(ci_ptr->mnc_fd, ncvarid(ci_ptr->mnc_fd, _dimnames[DIM_X]),
MIstart, 0.0);
miattputdbl(ci_ptr->mnc_fd, ncvarid(ci_ptr->mnc_fd, _dimnames[DIM_Y]),
MIstart, ci_ptr->dim_steps[DIM_Y] * ci_ptr->dim_lengths[DIM_Y]);
miattputstr(ci_ptr->mnc_fd, ncvarid(ci_ptr->mnc_fd, _dimnames[DIM_Z]),
MIunits, "mm");
miattputstr(ci_ptr->mnc_fd, ncvarid(ci_ptr->mnc_fd, _dimnames[DIM_X]),
MIunits, "mm");
miattputstr(ci_ptr->mnc_fd, ncvarid(ci_ptr->mnc_fd, _dimnames[DIM_Y]),
MIunits, "mm");
miattputstr(ci_ptr->mnc_fd, ncvarid(ci_ptr->mnc_fd, _dimnames[DIM_T]),
MIunits, "s");
ncendef(ci_ptr->mnc_fd);
}
/* ----------------------------- MNI Header -----------------------------------@NAME : CreateImageVars
@INPUT : CDF - ID of the MINC file in which to create MIimagemax
and MIimagemin variables
NumDim - *total* number of image dimensions (2,3, or 4)
DimIDs - ID's of the NumDim image dimensions
NCType - type of the image variable
Signed - TRUE or FALSE, for the image variable
ValidRange - for the image variable
@OUTPUT :
@RETURNS : TRUE on success
FALSE if any error creating either variable
(sets ErrMsg on error)
@DESCRIPTION: Create the MIimagemax and MIimagemin variables in a newly
created MINC file (must be in definition mode!). The
variables will depend on the two lowest (slowest-varying)
image dimensions, ie. frames and slices in the full 4-D
case. If the file has no frames or no slices (or both),
that will be handled properly.
@METHOD :
@GLOBALS : ErrMsg
@CALLS : MINC library
@CREATED : 93-10-28, Greg Ward: code moved from main()
@MODIFIED : 93-11-10, GPW: renamed and modified from CreateMinMax
---------------------------------------------------------------------------- */
Boolean CreateImageVars (int CDF, int NumDim, int DimIDs[],
nc_type NCType, Boolean Signed, double ValidRange[])
{
int image_id;
int max_id, min_id; /* ID's of the newly-created variables */
#ifdef DEBUG
printf ("CreateImageVars:\n");
printf (" Creating MIimage variable with %d dimensions\n", NumDim);
#endif
image_id = micreate_std_variable (CDF, MIimage, NCType, NumDim, DimIDs);
(void) miattputstr (CDF, image_id, MIsigntype, MI_SIGN_STR(Signed));
(void) miattputstr (CDF, image_id, MIcomplete, MI_FALSE);
(void) ncattput (CDF, image_id, MIvalid_range, NC_DOUBLE, 2, ValidRange);
/*
* Create the image-max and image-min variables. They should be
* dependent on the "non-image" dimensions (ie. time and slices,
* if they exist), so pass NumDim-2 as the number of
* dimensions, and DimIDs as the list of dimension ID's --
* micreate_std_variable should then only look at the first one
* or two dimension IDs in the list.
*/
#ifdef DEBUG
printf (" creating MIimagemin and MIimagemax with %d dimensions\n",
NumDim-2);
#endif
max_id = micreate_std_variable (CDF, MIimagemax, NC_DOUBLE,
NumDim-2, DimIDs);
min_id = micreate_std_variable (CDF, MIimagemin, NC_DOUBLE,
NumDim-2, DimIDs);
if ((max_id == MI_ERROR) || (min_id == MI_ERROR))
{
sprintf (ErrMsg, "Error creating image max/min variables: %s\n",
NCErrMsg (ncerr, errno));
return (FALSE);
}
return (TRUE);
} /* CreateImageVars () */
int main()
{
char *file = "test.mnc";
int cdfid;
int dim[MAX_VAR_DIMS];
cdfid=nccreate(file, NC_CLOBBER);
dim[0]=ncdimdef(cdfid, MIzspace, 3L);
dim[1]=ncdimdef(cdfid, MIyspace, 5L);
dim[2]=ncdimdef(cdfid, MIxspace, 6L);
(void) micreate_std_variable(cdfid, MIimage, NC_SHORT, 3, dim);
(void) micreate_std_variable(cdfid, MIimagemax, NC_DOUBLE, 1, dim);
(void) micreate_std_variable(cdfid, MIimagemin, NC_DOUBLE, 1, dim);
(void) micreate_std_variable(cdfid, MIzspace, NC_DOUBLE, 1, NULL);
(void) micreate_std_variable(cdfid, MIzspace_width, NC_DOUBLE, 1, &dim[0]);
(void) micreate_group_variable(cdfid, MIpatient);
(void) micreate_group_variable(cdfid, MIstudy);
(void) micreate_group_variable(cdfid, MIacquisition);
(void) ncclose(cdfid);
return 0;
}
static void
create_dimension(struct conversion_info *ci_ptr, int index, int length)
{
ci_ptr->dim_lengths[index] = length;
if (index == DIM_W && length <= 1) {
return;
}
ci_ptr->dim_ids[index] = ncdimdef(ci_ptr->mnc_fd,
_dimnames[index],
length);
if (index == DIM_T) {
micreate_std_variable(ci_ptr->mnc_fd, _dimnames[index],
NC_DOUBLE, 1, &ci_ptr->dim_ids[index]);
micreate_std_variable(ci_ptr->mnc_fd, MItime_width,
NC_DOUBLE, 1, &ci_ptr->dim_ids[index]);
}
else if (index != DIM_W) {
micreate_std_variable(ci_ptr->mnc_fd, _dimnames[index],
NC_DOUBLE, 0, NULL);
}
}
int main(int argc, char **argv)
{
int icv, cdfid, img, max, min, dimvar;
static int dim[MAX_VAR_DIMS];
static struct { long len; char *name;} diminfo[] = {
{ 7, MIzspace },
{ 9, MIyspace },
{ 2, MIxspace }
};
/* static struct { long len; char *name;} diminfo[]=
{3, MIzspace, 4, MIyspace, 5, MIxspace}; */
static int numdims=sizeof(diminfo)/sizeof(diminfo[0]);
static long coord[]={0,0,0};
static long count[]={3,4,5};
double dvalue;
short int ivalue[]={
111, 113, 115, 117, 119,
121, 123, 125, 127, 129,
131, 133, 135, 137, 139,
141, 143, 145, 147, 149,
211, 213, 215, 217, 219,
221, 223, 225, 227, 229,
231, 233, 235, 237, 239,
241, 243, 245, 247, 249,
311, 313, 315, 317, 319,
321, 323, 325, 327, 329,
331, 333, 335, 337, 339,
341, 343, 345, 347, 349
};
int i, j, k;
int cflag = 0;
#if MINC2
if (argc == 2 && !strcmp(argv[1], "-2")) {
cflag = MI2_CREATE_V2;
}
#endif /* MINC2 */
icv=miicv_create();
miicv_setint(icv, MI_ICV_XDIM_DIR, MI_ICV_NEGATIVE);
miicv_setint(icv, MI_ICV_YDIM_DIR, MI_ICV_NEGATIVE);
miicv_setint(icv, MI_ICV_ZDIM_DIR, MI_ICV_NEGATIVE);
miicv_setint(icv, MI_ICV_NUM_IMGDIMS, 3);
miicv_setint(icv, MI_ICV_DIM_SIZE+0, 5);
miicv_setint(icv, MI_ICV_DIM_SIZE+1, 4);
miicv_setint(icv, MI_ICV_DIM_SIZE+2, 3);
miicv_setint(icv, MI_ICV_DO_DIM_CONV, TRUE);
miicv_setint(icv, MI_ICV_KEEP_ASPECT, FALSE);
miicv_setint(icv, MI_ICV_DO_NORM, TRUE);
cdfid=micreate("test.mnc", NC_CLOBBER | cflag);
for (i=0; i<numdims; i++) {
dim[i]=ncdimdef(cdfid, diminfo[i].name, diminfo[i].len);
dimvar=micreate_std_variable(cdfid, diminfo[i].name, NC_DOUBLE,
0, &dim[i]);
miattputdbl(cdfid, dimvar, MIstep, 0.8);
miattputdbl(cdfid, dimvar, MIstart, 22.0);
}
img=micreate_std_variable(cdfid, MIimage, NC_SHORT,
numdims, dim);
max=micreate_std_variable(cdfid, MIimagemax, NC_DOUBLE, 1, dim);
min=micreate_std_variable(cdfid, MIimagemin, NC_DOUBLE, 1, dim);
ncendef(cdfid);
for (i=0; i<diminfo[0].len; i++) {
dvalue = 200;
coord[0]=i;
ncvarput1(cdfid, max, coord, &dvalue);
dvalue = -dvalue;
ncvarput1(cdfid, min, coord, &dvalue);
}
coord[0]=0;
miicv_attach(icv, cdfid, img);
miicv_inqdbl(icv, MI_ICV_ADIM_START, &dvalue);
printf("adim start = %g", dvalue);
miicv_inqdbl(icv, MI_ICV_ADIM_STEP, &dvalue);
printf(", step = %g\n", dvalue);
miicv_inqdbl(icv, MI_ICV_BDIM_START, &dvalue);
printf("bdim start = %g", dvalue);
miicv_inqdbl(icv, MI_ICV_BDIM_STEP, &dvalue);
printf(", step = %g\n", dvalue);
miicv_inqdbl(icv, MI_ICV_NORM_MAX, &dvalue);
printf("norm : max = %g", dvalue);
miicv_inqdbl(icv, MI_ICV_NORM_MIN, &dvalue);
printf(", min = %g\n", dvalue);
miicv_put(icv, coord, count, ivalue);
miicv_get(icv, coord, count, ivalue);
for (i=0; i<3; i++) {
for (j=0; j<4; j++) {
for (k=0; k<5; k++) {
printf("%5d",ivalue[i*20+j*5+k]);
}
printf("\n");
}
}
miclose(cdfid);
miicv_free(icv);
return 0;
}
int main(int argc, char **argv)
{
int icv, cdfid, img, max, min;
static char *typenm[]={"short", "double"};
static char *boolnm[] = {"true", "false"};
static nc_type intypes[] = {NC_SHORT, NC_DOUBLE};
static int norms[] = {TRUE, FALSE};
static nc_type outtypes[] = {NC_SHORT, NC_DOUBLE};
static int maxpresent[] = {TRUE, FALSE};
static int valpresent[] = {TRUE, FALSE};
static int dim[MAX_VAR_DIMS];
static struct { long len; char *name;} diminfo[] = {
{ 3, MIzspace },
{ 1, MIyspace },
{ 1, MIxspace }
};
static int numdims=sizeof(diminfo)/sizeof(diminfo[0]);
static long coord[]={0,0,0};
static long count[]={1,1,1};
static double max_values[] = {0.4, 0.6, 0.8};
double dvalue;
short int ivalue;
int i, intype, inorm, outtype, imax, ival;
int cflag = 0;
char filename[256];
#if MINC2
if (argc == 2 && !strcmp(argv[1], "-2")) {
cflag = MI2_CREATE_V2;
}
#endif /* MINC2 */
snprintf(filename, sizeof(filename), "test_icv_range-%d.mnc", getpid());
for (intype=0; intype<MAX_IN_TYPES; intype++) {
for (inorm=0; inorm<MAX_NORM; inorm++) {
icv=miicv_create();
miicv_setint(icv, MI_ICV_TYPE, intypes[intype]);
miicv_setint(icv, MI_ICV_DO_NORM, norms[inorm]);
miicv_setdbl(icv, MI_ICV_VALID_MAX, 20000.0);
miicv_setdbl(icv, MI_ICV_VALID_MIN, 0.0);
for (outtype=0; outtype<MAX_OUT_TYPES; outtype++) {
for (imax=0; imax<MAX_MAX; imax++) {
for (ival=0; ival<MAX_VAL; ival++) {
printf(
"in : %s, out : %s, norm : %s, imgmax : %s, valid : %s\n",
typenm[intype], typenm[outtype], boolnm[inorm],
boolnm[imax], boolnm[ival]);
cdfid=micreate(filename, NC_CLOBBER | cflag);
for (i=0; i<numdims; i++)
dim[i]=ncdimdef(cdfid, diminfo[i].name, diminfo[i].len);
img=micreate_std_variable(cdfid, MIimage, outtypes[outtype],
numdims, dim);
if (maxpresent[imax]) {
max=micreate_std_variable(cdfid, MIimagemax, NC_DOUBLE,
1, dim);
min=micreate_std_variable(cdfid, MIimagemin, NC_DOUBLE,
1, dim);
}
if (valpresent[ival]) {
dvalue = 32000;
ncattput(cdfid, img, MIvalid_max, NC_DOUBLE, 1, &dvalue);
dvalue = 0;
ncattput(cdfid, img, MIvalid_min, NC_DOUBLE, 1, &dvalue);
}
ncendef(cdfid);
if (maxpresent[imax]) {
for (i=0; i<3; i++) {
dvalue = max_values[i];
coord[0]=i;
ncvarput1(cdfid, max, coord, &dvalue);
dvalue = -dvalue;
ncvarput1(cdfid, min, coord, &dvalue);
}
coord[0]=0;
}
miicv_attach(icv, cdfid, img);
if (intypes[intype]==NC_DOUBLE) {
dvalue = 0.2;
miicv_put(icv, coord, count, &dvalue);
}
else {
ivalue = 12500;
miicv_put(icv, coord, count, &ivalue);
}
dvalue = 0;
mivarget1(cdfid, img, coord, NC_DOUBLE, MI_SIGNED, &dvalue);
printf(" file value = %g\n", dvalue);
if (intypes[intype]==NC_DOUBLE) {
miicv_get(icv, coord, count, &dvalue);
}
else {
miicv_get(icv, coord, count, &ivalue);
dvalue=ivalue;
}
printf(" icv value = %g\n", dvalue);
miclose(cdfid);
}
}
}
miicv_free(icv);
}
}
unlink(filename);
return 0;
}
MNCAPI int
minc_save_start(char *path, /* Path to the file */
int filetype, /* Date type as stored in the file */
long ct, /* Total length of time axis, in voxels */
long cz, /* Total length of Z axis, in voxels */
long cy, /* Total length of Y axis, in voxels */
long cx, /* Total length of X axis, in voxels */
double dt, /* Sample width along time axis, in seconds */
double dz, /* Sample width along Z axis, in mm */
double dy, /* Sample width along Y axis, in mm */
double dx, /* Sample width along X axis, in mm */
void *infoptr, /* Opaque file structure information */
const char *history) /* New history information */
{
int fd; /* MINC file descriptor */
int dim_id[MI_S_NDIMS]; /* netCDF dimension ID array */
int var_ndims; /* Number of dimensions per variable */
int var_dims[MI_S_NDIMS]; /* Dimension ID's per variable */
int i, j; /* Generic loop counters */
int old_ncopts; /* For supressing fatal error messages */
struct file_info *p_file; /* For accessing the file structure */
struct var_info *p_var;
struct att_info *p_att;
int var_id; /* netCDF ID for variable */
char *signstr;
nc_type nctype;
old_ncopts =get_ncopts();
set_ncopts(0);
fd = micreate(path, NC_CLOBBER);
set_ncopts(old_ncopts);
if (fd < 0) {
return (MINC_STATUS_ERROR);
}
if (ct > 0) {
dim_id[MI_S_T] = ncdimdef(fd, MItime, ct);
micreate_std_variable(fd, MItime, NC_INT, 0, NULL);
if (dt > 0.0) {
miattputdbl(fd, ncvarid(fd, MItime), MIstep, dt);
}
}
else {
dim_id[MI_S_T] = -1;
}
if (cz > 0) {
dim_id[MI_S_Z] = ncdimdef(fd, MIzspace, cz);
micreate_std_variable(fd, MIzspace, NC_INT, 0, NULL);
if (dz > 0.0) {
miattputdbl(fd, ncvarid(fd, MIzspace), MIstep, dz);
}
}
else {
dim_id[MI_S_Z] = -1;
}
if (cy > 0) {
dim_id[MI_S_Y] = ncdimdef(fd, MIyspace, cy);
micreate_std_variable(fd, MIyspace, NC_INT, 0, NULL);
if (dy > 0.0) {
miattputdbl(fd, ncvarid(fd, MIyspace), MIstep, dy);
}
}
else {
return (MINC_STATUS_ERROR); /* Must define Y */
}
if (cx > 0) {
dim_id[MI_S_X] = ncdimdef(fd, MIxspace, cx);
micreate_std_variable(fd, MIxspace, NC_INT, 0, NULL);
if (dx > 0.0) {
miattputdbl(fd, ncvarid(fd, MIxspace), MIstep, dx);
}
}
else {
return (MINC_STATUS_ERROR); /* Must define X */
}
/* The var_dims[] array is the array of actual dimension ID's to
* be used in defining the image variables. Here I set it up by
* copying all valid dimension ID's from the dim_id[] array.
*/
var_ndims = 0;
for (i = 0; i < MI_S_NDIMS; i++) {
if (dim_id[i] >= 0) {
var_dims[var_ndims] = dim_id[i];
var_ndims++;
}
}
minc_simple_to_nc_type(filetype, &nctype, &signstr);
/* Create the image variable with the standard
* dimension order, and the same type as the template
* file.
*/
micreate_std_variable(fd, MIimage, nctype, var_ndims, var_dims);
micreate_std_variable(fd, MIimagemin, NC_DOUBLE, 1, var_dims);
micreate_std_variable(fd, MIimagemax, NC_DOUBLE, 1, var_dims);
/* Copy information from the infoptr to the output.
*/
if ((p_file = infoptr) != NULL) {
old_ncopts =get_ncopts();
set_ncopts(0);
for (i = 0; i < p_file->file_natts; i++) {
p_att = &p_file->file_atts[i];
if (strcmp(p_att->att_name, "ident") != 0) {
ncattput(fd, NC_GLOBAL, p_att->att_name, p_att->att_type,
p_att->att_len, p_att->att_val);
}
}
for (i = 0; i < p_file->file_nvars; i++) {
p_var = &p_file->file_vars[i];
if ((var_id = ncvarid(fd, p_var->var_name)) < 0) {
var_id = ncvardef(fd, p_var->var_name, p_var->var_type,
p_var->var_ndims, p_var->var_dims);
}
for (j = 0; j < p_var->var_natts; j++) {
p_att = &p_var->var_atts[j];
ncattput(fd, var_id, p_att->att_name, p_att->att_type,
p_att->att_len, p_att->att_val);
}
}
set_ncopts(old_ncopts);
}
miattputstr(fd, ncvarid(fd, MIimage), MIcomplete, MI_FALSE);
miattputstr(fd, ncvarid(fd, MIimage), MIsigntype, signstr);
miappend_history(fd, history);
ncendef(fd);
return fd;
}
/* ----------------------------- MNI Header -----------------------------------
@NAME : setup_variables
@INPUT : inmincid - id of input minc file (MI_ERROR if no file)
mincid - id of output minc file
volume_info - volume information
arg_string - string giving argument list
@OUTPUT : (nothing)
@RETURNS : (nothing)
@DESCRIPTION: Routine to set up variables in the output minc file
@METHOD :
@GLOBALS :
@CALLS :
@CREATED : August 26, 1993 (Peter Neelin)
@MODIFIED :
---------------------------------------------------------------------------- */
static void setup_variables(int inmincid, int mincid,
Volume_Info *volume_info,
char *arg_string)
{
int dim[MAX_VAR_DIMS], ndims, idim, varid;
int excluded_vars[10], nexcluded;
/* Create the dimensions */
ndims = NUMBER_OF_DIMENSIONS;
dim[0] = ncdimdef(mincid, volume_info->dimension_names[0],
volume_info->nslices);
dim[1] = ncdimdef(mincid, volume_info->dimension_names[1],
volume_info->nrows);
dim[2] = ncdimdef(mincid, volume_info->dimension_names[2],
volume_info->ncolumns);
/* If an input file is provided, copy all header info from that file except
image, image-max, image-min */
if (inmincid != MI_ERROR) {
/* Look for the image variable and the image-max/min variables so that
we can exclude them from the copy */
nexcluded = 0;
excluded_vars[nexcluded] = ncvarid(inmincid, MIimage);
if (excluded_vars[nexcluded] != MI_ERROR) nexcluded++;
excluded_vars[nexcluded] = ncvarid(inmincid, MIimagemax);
if (excluded_vars[nexcluded] != MI_ERROR) nexcluded++;
excluded_vars[nexcluded] = ncvarid(inmincid, MIimagemin);
if (excluded_vars[nexcluded] != MI_ERROR) nexcluded++;
/* Copy the variable definitions */
(void) micopy_all_var_defs(inmincid, mincid, nexcluded, excluded_vars);
}
/* Set up the dimension variables. If the variable doesn't exist, create
it (either no input file or variable did not exist in it). If the
dimensions are not standard, then no variable is created. */
for (idim=0; idim < ndims; idim++) {
ncopts = 0;
varid = ncvarid(mincid, volume_info->dimension_names[idim]);
if (varid == MI_ERROR) {
varid = micreate_std_variable(mincid,
volume_info->dimension_names[idim],
NC_INT, 0, NULL);
}
ncopts = NC_OPTS_VAL;
if (varid != MI_ERROR) {
(void) miattputdbl(mincid, varid, MIstep,
volume_info->step[idim]);
(void) miattputdbl(mincid, varid, MIstart,
volume_info->start[idim]);
}
}
/* Create the image, image-max and image-min variables */
setup_image_variables(inmincid, mincid, ndims, dim);
/* Add the time stamp to the history */
update_history(mincid, arg_string);
/* Put the file in data mode */
(void) ncendef(mincid);
/* Copy over variable values */
if (inmincid != MI_ERROR) {
(void) micopy_all_var_values(inmincid, mincid,
nexcluded, excluded_vars);
}
}
/* Test case 1 - file creation & definition.
*/
static int test1(struct testinfo *ip, struct dimdef *dims, int ndims)
{
int varid;
int stat;
int i;
/* Test case #1 - file creation
*/
ip->name = micreate_tempfile();
if (ip->name == NULL) {
FUNC_ERROR("micreate_tempfile\n");
}
ip->fd = micreate(ip->name, NC_CLOBBER);
if (ip->fd < 0) {
FUNC_ERROR("micreate");
}
/* Have to use ncdimdef() here since there is no MINC equivalent. Sigh.
*/
for (i = 0; i < ndims; i++) {
/* Define the dimension
*/
ip->dim[i] = ncdimdef(ip->fd, dims[i].name, dims[i].length);
if (ip->dim[i] < 0) {
FUNC_ERROR("ncdimdef");
}
/* Create the dimension variable.
*/
varid = micreate_std_variable(ip->fd, dims[i].name, NC_DOUBLE, 0,
&ip->dim[i]);
if (varid < 0) {
FUNC_ERROR("micreate_std_variable");
}
stat = miattputdbl(ip->fd, varid, MIstep, 0.8);
if (stat < 0) {
FUNC_ERROR("miattputdbl");
}
stat = miattputdbl(ip->fd, varid, MIstart, 22.0);
if (stat < 0) {
FUNC_ERROR("miattputdbl");
}
}
/* Create the image-max variable.
*/
ip->maxid = micreate_std_variable(ip->fd, MIimagemax, NC_FLOAT, 0, NULL);
if (ip->maxid < 0) {
FUNC_ERROR("micreate_std_variable");
}
/* Create the image-min variable.
*/
ip->minid = micreate_std_variable(ip->fd, MIimagemin, NC_FLOAT, 0, NULL);
if (ip->minid < 0) {
FUNC_ERROR("micreate_std_variable");
}
ip->imgid = micreate_std_variable(ip->fd, MIimage, NC_INT, ndims, ip->dim);
if (ip->imgid < 0) {
FUNC_ERROR("micreate_std_variable");
}
ip->test_group = ncvardef(ip->fd,(char*)"test",NC_INT,0,0);/* micreate_group_variable(ip->fd,(char*)"test");*/
if(ip->test_group<0)
{
FUNC_ERROR("micreate_group_variable");
}
ip->large_attribute=calloc(ip->attribute_size,sizeof(char));
memset(ip->large_attribute,'X',ip->attribute_size-1);
ip->test_attribute = ncattput(ip->fd, ip->test_group, "test", NC_CHAR, ip->attribute_size, ip->large_attribute);
return (0);
}
/* Test case 1 - file creation & definition.
*/
int
test1(struct testinfo *ip, struct dimdef *dims, int ndims)
{
int fd2;
int varid;
int stat;
int i;
/* Test case #1 - file creation
*/
ip->name = micreate_tempfile();
if (ip->name == NULL) {
FUNC_ERROR("micreate_tempfile\n");
}
ip->fd = micreate(ip->name, NC_CLOBBER);
if (ip->fd < 0) {
FUNC_ERROR("micreate");
}
/* Try to create another file of the same name - should fail.
*/
fd2 = micreate(ip->name, NC_NOCLOBBER);
if (fd2 >= 0) {
FUNC_ERROR("micreate");
}
/* Try to open the file for write - should fail.
*/
/* VF: it doesn't fail!
fd2 = miopen(ip->name, NC_WRITE);
if (fd2 >= 0) {
FUNC_ERROR("miopen");
} */
/* Have to use ncdimdef() here since there is no MINC equivalent. Sigh.
*/
for (i = 0; i < ndims; i++) {
/* Define the dimension
*/
ip->dim[i] = ncdimdef(ip->fd, dims[i].name, dims[i].length);
if (ip->dim[i] < 0) {
FUNC_ERROR("ncdimdef");
}
/* Create the dimension variable.
*/
varid = micreate_std_variable(ip->fd, dims[i].name, NC_DOUBLE, 0,
&ip->dim[i]);
if (varid < 0) {
FUNC_ERROR("micreate_std_variable");
}
stat = miattputdbl(ip->fd, varid, MIstep, 0.8);
if (stat < 0) {
FUNC_ERROR("miattputdbl");
}
stat = miattputdbl(ip->fd, varid, MIstart, 22.0);
if (stat < 0) {
FUNC_ERROR("miattputdbl");
}
}
/* Try to create a bogus variable. This should trigger an error.
*/
varid = micreate_std_variable(ip->fd, "xyzzy", NC_DOUBLE, 0, NULL);
if (varid >= 0) {
FUNC_ERROR("micreate_std_variable");
}
/* Create the image-max variable.
*/
ip->maxid = micreate_std_variable(ip->fd, MIimagemax, NC_FLOAT, 0, NULL);
if (ip->maxid < 0) {
FUNC_ERROR("micreate_std_variable");
}
/* Create the image-min variable.
*/
ip->minid = micreate_std_variable(ip->fd, MIimagemin, NC_FLOAT, 0, NULL);
if (ip->minid < 0) {
FUNC_ERROR("micreate_std_variable");
}
ip->imgid = micreate_std_variable(ip->fd, MIimage, NC_INT, ndims, ip->dim);
if (ip->imgid < 0) {
FUNC_ERROR("micreate_std_variable");
}
return (0);
}
