RasterMap Class Reference

#include <rastermap.h>

Public Member Functions
int	MapWidth (void)
int	MapHeight (void)
int	Get (int a_x, int a_y)
char *	GetID (void)
int *	GetMagicP (int a_x, int a_y)
void	Put (int a_x, int a_y, int a_elem)
	RasterMap (const char *a_mapfile, Landscape *m_landscape)
	~RasterMap (void)
void	Manipulation1 ()
int	CellReplacementNeighbour (int a_x, int a_y, int a_polyref)
	A method for helping remove tiny polygons. More...
bool	MissingCellReplace (int a_x, int a_y, bool a_fieldsonly)
	A method for removing missing polygons. More...
bool	MissingCellReplaceWrap (int a_x, int a_y, bool a_fieldsonly)
	A method for removing missing polygons - tests for edge conditions. More...

Protected Member Functions
void	Init1 (const char *a_mapfile, Landscape *m_landscape)
void	Init2 (const char *a_mapfile, Landscape *m_landscape)

Private Attributes
char	m_id [12]
int *	m_map
int	m_width
int	m_height
int	m_x
int	m_y
Landscape *	m_landscape

Detailed Description

Definition at line 45 of file rastermap.h.

Constructor & Destructor Documentation

RasterMap::RasterMap	(	const char *	a_mapfile,
		Landscape *	m_landscape
	)

Definition at line 46 of file rastermap.cpp.

References cfg_manipulation1_on, g_msg, CfgBool::value(), MapErrorMsg::Warn(), and WARN_FILE.

{
    m_landscape = a_landscape;
    FILE * IFile;

    IFile = fopen(a_mapfile, "rb");
    if (!IFile) {
        g_msg->Warn(WARN_FILE, "RasterMap::RasterMap(): Unable to open file",
            a_mapfile);
        exit(1);
    }
    fread(m_id, 1, 12, IFile);
    fclose(IFile);
    if (strcmp(m_id, "An LSB File") == 0 ) Init1(a_mapfile, a_landscape); else Init2(a_mapfile, a_landscape);
    if (cfg_manipulation1_on.value()) {
        Manipulation1();
        exit(10);
    }
}

RasterMap::~RasterMap

(

void

)

Definition at line 171 of file rastermap.cpp.

{
  free( m_map );
}

Member Function Documentation

int RasterMap::CellReplacementNeighbour	(	int	a_x,
		int	a_y,
		int	a_polyref
	)

A method for helping remove tiny polygons.

Replaces a cell value with the most common value from the surrounding 8 cells. If near the edge nothing is done. Return value is the new value of the cell or -1 if nothing is done. First tests to make sure that the polygon we are removing is at these coords.

Definition at line 216 of file rastermap.cpp.

References g_msg, and MapErrorMsg::Warn().

{
    if ((a_x < 1) || (a_x > m_width - 2) || (a_y < 1) || (a_y > m_height - 2)) return -1;
    int toreplace = Get(a_x, a_y);
    if (toreplace != a_polyref)
    {
        g_msg->Warn("RasterMap::CellReplacementNeighbour: x,y pair does not match polyref ", a_polyref);
        exit(0);
    }
    int surroundingcells[8];
    surroundingcells[0] = Get(a_x - 1, a_y - 1);
    surroundingcells[1] = Get(a_x    , a_y - 1);
    surroundingcells[2] = Get(a_x + 1, a_y - 1);
    surroundingcells[3] = Get(a_x - 1, a_y);
    surroundingcells[4] = Get(a_x + 1, a_y);
    surroundingcells[5] = Get(a_x - 1, a_y + 1);
    surroundingcells[6] = Get(a_x    , a_y + 1);
    surroundingcells[7] = Get(a_x + 1, a_y + 1);
    // Find out how many duplicates we have
    int count[8];
    for (int i = 0; i < 8; i++)
    {
        count[i] = 0;
        for (int j = 0; j < 8; j++)
        {
            if (surroundingcells[j] == surroundingcells[i]) count[i]++;
        }
    }
    int found = -1;
    int index = 0;
    for (int i = 0; i < 8; i++)
    {
        if (count[i]>found) index = i;
    }
    return surroundingcells[index];
}

int RasterMap::Get ( int  a_x, int  a_y  )

inline

Definition at line 86 of file rastermap.h.

References g_msg, MapErrorMsg::Warn(), and WARN_BUG.

