WebSVN - public iLand - Diff - Rev 1064 and 1068

 GrassCover::GrassCover()
+{
     mLayers = new GrassCoverLayers();
     mLayers->setGrid(mGrid, this);
     mEnabled = false;
+    mType = Invalid;
+}
 GrassCover::~GrassCover()
+{
     delete mLayers;
         return;
+    }
     // create the grid
     mGrid.setup(GlobalSettings::instance()->model()->grid()->metricRect(), GlobalSettings::instance()->model()->grid()->cellsize());
     mGrid.wipe();
-    QString formula = xml.value("model.settings.grass.grassPotential");
-    if (formula.isEmpty())
-        throw IException("setup of 'grass': required expression 'grassPotential' is missing.");
-    mGrassPotential.setExpression(formula);
-    mGrassPotential.linearize(0.,1., qMin(GRASSCOVERSTEPS, 1000));
-    formula = xml.value("model.settings.grass.grassEffect");
-    if (formula.isEmpty())
-        throw IException("setup of 'grass': required expression 'grassEffect' is missing.");
-    mGrassEffect.setExpression(formula);
-    mMaxTimeLag = xml.valueDouble("model.settings.grass.maxTimeLag");
-    if (mMaxTimeLag==0)
-        throw IException("setup of 'grass': value of 'maxTimeLag' is invalid or missing.");
-    mGrowthRate = GRASSCOVERSTEPS / mMaxTimeLag;
+    mType = Invalid;
+    QString type = xml.value("model.settings.grass.type");
+    if (type == QStringLiteral("pixel"))
+        mType = Pixel;
-    // set up the effect on regeneration in NSTEPS steps
-    for (int i=0;i<GRASSCOVERSTEPS;++i) {
-        double effect = mGrassEffect.calculate(i/double(GRASSCOVERSTEPS-1));
-        mEffect[i] = limit(effect, 0., 1.);
-    }
+    if (type==QStringLiteral("continuous"))
+        mType = Continuous;
-    mMaxState = limit(mGrassPotential.calculate(1.f), 0., 1.)*(GRASSCOVERSTEPS-1); // the max value of the potential function
+    if (mType == Invalid)
+        throw IException("GrassCover::setup: invalid 'grass.type'. Allowed: 'continous' and 'pixel'.");
+    if (mType == Pixel) {
+        // setup of pixel based / discrete approach
+        QString formula = xml.value("model.settings.grass.grassDuration");
+        if (formula.isEmpty())
+            throw IException("GrassCover::setup(): missing equation for 'grassDuration'.");
+        mPDF.setup(formula, 0., 100.);
+        mGrassLIFThreshold = xml.valueDouble("model.settings.grass.LIFThreshold", 0.2);
+        // clear array
+        for (int i=0;i<GRASSCOVERSTEPS;++i) {
+            mEffect[i] = 0.;
+        }
+    } else {
+        // setup of continuous grass concept
+        QString formula = xml.value("model.settings.grass.grassPotential");
+        if (formula.isEmpty())
+            throw IException("setup of 'grass': required expression 'grassPotential' is missing.");
+        mGrassPotential.setExpression(formula);
+        mGrassPotential.linearize(0.,1., qMin(GRASSCOVERSTEPS, 1000));
+        formula = xml.value("model.settings.grass.grassEffect");
+        if (formula.isEmpty())
+            throw IException("setup of 'grass': required expression 'grassEffect' is missing.");
+        mGrassEffect.setExpression(formula);
+        mMaxTimeLag = xml.valueDouble("model.settings.grass.maxTimeLag");
+        if (mMaxTimeLag==0)
+            throw IException("setup of 'grass': value of 'maxTimeLag' is invalid or missing.");
+        mGrowthRate = GRASSCOVERSTEPS / mMaxTimeLag;
+        // set up the effect on regeneration in NSTEPS steps
+        for (int i=0;i<GRASSCOVERSTEPS;++i) {
+            double effect = mGrassEffect.calculate(i/double(GRASSCOVERSTEPS-1));
+            mEffect[i] = limit(effect, 0., 1.);
+        }
+        mMaxState = limit(mGrassPotential.calculate(1.f), 0., 1.)*(GRASSCOVERSTEPS-1); // the max value of the potential function
+    }
     GlobalSettings::instance()->controller()->addLayers(mLayers, QStringLiteral("grass cover"));
     mEnabled = true;
     qDebug() << "setup of grass cover complete.";
 void GrassCover::setInitialValues(const QVector<float *> &LIFpixels, const int percent)
