}
} else {
- fprintf(stderr, "Not supported mode: %ld in lmrcImageYProjection\n", img->HeaderMode);
+ fprintf(stderr, "Not supported mode: %d in lmrcImageYProjection\n", img->HeaderMode);
exit(EXIT_FAILURE);
}
mrcStatDataSet(proj, 0);
sumForCylinderSurface = 0;
sdForCylinderSurface = 0;
- rmax = MIN(offsetx*img->HeaderLength.x,
- offsety*img->HeaderLength.y);
+ rmax = MIN(offsetx*img->HeaderLength.x, offsety*img->HeaderLength.y);
if (info->flagrmax){
if(info->rmax<rmax){
fprintf(stderr,"rmax from option:%f\n",info->rmax);
}
case meanOfSphereSurface:{
DEBUGPRINT("meanOfSphereSurface\n");
- offsetx = floor ((img->HeaderN.x-1)/2); /* PIXEL */
- offsety = floor ((img->HeaderN.y-1)/2); /* PIXEL */
- offsetz = floor ((img->HeaderN.z-1)/2); /* PIXEL */
count = 0;
sumForSphereSurface = 0;
- rmax = MIN3(offsetx*img->HeaderLength.x,
- offsety*img->HeaderLength.y,
- offsetz*img->HeaderLength.z);
-
+ offsetx = floor ((img->HeaderN.x-1)/2); /* PIXEL */
+ offsety = floor ((img->HeaderN.y-1)/2); /* PIXEL */
+
+ if(1<img->HeaderN.z) {
+ offsetz = floor ((img->HeaderN.z-1)/2); /* PIXEL */
+ rmax = MIN3(offsetx*img->HeaderLength.x,
+ offsety*img->HeaderLength.y,
+ offsetz*img->HeaderLength.z);
+ } else {
+ offsetz = 0;
+ rmax = MIN(offsetx*img->HeaderLength.x,
+ offsety*img->HeaderLength.y);
+ }
if (info->flagrmax){
if(info->rmax<rmax){
fprintf(stderr,"rmax from option:%f\n",info->rmax);
info->sdOfEdge = sqrt(sumForSd/count - SQR(info->meanOfEdge));
info->seOfEdge = info->sdOfEdge/sqrt(count);
} else {
- fprintf(stderr, "pixel number %d for edge\n", count);
+ fprintf(stderr, "pixel number %f for edge\n", count);
exit(EXIT_FAILURE);
}
break;
info->sdOfEdge = sqrt(sumForSd/count - SQR(info->meanOfEdge));
info->seOfEdge = info->sdOfEdge/sqrt(count);
} else {
- fprintf(stderr, "pixel number %d for edge\n", count);
+ fprintf(stderr, "pixel number %f for edge\n", count);
exit(EXIT_FAILURE);
}
break;
-/bin/bash: mrcImageInfo: そのようなファイルやディレクトリはありません
+<HTML>
+<HEAD>
+<TITLE>mrcImageInfo</TITLE>
+</HEAD>
+<BODY>
+<H1>mrcImageInfo</H1>
+<H2>Usage</H2>
+<PRE>
+Usage: mrcImageInfo
+Options:
+ [-i[nput] In (NULL ).as(inFile::mrcImage ) ] :Essential :Input
+ [-o[utput] Out (stdout ).as(outFile::ASCII ) ] :Optional :Output
+ [-X[Section] X (0.0 ).as(Real ) ] :Optional :X Section at y=X(value)
+ [-Y[Section] Y (0.0 ).as(Real ) ] :Optional :Y Section at y=Y(value)
+ [-x[Projection] ] :Optional :x Projection[OutPut ASCII FILE]
+ [-y[Projection] ] :Optional :y Projection[OutPut ASCII FILE]
+ [-H[istgram] H (256 ).as(Integer ) ] :Optional :Histgram: Class Number: Normalize
+ [-H[istgram]2 H2 (1 ).as(Real ) ] :Optional :Histgram2: Class Width: Integer: Normalize
+ [-H[istgram]3 H3 (256 ).as(Integer ) ] :Optional :Histgram: Class Number: Count
+ [-H[istgram]4 H4 (1 ).as(Real ) ] :Optional :Histgram2: Class Width: Integer: Count
+ [-H[istgram]5 H5 (1 ).as(Integer ) ] :Optional :Histgram5: Class Number: Count
+ [-H[istgram]6 H6 (1 ).as(Real ) ] :Optional :Histgram6: Class Width: Integer: Count
+ [-Low Low (1 ).as(Real ) ] :Optional :Low : thres
+ [-High High (1 ).as(Real ) ] :Optional :High: thres
+ [-I[nfo] ] :Optional :Information
+ [-c[onfig] configFile (NULL ).as(inFile ) ] :Optional :ConfigurationFile
+ [-rmax rmax (250 ).as(Real ) ] :Optional :rmax [A]
+ [-m[ode] mode (0 ).as(Integer ) ] :Optional :Mode
+----- Additional Usage -----
+The -I option uses -rmax or -m options
+-m 0: mean of all density
+ 1: mean of sphere surface density at rmax
+ 2: mean of cylinder surface density at rmax
+ 3: mean of edge pixels
+ 4: mean of Centre
+ 5: mean of Sparse
+ 6: mean of 2DEdge
+ 7: RMS of all pixels
+</PRE>
+</BODY>
+</HTML>