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)
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	Init (const char *a_mapfile, Landscape *m_landscape)

Private Attributes
char	m_id [12]
	A 12 character ID for the map. More...
int *	m_map
	The raster of the landscape. More...
int	m_width
	The width of the landscape. More...
int	m_height
	The height of the landscape. More...
Landscape *	m_landscape
	A pointer to the parent landscape class. More...

Constructor & Destructor Documentation

RasterMap()

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

{
    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);
    Init(a_mapfile, a_landscape);
}

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

~RasterMap()

RasterMap::~RasterMap

(

void

)

{
  delete[] m_map;
}

Member Function Documentation

CellReplacementNeighbour()

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.

{
    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;
        found = count[i];
    }
    // Here we just make sure that we have not chosen another illegal. If so quietly ignore this replacement
    if (m_landscape->SupplyPolygonArea(surroundingcells[index]) < 2) return -1;
    // All OK so replace
    Put(a_x, a_y, surroundingcells[index]);
    return 1; // Signals that replacement occured
}

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

Get()

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

inline

                                            {
#ifndef DNDEBUG
    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!", "" );
    g_msg->WarnAddInfo(WARN_BUG, "X Coord ", a_x);
    g_msg->WarnAddInfo(WARN_BUG, "Y Coord ", a_y);
    g_msg->WarnAddInfo(WARN_BUG, "X max ", m_width);
    g_msg->WarnAddInfo(WARN_BUG, "Y max ", m_height);
    exit(1);
  }
#endif
  return m_map[ a_y * m_width + a_x ];
}

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

Referenced by Landscape::GetActualGooseGrazingForage(), Landscape::GetOwner_tole(), PesticideMap::Spray(), Landscape::SupplyAttIsForest(), Landscape::SupplyAttIsHigh(), Landscape::SupplyAttIsUrbanNoVeg(), Landscape::SupplyAttIsVegGrass(), Landscape::SupplyAttIsVegPatchy(), Landscape::SupplyAttIsWater(), Landscape::SupplyAttIsWoody(), Landscape::SupplyAttUserDefinedBool(), Landscape::SupplyBirdMaizeForage(), Landscape::SupplyBirdSeedForage(), Landscape::SupplyCentroidX(), Landscape::SupplyCentroidY(), Landscape::SupplyCountryDesig(), Landscape::SupplyCropType(), Landscape::SupplyDeadBiomass(), Landscape::SupplyElementSubType(), Landscape::SupplyElementType(), Landscape::SupplyElementTypeCC(), Landscape::SupplyFarmOwner(), Landscape::SupplyFarmOwnerIndex(), Landscape::SupplyFarmRotFilename(), Landscape::SupplyFarmType(), Landscape::SupplyGrazingPressure(), Landscape::SupplyGreenBiomass(), Landscape::SupplyGreenBiomassProp(), Landscape::SupplyGreenBiomassTotal(), Landscape::SupplyHasTramlines(), Landscape::SupplyInsects(), Landscape::SupplyInterestedGreenBiomass(), Landscape::SupplyInterestedGreenBiomassTotal(), Landscape::SupplyIsVeg(), Landscape::SupplyJustSprayed(), Landscape::SupplyLAGreen(), Landscape::SupplyLastSownVeg(), Landscape::SupplyLastTreatment(), Landscape::SupplyLATotal(), Landscape::SupplyNectar(), Landscape::SupplyOpenness(), Landscape::SupplyOwner_tole(), Landscape::SupplyPollen(), Landscape::SupplyPolyLEptr(), Landscape::SupplyPolyRef(), Landscape::SupplyPolyRefCC(), Landscape::SupplyPolyRefIndex(), Landscape::SupplyRoadWidth(), Landscape::SupplySoilType(), Landscape::SupplySoilTypeR(), Landscape::SupplySugar(), Landscape::SupplyTotalNectar(), Landscape::SupplyTotalPollen(), Landscape::SupplyTrafficLoad(), Landscape::SupplyUMRef(), Landscape::SupplyVegAge(), Landscape::SupplyVegBiomass(), Landscape::SupplyVegCover(), Landscape::SupplyVegDensity(), Landscape::SupplyVegDigestibility(), Landscape::SupplyVegGrowthStage(), Landscape::SupplyVegHeight(), Landscape::SupplyVegType(), Landscape::SupplyWeedBiomass(), and Pesticide::TwinMapSpray().

GetID()

char* RasterMap::GetID ( void  )

inline

{ return m_id; }

GetMagicP()

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

inline

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

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

Init()

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

protected

{
    m_landscape = a_landscape;
    char error_num[20];
    int   mx;
    int   my;
 
    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); // Needs smth that is 12 characters 
 
    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));
 
    mx = m_width;
    my = m_height;
 
    int pref = -1;
    int lastref = -1; // the first "last reference" is not good 
 
    for (int x = 0; x < mx; x++)
    {
        for (int y = 0; y < my; y++)
        {
            pref = Get(x, y);
 
            if (pref < 0)
            {
                // This fills in holes caused by negative polyrefs, but if two consequetive holes are found, then it will exit
                sprintf(error_num, "%d", pref);
                g_msg->Warn(WARN_FILE, "RasterMap::RasterMap(): Unknown/Negative polygon map number (could be a hole):", error_num);
                g_msg->WarnAddInfo(WARN_FILE, "X:", x);
                g_msg->WarnAddInfo(WARN_FILE, "Y:", y);
 
                if (lastref >= 0)
                {
                    sprintf(error_num, "%d", lastref);
                    g_msg->Warn(WARN_FILE, "Last polygon map number was good though, using it:", error_num);
                    Put(x, y, lastref);
                }
                else
                {
                    g_msg->Warn(WARN_FILE, "Last polygon map number was bad also, exiting","");
                    exit(1);
                }
            }
 
            lastref = pref;
 
            if (m_landscape->SupplyPolymapping(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);
            }
 
            if (m_landscape->SupplyPolymapping(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();
 
}

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

MapHeight()

int RasterMap::MapHeight ( void  )

inline

{ return m_height; }

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

MapWidth()

int RasterMap::MapWidth ( void  )

inline

{ return m_width; }

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

MissingCellReplace()

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

A method for removing missing polygons.

{
    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(g_random_fnc(8) + 8); // [8,9,..,15)
        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;
}

References g_random_fnc().

MissingCellReplaceWrap()

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

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

{
    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(g_random_fnc(8) + 8); // [8,9,..,15)
        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(random(8) + 8); // [8,9,..,15)
        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;
}

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

Put()

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

inline

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

Member Data Documentation

m_height

int RasterMap::m_height

private

The height of the landscape.

m_id

char RasterMap::m_id[12]

private

A 12 character ID for the map.

m_landscape

Landscape* RasterMap::m_landscape

private

A pointer to the parent landscape class.

m_map

int* RasterMap::m_map

private

The raster of the landscape.

m_width

int RasterMap::m_width

private

The width of the landscape.

---

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

• /builds/ALMaSS/ALMaSS_MIDox/Landscape/Rastermap.h
• /builds/ALMaSS/ALMaSS_MIDox/Landscape/Rastermap.cpp