DEBUGPRINT("Start lmrcImage2Dto3DSingle\n")
if(NULL!=in) {
tmp.Header = prj.Header = in->Header;
+ tmp.HeaderMode = prj.HeaderMode = mrcFloatImage;
tmp.HeaderN.z = prj.HeaderN.z = 1;
mrcInit(&tmp, NULL); mrcInit(&prj, NULL);
if(linfo->flagDouble) {
/* Preparing Rho Filtering */
DEBUGPRINT("Start RhoFiltering2\n");
+ linfo->rhoInfo.Rmax = linfo->Rmax;
lmrcImageRhoFiltering2(out, counter, counterNum, Matrix, &linfo->rhoInfo, 0);
DEBUGPRINT("End RhoFiltering2\n");
prj = tmp;
if( -0.5<=prjx && prjx<prj->HeaderN.x-0.5
&&-0.5<=prjy && prjy<prj->HeaderN.y-0.5) {
data = out->ImageFloatImage[index];
- mrcPixelDataGet(prj, prjx, prjy, 0.0, &projdata, mrcPixelRePart, linfo->InterpolationMode);
+ //mrcPixelDataGet(prj, prjx, prjy, 0.0, &projdata, mrcPixelRePart, linfo->InterpolationMode);
if(linfo->flagDouble && linfo->rhoInfo.flagThicknessWeight) {
normz = fabs(tmpz/g3z);
if(normz < 0.9) {
}
sum+=data;
}
- out2D->ImageFloatImage[(int)(prjx+prjy*out2D->HeaderN.x+llinfo->section*out2D->HeaderN.x*out2D->HeaderN.y)] = sum/max;
+ out2D->ImageFloatImage[(int)(prjx+prjy*out2D->HeaderN.x)] = sum/max;
//mrcPixelDataSet(out2D, prjx, prjy, llinfo->section, sum/max, mrcPixelRePart);
+ //mrcPixelDataSet(out2D, prjx, prjy, 0, sum/max, mrcPixelRePart);
}
}
}