Referenced by Landscape::GetActualGooseGrazingForage(), Landscape::SupplyBirdMaizeForage(), Landscape::SupplyBirdSeedForage(), Landscape::SupplyCentroidX(), Landscape::SupplyCentroidY(), Landscape::SupplyLEHigh(), Landscape::SupplyOpenness(), Landscape::SupplySoilType(), Landscape::SupplySoilTypeR(), and Pesticide::TwinMapSpray().

                                            {
    if (a_x<0 || a_x>=m_width ||
      a_y<0 || a_y>=m_height ) {
      g_msg->Warn( WARN_BUG,
          "RasterMap::Get(): Coordinates out of range!", "" );
        exit(1);
  }
  return m_map[ a_y * m_width + a_x ];
}

char* RasterMap::GetID ( void  )

inline

Definition at line 60 of file rastermap.h.

{ return m_id; }

int * RasterMap::GetMagicP ( int  a_x, int  a_y  )

inline

Definition at line 79 of file rastermap.h.

Referenced by Landscape::hb_Add().

{
  return &m_map[ a_y * m_width + a_x ];
}

void RasterMap::Init1 ( const char *  a_mapfile, Landscape *  m_landscape  )

protected

Definition at line 66 of file rastermap.cpp.

References cfg_rectangularmaps_on, g_msg, m_polymapping, CfgBool::value(), MapErrorMsg::Warn(), WARN_FATAL, and WARN_FILE.

{
    m_landscape = a_landscape;
    FILE * IFile;
    char error_num[20];
    int lastref = -9999;

    IFile = fopen(a_mapfile, "rb");
    if (!IFile) {
        g_msg->Warn(WARN_FILE, "RasterMap::RasterMap(): Unable to open file",
            a_mapfile);
        exit(1);
    }
    fread(m_id, 1, 12, IFile);
    fread(&m_width, 1, sizeof(int), IFile);
    if (!cfg_rectangularmaps_on.value())m_height = m_width; else fread(&m_height, 1, sizeof(int), IFile);

    m_map = (int*)malloc(sizeof(int)* m_width*m_height);
    if (m_map == NULL) {
        g_msg->Warn(WARN_FATAL,
            "RasterMap::RasterMap(): Out of memory.", "");
        exit(1);
    }

    fread(m_map, 1, sizeof(int)*m_width*m_height, IFile);
    m_x = m_width;
    m_y = m_height;
    for (int x = 0; x<m_x; x++) {
        for (int y = 0; y<m_y; y++) {
            int i = Get(x, y);
            if (i<0) {
                //if ( i<0 || i> 250000 ) { REMOVED 5/7/2013 due to problems with large maps overstepping this limit
                // This fills in holes caused by -9999, but if two consequetive holes are found, then it will exit
                i = lastref;
                if (i<0) {
                    //        if ( i<0 || i> 250000 ) {
                    sprintf(error_num, "%d", Get(x, y));
                    g_msg->Warn(WARN_FILE, "RasterMap::RasterMap(): Unknown polygon map number (could be a hole):", error_num);
                    exit(1);
                }
                Put(x, y, i);
            }
            lastref = i;
            if (m_polymapping[Get(x, y)] == -1)
            {
                // Tripping out here means that the lsb file has a polygon that is not
                // listed in the polygonrefs file
                sprintf(error_num, "%d", Get(x, y));
                g_msg->Warn(WARN_FILE, "RasterMap::RasterMap():" " Unknown polygon ref number:", error_num);
                exit(1);
            }
        }
    }
    fclose(IFile);
}

void RasterMap::Init2 ( const char *  a_mapfile, Landscape *  m_landscape  )

protected

Definition at line 122 of file rastermap.cpp.

References g_msg, m_polymapping, MapErrorMsg::Warn(), WARN_FATAL, and WARN_FILE.