+{
     if (!enabled())
         return;
-    grid_type cval = limit(percent / 100., 0., 1.)*(GRASSCOVERSTEPS-1);
-    if (cval > mMaxState)
-        cval = mMaxState;
+    if (mType == Continuous) {
+        grass_grid_type cval = limit(percent / 100., 0., 1.)*(GRASSCOVERSTEPS-1);
+        if (cval > mMaxState)
+            cval = mMaxState;
-    Grid<float> *lif_grid = GlobalSettings::instance()->model()->grid();
-    for (QVector<float *>::const_iterator it = LIFpixels.constBegin(); it!=LIFpixels.constEnd(); ++it)
-        mGrid.valueAtIndex(lif_grid->indexOf(*it)) = cval;
+        Grid<float> *lif_grid = GlobalSettings::instance()->model()->grid();
+        for (QVector<float *>::const_iterator it = LIFpixels.constBegin(); it!=LIFpixels.constEnd(); ++it)
+            mGrid.valueAtIndex(lif_grid->indexOf(*it)) = cval;
+    } else {
+        // mType == Pixel
+        Grid<float> *lif_grid = GlobalSettings::instance()->model()->grid();
+        for (QVector<float *>::const_iterator it = LIFpixels.constBegin(); it!=LIFpixels.constEnd(); ++it) {
+            if (percent > irandom(0,100))
+                mGrid.valueAtIndex(lif_grid->indexOf(*it)) = mPDF.get();
+            else
+                mGrid.valueAtIndex(lif_grid->indexOf(*it)) = 0;
+        }
+    }
+}
 void GrassCover::execute()
+{
     if (!enabled())
         return;
     // Main function of the grass submodule
     float *lif = GlobalSettings::instance()->model()->grid()->begin();
     float *end_lif = GlobalSettings::instance()->model()->grid()->end();
-    grid_type *gr = mGrid.begin();
+    grass_grid_type *gr = mGrid.begin();
-    // loop over every LIF pixel
-    int skipped=0;
-    for (; lif!=end_lif;++lif, ++gr) {
-        // calculate potential grass vegetation cover
-        if (*lif == 1.f && *gr==mMaxState) {
-            ++skipped;
-            continue;
-        }
+    if (mType == Continuous) {
+        // loop over every LIF pixel
+        int skipped=0;
+        for (; lif!=end_lif;++lif, ++gr) {
+            // calculate potential grass vegetation cover
+            if (*lif == 1.f && *gr==mMaxState) {
+                ++skipped;
+                continue;
+            }
-        int potential = limit(mGrassPotential.calculate(*lif), 0., 1.)*(GRASSCOVERSTEPS-1);
-        *gr = qMin( int(*gr) + mGrowthRate, potential);
+            int potential = limit(mGrassPotential.calculate(*lif), 0., 1.)*(GRASSCOVERSTEPS-1);
+            *gr = qMin( int(*gr) + mGrowthRate, potential);
+        }
+        //qDebug() << "skipped" << skipped;
+    } else {
+        // type = Pixel
+        for (; lif!=end_lif;++lif, ++gr) {
+            if (*gr>1)
+                (*gr)--; // count down the years (until gr=1)
+            if (*gr==0 && *lif>mGrassLIFThreshold) {
+                // enable grass cover
+                double v = qMax(mPDF.get(), 0.);
+                *gr = v + 1;
+            }
+            if (*gr==1 && *lif<mGrassLIFThreshold) {
+                // now LIF is below the threshold - this enables the pixel get grassy again
+                *gr = 0;
+            }
+        }
+    }
-    qDebug() << "skipped" << skipped;
+}
-double GrassCoverLayers::value(const grid_type &data, const int index) const
+double GrassCoverLayers::value(const grass_grid_type &data, const int index) const
+{
     if (!mGrassCover->enabled()) return 0.;
     switch(index){
     case 0: return mGrassCover->effect(data); //effect
-    case 1: return data/double(GRASSCOVERSTEPS-1); // state
+    case 1: return mGrassCover->cover(data); // cover
     default: throw IException(QString("invalid variable index for a GrassCoverLayers: %1").arg(index));
+    }
+}
 const QVector<LayeredGridBase::LayerElement> &GrassCoverLayers::names()
+{
     if (mNames.isEmpty())
         mNames = QVector<LayeredGridBase::LayerElement>()
                 << LayeredGridBase::LayerElement(QLatin1Literal("effect"), QLatin1Literal("prohibiting effect on regeneration [0..1]"), GridViewGreens)
-                << LayeredGridBase::LayerElement(QLatin1Literal("cover"), QLatin1Literal("current grass cover on pixels [0..1]"), GridViewGreens);
+                << LayeredGridBase::LayerElement(QLatin1Literal("cover"), QLatin1Literal("current grass cover on pixels [0..1 for continuous, or #years for pixel mode]"), GridViewGreens);
     return mNames;
+}

Subversion Repositories public iLand

(root)/src/core/grasscover.cpp @ 1222 - Rev 1064 → 1068