Layers Class Reference

Define layer objects at the regional level. More...

#include <Layers.h>

Inheritance diagram for Layers:

Collaboration diagram for Layers:

Public Member Functions
	Layers (ThreadManager *MTHREAD_h, string name_h, string label_h, bool isInteger_h, bool dynamicContent_h, string fullFilename_h, bool display_h=true)
	In the constructor we set the main layer properties.
	~Layers ()
void	addLegendItem (int ID_h, string label_h, int rColor_h, int gColor_h, int bColor_h, double minValue_h, double maxValue_h)
	Add a legend item.
void	addLegendItems (vector< LegendItems > legendItems_h)
vector< LegendItems >	getLegendItems ()
QColor	getColor (double ID_h)
	Evaluates all the legend items to find the one that match the input code, and return its color as a QColor.
string	getCategory (double ID_h)
	Evaluates all the legend items to find the one that match the input code, and return its label.
double	filterExogenousDataset (double code_h)
	Used to reclassify the land use map for "generic" categories.
void	countMyPixels (bool debug=false)
	Count the pixels going to each legend item and print them if debug==true.
void	randomShuffle ()
	For some sensitivity analisys, random the values for this layer for not-empty values (only integer layers)
bool	getIsInteger ()
	Return if the layer is integer or not (If integer on each legend item: minValue==maxValue==ID)
void	print ()
	Print the layer content as an ASCII grid map with its companion files (classification and colors). It always print the whole region, even when subregion is actived.
void	printBinMap ()
	Print a binary reppresentation of the data (a standard image, e.g. a .png file). It prints only the subregion if this is active.
string	getName () const
string	getFilename ()
	Return the filename of the associated dataset.
bool	getDynamicContent ()
	Return true if the content may change during simulation period.
bool	getDisplay ()
Public Member Functions inherited from BaseClass
	BaseClass ()
	~BaseClass ()
void	msgOut (const int &msgCode_h, const string &msg_h, const bool &refreshGUI_h=true) const
	Overloaded function to print the output log.
void	msgOut (const int &msgCode_h, const int &msg_h, const bool &refreshGUI_h=true) const
	Overloaded function to print the output log.
void	msgOut (const int &msgCode_h, const double &msg_h, const bool &refreshGUI_h=true) const
	Overloaded function to print the output log.
int	s2i (const string &string_h) const
	string to integer conversion
double	s2d (const string &string_h) const
	string to double conversion
double	s2d (const string &string_h, const bool &replaceComma) const
	string to double conversion
bool	s2b (const string &string_h) const
	string to bool conversion
string	i2s (const int &int_h) const
	integer to string conversion
string	d2s (const double &double_h) const
	double to string conversion
string	b2s (const bool &bool_h) const
	bool to string conversion
vector< int >	s2i (const vector< string > &string_h) const
	string to integer conversion (vector)
vector< double >	s2d (const vector< string > &string_h, const bool &replaceComma=false) const
	string to double conversion (vector)
vector< bool >	s2b (const vector< string > &string_h) const
	string to bool conversion (vector)
vector< string >	i2s (const vector< int > &int_h) const
	integer to string conversion (vector)
vector< string >	d2s (const vector< double > &double_h) const
	double to string conversion (vector)
vector< string >	b2s (const vector< bool > &bool_h) const
	bool to string conversion (vector)
int	getType (const string &type_h) const
	Return a type according to enum TYPE_* from a string (eg: "string" -> TYPE_STRING (2))
void	refreshGUI () const
	Ping to periodically return the control to the GUI.
template<typename T >
string	toString (const T &x) const
template<typename T >
T	stringTo (const std::string &s) const
int	vSum (const vector< int > &vector_h) const
double	vSum (const vector< double > &vector_h) const
int	vSum (const vector< vector< int > > &vector_h) const
double	vSum (const vector< vector< double > > &vector_h) const
void	tokenize (const string &str, vector< string > &tokens, const string &delimiter=" ") const
	Tokenize a string using a delimiter (default is space)
