F5image.c

Go to the documentation of this file.

#include "F5image.h"
#include "F5defs.h"
#include "F5X.h"
#include "F5B.h"
#include "F5R.h"
#include "F5F.h"
#include <string.h>
 
void initF5image(F5image*This)
{
        This->height = 0;
        This->width  = 0;
 
        This->time      = 0.0;
        This->viewangle = 0.0;
 
        This->viewpoint.x =
        This->viewpoint.y =
        This->viewpoint.z = 0.0;
 
        This->origin.x =
        This->origin.y =
        This->origin.z = 0.0;
 
        This->x_direction.x =
        This->x_direction.y =
        This->x_direction.z = 0.0;
 
        This->y_direction.x =
        This->y_direction.y =
        This->y_direction.z = 0.0;
 
        This->directions   = NULL;
 
        This->intensity    = NULL;
        This->transparency = NULL;
        This->zvalues      = NULL;
 
        This->intensity16    = NULL;
        This->transparency16 = NULL;
        This->zvalues16      = NULL;
 
        This->intensity8    = NULL;
        This->transparency8 = NULL;
        This->zvalues8      = NULL;
 
        This->alpha         = NULL;
 
        This->spectrum      = NULL;
}
 
 
void    saveF5image(F5image*This, const char*filename,
                    const char*vp_txt, const char*scenery)
{
char    buf[4096];
hid_t   File_id = F5append(filename);
 
        printf("saveF5image(, %s, %s, %s)\n", filename, vp_txt, scenery);
 
        if (!vp_txt || !*vp_txt)
                vp_txt = "StandardView";
 
        if (File_id<0)
        {
                printf("saveF5image(): Cannot open %s\n", filename);
                return;
        }
 
        if (!scenery)
        {
        char*s;
                scenery = strrchr(filename, '/');
                if (scenery)
                        strcpy(buf, scenery+1);
                else
                        strcpy(buf, filename);
                s = strrchr(buf, '.');
                if (s) *s = 0;
                scenery = buf;
                printf("saveF5image(): scenery set to '%s'\n", scenery);
        }
/*
        if (This->intensity)
                puts("Have image float data");
        if (This->transparency)
                puts("Have image transparency data");
        if (This->alpha)
                puts("Have image alpha data");
*/
        {
        F5Path*myPath = F5Rcreate_cartesian_nD(File_id, This->time,
                                               scenery, /* gridname, */
                                               2,
                                               vp_txt); /* coordinate_system); */
 
/*
        hid_t   Slice_hid = F5Bappend_slice(File_id, &This->time);
                {
                hid_t   Grid_hid    = F5Gappend(Slice_hid,  scenery, 0) ;
                hid_t   Charts_hid  = F5Gappend(Grid_hid,  FIBER_HDF5_CHARTS_PER_GRID, 0) ;
                hid_t   VPChart_hid = F5Gappend(Charts_hid, vp_txt, 0)  ;
                hid_t   Transf_hid  = F5Gappend(VPChart_hid, FIBER_HDF5_DEFAULT_CHART, 0)  ;
                hid_t   Top_hid     = F5Gappend(Grid_hid,  FIBER_HDF5_POINTS, 0)  ;
                hid_t   Rep_hid     = F5Gappend(Top_hid,  vp_txt, 0)     ;
*/
                hsize_t imgdims[2];
                        imgdims[0] = This->width ;
                        imgdims[1] = This->height;
 
                        { F5_vec3_double_t delta; delta.x = delta.y = delta.z = 1.0;
 
                        F5Fwrite_linear(myPath, FIBER_HDF5_POSITIONS_STRING, 2, imgdims,
                                        F5T_VEC3_DOUBLE, (void*)&(This->origin), &delta);
                        }
 
                        /*
                        //
                        // Storing image metadata information
                        //
                        {
                        hid_t   Pos_hid    = F5Gappend(Rep_hid,   FIBER_HDF5_POSITIONS_STRING, 0) ;
                        hid_t   space_hid = H5Screate(H5S_SCALAR);
                        hid_t   attr_hid  = F5Aappend(Pos_hid, "origin", F5T_VEC3_DOUBLE, space_hid, H5P_DEFAULT);
                                H5Awrite(attr_hid, F5T_VEC3_DOUBLE, (void*)&(This->origin));
                                H5Aclose(attr_hid);
 
                                attr_hid  = F5Aappend(Pos_hid, "x", F5T_VEC3_DOUBLE, space_hid, H5P_DEFAULT);
                                H5Awrite(attr_hid, F5T_VEC3_DOUBLE, (void*)&(This->x_direction));
                                H5Aclose(attr_hid);
 
                                attr_hid  = F5Aappend(Pos_hid, "y", F5T_VEC3_DOUBLE, space_hid, H5P_DEFAULT);
                                H5Awrite(attr_hid, F5T_VEC3_DOUBLE, (void*)&(This->y_direction));
                                H5Aclose(attr_hid);
 
                                attr_hid  = F5Aappend(Transf_hid, "shift", F5T_VEC3_DOUBLE, space_hid, H5P_DEFAULT);
                                H5Awrite(attr_hid, F5T_VEC3_DOUBLE, (void*)&(This->viewpoint));
                                H5Aclose(attr_hid);
 
                                H5Sclose(space_hid);
                                H5Gclose(Pos_hid);
                        }
                        */
 
                        //
                        // Storing image data information
                        //
                        {
//                      hid_t img_ds = H5Screate_simple(2, imgdims, NULL);
                        hid_t creator_hid = H5Pcreate(H5P_DATASET_CREATE);
 
                                if (This->directions)
                                {
                                        F5Fwrite( myPath, F5_IMAGE_FINAL_VIEW, 2, imgdims, F5T_RGB_REAL, F5T_RGB_REAL, This->directions, creator_hid);
                        //      hid_t   id = F5Dappend(Rep_hid, F5_IMAGE_FINAL_VIEW, F5T_VEC3_FLOAT, img_ds, creator_hid);
                        //              H5Dwrite(id, F5T_VEC3_FLOAT, H5S_ALL, H5S_ALL, H5P_DEFAULT, This->directions);
                        //              H5Dclose(id);
                                }
 
                                if (This->intensity)
                                {
                                //                      puts("Saving intensity...");
                                //      F5FwriteIMAGE( myPath, F5_IMAGE_RGB_FLOAT_INTENSITY, 2, imgdims, F5T_RGB_REAL, F5T_RGB_REAL, This->intensity, creator_hid);
                                        F5Fwrite     ( myPath, F5_IMAGE_RGB_FLOAT_INTENSITY, 2, imgdims, F5T_RGB_REAL, F5T_RGB_REAL, This->intensity, creator_hid);
 
//                              hid_t   id = F5Dappend(Rep_hid, F5_IMAGE_RGB_FLOAT_INTENSITY, F5T_RGB_REAL, img_ds, creator_hid);
//                                      H5Dwrite(id, F5T_RGB_REAL, H5S_ALL, H5S_ALL, H5P_DEFAULT, This->intensity);
//                                      H5Dclose(id);
                                }
 
                                if (This->transparency)
                                {
                                        F5Fwrite( myPath, F5_IMAGE_RGB_FLOAT_TRANSPARENCY, 2, imgdims, F5T_RGB_REAL, F5T_RGB_REAL, This->transparency, creator_hid);
 
//                              hid_t   id = F5Dappend(Rep_hid, F5_IMAGE_RGB_FLOAT_TRANSPARENCY, F5T_RGB_REAL, img_ds, creator_hid);
//                                      H5Dwrite(id, F5T_RGB_REAL, H5S_ALL, H5S_ALL, H5P_DEFAULT, This->transparency);
//                                      H5Dclose(id);
                                }
                                /*
                                if (This->zvalues)
                                {
                                hid_t   id = F5Dappend(Rep_hid, F5_IMAGE_ZBUFFER, H5T_NATIVE_FLOAT, img_ds, creator_hid);
                                        H5Dwrite(id, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, H5P_DEFAULT, This->zvalues);
                                        H5Dclose(id);
                                }
 
                                if (This->intensity16)
                                {
                                hid_t   id = F5Dappend(Rep_hid, F5_IMAGE_RGB16_INTENSITY, F5T_RGB16, img_ds, creator_hid);
                                        H5Dwrite(id, F5T_RGB16, H5S_ALL, H5S_ALL, H5P_DEFAULT, This->intensity16);
                                        H5Dclose(id);
                                }
 
                                if (This->transparency16)
                                {
                                hid_t   id = F5Dappend(Rep_hid, F5_IMAGE_RGB16_TRANSPARENCY, F5T_RGB16, img_ds, creator_hid);
                                        H5Dwrite(id, F5T_RGB16, H5S_ALL, H5S_ALL, H5P_DEFAULT, This->transparency16);
                                        H5Dclose(id);
                                }
 
                                if (This->intensity8)
                                {
                                hid_t   id = F5Dappend(Rep_hid,  F5_IMAGE_RGB8_INTENSITY, F5T_RGB, img_ds, creator_hid);
                                        H5Dwrite(id, F5T_RGB, H5S_ALL, H5S_ALL, H5P_DEFAULT, This->intensity8);
                                        H5Dclose(id);
                                }
 
                                if (This->transparency8)
                                {
                                hid_t   id = F5Dappend(Rep_hid, F5_IMAGE_RGB8_TRANSPARENCY, F5T_RGB, img_ds, creator_hid);
                                        H5Dwrite(id, F5T_RGB, H5S_ALL, H5S_ALL, H5P_DEFAULT, This->transparency8);
                                        H5Dclose(id);
                                }
 
                                if (This->alpha)
                                {
                                hid_t   id = F5Dappend(Rep_hid, F5_IMAGE_ALPHA_CHANNEL, H5T_NATIVE_UCHAR, img_ds, creator_hid);
                                        H5Dwrite(id,  H5T_NATIVE_UCHAR, H5S_ALL, H5S_ALL, H5P_DEFAULT, This->alpha);
                                        H5Dclose(id);
                                }
 
                                H5Sclose(img_ds);
                                */
                                H5Pclose(creator_hid);
                        }
/*
                        H5Gclose(Transf_hid);
                        H5Gclose(VPChart_hid);
                        H5Gclose(Charts_hid);
                        H5Gclose(Rep_hid);
                        H5Gclose(Top_hid);
                        H5Gclose(Grid_hid);
                }
 
                H5Gclose( Slice_hid );
*/
                F5close( myPath );
        }
 
        H5Fclose(File_id);
}