SpatialAnalysis/mainMfEBA.cpp at master · lsaravia/SpatialAnalysis · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
#include "mf.h"
#include "RWFile.h"


int main(int argc, char * argv[])
{
	RWFile file;
	simplmat <double> data;
	simplmat <double> q;


	if( argc < 6)
	{
		cerr << "Usage: mfEBA inputFile qFile minBox maxBox deltaBox option" << endl
        	 << "       option N: Not normalize measure" << endl
        	 << "              S: Normalize measure" << endl
        	 << "              A: Normalize measure and save" << endl;

		exit(1);
	}
   	string fname = argv[1];

	if( fname.find(".rst") != string::npos )
    {
		if(!file.ReadIdrisi(argv[1], data))
			exit(1);
    }
	else if( fname.find(".tif")!=string::npos )
	{
		if(!file.ReadTiff(argv[1], data))
			exit(1);
	}
	else
    {
		if(!file.ReadSeed(argv[1], data))
			exit(1);
    }


	if(!file.ReadSeed(argv[2], q))
		exit(1);

	int minBox = atoi(argv[3]);
	int maxBox = atoi(argv[4]);
	int deltaBox = atoi(argv[5]);

    char norm=argv[6][0];
    if(norm=='A') norm='S';

	MultifractalEBA(data, q,argv[1] ,minBox, maxBox, deltaBox, norm);

    if(argv[6][0]=='A')
    {
	    string::size_type pos=0;
		if( (pos=fname.find(".")) == string::npos )
	    {
			fname += "Norm.sed";
	    }
		else
		{
	    	fname = fname.substr(0,pos) + "Norm.sed";
		}

    	if(!file.WriteSeed(fname.c_str(),data))
			exit(1);
    }
}