void	untokenize (string &str, vector< string > &tokens, const string &delimiter=" ") const
template<typename K , typename V >
V	findMap (const map< K, V > &mymap, const K &key, const int &error_level=MSG_CRITICAL_ERROR, const V &notFoundValue=numeric_limits< V >::min()) const
	Lookup a map for a value. Return the value starting from the key.
template<typename K , typename V >
void	changeMapValue (map< K, V > &mymap, const K &key, const V &value, const int &error_level=MSG_CRITICAL_ERROR)
	Change the value stored in a map given the key and the new value.
template<typename K , typename V >
void	incrMapValue (map< K, V > &mymap, const K &key, const V &value, const int &error_level=MSG_CRITICAL_ERROR)
	Increments a value stored in a map of the specified value, given the key.
template<typename K , typename V >
void	incrOrAddMapValue (map< K, V > &mymap, const K &key, const V &value)
	Increments a value stored in a map of the specified value, given the key.
template<typename K , typename V >
void	resetMapValues (map< K, V > &mymap, const V &value)
	Reset all values stored in a map to the specified one.
template<typename K , typename V >
map< K, V >	vectorToMap (const vector< K > &keys, const V &value=0.0)
	Returns a map built using the given vector and the given (scalar) value as keys/values pairs.
template<typename T >
vector< T >	positionsToContent (const vector< T > &vector_h, const vector< int > &positions)
	Return a vector of content from a vector and a vector of positions (int)
template<typename V >
void	debugMap (const map< iisskey, V > &mymap)
	Debug a map.
template<typename K , typename V >
void	debugMap (const map< K, V > &mymap, const K &key)
template<typename K >
int	getMaxPos (const vector< K > &v)
	Returns the position of the maximum element in the vector (the last one in case of multiple equivalent maxima)
template<typename K >
int	getMinPos (const vector< K > &v)
	Returns the position of the minimum element in the vector (the first one in case of multiple equivalent minima)
template<typename K >
K	getMax (const vector< K > &v)
	Returns the value of the maximum element in the vector (the last one in case of multiple equivalent maxima)
template<typename K >
K	getMin (const vector< K > &v)
	Returns the value of the minimum element in the vector (the first one in case of multiple equivalent minima)
template<typename K >
double	getAvg (const vector< K > &v)
	Returns the average of the elements in the vector.
template<typename K >
double	getSd (const vector< K > &v, bool sample=true)
template<typename K >
int	getPos (const K &element, const vector< K > &v, const int &msgCode_h=MSG_CRITICAL_ERROR)
template<typename K >
bool	inVector (const K &element, const vector< K > &v)
double	normSample (const double &avg, const double &stdev, const double &minval=NULL, const double &maxval=NULL) const
	Sample from a normal distribution with bounds. Slower (double time, but still you see the diff only after milion of loops).
template<typename K >
K	normSample (normal_distribution< K > &d, std::mt19937 &gen, const K &minval=NULL, const K &maxval=NULL) const
	Sample from a normal distribution with bounds. Faster (half time) as the normal_distribution is made only once.
template<typename T >
std::string	toString (const T &x) const

Private Attributes
string	name
	ID of the layer (no spaces allowed)
string	label
	Label of the layer (spaces allowed)
bool	isInteger
	Type of the layer (true==integer, false==double. If true, on each legend item: minValue==maxValue==ID)
bool	dynamicContent
	True if the content may change during simulation year.
bool	display
	Normally true, but some layers used to just keep data shoudn't be normally processed.
string	fullFileName
	Filename of the associated dataset (map)
vector< LegendItems >	legendItems
	Vector of legend items.
vector< ReclassRules >	reclassRulesVector
	Vector of initial reclassification rules.

Additional Inherited Members
Protected Attributes inherited from BaseClass
ThreadManager *	MTHREAD
	Pointer to the Thread manager.

Detailed Description

Define layer objects at the regional level.

Layer class (setting, legend...)
This class define layer objects, including:

