//--------------------------------------------------------------------- //--------------------------------------------------------------------- /// <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); } } } }