//---------------------------------------------------------------------
//---------------------------------------------------------------------
/// <summary>
/// Reads the initial communities map, finds all unique site keys, and sets aside sites to process first and second
/// </summary>
private void ProcessInitialCommunitiesMap(string initialCommunitiesMap,
Landis.Library.InitialCommunities.IDataset communities, ref List <ActiveSite> processFirst,
ref List <ActiveSite> processSecond)
{
IInputRaster <uintPixel> map = ModelCore.OpenRaster <uintPixel>(initialCommunitiesMap);
Dictionary <uint, ActiveSite> uniqueKeys = new Dictionary <uint, ActiveSite>();
using (map)
{
uintPixel pixel = map.BufferPixel;
foreach (Site site in ModelCore.Landscape.AllSites)
{
map.ReadBufferPixel();
uint mapCode = pixel.MapCode.Value;
if (!site.IsActive)
{
continue;
}
ActiveSite activeSite = (ActiveSite)site;
var initialCommunity = communities.Find(mapCode);
if (initialCommunity == null)
{
throw new ApplicationException(string.Format("Unknown map code for initial community: {0}", mapCode));
}
sitesAndCommunities.Add(activeSite, initialCommunity);
uint key = SiteCohorts.ComputeKey((ushort)initialCommunity.MapCode, Globals.ModelCore.Ecoregion[site].MapCode);
if (!uniqueKeys.ContainsKey(key))
{
uniqueKeys.Add(key, activeSite);
processFirst.Add(activeSite);
}
else
{
processSecond.Add(activeSite);
}
if (!allKeys.ContainsKey(activeSite))
{
allKeys.Add(activeSite, key);
}
}
}
}