• a set of layer proprieties (name(ID), label, associated dataset, typology (integer or double)
• a vector of legend items, associating one color to each value or interval
• a vector of reclassification rule, when we need to work with a level of depth different of those coming with the dataset

  Author

  Antonello Lobianco anton.nosp@m.ello.nosp@m.@regm.nosp@m.as.o.nosp@m.rg

Definition at line 49 of file Layers.h.

Constructor & Destructor Documentation

Layers()

Layers	(	ThreadManager *	MTHREAD_h,
		string	name_h,
		string	label_h,
		bool	isInteger_h,
		bool	dynamicContent_h,
		string	fullFilename_h,
		bool	display_h = true
	)

In the constructor we set the main layer properties.

Definition at line 32 of file Layers.cpp.

{
  MTHREAD=MTHREAD_h;
  name = name_h;
  label = label_h;
  isInteger = isInteger_h;
  dynamicContent = dynamicContent_h;
  fullFileName = fullFilename_h;
  display = display_h;
}

~Layers()

~Layers

(

)

Definition at line 43 of file Layers.cpp.

{
}

Member Function Documentation

addLegendItem()

void addLegendItem	(	int	ID_h,
		string	label_h,
		int	rColor_h,
		int	gColor_h,
		int	bColor_h,
		double	minValue_h,
		double	maxValue_h
	)

Add a legend item.

See also

LegendItems

Definition at line 48 of file Layers.cpp.

                                                                                                                             {
 
  for (uint i=0;i<legendItems.size();i++){
    if (legendItems.at(i).ID == ID_h){
      msgOut(MSG_ERROR, "Trying to add a legend item that already exist on this layer (layer: "+label+" - legend label: "+label_h+")");
      //cout << "ID: "<<ID_h<<" Label: "<<label_h<<" minValue: "<<minValue_h << " maxValue: "<<maxValue_h<<endl;
      return;
    }
  }
 
  LegendItems ITEM;
  ITEM.ID = ID_h;
  ITEM.label = label_h;
  ITEM.rColor = rColor_h;
  ITEM.gColor = gColor_h;
  ITEM.bColor = bColor_h;
  ITEM.minValue = minValue_h;
  ITEM.maxValue = maxValue_h;
  ITEM.cashedCount=0;
  legendItems.push_back(ITEM);  
 
}

Here is the call graph for this function:

addLegendItems()

void addLegendItems

(

vector< LegendItems >

legendItems_h

)

Definition at line 72 of file Layers.cpp.

                                                       {
  vector <LegendItems> toAdd;
  for(uint i=0; i<legendItems_h.size();i++){
    bool existing = false;
    for (uint j=0;j<legendItems.size();j++){
      if(legendItems_h[i].ID == legendItems[j].ID){
        existing = true;
        break;
      }
    }
    if(existing){
      msgOut(MSG_WARNING, "Legend item "+i2s(legendItems_h[i].ID)+" non added on layer "+this->name+" as already existing.");
    } else {
      toAdd.push_back(legendItems_h[i]);
    }
  }
  legendItems.insert( legendItems.end(), toAdd.begin(), toAdd.end() );
}

Referenced by Gis::applyForestReclassification().

Here is the call graph for this function:

Here is the caller graph for this function:

countMyPixels()

void countMyPixels

(

bool

debug = false

)

Count the pixels going to each legend item and print them if debug==true.

Definition at line 188 of file Layers.cpp.

                               {
 
  for (uint i=0; i<legendItems.size(); i++){
    legendItems.at(i).cashedCount=0; //initialized with 0 values...
  }
  double totPixels = MTHREAD->GIS->getXyNPixels();
  double pixelValue;
  for (uint j=0;j<totPixels;j++){
    pixelValue = MTHREAD->GIS->getPixel(j)->getDoubleValue(name);
    if (isInteger){
      for(uint i=0; i<legendItems.size(); i++){
        if (legendItems.at(i).ID == ((int)pixelValue)){
          legendItems.at(i).cashedCount++;
          break;
        }
      }
    }
    else {
      for(uint i=0; i<legendItems.size(); i++){
        if (pixelValue < legendItems.at(i).maxValue &&  pixelValue >= legendItems.at(i).minValue){
          legendItems.at(i).cashedCount++;
          break;
        }
      }
    }
  }
  if (debug){
    msgOut(MSG_INFO, "Layer statistics - Count by Legend items");
    msgOut(MSG_INFO, "Layer name: "+label);
    msgOut(MSG_INFO, "Total plots: "+ d2s(totPixels));
    for(uint i=0;i<legendItems.size();i++){
      msgOut(MSG_INFO, legendItems.at(i).label+": "+i2s(legendItems.at(i).cashedCount));
    }
  }
}

Here is the call graph for this function:

filterExogenousDataset()

double filterExogenousDataset

(

double

code_h

)

Used to reclassify the land use map for "generic" categories.

Used in the init stage, this function take as input the real map code as just read from the map file, and filter it according to the reclassification rules.

See also

ReclassRules

Definition at line 97 of file Layers.cpp.

                                           {
  bool check =false;
  std::vector <double> cumPVector;
  std::vector <double> outCodesVector;
  double cumP = 0;
  double returnCode=0;
  
  for(uint i=0; i<reclassRulesVector.size(); i++){
    if (reclassRulesVector.at(i).inCode == code_h){
      check = true;
      cumP += reclassRulesVector.at(i).p;
      cumPVector.push_back(cumP);
      outCodesVector.push_back(reclassRulesVector.at(i).outCode);
    }
  }
  if (!check) {return code_h;}
  if (cumP <= 0.99999999 || cumP >= 1.00000001){msgOut(MSG_CRITICAL_ERROR,"the sum of land use reclassification rules is not 1 for at least one input code (input code:  "+d2s(code_h)+"; cumP: "+d2s(cumP)+")");}
  double random;
  //srand(time(NULL)); // this would re-initialise the random seed
  random = ((double)rand() / ((double)(RAND_MAX)+(double)(1)) );
  for(uint i=0; i<cumPVector.size(); i++){
    if (random <= cumPVector.at(i)){
      returnCode = outCodesVector.at(i);
      break;  
    }
  }  
  return returnCode;
}

Here is the call graph for this function:

getCategory()

string getCategory

(

double

ID_h

)

Evaluates all the legend items to find the one that match the input code, and return its label.

This function take as input the value stored in the pixel for the specific layer, loops over the legend item and find the one that match it, returning its label.
If the layer is of type integer, the match is agains legendItem IDs, otherwise we compare the legendItem ranges.

See also

LegendItems

Definition at line 162 of file Layers.cpp.

                              {
  if (ID_h == MTHREAD->GIS->getNoValue()){
    return "";
  }
  if (isInteger){
    for(uint i=0; i<legendItems.size(); i++){
      if (legendItems.at(i).ID == ((int)ID_h)){
        return legendItems.at(i).label;
      }
    }
    return "";
  }
  else {
    for(uint i=0; i<legendItems.size(); i++){
      if (ID_h < legendItems.at(i).maxValue &&  ID_h >= legendItems.at(i).minValue){
        return legendItems.at(i).label;
      }
    }
    return "";  
  }
}

Here is the call graph for this function:

getColor()

QColor getColor

(

double

ID_h

)

Evaluates all the legend items to find the one that match the input code, and return its color as a QColor.

This function take as input the value stored in the pixel for the specific layer, loops over the legend item and find the one that match it, returning its color.
If the layer is of type integer, the match is agains legendItem IDs, otherwise we compare the legendItem ranges.

See also

LegendItems

Definition at line 132 of file Layers.cpp.

                           {
  QColor nocolor(255,255,255);
  if (ID_h == MTHREAD->GIS->getNoValue()){
    return nocolor;
  }
  if (isInteger){
    for(uint i=0; i<legendItems.size(); i++){
      if (legendItems.at(i).ID == ((int)ID_h)){
        QColor color(legendItems.at(i).rColor, legendItems.at(i).gColor, legendItems.at(i).bColor);
        return color;
      }
    }
    return nocolor;
  }
  else {
    for(uint i=0; i<legendItems.size(); i++){
      if (ID_h < legendItems.at(i).maxValue &&  ID_h >= legendItems.at(i).minValue){
        QColor color(legendItems.at(i).rColor, legendItems.at(i).gColor, legendItems.at(i).bColor);
        return color;
      }
    }
    return nocolor;  
  }
}

Referenced by printBinMap().

Here is the call graph for this function:

Here is the caller graph for this function:

getDisplay()

bool getDisplay ( )

inline

Definition at line 94 of file Layers.h.

{return display;}

getDynamicContent()

bool getDynamicContent ( )

inline

Return true if the content may change during simulation period.

Definition at line 93 of file Layers.h.

{return dynamicContent;}

getFilename()

string getFilename ( )

inline

Return the filename of the associated dataset.

Definition at line 91 of file Layers.h.

{return fullFileName;}

getIsInteger()

bool getIsInteger ( )

inline

Return if the layer is integer or not (If integer on each legend item: minValue==maxValue==ID)

Definition at line 83 of file Layers.h.

{return isInteger;}

getLegendItems()

vector< LegendItems > getLegendItems ( )

inline

Definition at line 70 of file Layers.h.

{return legendItems;};

Referenced by Gis::applyForestReclassification().

Here is the caller graph for this function:

getName()

string getName ( ) const

inline

Definition at line 89 of file Layers.h.

{return name;}

print()

void print

(

)

Print the layer content as an ASCII grid map with its companion files (classification and colors). It always print the whole region, even when subregion is actived.

Definition at line 251 of file Layers.cpp.

             {
 
  if(MTHREAD->MD->getIntSetting("outputLevel",DATA_NOW)<OUTVL_MAPS) return;
  //if(!display || !dynamicContent) return; // already checked in Gis::printLayers() and allowed only for the first year
  string mapBaseDirectory = MTHREAD->MD->getBaseDirectory()+MTHREAD->MD->getOutputDirectory()+"maps/";
  string mapGridOutputDirectory = mapBaseDirectory+"asciiGrids/";
  string catsOutputDirectory = mapBaseDirectory+"cats/";
  string coloursOutputDirectory = mapBaseDirectory+"colr/";
 
  string mapFilename      = mapGridOutputDirectory +name+ "_" +i2s(MTHREAD->SCD->getYear()) +"_" +MTHREAD->getScenarioName();
  string catsFilename     = catsOutputDirectory    +name+ "_" +i2s(MTHREAD->SCD->getYear()) +"_" +MTHREAD->getScenarioName();
  string coloursFilename  = coloursOutputDirectory +name+ "_" +i2s(MTHREAD->SCD->getYear()) +"_" +MTHREAD->getScenarioName();
    string filenameListIntLayers = mapBaseDirectory+"integerListLayers/"+MTHREAD->getScenarioName();
    string filenameListFloatLayers = mapBaseDirectory+"floatListLayers/"+MTHREAD->getScenarioName();
 
  // printing the map...
  string header;
  if(MTHREAD->MD->getIntSetting("mapOutputFormat",DATA_NOW) == 1){ // GRASS ASCII Grid
    header =      "north: " + d2s(MTHREAD->GIS->getGeoTopY()) + "\n"
          + "south: " + d2s(MTHREAD->GIS->getGeoBottomY()) + "\n"
          + "east: " + d2s(MTHREAD->GIS->getGeoRightX()) + "\n"
          + "west: " + d2s(MTHREAD->GIS->getGeoLeftX()) + "\n"
          + "rows: " + i2s(MTHREAD->GIS->getYNPixels()) + "\n"
          + "cols: " + i2s(MTHREAD->GIS->getXNPixels()) + "\n"
          + "null: " + d2s(MTHREAD->GIS->getNoValue()) + "\n";
 
  } else if(MTHREAD->MD->getIntSetting("mapOutputFormat",DATA_NOW) == 2){
    header =      "ncols: " + i2s(MTHREAD->GIS->getXNPixels()) + "\n"
          + "lrows: " + i2s(MTHREAD->GIS->getYNPixels()) + "\n"
          + "xllcornel: " + d2s(MTHREAD->GIS->getGeoLeftX()) + "\n"
          + "yllcorner: " + d2s(MTHREAD->GIS->getGeoBottomY()) + "\n"
          + "cellsize: "  + d2s(MTHREAD->GIS->getXMetersByPixel()) + "\n"
          + "nodata_value: " + d2s(MTHREAD->GIS->getNoValue()) + "\n";
    if(MTHREAD->GIS->getXMetersByPixel() != MTHREAD->GIS->getYMetersByPixel()){
      msgOut(MSG_ERROR, "The X resolution is different to the Y resolution. I am exporting the map in ArcInfo ASCII Grid format using the X resolution, but be aware that it is incorrect, as this format doesn't support different X-Y resolutions.");
    }
  
  } else {
    msgOut(MSG_ERROR,"Map not print for unknow output type.");
  }
 
  ofstream outm; //out map
  outm.open(mapFilename.c_str(), ios::out); //ios::app to append..
  if (!outm){ msgOut(MSG_ERROR,"Error in opening the file "+mapFilename+".");}
  outm << header << "\n";
 
  for (int i=0;i<MTHREAD->GIS->getYNPixels();i++){
    for (int j=0;j<MTHREAD->GIS->getXNPixels();j++){
      outm << MTHREAD->GIS->getPixel(j, i)->getDoubleValue(name) << " ";
    }
    outm << "\n";
  }
  outm.close();
 
  //printing the cat file
  ofstream outc; //out category file
  outc.open(catsFilename.c_str(), ios::out); //ios::app to append..
  if (!outc){ msgOut(MSG_ERROR,"Error in opening the file "+catsFilename+".");}
  outc << "# " << name  << "_-_" << i2s(MTHREAD->SCD->getYear()) << "\n\n\n";
  outc << "0.00 0.00 0.00 0.00"<<"\n";
 
  if (isInteger){
    for(uint i=0;i<legendItems.size();i++){
      outc << legendItems[i].ID << ":"<< legendItems[i].label << "\n";
    }
  }
  else {
    for(uint i=0;i<legendItems.size();i++){
      outc << legendItems[i].minValue << ":"<< legendItems[i].maxValue << ":"<< legendItems[i].label << "\n";
    }
  }
 
  //printing the colour legend file
  ofstream outcl; //out colour file
  outcl.open(coloursFilename.c_str(), ios::out); //ios::app to append..
  if (!outcl){ msgOut(MSG_ERROR,"Error in opening the file "+coloursFilename+".");}
  outcl << "% " << name  << "_-_" << i2s(MTHREAD->SCD->getYear()) << "\n\n\n";
  
  if (isInteger){
    for(uint i=0;i<legendItems.size();i++){
      outcl << legendItems[i].ID << ":"<< legendItems[i].rColor << ":" << legendItems[i].gColor << ":" << legendItems[i].bColor << "\n";
    }
  }
  else {
    for(uint i=0;i<legendItems.size();i++){
      outcl << legendItems[i].minValue << ":"<< legendItems[i].rColor << ":" << legendItems[i].gColor << ":" << legendItems[i].bColor << " "<< legendItems[i].maxValue << ":"<< legendItems[i].rColor << ":" << legendItems[i].gColor << ":" << legendItems[i].bColor << "\n";
    }
  }
 
  // adding the layer to the list of saved layers..
  ofstream outList;
  if (isInteger){
    outList.open(filenameListIntLayers.c_str(), ios::app); // append !!!
    outList << name << "_" << MTHREAD->SCD->getYear() << "_" << MTHREAD->getScenarioName() << "\n";
  }
  else {
    outList.open(filenameListFloatLayers.c_str(), ios::app); // append !!!
    outList << name << "_" << MTHREAD->SCD->getYear() << "_" << MTHREAD->getScenarioName() << "\n";
  }
  outList.close();
}

Here is the call graph for this function:

printBinMap()

void printBinMap

(

)

Print a binary reppresentation of the data (a standard image, e.g. a .png file). It prints only the subregion if this is active.

Definition at line 354 of file Layers.cpp.

                   {
 
  if(!display || !dynamicContent) return;
 
  int xNPixels            = MTHREAD->GIS->getXNPixels();
  int subXR               = MTHREAD->GIS->getSubXR();
  int subXL               = MTHREAD->GIS->getSubXL();
  int subYT               = MTHREAD->GIS->getSubYT();
  int subYB               = MTHREAD->GIS->getSubYB();
 
  string mapBaseDirectory = MTHREAD->MD->getBaseDirectory()+MTHREAD->MD->getOutputDirectory()+"maps/bitmaps/";
  string mapFilename      = mapBaseDirectory +name+ "_" +i2s(MTHREAD->SCD->getYear()) +"_" +MTHREAD->getScenarioName()+".png";
 
  QImage image = QImage(subXR-subXL+1, subYB-subYT+1, QImage::Format_RGB32);
  image.fill(qRgb(255, 255, 255));
  for (int countRow=subYT;countRow<subYB;countRow++){
    for (int countColumn=subXL;countColumn<subXR;countColumn++){
      double value = MTHREAD->GIS->getPixel(countRow*xNPixels+countColumn)->getDoubleValue(name);
      QColor color = this->getColor(value);
      image.setPixel(countColumn-subXL,countRow-subYT,color.rgb());
    }
  }
  image.save(mapFilename.c_str());
}

Here is the call graph for this function:

randomShuffle()

void randomShuffle

(

)

For some sensitivity analisys, random the values for this layer for not-empty values (only integer layers)

Definition at line 224 of file Layers.cpp.

                     {
 
 
  vector <double> origValues;
  int maskValue = -MTHREAD->GIS->getNoValue();
  double totPixels = MTHREAD->GIS->getXyNPixels();
  for (uint i=0;i<totPixels;i++){
    double pxValue= MTHREAD->GIS->getPixel(i)->getDoubleValue(name);
    if(pxValue != MTHREAD->GIS->getNoValue()){
      origValues.push_back(pxValue);
      MTHREAD->GIS->getPixel(i)->changeValue(name,maskValue);
    }
  }
  random_shuffle(origValues.begin(), origValues.end()); // randomize the elements of the array.
 
  for (uint i=0;i<totPixels;i++){
    double pxValue= MTHREAD->GIS->getPixel(i)->getDoubleValue(name);
    if(pxValue != MTHREAD->GIS->getNoValue()){
      double toChangeValue = origValues.at(origValues.size()-1);
      //cout << toChangeValue << endl;
      origValues.pop_back();
      MTHREAD->GIS->getPixel(i)->changeValue(name,toChangeValue);
    }
  }
 
}

Here is the call graph for this function:

Member Data Documentation

display

bool display

private

Normally true, but some layers used to just keep data shoudn't be normally processed.

Definition at line 102 of file Layers.h.

Referenced by getDisplay(), Layers(), and printBinMap().

dynamicContent

bool dynamicContent

private

True if the content may change during simulation year.

Definition at line 101 of file Layers.h.

Referenced by getDynamicContent(), Layers(), and printBinMap().

fullFileName

string fullFileName

private

Filename of the associated dataset (map)

Definition at line 103 of file Layers.h.

Referenced by getFilename(), and Layers().

isInteger

bool isInteger

private

Type of the layer (true==integer, false==double. If true, on each legend item: minValue==maxValue==ID)

Definition at line 100 of file Layers.h.

Referenced by countMyPixels(), getCategory(), getColor(), getIsInteger(), Layers(), and print().

label

string label

private

Label of the layer (spaces allowed)

Definition at line 99 of file Layers.h.

Referenced by addLegendItem(), countMyPixels(), and Layers().

legendItems

vector<LegendItems> legendItems

private

Vector of legend items.

See also

LegendItems

Definition at line 104 of file Layers.h.

Referenced by addLegendItem(), addLegendItems(), countMyPixels(), getCategory(), getColor(), getLegendItems(), and print().

name

string name

private

ID of the layer (no spaces allowed)

Definition at line 98 of file Layers.h.

Referenced by addLegendItems(), countMyPixels(), getName(), Layers(), print(), printBinMap(), and randomShuffle().

reclassRulesVector

vector<ReclassRules> reclassRulesVector

private

Vector of initial reclassification rules.

See also

ReclassRules

Definition at line 105 of file Layers.h.

Referenced by filterExogenousDataset().

---

The documentation for this class was generated from the following files:

• /exports/ffsm/src_latest/src/Layers.h
• /exports/ffsm/src_latest/src/Layers.cpp