{
    m_landscape = a_landscape;
    char error_num[20];

    ifstream IFile(a_mapfile, ios::binary);
    if (!IFile.is_open()) {
        g_msg->Warn(WARN_FILE, "RasterMap::RasterMap(): Unable to open file",
            a_mapfile);
        exit(1);
    }
    IFile.read(m_id,12);
    IFile.read((char*)&m_width, sizeof (int));
    IFile.read((char*)&m_height, sizeof (int));
    m_map = new int[ m_width*m_height];
    if (m_map == NULL) {
        g_msg->Warn(WARN_FATAL,
            "RasterMap::RasterMap(): Out of memory.", "");
        exit(1);
    }
    IFile.read((char*)m_map, m_width*m_height * sizeof(int));

    m_x = m_width;
    m_y = m_height;
    int pref = -1;
    for (int x = 0; x<m_x; x++) 
    {
        for (int y = 0; y<m_y; y++) 
        {
            pref = Get(x, y);
            if (pref<0)
            {
                sprintf(error_num, "%d", Get(x, y));
                g_msg->Warn(WARN_FILE, "RasterMap::RasterMap(): Unknown polygon map number (could be a hole):", error_num);
                exit(1);
            }
            if (m_polymapping[Get(x, y)] == -1)
            {
                // Tripping out here means that the lsb file has a polygon that is not
                // listed in the polygonrefs file
                sprintf(error_num, "%d", Get(x, y));
                g_msg->Warn(WARN_FILE, "RasterMap::RasterMap():" " Unknown polygon ref number:", error_num);
                exit(1);
            }
        }
    }
    IFile.close();
}

void RasterMap::Manipulation1

(

)

Definition at line 176 of file rastermap.cpp.

References cfg_rastermapmanipmultiply, and CfgInt::value().

                              {

    FILE* o1;
    FILE* o2;
    o1=fopen("RasterDump1.lsb","wb");
    o2=fopen("RasterDump2.lsb","wb");
    char * l_id = m_id;
    int new_height = m_height*2;
    fwrite( l_id, 1, 12, o1 );
    fwrite( & new_height, 1, sizeof( int ), o1 );
        for (int y=0; y< m_width; y++ ) {
    for ( int x = 0; x < m_width ; x++ ) {
            int l_poly = Get(x,y);
            fwrite( & l_poly, 1, sizeof( int ), o1 );
            fwrite( & l_poly, 1, sizeof( int ), o1 );
        }
    for ( int x = 0; x < m_width ; x++ ) {
            int l_poly = Get(x,y);
            fwrite( & l_poly, 1, sizeof( int ), o1 );
            fwrite( & l_poly, 1, sizeof( int ), o1 );
        }
    }
    int multiply=cfg_rastermapmanipmultiply.value();
    new_height = m_height*multiply;
    fwrite( l_id, 1, 12, o2 );
    fwrite( &new_height, 1, sizeof( int ), o2 );
    for (int m=0; m<multiply; m++) {
        for (int y=0; y< m_width; y++ ) {
            for (int mm=0; mm<multiply; mm++) {
                for ( int x = 0; x < m_width ; x++ ) {
                    int l_poly = Get(x,y);
                    fwrite( & l_poly, 1, sizeof( int ), o2 );
                }
            }
        }
    }
    fclose( o1 );
    fclose( o2 );
}

int RasterMap::MapHeight ( void  )

inline

Definition at line 58 of file rastermap.h.

Referenced by Landscape::hb_Add(), and Pesticide::Pesticide().

{ return m_height; }

int RasterMap::MapWidth ( void  )

inline

Definition at line 57 of file rastermap.h.

Referenced by Landscape::hb_Add(), Landscape::Landscape(), and Pesticide::Pesticide().

{ return m_width; }

bool RasterMap::MissingCellReplace	(	int	a_x,
		int	a_y,
		bool	a_fieldsonly
	)

A method for removing missing polygons.

Definition at line 259 of file rastermap.cpp.

{
    int surroundingcells[8];
    surroundingcells[0] = Get(a_x - 1, a_y - 1);
    surroundingcells[1] = Get(a_x, a_y - 1);
    surroundingcells[2] = Get(a_x + 1, a_y - 1);
    surroundingcells[3] = Get(a_x - 1, a_y);
    surroundingcells[4] = Get(a_x + 1, a_y);
    surroundingcells[5] = Get(a_x - 1, a_y + 1);
    surroundingcells[6] = Get(a_x, a_y + 1);
    surroundingcells[7] = Get(a_x + 1, a_y + 1);
    // If we prefer fields then check if we have any
    if (a_fieldsonly)
    {
        unsigned offset = (unsigned)(((double)rand() / RAND_MAX + 1) * 8);
        for (unsigned i = 0; i < 8; i++)
        {
            if (m_landscape->IsFieldType(m_landscape->SupplyElementType(surroundingcells[(i + offset) % 8])))
            {
                // We have a next door field, so use this
                Put(a_x, a_y, surroundingcells[(i + offset) % 8]);
                return true;
            }
        }

    }
    return false;
}

