Add the model using LERP instead of intervals.

Antoine Goutenoir
1 parent c20659f0
Showing 2 changed files with 74 additions and 0 deletions Show diff stats
content.yml
flaskr/laws/travel_emission_lerp_fit.py
@@ -62,6 +62,34 @@ models:
           # b is optional and defaults to 0
 #          - "0 500 4"
 #          - 500 1000 8
+  - name: MyClimate
+    # Slugged version of the display name above
+    # Only lowercase alphanumeric, starting with a letter, using - as liaison
+    # In other words, MUST match ^[a-z]([a-z0-9-]*[a-z0-9])?$ eg: icao-with-rfi
+    # MUST be unique, two models MUST NOT have the same slug.
+    slug: my-climate
+    # There MUST exist a python file named like this in `flaskr/laws`
+    # And it MUST contain a class named EmissionModel
+    # Please keep this lowercased, letters-only (or I will breathe fire)
+    file: travel_emission_lerp_fit
+    # The configuration that will be fed to the model.
+    # May be anything, really.  Go bonkers!
+    config:
+      plane_emission_linear_fit:
+        # A coefficient applied to the distance
+        connecting_flights_scale: 1.05
+        # Radiative Forcing Index
+        # Multiplier after scaling
+        rfi: 1.9
+        # Flat scalar to add before scaling with intervals
+        offset_before: 0
+        # Flat scalar to multiply before scaling with intervals
+        scale_before: 1
+        points:
+          - [   300.0,  120.42928309 ]
+          - [  1500.0,  278.88301355 ]
+          - [  2500.0,  438.31863895 ]
+          - [ 20000.0, 3335.62849772 ]
  
  
 # The content is Markdown.  HTML is also allowed.
@@ -0,0 +1,46 @@
+import numpy as np
+from .travel_emission_linear_fit import EmissionModel as BaseEmissionModel
+from flaskr.content import Struct
+
+
+class EmissionModel(BaseEmissionModel):
+
+    def apply_scaling_law(self, distance, config):
+
+        assert config.points
+        assert len(config.points) > 0
+
+        footprint = None
+
+        sample_points = sorted(config.points, key=lambda p: float(p[0]))
+
+        previous_point = [0, 0]
+        for i, sample_point in enumerate(sample_points):
+            if distance <= sample_point[0]:
+                t = 0
+                if sample_point[0] != previous_point[0]:
+                    t = (distance - previous_point[0]) * 1.0 \
+                                       / \
+                        (sample_point[0] - previous_point[0])
+
+                footprint = previous_point[1] + t * (sample_point[1] - previous_point[1])
+                previous_point = sample_point
+                break
+
+            previous_point = sample_point
+
+        if footprint is None:
+            if len(config.points) == 1:
+                last = sample_points[0]
+                penu = [0, 0]
+            else:
+                last = sample_points[-1]
+                penu = sample_points[-2]
+
+            t = 0
+            if last[0] != penu[0]:
+                t = (distance - last[0]) * 1.0 / (last[0] - penu[0])
+
+            footprint = last[1] + t * (last[1] - penu[1])
+
+        return footprint