summaryrefslogtreecommitdiffstats
path: root/tools/h4toh5image.c
diff options
context:
space:
mode:
Diffstat (limited to 'tools/h4toh5image.c')
-rw-r--r--tools/h4toh5image.c861
1 files changed, 0 insertions, 861 deletions
diff --git a/tools/h4toh5image.c b/tools/h4toh5image.c
deleted file mode 100644
index 8fbbdc0..0000000
--- a/tools/h4toh5image.c
+++ /dev/null
@@ -1,861 +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) {
-
- 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;
- int32 chunk_dims[2];
- int32 c_flags;
-
- /* 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];
- herr_t ret;
- hid_t create_plist;
-
- /* 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,
- NULL, 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;
- }
-
- 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. */
-
-
-
- h5dims[0] = edges[1]-start[1];
- h5dims[1] = edges[0]-start[0];
-
- 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);
-
- if(c_flags == HDF_CHUNK || c_flags == (HDF_CHUNK | HDF_COMP)
- || c_flags == (HDF_CHUNK | HDF_NBIT) ){
-
- 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);
- 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);
- 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 { /* compound datatype. */
-
- h5_ctype = H5Tcreate(H5T_COMPOUND,ncomp*h4size);
- if (h5_ctype < 0) {
- printf("error in generating hdf5 compound data type. \n");
- free(image_data);
- free(h5cimage_name);
- H5Pclose(create_plist);
- return FAIL;
- }
-
- h5_cmemtype = H5Tcreate(H5T_COMPOUND,ncomp*h4memsize);
- if (h5_cmemtype < 0) {
- printf("error in generating hdf5 memory compound data type. \n");
- free(image_data);
- free(h5cimage_name);
- H5Pclose(create_plist);
- return FAIL;
- }
-
-
- fielddim[0] = ncomp;
-
- {
- hid_t arr_type; /* Array datatype for inserting fields */
-
- /* Create array datatype */
- if((arr_type=H5Tarray_create(h5ty_id,1,fielddim,NULL))<0) {
- printf("error creating array datatype.\n");
- free(image_data);
- free(h5cimage_name);
- H5Pclose(create_plist);
- return FAIL;
- }
-
- ret = H5Tinsert(h5_ctype,"HDF4Image_data",0,arr_type);
- if(ret < 0) {
- printf("error in inserting array of compound datatype. \n");
- free(image_data);
- free(h5cimage_name);
- H5Pclose(create_plist);
- return FAIL;
- }
-
- /* Close array datatype */
- if(H5Tclose(arr_type)<0) {
- printf("error closing array datatype.\n");
- free(image_data);
- free(h5cimage_name);
- H5Pclose(create_plist);
- return FAIL;
- }
-
- /* Create array datatype */
- if((arr_type=H5Tarray_create(h5memtype,1,fielddim,NULL))<0) {
- printf("error creating array datatype.\n");
- free(image_data);
- free(h5cimage_name);
- H5Pclose(create_plist);
- return FAIL;
- }
-
- ret = H5Tinsert(h5_cmemtype,"HDF4Image_data",0,arr_type);
- if(ret < 0) {
- printf("error in inserting array of compound datatype at memory. \n");
- free(image_data);
- free(h5cimage_name);
- H5Pclose(create_plist);
- return FAIL;
- }
-
- /* Close array datatype */
- if(H5Tclose(arr_type)<0) {
- printf("error closing array datatype.\n");
- free(image_data);
- free(h5cimage_name);
- H5Pclose(create_plist);
- return FAIL;
- }
- }
-
- h5d_sid = H5Screate_simple(2,h5dims,NULL);
- if(h5d_sid < 0) {
- printf("error in creating space. \n");
- free(image_data);
- free(h5cimage_name);
- H5Pclose(create_plist);
- return FAIL;
- }
-
- h5dset = H5Dcreate(h5_group,h5cimage_name,h5_ctype,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,h5_cmemtype,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;
- }
- ret = H5Tclose(h5_ctype);
- if(ret < 0) {
- printf("error in closing h5_ctype. \n");
- }
- ret = H5Tclose(h5_cmemtype);
- if(ret <0) {
- printf("error in closing h5_cmemtype. \n");
- }
- }
-
-/* 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(image_data);
- 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);
- free(image_data);
- 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);
- free(image_data);
- 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");
- free(image_data);
- 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);
- free(image_data);
- return FAIL;
- }
-
- 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);
- free(image_data);
- 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);
- free(image_data);
- return FAIL;
- }
-
- 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);
- free(image_data);
- 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);
- free(image_data);
- return FAIL;
- }
-
- /* deal with palette. */
-
- if(gr_palette(file_id,ri_id,h5dset,h5_palgroup)== FAIL) {
- printf("error in translating palette into h5 dataset.\n");
- H5Pclose(create_plist);
- H5Sclose(h5d_sid);
- H5Dclose(h5dset);
- free(h5cimage_name);
- free(image_data);
- return FAIL;
- }
-
- ret = H5Pclose(create_plist);
- ret = H5Sclose(h5d_sid);
- ret = H5Dclose(h5dset);
- istat = GRendaccess(ri_id);
- free(image_data);
- 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) {
-
- 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);
-
- 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;
-}