bool RasterMap::MissingCellReplaceWrap	(	int	a_x,
		int	a_y,
		bool	a_fieldsonly
	)

A method for removing missing polygons - tests for edge conditions.

Definition at line 288 of file rastermap.cpp.

References tole_Field, tole_Orchard, tole_PermanentSetaside, tole_PermPasture, tole_PermPastureLowYield, tole_PermPastureTussocky, tole_PermPastureTussockyWet, and tole_Vildtager.

{
    int surroundingcells[8];
    if (a_x >= 1 && a_y >= 1) surroundingcells[0] = Get(a_x - 1, a_y - 1); else surroundingcells[0] = Get(a_x, a_y);
    if (a_y >= 1) surroundingcells[1] = Get(a_x, a_y - 1); else surroundingcells[1] = Get(a_x, a_y);
    if (a_x < m_width-1 && a_y >= 1) surroundingcells[2] = Get(a_x + 1, a_y - 1); else surroundingcells[2] = Get(a_x, a_y);
    if (a_x >= 1) surroundingcells[3] = Get(a_x - 1, a_y); else surroundingcells[3] = Get(a_x, a_y);
    if (a_x < m_width-1) surroundingcells[4] = Get(a_x + 1, a_y); else surroundingcells[4] = Get(a_x, a_y);
    if (a_x >= 1 && a_y < m_height-1) surroundingcells[5] = Get(a_x - 1, a_y + 1); else surroundingcells[5] = Get(a_x, a_y);
    if (a_y < m_height-1) surroundingcells[6] = Get(a_x, a_y + 1); else surroundingcells[6] = Get(a_x, a_y);
    if (a_x < m_width-1 && a_y < m_height-1) surroundingcells[7] = Get(a_x + 1, a_y + 1); else surroundingcells[7] = Get(a_x, a_y);
    // If we prefer fields then check if we have any
    if (a_fieldsonly)
    {
        unsigned offset = (unsigned)(((double)rand() / RAND_MAX + 1) * 8);
        for (unsigned i = 0; i < 8; i++)
        {
            TTypesOfLandscapeElement tole = m_landscape->SupplyElementType(surroundingcells[(i + offset) % 8]);
            if ((tole == tole_Field) || (tole == tole_Orchard) || (tole == tole_PermanentSetaside) || (tole == tole_PermPasture) || (tole == tole_PermPastureLowYield)
                || (tole == tole_PermPastureTussocky) || (tole == tole_PermPastureTussockyWet) || (tole == tole_Vildtager))
            {
                // We have a next door field, so use this
                Put(a_x, a_y, surroundingcells[(i + offset) % 8]);
                return true;
            }
        }

    }
/*  else
    {
        unsigned offset = (unsigned)(((double)rand() / RAND_MAX + 1) * 8);
        for (unsigned i = 0; i < 8; i++)
        {
            TTypesOfLandscapeElement tole = m_landscape->SupplyElementType(surroundingcells[(i + offset) % 8]);
            if ((tole != tole_Missing))
            {
                // We have a next door useful polygon, so use this
                Put(a_x, a_y, surroundingcells[(i + offset) % 8]);
                return true;
            }
        }
    }
*/
    return false;
}

void RasterMap::Put ( int  a_x, int  a_y, int  a_elem  )

inline

Definition at line 62 of file rastermap.h.

{ m_map[ a_y * m_width + a_x ] = a_elem; }

Member Data Documentation

int RasterMap::m_height

private

Definition at line 50 of file rastermap.h.

char RasterMap::m_id[12]

private

Definition at line 47 of file rastermap.h.

Landscape* RasterMap::m_landscape

private

Definition at line 53 of file rastermap.h.

int* RasterMap::m_map

private

Definition at line 48 of file rastermap.h.

int RasterMap::m_width

private

Definition at line 49 of file rastermap.h.

int RasterMap::m_x

private

Definition at line 51 of file rastermap.h.

int RasterMap::m_y

private

Definition at line 52 of file rastermap.h.

---

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

• D:/CJT/Repos/ALMaSS_all/Landscape/rastermap.h
• D:/CJT/Repos/ALMaSS_all/Landscape/rastermap.cpp