CT图像重建SART算法matlab代码

#include <stdio.h>
#include <math.h>
#include <time.h>
#include <stdlib.h>

#define pi 3.1415926

#define Resolution 256
#define R Resolution
#define TubeNum 256
#define TN TubeNum
#define ProjNum 200
#define PN ProjNum
#define PixelSize 2.0/256
#define PS PixelSize
#define RelativeTubeSpace 1  /* Compared to PixelSize */
#define RTS RelativeTubeSpace
#define AbsoluteTubeSpace PS*RTS
#define ATS AbsoluteTubeSpace
#define IterationNum 999
#define IN IterationNum

int main ( )
{
    double begin=time(NULL);

    typedef struct {
        double l;
        int index;

    }Coef;

    Coef CoefArray[2*R];

    FILE * sino;    //Read sinogram in
    sino=fopen("Sinogram","rb");
    if (sino==NULL)
    {
                printf("Open file failure!\n");
    }
    double * Sinogram=new double[PN*TN];
    fread(Sinogram,sizeof(double),PN*TN,sino);
    if (feof(sino)==0)
    {
        printf("Error!Reading sinogram has not been finished.\n");
    }
    fclose(sino);

    FILE * ei;
    ei=fopen("EI","rb");
    if (ei==NULL)
    {
                printf("Open file failure!\n");
    }
    double * EI=new double[PN*TN];
    fread(EI,sizeof(double),PN*TN,ei);
    if (feof(ei)==0)
    {
        printf("Error!Reading EI has not been finished.\n");
    }