layer2mod.c

射线追踪程序

/*
 * This file is part of tomo3d
 *
 * Copyright (C) 2002, 2003, Sebastien Judenherc <[e-mail]>
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307
 * USA
 *
 */
/* transformation d'un fichier "layer" en fichier "block"
 * nlay ndir:
 *    nlay=couche a shifter
 *    ndir=0:shift en X, ndir=1:shift en X+Y, ndir=2:shift en Y
 *    ndir=10..17 -> shifts par 1/3 en X et Y
 */

#include <stdio.h>
#include <stdlib.h>
#include "tomo3d.h"

int main(int argc, char **argv)
{
        struct block_t * blk=NULL;
        struct layer_t *layers;
        int nbl,nbb,nc=-1,nd=-1,j,k,i;
        char layname[255],blkname[255];
        FILE *f;
        if (argc<2) {
                fprintf(stderr,"%s <name> [nlay ndir]\n",argv[0]);
                exit(0);
        }
        if (argc>3) {
                sscanf(argv[2],"%d",&nc);
                sscanf(argv[3],"%d",&nd);
        }
        sprintf(blkname,"%s.blk",argv[1]);
        sprintf(layname,"%s.lay",argv[1]);
        fprintf(stderr,"================ layer2mod ================\n");
        if ((blk=readlayers(layname, &nbb, &layers, &nbl))==NULL) {
                exit(0);
        }
        if (nc>-1) {
                fprintf(stderr,">>>> SHIFT: layer %d, dir %d\n",nc,nd);
        }
        printlayers(stderr,layers,nbl);
        if (nc>-1) {
                for (k=i=0; (i<nbl)&&(i<nc); i++) k+=layers[i].Nx*layers[i].Ny;
                for (j=k; j<k+layers[i].Nx*layers[i].Ny; j++) {
                        switch (nd) {
                                case 0 : blk[j].x+=(0.5*layers[i].dx);
                                         break;
                                case 1 :
                                         blk[j].x+=(0.5*layers[i].dx);
                                         blk[j].y+=(0.5*layers[i].dy);
                                         break;
                                case 2 :
                                         blk[j].y+=(0.5*layers[i].dy);
                                         break;
                                case 10:
                                         blk[j].x+=(1.*layers[i].dx/3.);
                                         break;
                                case 11:
                                         blk[j].x+=(2.*layers[i].dx/3.);
                                         break;
                                case 12:
                                         blk[j].y+=(1.*layers[i].dy/3.);
                                         break;
                                case 13:
                                         blk[j].y+=(1.*layers[i].dy/3.);
                                         blk[j].x+=(1.*layers[i].dx/3.);
                                         break;
                                case 14:
                                         blk[j].y+=(1.*layers[i].dy/3.);
                                         blk[j].x+=(2.*layers[i].dx/3.);
                                         break;
                                case 15:
                                         blk[j].y+=(2.*layers[i].dy/3.);
                                         break;
                                case 16:
                                         blk[j].y+=(2.*layers[i].dy/3.);
                                         blk[j].x+=(1.*layers[i].dx/3.);
                                         break;
                                case 17:
                                         blk[j].y+=(2.*layers[i].dy/3.);
                                         blk[j].x+=(2.*layers[i].dx/3.);
                                         break;
                                default:
                                         fprintf(stderr,"bad nd: %d\n",nd);
                                         exit(1);
                        }
                }
        }
        if (NULL==(f=fopen(blkname,"wt"))) {
                perror(blkname);
                exit(1);
        }
        printblocks(f,blk,nbb);
        fclose(f);
        return(0);
}