[svn-r5088] Purpose:

	Code removal
Description:
	Removed the HDF4 source files from the HDF5 tree. The directories
	will remain. Use the "-P" option when doing a cvs checkout or update
	to "prune" the empty directories from your personal tree.

1 files changed, 0 insertions, 1069 deletions

diff --git a/tools/h4toh5/h4toh5image.c b/tools/h4toh5/h4toh5image.c
deleted file mode 100644
index 7523974..0000000
--- a/tools/h4toh5/h4toh5image.c
+++ /dev/null
@@ -1,1069 +0,0 @@
-/*-------------------------------------------------------------------------
- *
- * Copyright (C) 2000   National Center for Supercomputing Applications.
- *                      All rights reserved.
- *
- *-------------------------------------------------------------------------
- */
-
-/******************************************************************************
-
-  Description: 
-
-1. converter
-
-See HDF4 to HDF5 mapping specification at
-(http://hdf.ncsa.uiuc.edu/HDF5/papers/h4toh5) for the default mapping 
-from HDF4 object to HDF5 object.
- 
-The whole converter includes 10 files, h4toh5util.h, h4toh5main.h, h4toh5util.c, h4toh5main.c, h4toh5sds.c, h4toh5image.c,h4toh5vdata.c,h4toh5vgroup.c,h4toh5pal.c and h4toh5anno.c.
-
-2. this file 
-
-converting an hdf4 image object into an hdf5 dataset, for three component image, this object will be converted into an hdf5 dataset with compound data type.
-
-Author:  Kent Yang(ymuqun@ncsa.uiuc.edu)
- 
-
-*****************************************************************************/
-
-#include "h4toh5main.h"
-
-/*-------------------------------------------------------------------------
- * Function:	Image_h4_to_h5
- *
- * Purpose:     translate Image object into hdf5 dataset
- *              
- * Return:	FAIL if failed, SUCCEED if successful.
- *
- * In :	        
-                ri_id: RI identifier
-		h5_group: hdf5 group id
-		h5_palgroup: hdf5 palette group id
-
- *-------------------------------------------------------------------------
- */	
-
-int Image_h4_to_h5(int32 file_id,int32 ri_id,hid_t h5_group,hid_t h5_palgroup,int h4_attr) {
-
-  int32    istat;
-  int32    ngrattrs;
-  int32    ncomp;
-  int      check_gloattr;
-  int32    start[2];
-  int32    edges[2];
-  int32    dimsizes[2];
-  uint16   gr_ref;
-  int32    image_dtype;
-
-  int      check_imagename;
-  int      i;
-  char     image_name[MAX_GR_NAME];
-  char     grlabel[MAX_GR_NAME];
-  char     image_class[MAX_GR_NAME];
-  char*    h5cimage_name;
-  void*    image_data;
-  HDF_CHUNK_DEF c_def_out;
-  hsize_t    chunk_dims[2];
-  hsize_t    chunk_dims24[3];
-  int32    c_flags;
-  int32    interlace_mode;
-
-  /* for checking compression */
-
-  sp_info_block_t info_block;
-  int16           special_code;
-  int32           access_id;
-  uint16          ri_ref;
-  int             gzip_level;
-  uint16          temp_tag;
-  uint16*          ptag_out;
-
-  /* define varibles for hdf5. */
-
-  hid_t    h5ty_id;
-  hid_t    h5memtype;
-
-  hid_t    h5_ctype;
-  hid_t    h5_cmemtype;
-
-  hid_t    h5d_sid;
-  hid_t    h5dset;
-
-  size_t   h4size;
-  size_t   h4memsize;
-  hsize_t   fielddim[1];
-  hsize_t  h5dims[2];
-  hsize_t  h5dims24[3];
-#ifdef H5_WANT_H5_V1_4_COMPAT
-  hsize_t bufsize;
-#else /* H5_WANT_H5_V1_4_COMPAT */
-  size_t bufsize;
-#endif /* H5_WANT_H5_V1_4_COMPAT */
-  herr_t   ret;
-  hid_t    create_plist;
-  hid_t    write_plist;
-
-
-  temp_tag = DFTAG_NULL;
-  ptag_out = &temp_tag;
-
-  /* zeroing out memory.*/
-
-  h4toh5_ZeroMemory(image_name,MAX_GR_NAME);
-  h4toh5_ZeroMemory(image_class,MAX_GR_NAME);
-  h4toh5_ZeroMemory(grlabel,MAX_GR_NAME);
-
-  /* Obtain information of the image.*/  
-
-  if(GRgetchunkinfo(ri_id,&c_def_out,&c_flags)==FAIL){
-    printf("error in getting chunking information. \n");
-    return FAIL;
-  }
-
-  istat = GRgetiminfo(ri_id, image_name, &ncomp, &image_dtype, 
-		      &interlace_mode, dimsizes, &ngrattrs);
-
-  if(istat == FAIL) {
-    printf("Cannot obtain GR info. at Image routine.\n");
-    return FAIL;
-  }
-
-  /* data type transferring from hdf4 to hdf5. */
-  if(h4type_to_h5type(image_dtype,&h5memtype,&h4memsize,
-		      &h4size,&h5ty_id)== FAIL) {
-    printf("failed to translate image datatype. \n");
-    return FAIL;
-  }
- 
-  /* check whether the datatype is string. */
-  if (h5ty_id == H5T_STRING) {
-    /* rechange string datatype into numerical datatype.*/
-
-    if(h5string_to_int(image_dtype,&h5memtype,h4memsize,
-		       &h5ty_id)== FAIL) {
-      printf("error in translating H5T_STRING to int.\n");
-      return FAIL;
-    }
-  }     
-
-  start[0]   = 0;
-  start[1]   = 0;
-  edges[0]   = dimsizes[0];
-  edges[1]   = dimsizes[1];
-  
-  image_data = malloc(h4memsize*dimsizes[0]*dimsizes[1]*ncomp);
-  
-  if(image_data == NULL) {
-    printf("error in allocating memory for image data. \n");
-    return FAIL;
-  }
-  /* GRreqimageil or GRreadimage are not working, comment this out for the time being.
-  if(interlace_mode == MFGR_INTERLACE_LINE){
-    istat = GRreqimageil(ri_id,MFGR_INTERLACE_PIXEL);
-    if(istat == FAIL){
-      printf("error in setting interlace_mode.\n");
-      free(image_data);
-      return FAIL;
-    }
-  }
-  */
-  istat = GRreadimage(ri_id, start, NULL, edges, (VOIDP)image_data);
-
-  if (istat == FAIL) {
-    printf("error in reading images.\n");
-    free(image_data);
-    return FAIL;
-  }
-
-  /* change the order of image dimension:
-     due to the difference of hdf4 image specification and
-     hdf5 image specification. Here we should separate 8-bit from
-     24-bit according to H4TOH5 mapping specification.*/
-  
-  if(ncomp == 1) {
-    h5dims[0] = edges[1]-start[1];
-    h5dims[1] = edges[0]-start[0];
-  }
-  
-  else {
-    if(interlace_mode == MFGR_INTERLACE_PIXEL){
-      h5dims24[0] = edges[1]-start[1];
-      h5dims24[1] = edges[0]-start[0];
-      h5dims24[2] = 3;
-    }
-    /* currently scan-line is not supported. */
-    else if (interlace_mode == MFGR_INTERLACE_LINE){
-      printf("currently line interleaving is not supported.\n");
-      printf("the image %s will not be converted.\n",image_name);
-       free(image_data);
-       return SUCCEED;     
-      /*h5dims24[0] = 3;
-      h5dims24[1] = edges[1]-start[1];
-      h5dims24[2] = edges[0]-start[0];*/
-    }
-    else if (interlace_mode == MFGR_INTERLACE_COMPONENT){
-      h5dims24[0] = 3;
-      h5dims24[1] = edges[1]-start[1];
-      h5dims24[2] = edges[0]-start[0];
-    }
-
-    else {/* treat as pixel */
-      h5dims24[0] = edges[1]-start[1];
-      h5dims24[1] = edges[0]-start[0];
-      h5dims24[2] = 3;
-    }
-  }
-
-  gr_ref = GRidtoref(ri_id);
-  if(gr_ref == 0) {
-    printf("error in obtaining gr reference number. \n");
-    free(image_data);
-    return FAIL;
-  }
-  
-  /* obtaining absolute path of image name.*/
-  check_imagename = -10;
-  h5cimage_name = get_name(gr_ref,2*num_images,gr_hashtab,&check_imagename);
-
-  if (h5cimage_name == NULL && check_imagename == 0 ) {
-     printf("error,cannot find image name.\n");
-     free(image_data);
-     return FAIL;
-  }
-
-  if (h5cimage_name == NULL && check_imagename == -1) {
-     printf("error,image name is not defined.\n");
-     free(image_data);
-     return FAIL;
-  }
-
-  if (h5cimage_name == NULL && check_imagename == -2) {
-     printf("error,not enough memory for get_name. \n");
-     free(image_data);
-     return FAIL;
-  }
-
-  /****  check number of component of the image object,
-	 and transfer HDF4 object into HDF5 object. ****/
-
-  if (ncomp <= 0) {
-     printf("error in obtaining image component\n");
-     free(image_data);
-     free(h5cimage_name);
-     return FAIL;
-  }
-
-  /* create property list. */
-
-  create_plist = H5Pcreate(H5P_DATASET_CREATE);
-  /* wait until the compression information can be obtained for image,
-     4/28/2001, Kent Yang.*/
-
-  /* the following code deals with compression. Has to use 
-     some middle-level APIs. */
-  ri_ref = 0;
-  /*
-   ri_ref = get_RIref(file_id,DFTAG_VG,gr_ref,ptag_out);
-
-  if(ri_ref >0 ){
-    if(*ptag_out == DFTAG_RI) printf("okay\n");
-    access_id = Hstartread(file_id,*ptag_out,ri_ref);
-  }
-  */
-  if(ri_ref == 0) 
-    access_id = FAIL;
-  
-
-  /*  if(access_id == FAIL){
-    printf("cannot find RI tag,the compression mode will not be checked");
-    printf(" when RIG is not chunked.\n");
-    }*/
-  
-  if(access_id != FAIL) {
-    istat = Hinquire(access_id,NULL,NULL,NULL,NULL,NULL,NULL,NULL,&special_code);
-    if(istat == FAIL) {
-      printf("failed to inquire information \n ");
-      free(image_data);
-       free(h5cimage_name);
-      H5Pclose(create_plist);
-      return FAIL;	
-    }
-
-    if(special_code >0){
-      
-      if(HDget_special_info(access_id,&info_block)==FAIL){
-	printf("fail to get special info.\n");
-	free(image_data);
-	free(h5cimage_name);
-	H5Pclose(create_plist);
-	return FAIL;	
-      }
-
-      if(info_block.key == SPECIAL_COMP) {
-   
-	if(c_flags == HDF_NONE){
-	  /* 1. the current HDF5 will not handle compression case itself, 
-	     in order that the converted HDF5 is compressed, we have to
-	     provide a chunking size. currently it is set to h5dim[i].*/
-
-	  if(ncomp == 1) {
-	    chunk_dims[0] = (hsize_t)(h5dims[0]);
-            chunk_dims[1] =(hsize_t)(h5dims[1]);
-	  
-	    if(H5Pset_chunk(create_plist, 2, chunk_dims)<0) {
-	      printf("failed to set up chunking information for ");
-	      printf("property list.\n");
-	      free(image_data);
-	      free(h5cimage_name);
-	      H5Pclose(create_plist);
-	      return FAIL;	
-	    }
-	  }
-
-	  else if(ncomp ==3) {/*24-bit chunk dimension is set to the current
-				HDF5 dimension.*/
-	    if(interlace_mode == MFGR_INTERLACE_PIXEL){
-	      chunk_dims24[0] = edges[1]-start[1]; 
-	      chunk_dims24[1] = edges[0]-start[0];
-	      chunk_dims24[2] = 3;
-	    }
-	    /* currently scan-line is not supported.
-	       else if (interlace_mode == MFGR_INTERLACE_LINE){
-	       chunk_dims24[0] = 3; 
-	       chunk_dims24[1] = edges[1]-start[1];
-	       chunk_dims24[2] = edges[0]-start[0];
-	       }
-	    */
-	    else if (interlace_mode == MFGR_INTERLACE_COMPONENT){
-	      chunk_dims24[1] = edges[1]-start[1]; 
-	      chunk_dims24[2] = edges[0]-start[0];
-	      chunk_dims24[0] = 3;
-	    }
-
-	    else {/* treat as pixel */
-	      chunk_dims24[0] = edges[1]-start[1]; 
-	      chunk_dims24[1] = edges[0]-start[0];
-	      chunk_dims24[2] = 3;
-	    }
-	    if(H5Pset_chunk(create_plist, 3, chunk_dims24)<0) {
-	      printf("failed to set up chunking information for ");
-	      printf("property list.\n");
-	      free(image_data);
-	      free(h5cimage_name);
-	      H5Pclose(create_plist);
-	      return FAIL;	
-	    }
-	  }
-	  if(H5Pset_deflate(create_plist,GZIP_COMLEVEL)<0){
-	    /*    if(H5Pset_deflate(create_plist,2)<0){*/
-	    printf("fail to set compression method for HDF5 file.\n"); 
-	     free(image_data);
-	    free(h5cimage_name);
-	    H5Pclose(create_plist);
-	    return FAIL;
-	  }
-	}
-
-      }
-    }
-
-  }
-  if(c_flags == HDF_CHUNK || c_flags == (HDF_CHUNK | HDF_COMP)
-     || c_flags == (HDF_CHUNK | HDF_NBIT)  ){
-
-    if(c_def_out.comp.comp_type == COMP_CODE_RLE || c_def_out.comp.comp_type == COMP_CODE_NBIT || c_def_out.comp.comp_type == COMP_CODE_SKPHUFF || c_def_out.comp.comp_type == COMP_CODE_DEFLATE || c_def_out.comp.comp_type == COMP_CODE_JPEG) {
-      if(ncomp ==1) {   
-	chunk_dims[0] = c_def_out.chunk_lengths[0]; 
-	chunk_dims[1] = c_def_out.chunk_lengths[1];
-       
-	if(H5Pset_chunk(create_plist, 2, (hsize_t *)chunk_dims)<0) {
-	  printf("failed to set up chunking information for ");
-	  printf("property list.\n");
-	  free(image_data);
-	  free(h5cimage_name);
-	  H5Pclose(create_plist);
-	  return FAIL;	
-	}
-      }
-
-      else if(ncomp ==3) {
-	if(interlace_mode == MFGR_INTERLACE_PIXEL){
-      chunk_dims24[0] = c_def_out.chunk_lengths[1]; 
-      chunk_dims24[1] = c_def_out.chunk_lengths[0];
-      chunk_dims24[2] = 3;
-    }
-    /* currently scan-line is not supported.
-    else if (interlace_mode == MFGR_INTERLACE_LINE){
-      chunk_dims24[0] = c_def_out.chunk_lengths[1]; 
-      chunk_dims24[2] = c_def_out.chunk_lengths[0];
-      chunk_dims24[1] = 3;
-    }
-    */
-    else if (interlace_mode == MFGR_INTERLACE_COMPONENT){
-      chunk_dims24[1] = c_def_out.chunk_lengths[1]; 
-      chunk_dims24[2] = c_def_out.chunk_lengths[0];
-      chunk_dims24[0] = 3;
-    }
-
-    else {/* treat as pixel */
-      chunk_dims24[0] = c_def_out.chunk_lengths[1]; 
-      chunk_dims24[1] = c_def_out.chunk_lengths[0];
-      chunk_dims24[2] = 3;
-    }
-      }
-	if(c_def_out.comp.comp_type == COMP_CODE_DEFLATE)
-	  gzip_level = c_def_out.comp.cinfo.deflate.level;
-	else gzip_level = GZIP_COMLEVEL;
-	if(H5Pset_deflate(create_plist,gzip_level)<0){
-	  printf("fail to set compression method for HDF5 file.\n"); 
-	   free(image_data);
-	   free(h5cimage_name);
-	   H5Pclose(create_plist);
-	}
-    }
-  }
-
-  /* HDF4 can support various compression methods including simple RLE, NBIT, Skip Huffman, gzip,Jpeg , HDF5 currently only supports gzip compression. 
- By default, we will compress HDF5 dataset by using gzip compression if HDF5 file is compressed. */
-
-  /* we don't use data transfer property list. 
- write_plist = H5Pcreate(H5P_DATASET_XFER);
- bufsize = h4memsize *h5dims[1]*ncomp;
-
-  if(H5Pset_buffer(write_plist,bufsize,NULL,NULL)<0) {
-    printf("fail to create data transfer property list.\n");
-    free(image_data);
-    free(h5cimage_name);
-    H5Pclose(create_plist);
-    return FAIL;		
-  }
-  */
- if (ncomp == 1) {
-
-     h5d_sid = H5Screate_simple(2,h5dims,NULL);
-     
-     if(h5d_sid <0) {
-       printf("error in creating space for dataset. \n");
-       free(image_data);
-       free(h5cimage_name);
-       H5Pclose(create_plist);
-       return FAIL;
-     }
-
-     h5dset  = H5Dcreate(h5_group,h5cimage_name,h5ty_id,h5d_sid,create_plist);
-
-     if(h5dset < 0) {
-       printf("error in creating hdf5 dataset converted from images. \n");
-       free(image_data);
-       free(h5cimage_name);
-       H5Pclose(create_plist);
-       return FAIL;
-     }
-
-     if (H5Dwrite(h5dset,h5memtype,h5d_sid,h5d_sid,H5P_DEFAULT,
-		  image_data)<0) {
-        printf("error writing data for hdf5 dataset converted from images.\n");
-	free(image_data);
-	free(h5cimage_name);
-        H5Pclose(create_plist);
-        return FAIL;
-     }
-	
- }
-
-   else { /* 24-bit image */
-
-     h5d_sid     = H5Screate_simple(3,h5dims24,NULL);
-
-     h5dset      = H5Dcreate(h5_group,h5cimage_name,h5ty_id,h5d_sid,
-			     create_plist);
-     if(h5dset < 0) {
-       printf("error in creating dataset. \n");
-       free(image_data);
-       free(h5cimage_name);
-       H5Pclose(create_plist);
-       return FAIL;
-     }
-
-     if (H5Dwrite(h5dset,h5memtype,h5d_sid,h5d_sid,H5P_DEFAULT,
-		  (void *)image_data)<0) {
-        printf("error writing data\n");
-	free(image_data);
-	free(h5cimage_name);
-        H5Pclose(create_plist);
-	return FAIL;
-     } 
-   }
- 
- free(image_data);
-/* convert image annotation into attribute of image dataset.
-     Since there is no routines to find the exact tag of image object,
-     we will check three possible object tags of image objects, that is:
-     DFTAG_RIG,DFTAG_RI,DFTAG_RI8. If the object tag of image object is 
-     falling out of this scope, we will not convert annotations into
-     hdf5 attributes; it is user's responsibility to make sure object tags
-     for image objects are only one of the above three tags.*/
-
-  if(Annoobj_h4_to_h5(file_id,gr_ref,DFTAG_RIG,h5dset)== FAIL){
-    printf("failed to convert image annotation into hdf5 attribute.\n");
-    free(h5cimage_name);
-    H5Pclose(create_plist);
-    H5Sclose(h5d_sid); 
-    H5Dclose(h5dset);
-    return FAIL;
-  }
-
-  if(Annoobj_h4_to_h5(file_id,gr_ref,DFTAG_RI,h5dset)== FAIL){
-    printf("failed to convert image annotation into hdf5 attribute.\n");
-    free(h5cimage_name);
-    H5Pclose(create_plist);
-    H5Sclose(h5d_sid); 
-    H5Dclose(h5dset);
-    return FAIL;
-  }
-
-  if(Annoobj_h4_to_h5(file_id,gr_ref,DFTAG_RI8,h5dset)== FAIL){
-    printf("failed to convert image annotation into hdf5 attribute.\n");
-    free(h5cimage_name);
-    H5Pclose(create_plist);
-    H5Sclose(h5d_sid); 
-    H5Dclose(h5dset);
-    return FAIL;
-  }
-  
-  
- /************************************/
-  /* translate GR attributes into HDF5 dataset attribute.*/
-
-  check_gloattr = 0;
-  if(gr_tranattrs(ri_id,h5dset,ngrattrs,check_gloattr)==FAIL){ 
-    printf(" cannot obtain attributes. \n");
-    H5Pclose(create_plist);
-    H5Sclose(h5d_sid); 
-    H5Dclose(h5dset);
-    return FAIL;
-  }
-
-  /*  deal with h5dset predefined and user-defined attributes. 
-      Obtain the name and data type and the total number of attributes.  
-      Data attribute at hdf4 is only one-dimensional array. */
-
- 
-  if (ncomp == 1 && h4size == 1)
-     strcpy(grlabel,RAST8LABEL);
-  else if(ncomp == 3 && h4size == 1)
-     strcpy(grlabel,RAST24LABEL);
-  else
-     strcpy(grlabel,GRLABEL);
-
-  strcpy(image_class,IM_CLASS);
-
-  /* transfer hdf4 predefined attributes into hdf5 dataset.*/
-  if(h4_transpredattrs(h5dset,HDF4_OBJECT_TYPE,grlabel)==FAIL){
-    printf("error in getting hdf4 image type attribute \n");
-    H5Pclose(create_plist);
-    H5Sclose(h5d_sid); 
-    H5Dclose(h5dset);
-    free(h5cimage_name);
-    return FAIL;
-  }
-
-  if(h4_attr!=0) {
-  if(h4_transpredattrs(h5dset,HDF4_OBJECT_NAME,image_name)==FAIL){
-    printf("error in getting hdf4 image name attribute. \n");
-    H5Pclose(create_plist);
-    H5Sclose(h5d_sid); 
-    H5Dclose(h5dset);
-    free(h5cimage_name);
-    return FAIL;
-  }
-  }
-  if(h4_transpredattrs(h5dset,HDF4_IMAGE_CLASS,image_class)==FAIL){
-    printf("error in getting hdf4 image class attribute. \n");
-    H5Pclose(create_plist);
-    H5Sclose(h5d_sid); 
-    H5Dclose(h5dset);
-    free(h5cimage_name);
-    return FAIL;
-  }
-
-  if(ncomp >1) {
-    if(interlace_mode == MFGR_INTERLACE_PIXEL){
-      if(h4_transpredattrs(h5dset,INTERLACE_MODE,PIXEL_INTERLACE)==FAIL){
-	printf("unable to generate image pixel attribute.\n");
-      H5Pclose(create_plist);
-      H5Sclose(h5d_sid); 
-      H5Dclose(h5dset);
-      free(h5cimage_name);
-      return FAIL;
-      }
-    }
-    /* currently scan-line is not supported.
-    else if (interlace_mode == MFGR_INTERLACE_LINE){
-     if(h4_transpredattrs(h5dset,INTERLACE_MODE,LINE_INTERLACE)==FAIL){
-      printf("unable to generate image line attribute.\n");
-      H5Pclose(create_plist);
-      H5Sclose(h5d_sid); 
-      H5Dclose(h5dset);
-      free(h5cimage_name);
-      return FAIL;
-      }
-      }*/
-    else if (interlace_mode == MFGR_INTERLACE_COMPONENT){
-      if(h4_transpredattrs(h5dset,INTERLACE_MODE,PLANE_INTERLACE)==FAIL){
-       printf("unable to generate image component attribute.\n");
-      H5Pclose(create_plist);
-      H5Sclose(h5d_sid); 
-      H5Dclose(h5dset);
-      free(h5cimage_name);
-      return FAIL;
-      }
-    }
-   
-    else {/* treat as pixel interlace mode. */
-     if(h4_transpredattrs(h5dset,INTERLACE_MODE,PIXEL_INTERLACE)==FAIL){
-      printf("unable to generate image pixel attribute.\n");
-      H5Pclose(create_plist);
-      H5Sclose(h5d_sid); 
-      H5Dclose(h5dset);
-      free(h5cimage_name);
-      return FAIL;
-      }
-    } 
-      
-  }
-  
-  if(h4_attr !=0){
- 
-  gr_ref = GRidtoref(ri_id);
-
-  if(gr_ref == 0) {
-    printf("error in obtaining reference number of GR.\n");
-    H5Pclose(create_plist);
-    H5Sclose(h5d_sid); 
-    H5Dclose(h5dset);
-    free(h5cimage_name);
-    return FAIL;
-  }
-
-  if(h4_transnumattr(h5dset,HDF4_REF_NUM,gr_ref)==FAIL) {
-    printf("error in getting hdf4 image number attribute.\n");
-    H5Pclose(create_plist);
-    H5Sclose(h5d_sid); 
-    H5Dclose(h5dset);
-    free(h5cimage_name);
-    return FAIL;
-  }
-  }
-  /* deal with palette. */
-
-  if(gr_palette(file_id,ri_id,h5dset,h5_palgroup,h4_attr)== FAIL) {
-    printf("error in translating palette into h5 dataset.\n");
-    H5Pclose(create_plist);
-    H5Sclose(h5d_sid); 
-    H5Dclose(h5dset);
-    free(h5cimage_name);
-    return FAIL;
-  }
- 
-  ret   = H5Pclose(create_plist);
-  ret   = H5Sclose(h5d_sid); 
-  ret   = H5Dclose(h5dset);
-  istat = GRendaccess(ri_id);
-  free(h5cimage_name);
-  return SUCCEED;
-  }
-
-/**** palette routine. ****/
-/*-------------------------------------------------------------------------
- * Function:	gr_palette
- *
- * Purpose:     translate palette  into hdf5 dataset
- *              
- * Return:	FAIL if failed, SUCCEED if successful.
- *
- * In :	        
-                file_id: HDF4 identifier
-		ri: raster image id
-		h5dset: hdf5 dataset
-		h5_palgroup: hdf5 palette group
-
-   Out:         
- *-------------------------------------------------------------------------
- */	
-
-int gr_palette(int32 file_id,int32 ri_id,hid_t h5dset,hid_t h5_palgroup,int h4_attr) {
-
-  int32      pal_id;
-  uint16     pal_ref;
-  char       palref_str[MAXREF_LENGTH];
-  char       palg_name[MAX_GR_NAME];
-  char       image_index[MAX_GR_NAME];
-  int        check_pal;
-  int        check_palname;
-  int        pal_stat;
-  char*      h5pal_name=NULL;
-
-
-  /* get palette id */
-  pal_id  = GRgetlutid(ri_id,0);
-  if(pal_id == FAIL) {
-    printf("error in obtaining palette id. \n");
-    return FAIL;
-  }
-
-  pal_ref = GRluttoref(pal_id);
- 
-  if(pal_ref >0) {
-
-    /* convert reference number into string format. */
-    if(conv_int_str(pal_ref,palref_str)==FAIL) {
-      printf("error in converting palette reference number into string.\n");
-      return FAIL;
-    }
-
-    /* check whether this palette has been looked up already. */
-    check_pal = lookup(pal_ref,PAL_HASHSIZE,pal_hashtab);
-
-    if( check_pal < 0) {
-      printf("error at looking up palette table. \n");
-      return FAIL;
-    }
-
-    /* if check_pal equals to 1, this palette has already been 
-       converted into hdf5 dataset, just obtain the palette name.
-       if check_pal equals to 0, we will do the converting. */
-
-    if(check_pal == 1) {
-
-      h5pal_name = get_name(pal_ref,PAL_HASHSIZE,pal_hashtab,
-			    &check_palname);
-    
-      if (h5pal_name == NULL && check_palname == 0 ) {		      
-	 printf("error,cannot find group\n");			      
-	 return FAIL;							      
-      }							      
-									      
-      if (h5pal_name == NULL && check_palname == -1 ) {	       
-         printf("error,group name is not defined.\n");		      
-	 return FAIL;							      
-      }	
-						      
-    }							      
-    
-    if(check_pal == 0) {
-      /* do converting. */
-      strcpy(palg_name,HDF4_PALG);
-
-      /* obtain hdf5 dataset name converted from palette,
-	 no name for hdf4 palette.*/
-      h5pal_name = get_obj_aboname(NULL,palref_str,palg_name,HDF4_PALETTE);
-      if(h5pal_name == NULL) {
-	printf("error in getting hdf5 palette name.\n");
-	return FAIL;
-      }
-
-      if(set_name(pal_ref,PAL_HASHSIZE,pal_hashtab,h5pal_name)==FAIL) {
-	printf("error in setting object name.\n");	
-	free(h5pal_name);
-	return FAIL;
-      }
-
-      pal_stat = Palette_h4_to_h5(file_id,pal_id,h5_palgroup,h5pal_name,h4_attr);
-
-      if(pal_stat == FAIL) {
-        printf("error occurring in transferring palette into dataset. \n");
-	free(h5pal_name);
-	return FAIL;
-      }
-
-    }
-
-    if(create_pal_objref(h5dset,h5_palgroup,h5pal_name)== FAIL) {
-      printf("error in creating palette object reference.\n");
-      free(h5pal_name);
-      return FAIL;
-    }
-
-    if(h5pal_name != NULL) free(h5pal_name);
-
-    strcpy(image_index,HDF4_IMAGE_INDEXED);
-    if(h4_transpredattrs(h5dset,HDF4_IMAGE_SUBCLASS,image_index)== FAIL) {
-      printf("failed to transfer hdf4 image indexed.\n");
-      return FAIL;
-    }
-  }
-  return SUCCEED;
-}
-/***** end of palette application. *****/
-/*-------------------------------------------------------------------------
- * Function:	gr_tranattrs
- *
- * Purpose:     translate attributes of Image object into hdf5 dataset
- *              
- * Return:	FAIL if failed, SUCCEED if successful.
- *
- * In :	        
-                sri_id: RI identifier
-		sh5_dset: hdf5 dataset
-		snum_grattrs: number of attribute
-		check_gloflag: flag to check whether this attribute belongs
-		to gr interface.
-
-   Out:         
- *-------------------------------------------------------------------------
- */	
-int gr_tranattrs(int32 sri_id, hid_t sh5_dset,int snum_grattrs,
-		 int check_gloflag) {
-    
-  char      sgratrr_name[2*MAX_NC_NAME];
-  char      grglo[MAX_NC_NAME];
-  char*     grrepattr_name;
-  int32     count_sgradata;
-  int32     sgr_atype;
-  size_t    sh4_amemsize;
-  size_t    sh4_asize;
-
-  hid_t     sh5a_sid;
-  hid_t     sh5a_id;
-  hid_t     sh5_atype;
-  hid_t     sh5_amemtype;
-  hid_t     sh5str_type;
-  hid_t     sh5str_memtype;
-  hsize_t   sh5dims[MAX_VAR_DIMS];
-  void*     sgr_adata;
-  herr_t    sret;
-  int       i;
-
- 
-  for (i =0;i <snum_grattrs;i++) {
-      
-     if (GRattrinfo(sri_id,i,sgratrr_name,&sgr_atype,&count_sgradata)==FAIL){  
-        printf("unable to obtain attribute information. \n"); 
-        return FAIL;
-     }
-     
-     /*convert datatype for attribute. */
-
-     if(h4type_to_h5type(sgr_atype,&sh5_amemtype,&sh4_amemsize,
-			 &sh4_asize,&sh5_atype)==FAIL){
-       printf("unable to do type transferring.\n");
-       return FAIL;
-     }
-
-     sgr_adata = malloc(sh4_amemsize*count_sgradata);
-
-     if(GRgetattr(sri_id,i,(VOIDP)sgr_adata)==FAIL){
-       printf("unable to get GR attributes. \n");
-       return FAIL;
-     }
-	
-     /* if attribute doesn't have name, a default name is set. */
-     if(sgratrr_name[0] == '\0') {
-       grrepattr_name = trans_obj_name(DFTAG_RIG,i);
-       strcpy(sgratrr_name,grrepattr_name);
-       free(grrepattr_name);
-     }
-
-     /* if the sds attribute is a file attribute. */
-     if(check_gloflag == 1){
-       strcpy(grglo,GLOIMAGE);
-       strcat(sgratrr_name,"_");
-       strcat(sgratrr_name,grglo);
-     }
-     /* now do attribute-transferring.
-       1. deal with string data type
-       2. set attribute space.
-       3. get attribute name, set property list. */
-           
-     if (sh5_atype == H5T_STRING) {
-
-	sh5a_sid = H5Screate(H5S_SCALAR);
-
-	if (sh5a_sid < 0) {
-	   printf("failed to create attribute space for IMAGE. \n"); 
-	   return FAIL;
-	}
-
-	if ((sh5str_type = mkstr(count_sgradata*sh4_asize,H5T_STR_SPACEPAD))<0){
-           printf("error in making string for image attribute \n");
-	   return FAIL;
-	}
-
-        /* check this line later. */
-       	if ((sh5str_memtype = mkstr(count_sgradata*sh4_amemsize,
-				    H5T_STR_SPACEPAD))<0){
-	  printf("error in making memory string. \n");
-	  return FAIL;
-	}
-
-        sh5a_id = H5Acreate(sh5_dset,sgratrr_name,sh5str_type,sh5a_sid,
-			    H5P_DEFAULT);
-        if (sh5a_id <0) {
-	   printf("failed to obtain attribute id for IMAGE. \n");
-	   return FAIL;
-	}
-
-        sret = H5Awrite(sh5a_id,sh5str_memtype,(void *)sgr_adata);
-
-	if (sret <0) {
-	  printf("failed to obtain attribute of IMAGE.\n ");
-	  return FAIL;
-	}
-
-        sret = H5Sclose(sh5a_sid);
-        sret = H5Aclose(sh5a_id);
-     }
-	 
-     else {
-      
-       if (count_sgradata == 1) {
-
-	 sh5a_sid = H5Screate(H5S_SCALAR);
-
-	 if (sh5a_sid < 0) {
-	   printf("failed to create space id. \n");
-	   return FAIL;
-	 }
-       }
-       else {
-
-          sh5dims[0] = count_sgradata;
-	  sh5a_sid =  H5Screate_simple(1,sh5dims,NULL);
-          
-	  if (sh5a_sid < 0) {
-	    printf("failed to create attribute space. \n");
-	    return FAIL;
-	  }
-       }
-
-       sh5a_id = H5Acreate(sh5_dset,sgratrr_name,sh5_atype,sh5a_sid,
-			   H5P_DEFAULT);
-
-       if(sh5a_id <0) {
-	 printf("failed to obtain attribute id. \n");
-         return FAIL;
-       }
-
-       sret = H5Awrite(sh5a_id,sh5_amemtype,(void *)sgr_adata);
-
-       if(sret <0) {
-	 printf("failed to obtain attribute.\n ");
-         return FAIL;
-       }
-       sret = H5Aclose(sh5a_id);
-       sret = H5Sclose(sh5a_sid);
-
-     }
-
-     free(sgr_adata);
-
-  }
-
-  return SUCCEED;
-}
-	
-/*-------------------------------------------------------------------------
- * Function:	create_pal_objref
- *
- * Purpose:     create object reference for palette
- *              
- * Return:	FAIL if failed, SUCCEED if successful.
- *
- * In :	        
-		h5dset: hdf5 dataset
-		h5_palgroup: hdf5 palette group
-		h5pal_name: hdf5 palette name
-
-   Out:         
- *-------------------------------------------------------------------------
- */	
-
-int create_pal_objref(hid_t h5dset,hid_t h5_palgroup,char *h5pal_name){
-
-    hobj_ref_t pal_refdat;
-    hsize_t    pal_refDims[1];
-    hid_t      pal_refSpace;
-    hid_t      pal_refType;
-    hid_t      attribID;
-    herr_t     ret;
-
-    pal_refDims[0] = 1;
-    pal_refSpace   = H5Screate_simple(1,pal_refDims,NULL);
-
-    if(pal_refSpace < 0) {
-      printf("error in obtaining reference space. \n");
-      return FAIL;
-    }
-
-    pal_refType    = H5Tcopy(H5T_STD_REF_OBJ);
-    if(pal_refType < 0) {
-      printf("error in obtaining reference type. \n");
-      H5Sclose(pal_refSpace);
-      return FAIL;
-    }
-
-    ret            = H5Rcreate(&pal_refdat,h5_palgroup,h5pal_name,
-				H5R_OBJECT,-1);
-    if(ret < 0) {
-      printf("error in creating reference space. \n");
-      H5Sclose(pal_refSpace);
-      H5Tclose(pal_refType);
-      return FAIL;
-    }
-
-    attribID       = H5Acreate(h5dset,PALETTE,pal_refType,pal_refSpace,
-				H5P_DEFAULT);
-
-    if(attribID < 0) {
-      printf("error in obtaining attribute ID. \n");
-      H5Sclose(pal_refSpace);
-      H5Tclose(pal_refType);
-      return FAIL;
-    }
-
-    ret            = H5Awrite(attribID,pal_refType,(void *)&pal_refdat);
-
-      
-    H5Sclose(pal_refSpace);
-    if(H5Tclose(pal_refType)<0) {
-      printf("error closing palette reference type.\n");
-      H5Aclose(attribID);
-    }   
-    H5Aclose(attribID);
-    return SUCCEED;
-}
-
-
-uint16 get_RIref(int32 file_id,uint16 tag,int32 gr_ref,uint16* ptag_out){
-
-  DFdi di;
-  int32 found,GroupID;
-  int ri_ref = 0;
-  
-  if((GroupID = DFdiread(file_id,tag,(uint16)gr_ref))<0){
-    printf("cannot find gr_ref\n");
-    return ri_ref;
-  }
-
-  found = 0;
-  di.tag = DFTAG_NULL;
-  di.ref = 0;
-  while((found == 0) &&(DFdiget(GroupID,&di.tag,&di.ref)==0)){
-    if(di.tag == DFTAG_CI ||di.tag == DFTAG_RI || di.tag == DFTAG_CI8 ||
-       di.tag == DFTAG_RI8)
-      found = 1;
-  }
-
-  *ptag_out = di.tag;
-  ri_ref = di.ref;
-  ri_ref = 1;
-  if(!found) {
-    printf("cannot find ri_ref\n");
-    return 0;
-  }
-  DFdifree(GroupID);
-  return ri_ref;
-
-}