// bbSetOutputOut( bbGetInputIn() );
// std::cout << "Output value = " <<bbGetOutputOut() << std::endl;
-
+printf("EED ReadMHDPlane::Process Start\n");
vtkImageData *newImage;
long int newHeaderSize;
char mystring[250];
char strTmp[20];
bool ok=true;
- FILE *ffIn = fopen(bbGetInputFileName().c_str(),"r+");
- if (ffIn!=NULL){
- FILE *ffOut = fopen(newFileName.c_str(),"w+");
+ FILE *ffIn = fopen(bbGetInputFileName().c_str(),"r");
+ if (ffIn!=NULL)
+ {
+printf("EED ReadMHDPlane::Process 1\n");
+ FILE *ffOut = fopen(newFileName.c_str(),"w");
+printf("EED ReadMHDPlane::Process 2 %p %s\n",ffOut, newFileName.c_str() );
while(!feof(ffIn))
{
strcpy(mystring,"\n");
fgets(mystring,250,ffIn);
+printf("EED ReadMHDPlane::Process MHD:%s\n", mystring);
if (strncmp("NDims",mystring,5)==0) {strcpy(mystring,"NDims = 2\n");}
if (strncmp("DimSize",mystring,6)==0)
{
newHeaderSize = sx*sy*bbGetInputSlice();
if (bbGetInputSlice()>=sz) {ok=false;}
} // if
+ if (strncmp("ElementType = MET_CHAR",mystring,24)==0) { newHeaderSize=newHeaderSize*1; }
+ if (strncmp("ElementType = MET_UCHAR",mystring,24)==0) { newHeaderSize=newHeaderSize*1; }
+ if (strncmp("ElementType = MET_USHORT",mystring,24)==0) { newHeaderSize=newHeaderSize*2; }
if (strncmp("ElementType = MET_USHORT",mystring,24)==0) { newHeaderSize=newHeaderSize*2; }
- if (strncmp("ElementType = MET_SHORT",mystring,23)==0) { newHeaderSize=newHeaderSize*2; }
- if (strncmp("ElementType = MET_FLOAT",mystring,23)==0) { newHeaderSize=newHeaderSize*2; }
+ if (strncmp("ElementType = MET_SHORT",mystring,23)==0) { newHeaderSize=newHeaderSize*2; }
+ if (strncmp("ElementType = MET_FLOAT",mystring,23)==0) { newHeaderSize=newHeaderSize*4; }
+ if (strncmp("ElementType = VTK_FLOAT",mystring,23)==0) { newHeaderSize=newHeaderSize*4; }
if (strncmp("Offset",mystring,6)==0) {strcpy(mystring,"Offset = 0 0 0\n");}
if (strncmp("HeaderSize",mystring,10)==0) {strcpy(mystring,"");}
if (strncmp("ElementDataFile",mystring,15)==0)
{
fprintf(ffOut,"HeaderSize = %ld\n\n", newHeaderSize );
} // if
+printf("EED ReadMHDPlane::Process 3\n");
fprintf(ffOut,mystring);
+printf("EED ReadMHDPlane::Process 4\n");
} // while
fclose(ffIn);
fclose(ffOut);
} // if ok
} else {
ok=false;
+ printf("EED ERROR: Problem openin:%s\n", bbGetInputFileName().c_str() );
}
if (ok==false)
{
- int i,j;
+ int i;
int sizeX, sizeY, sizeZ;
sizeX = 200;
sizeY = sizeX;
} // if ok
bbSetOutputOut( newImage );
+printf("EED ReadMHDPlane::Process End\n");
}
//=====