apps/pnm2vti/pnm2vti.cpp

Go to the documentation of this file.

#include "vtkImageData.h"
#include "vtkPNMReader.h"
#include "vtkXMLImageDataWriter.h"

#include "DGVAliases.h"

void processArgs(int argcount, char *argvars[], int &binary, string &prefix, int &startSlice, int &endSlice);

int main(int argc, char *argv[])
{
    int start = 0, end = 10, tmp = 0;
    string filePrefix, tmpStr1, tmpStr2;
    bool writeBinary = false, isPrefix = true;

    cout << ">| Portable aNy Map (PNM) Format to VTK Image Format (VTI)" << endl;
    cout << ">| PNM formats have the extension *.pbm *.pgm *.pnm *.ppm" << endl;
    cout << ">| It may read in a single file or a stack like: <FileName>.pgm.[Number]" << endl;
    cout << ">| Copyright Shekhar Chandra, 2007-2009." << endl;

        processArgs(argc,argv,tmp,filePrefix,start,end);

        if(start == 0 && end == 0)
        isPrefix = false;
    if(tmp == 1)
        writeBinary = true;

    cout << ">| File to be read";
    if(isPrefix)
    {
        cout << " with prefix: " << filePrefix << endl;
        cout << ">| Slices to be read: " << start << "-" << end << endl;
    }
    else
        cout << " with name: " << filePrefix << endl;
    cout << ">| VTI file written to : " << filePrefix+".vti";
    if(writeBinary)
        cout << " (Binary)" << endl;
    else
        cout << " (ASCII)" << endl;
    cout << ">| You may wish to rename this" << endl;

    vtkPNMReader* reader = vtkPNMReader::New();
    tmpStr1 = filePrefix + "." + numToString(start);
    tmpStr2 = filePrefix + "." + numToString(end);
    if( (reader->CanReadFile(tmpStr1.c_str())
        && reader->CanReadFile(tmpStr2.c_str()))
        || reader->CanReadFile(filePrefix.c_str()) )
    {
        reader->SetDebug(true);
        if(isPrefix)
        {
            int bounds[6] = {0};

            bounds[4] = start;
            bounds[5] = end;

            reader->SetDataExtent(bounds);
            reader->SetFilePrefix(filePrefix.c_str());
            reader->Update();
        }
        else
        {
            reader->SetFileName(filePrefix.c_str());
            reader->Update();
        }

        vtkImageData *data = reader->GetOutput();

        cout << ">| Writing VTI File ";
        vtkXMLImageDataWriter *XMLwriter = vtkXMLImageDataWriter::New();
            filePrefix += ".vti";
            XMLwriter->SetFileName(filePrefix.c_str());
            if(writeBinary)
                XMLwriter->SetDataModeToBinary();
            else
                XMLwriter->SetDataModeToAscii();
            XMLwriter->SetInput(data);
            XMLwriter->Write();
        cout << "Done." << endl;

        XMLwriter->Delete();
    }
    reader->Delete();

    return EXIT_SUCCESS;
}

void processArgs(int argcount, char *argvars[], int &binary, string &prefix, int &startSlice, int &endSlice)
{
    switch(argcount)
    {
        case 2:
            prefix = argvars[1];
            binary = 0;
            startSlice = 0;
            endSlice = 0;
            break;
        case 3:
            prefix = argvars[1];
            binary = atoi(argvars[2]);
            startSlice = 0;
            endSlice = 0;
            break;
        case 5:
            prefix = argvars[1];
            binary = atoi(argvars[2]);
            startSlice = atoi(argvars[3]);
            endSlice = atoi(argvars[4]);
            break;
        default:
            printf(">| Usage: %s <File Name/Prefix> [Write Binary?] [Start Bound] [End Bound]\n",argvars[0]);
            printf(">| <File Name/Prefix> is required and the rest are optional.\n");
            printf(">| [Write Binary] should be a Boolean.");
            exit(EXIT_FAILURE);
    }
}