gdata.c

Go to the documentation of this file.

 
#include "gdata.h"
#include "jbxl_state.h"
 
 
int     ZeroBase = 0;
double  RZxy     = 1.0;     // 3D 画像での x,y 対 z方向比　(z方向のピクセル間隔)/(x,y方向のピクセル間隔）
int     ChkRZxy  = OFF;     // RZxy　が設定されているなら ON, 設定されていないなら OFF 
 
 
void  free_CmnHead(CmnHead* hd)
{
    if (hd!=NULL) {
        if (hd->buf!=NULL)   free(hd->buf);
        if (hd->grptr!=NULL) free(hd->grptr);
        init_CmnHead(hd);
    }
}
 
 
void  init_CmnHead(CmnHead* hd)
{
    if (hd!=NULL) {
        memset(hd, 0, sizeof(CmnHead));
        hd->kind = HEADER_NONE;
    }
}
 
 
BSGraph  make_BSGraph(int xs, int ys, int zs)
{
    int  i;
    BSGraph  vp;
 
    memset(&vp, 0, sizeof(BSGraph));
    if (xs==0 || ys==0) {
        vp.state = JBXL_GRAPH_IVDARG_ERROR;
        return vp;
    }
 
    vp.xs = xs;
    vp.ys = ys;
    if (zs>0) vp.zs = zs;
    else      vp.zs = 1;
    vp.state = JBXL_NORMAL;
 
    vp.gp = (uByte*)malloc(vp.xs*vp.ys*vp.zs*sizeof(uByte));
    if (vp.gp==NULL) {
        memset(&vp, 0, sizeof(BSGraph));
        vp.state = JBXL_GRAPH_MEMORY_ERROR;
        return vp;
    }
    for (i=0; i<vp.xs*vp.ys*vp.zs; i++) vp.gp[i] = 0;
 
    return vp;
}
 
 
WSGraph  make_WSGraph(int xs, int ys, int zs)
{
    int  i;
    WSGraph  vp;
 
    memset(&vp, 0, sizeof(WSGraph));
    if (xs==0 || ys==0) {
        vp.state = JBXL_GRAPH_IVDARG_ERROR;
        return vp;
    }
 
    vp.xs = xs;
    vp.ys = ys;
    if (zs>0) vp.zs = zs;
    else      vp.zs = 1;
    vp.state = JBXL_NORMAL;
 
    vp.gp = (sWord*)malloc(vp.xs*vp.ys*vp.zs*sizeof(sWord));
    if (vp.gp==NULL) {
        memset(&vp, 0, sizeof(WSGraph));
        vp.state = JBXL_GRAPH_MEMORY_ERROR;
        return vp;
    }
    for (i=0; i<vp.xs*vp.ys*vp.zs; i++) vp.gp[i] = 0;
 
    return vp;
}
 
 
FSGraph  make_FSGraph(int xs, int ys, int zs)
{
    int  i;
    FSGraph  vp;
 
    memset(&vp, 0, sizeof(FSGraph));
    if (xs==0 || ys==0) {
        vp.state = JBXL_GRAPH_IVDARG_ERROR;
        return vp;
    }
 
    vp.xs = xs;
    vp.ys = ys;
    if (zs>0) vp.zs = zs;
    else      vp.zs = 1;
    vp.state = JBXL_NORMAL;
 
    vp.gp = (double*)malloc(vp.xs*vp.ys*vp.zs*sizeof(double));
    if (vp.gp==NULL) {
        memset(&vp, 0, sizeof(FSGraph));
        vp.state = JBXL_GRAPH_MEMORY_ERROR;
        return vp;
    }
    for (i=0; i<vp.xs*vp.ys*vp.zs; i++) vp.gp[i] = 0.0;
 
    return vp;
}
 
 
ISGraph  make_ISGraph(int xs, int ys, int zs)
{
    int  i;
    ISGraph  vp;
 
    memset(&vp, 0, sizeof(ISGraph));
    if (xs==0 || ys==0) {
        vp.state = JBXL_GRAPH_IVDARG_ERROR;
        return vp;
    }
 
    vp.xs = xs;
    vp.ys = ys;
    if (zs>0) vp.zs = zs;
    else      vp.zs = 1;
    vp.state = JBXL_NORMAL;
 
    vp.gp = (int*)malloc(vp.xs*vp.ys*vp.zs*sizeof(int));
    if (vp.gp==NULL) {
        memset(&vp, 0, sizeof(ISGraph));
        vp.state = JBXL_GRAPH_MEMORY_ERROR;
        return vp;
    }
    for (i=0; i<vp.xs*vp.ys*vp.zs; i++) vp.gp[i] = 0;
 
    return vp;
}
 
 
VSGraph  make_VSGraph(int xs, int ys, int zs)
{
    int  i;
    VSGraph  vp;
 
    memset(&vp, 0, sizeof(VSGraph));
    if (xs==0 || ys==0) {
        vp.state = JBXL_GRAPH_IVDARG_ERROR;
        return vp;
    }
 
    vp.xs = xs;
    vp.ys = ys;
    if (zs>0) vp.zs = zs;
    else      vp.zs = 1;
    vp.state = JBXL_NORMAL;
 
    vp.gp = (vector*)malloc(xs*ys*zs*sizeof(vector));
    if (vp.gp==NULL) {
        memset(&vp, 0, sizeof(VSGraph));
        vp.state = JBXL_GRAPH_MEMORY_ERROR;
        return vp;
    }
    for (i=0; i<xs*ys*zs; i++) vp.gp[i] = set_vector(0.0, 0.0, 0.0);
 
    return vp;
}
 
 
MSGraph  make_MSGraph(int xs, int ys, int zs, int depth)
{
    MSGraph  vp;
 
    memset(&vp, 0, sizeof(MSGraph));
    if (xs==0 || ys==0) {
        vp.state = JBXL_GRAPH_IVDARG_ERROR;
        return vp;
    }
 
    vp.xs = xs;
    vp.ys = ys;
    if (zs>0) vp.zs = zs;
    else      vp.zs = 1;
    if (depth>0) vp.depth = depth;
    else         vp.depth = 1;
    vp.state = JBXL_NORMAL;
 
    vp.gp = (unsigned char*)malloc(xs*ys*zs*((depth+7)/8));
    if (vp.gp==NULL) {
        memset(&vp, 0, sizeof(MSGraph));
        vp.state = JBXL_GRAPH_MEMORY_ERROR;
        return vp;
    }
    memset(vp.gp, 0, xs*ys*zs*((depth+7)/8));
    return vp;
}
 
 
BSGraph* new_BSGraph(int xs, int ys, int zs)
{
    int  i;
    BSGraph* vp;
 
    vp = (BSGraph*)malloc(sizeof(BSGraph));
    if (vp==NULL) return NULL;
    memset(vp, 0, sizeof(BSGraph));
 
    if (xs==0 || ys==0) {
        vp->state = JBXL_GRAPH_IVDARG_ERROR;
        return vp;
    }
 
    vp->xs = xs;
    vp->ys = ys;
    if (zs>0) vp->zs = zs;
    else      vp->zs = 1;
    vp->state = JBXL_NORMAL;
 
    vp->gp = (uByte*)malloc(vp->xs*vp->ys*vp->zs*sizeof(uByte));
    if (vp->gp==NULL) {
        memset(vp, 0, sizeof(BSGraph));
        vp->state = JBXL_GRAPH_MEMORY_ERROR;
        return vp;
    }
    for (i=0; i<vp->xs*vp->ys*vp->zs; i++) vp->gp[i] = 0;
 
    return vp;
}
 
 
WSGraph* new_WSGraph(int xs, int ys, int zs)
{
    int  i;
    WSGraph* vp;
 
    vp = (WSGraph*)malloc(sizeof(WSGraph));
    if (vp==NULL) return NULL;
    memset(vp, 0, sizeof(WSGraph));
 
    if (xs==0 || ys==0) {
        vp->state = JBXL_GRAPH_IVDARG_ERROR;
        return vp;
    }
 
    vp->xs = xs;
    vp->ys = ys;
    if (zs>0) vp->zs = zs;
    else      vp->zs = 1;
    vp->state = JBXL_NORMAL;
 
    vp->gp = (sWord*)malloc(vp->xs*vp->ys*vp->zs*sizeof(sWord));
    if (vp->gp==NULL) {
        memset(vp, 0, sizeof(WSGraph));
        vp->state = JBXL_GRAPH_MEMORY_ERROR;
        return vp;
    }
    for (i=0; i<vp->xs*vp->ys*vp->zs; i++) vp->gp[i] = 0;
 
    return vp;
}
 
 
FSGraph* new_FSGraph(int xs, int ys, int zs)
{
    int  i;
    FSGraph* vp;
 
    vp = (FSGraph*)malloc(sizeof(FSGraph));
    if (vp==NULL) return NULL;
    memset(vp, 0, sizeof(FSGraph));
 
    if (xs==0 || ys==0) {
        vp->state = JBXL_GRAPH_IVDARG_ERROR;
        return vp;
    }
 
    vp->xs = xs;
    vp->ys = ys;
    if (zs>0) vp->zs = zs;
    else      vp->zs = 1;
    vp->state = JBXL_NORMAL;
 
    vp->gp = (double*)malloc(vp->xs*vp->ys*vp->zs*sizeof(double));
    if (vp->gp==NULL) {
        memset(vp, 0, sizeof(FSGraph));
        vp->state = JBXL_GRAPH_MEMORY_ERROR;
        return vp;
    }
    for (i=0; i<vp->xs*vp->ys*vp->zs; i++) vp->gp[i] = 0.0;
 
    return vp;
}
 
 
ISGraph* new_ISGraph(int xs, int ys, int zs)
{
    int  i;
    ISGraph* vp;
 
    vp = (ISGraph*)malloc(sizeof(ISGraph));
    if (vp==NULL) return NULL;
    memset(vp, 0, sizeof(ISGraph));
 
    if (xs==0 || ys==0) {
        vp->state = JBXL_GRAPH_IVDARG_ERROR;
        return vp;
    }
 
    vp->xs = xs;
    vp->ys = ys;
    if (zs>0) vp->zs = zs;
    else      vp->zs = 1;
    vp->state = JBXL_NORMAL;
 
    vp->gp = (int*)malloc(vp->xs*vp->ys*vp->zs*sizeof(int));
    if (vp->gp==NULL) {
        memset(vp, 0, sizeof(ISGraph));
        vp->state = JBXL_GRAPH_MEMORY_ERROR;
        return vp;
    }
    for (i=0; i<vp->xs*vp->ys*vp->zs; i++) vp->gp[i] = 0;
 
    return vp;
}
 
 
VSGraph* new_VSGraph(int xs, int ys, int zs)
{
    int  i;
    VSGraph* vp;
 
    vp = (VSGraph*)malloc(sizeof(VSGraph));
    if (vp==NULL) return NULL;
    memset(vp, 0, sizeof(VSGraph));
 
    if (xs==0 || ys==0) {
        vp->state = JBXL_GRAPH_IVDARG_ERROR;
        return vp;
    }
 
    vp->xs = xs;
    vp->ys = ys;
    if (zs>0) vp->zs = zs;
    else      vp->zs = 1;
    vp->state = JBXL_NORMAL;
 
    vp->gp = (vector*)malloc(vp->xs*vp->ys*vp->zs*sizeof(vector));
    if (vp->gp==NULL) {
        memset(vp, 0, sizeof(VSGraph));
        vp->state = JBXL_GRAPH_MEMORY_ERROR;
        return vp;
    }
    for (i=0; i<vp->xs*vp->ys*vp->zs; i++) vp->gp[i] = set_vector(0.0, 0.0, 0.0);
 
    return vp;
}
 
 
MSGraph* new_MSGraph(int xs, int ys, int zs, int depth)
{
    MSGraph* vp;
 
    vp = (MSGraph*)malloc(sizeof(MSGraph));
    if (vp==NULL) return NULL;
    memset(vp, 0, sizeof(MSGraph));
 
    if (xs==0 || ys==0) {
        vp->state = JBXL_GRAPH_IVDARG_ERROR;
        return vp;
    }
 
    vp->xs = xs;
    vp->ys = ys;
    if (zs>0) vp->zs = zs;
    else      vp->zs = 1;
    if (depth>0) vp->depth = depth;
    else         vp->depth = 1;
    vp->state = JBXL_NORMAL;
 
    vp->gp = (unsigned char*)malloc(vp->xs*vp->ys*vp->zs*((depth+7)/8));
    if (vp->gp==NULL) {
        memset(vp, 0, sizeof(MSGraph));
        vp->state = JBXL_GRAPH_MEMORY_ERROR;
        return vp;
    }
    memset(vp->gp, 0, vp->xs*vp->ys*vp->zs*((depth+7)/8));
 
    return vp;
}
 
 
ISGraph  W2ISGraph(WSGraph vp)
{
    int  i;
    ISGraph  ix;
 
    ix.xs = vp.xs;
    ix.ys = vp.ys;
    ix.zs = vp.zs;
    ix.state = vp.state;
 
    ix.gp = (int*)malloc(ix.xs*ix.ys*ix.zs*sizeof(int));
    if (ix.gp==NULL) {
        memset(&ix, 0, sizeof(ISGraph));
        ix.state = JBXL_GRAPH_MEMORY_ERROR;
        return ix;
    }
    for (i=0; i<ix.xs*ix.ys*ix.zs; i++) ix.gp[i] = (int)vp.gp[i];
 
    return ix;
}
 
 
FSGraph  W2FSGraph(WSGraph vp)
{
    int  i;
    FSGraph  ix;
 
    ix.xs = vp.xs;
    ix.ys = vp.ys;
    ix.zs = vp.zs;
    ix.state = vp.state;
 
    ix.gp = (double*)malloc(ix.xs*ix.ys*ix.zs*sizeof(double));
    if (ix.gp==NULL) {
        memset(&ix, 0, sizeof(FSGraph));
        ix.state = JBXL_GRAPH_MEMORY_ERROR;
        return ix;
    }
    for (i=0; i<ix.xs*ix.ys*ix.zs; i++) ix.gp[i]=(double)vp.gp[i];
 
    return ix;
}
 
 
WSGraph  B2WSGraph(BSGraph vp)
{
    int  i;
    WSGraph  ix;
 
    ix.xs = vp.xs;
    ix.ys = vp.ys;
    ix.zs = vp.zs;
    ix.state = vp.state;
 
    ix.gp = (sWord*)malloc(ix.xs*ix.ys*ix.zs*sizeof(sWord));
    if (ix.gp==NULL) {
        memset(&ix, 0, sizeof(WSGraph));
        ix.state = JBXL_GRAPH_MEMORY_ERROR;
        return ix;
    }
    for (i=0; i<ix.xs*ix.ys*ix.zs; i++) ix.gp[i] = (sWord)vp.gp[i];
 
    return ix;
}
 
 
WSGraph  I2WSGraph(ISGraph vp)
{
    int  i;
    WSGraph  ix;
 
    ix.xs = vp.xs;
    ix.ys = vp.ys;
    ix.zs = vp.zs;
    ix.state = vp.state;
 
    ix.gp = (sWord*)malloc(ix.xs*ix.ys*ix.zs*sizeof(sWord));
    if (ix.gp==NULL) {
        memset(&ix, 0, sizeof(WSGraph));
        ix.state = JBXL_GRAPH_MEMORY_ERROR;
        return ix;
    }
    for (i=0; i<ix.xs*ix.ys*ix.zs; i++) ix.gp[i] = (sWord)vp.gp[i];
 
    return ix;
}
 
 
WSGraph  F2WSGraph(FSGraph vp)
{
    int  i;
    WSGraph ix;
 
    ix.xs = vp.xs;
    ix.ys = vp.ys;
    ix.zs = vp.zs;
    ix.state = vp.state;
 
    ix.gp = (sWord*)malloc(ix.xs*ix.ys*ix.zs*sizeof(sWord));
    if (ix.gp==NULL) {
        memset(&ix, 0, sizeof(WSGraph));
        ix.state = JBXL_GRAPH_MEMORY_ERROR;
        return ix;
    }
    for (i=0; i<ix.xs*ix.ys*ix.zs; i++) ix.gp[i] = (sWord)vp.gp[i];
 
    return ix;
}
 
 
FSGraph  V2FSGraph(VSGraph vp)
{
    int  i;
    FSGraph ix;
 
    ix = make_FSGraph(vp.xs, vp.ys, vp.zs);
    if (ix.gp==NULL) return ix;
 
    for (i=0; i<ix.xs*ix.ys*ix.zs; i++) ix.gp[i] = vp.gp[i].n;
    ix.state = vp.state;
 
    return ix;
}
 
 
WSGraph  V2WSGraph(VSGraph vp)
{
    int  i;
    WSGraph  ix;
 
    ix = make_WSGraph(vp.xs, vp.ys, vp.zs);
    if (ix.gp==NULL) return ix;
 
    for (i=0; i<ix.xs*ix.ys*ix.zs; i++) ix.gp[i] = (sWord)vp.gp[i].n;
    ix.state = vp.state;
 
    return ix;
}
 
 
void  init_IRBound(IRBound* rb)
{
    rb->xmax = rb->ymax = rb->zmax = 0;
    rb->xmin = rb->ymin = rb->zmin = SINTMAX;
    rb->misc = OFF;
}
 
 
void  init_DRBound(DRBound* rb)
{
    rb->xmax = rb->ymax = rb->zmax = 0.0;
    rb->xmin = rb->ymin = rb->zmin = DBL_MAX;
    rb->misc = OFF;
}
 
 
void  set_RZxy(double rzm)
{
    ChkRZxy = ON;
    RZxy = rzm;
    //fprintf(stderr,"set_RZxy: set RZxy to %f\n", rzm);
}
 
 
int  chk_RZxy(void)   
{
    if (ChkRZxy==OFF) {
        DEBUG_MODE PRINT_MESG("CHK_RZXY: WARNING: undefined RZxy!!!\n");
        //exit(1);
        return FALSE;
    }
    return TRUE;
}