コード例 #1
0
ファイル: reclass.cs プロジェクト: Leowisd/landisProFire2Core
        //This will perform a reclassification based upon the oldest cohort upon a landis stand.
        //The cohorts will be scaled into 16 age classes.
        public static void ageReclass(map8 m)
        {
            m.dim(snr, snc);

            m.rename("Age class representation");

            for (uint j = 1; j < map8.MapmaxValue; j++)
            {
                m.assignLeg(j, "");
            }


            string str;

            //J.Yang hard coding changing itr*sites.TimeStep to itr
            //J.Yang maxLeg is defined as 256 in map8.h, therefore, maximum age cohorts it can output is 254
            for (uint i = 1; i < map8.MaxValueforLegend - 4; i++)
            {
                str = string.Format("{0:   } - {1:   } yr", (i - 1) * time_step + 1, i * time_step);

                m.assignLeg(i, str);
            }


            m.assignLeg(0, "NoSpecies");

            m.assignLeg(map8.MaxValueforLegend - 1, "N/A");
            m.assignLeg(map8.MaxValueforLegend - 2, "Water");
            m.assignLeg(map8.MaxValueforLegend - 3, "NonForest");

            str = string.Format("	  >{0:   } yr", (map8.MaxValueforLegend - 4 - 1) * time_step);
            m.assignLeg(map8.MaxValueforLegend - 4, str);

            for (uint i = snr; i >= 1; i--)
            {
                for (uint j = 1; j <= snc; j++)
                {
                    if (PlugIn.gl_sites.locateLanduPt(i, j).active())
                    {
                        m[i, j] = 0;

                        uint myage = 0;

                        site local_site = PlugIn.gl_sites[i, j];

                        specie s = local_site.first();

                        while (s != null)
                        {
                            uint temp = s.oldest();

                            if (temp > myage)
                            {
                                myage = temp;
                            }

                            s = local_site.next();
                        }

                        m[i, j] = (ushort)(myage / time_step);
                    }

                    else if (PlugIn.gl_sites.locateLanduPt(i, j).lowland())
                    {
                        m[i, j] = (ushort)(map8.MaxValueforLegend - 3);
                    }

                    else if (PlugIn.gl_sites.locateLanduPt(i, j).water())
                    {
                        m[i, j] = (ushort)(map8.MaxValueforLegend - 2);
                    }

                    else
                    {
                        m[i, j] = (ushort)(map8.MaxValueforLegend - 1);
                    }
                }
            }
        }
コード例 #2
0
        //This will reclassify a singular site.  M is the number of possible output classes.
        private int reclassifySite(site site_in, int m)
        {
            float[] sval = new float[MAX_RECLASS];

            specie local_specie = site_in.first();

            int j = 1;

            while (local_specie != null)
            {
                float c = (float)local_specie.oldest() / maximum[j - 1];

                if (c > 1.0)
                {
                    c = 1.0f;
                }

                for (int i = 1; i <= m; i++)
                {
                    if (BOOL[i, j] != 0)
                    {
                        if (BOOL[i, j] > 0)
                        {
                            sval[i] += c;
                        }
                        else
                        {
                            sval[i] -= c;
                        }


                        if (sval[i] != 0)
                        {
                            if (sval[i] > 1.0)
                            {
                                sval[i] = 1.0f;
                            }

                            if (sval[i] < 0.0)
                            {
                                sval[i] = 0.0f;
                            }
                        }
                    }
                }

                j++;

                local_specie = site_in.next();
            }



            int mx = 0;

            float mxVal = 0.0f;

            for (int i = 1; i <= m; i++)
            {
                if (sval[i] > mxVal)
                {
                    mxVal = sval[i];

                    mx = i;
                }
            }


            if (mxVal > 0.0)
            {
                return(mx);
            }
            else
            {
                return(m + 1);
            }
        }