BugFix:

-Workarround for ScottPlot chrashes on NaNs ScottPlot/ScottPlot#4770
- Important signum issues for curvature and heading

Author: EMS90

TRA.Lib/Trasse.cs

@@ -5,6 +5,7 @@
 using System.Numerics;
 using System.Runtime.CompilerServices;
 using System.Collections.Specialized;
+using ScottPlot.Colormaps;
 
 #if USE_SCOTTPLOT
 using ScottPlot.Plottables;
@@ -539,10 +540,10 @@ public void Plot()
 
                 if (interpolation.H != null && interpolation.s != null)
                 {
-                    var scatterH = PlotG.Plot.Add.Scatter(interpolation.Y, interpolation.H, color);
+                    var scatterH = PlotG.Plot.Add.Scatter(interpolation.Y, interpolation.H.Where(i => !double.IsNaN(i)).ToArray(), color); //BugFix as ScottPlot Crashes on NaNs should be fixed in future release https://github.com/ScottPlot/ScottPlot/pull/4770
                     //scatterH.LegendText = "Elevation";
                     scatterH.Axes.YAxis = PlotG.Plot.Axes.Left;
-                    var scatterSlope = PlotG.Plot.Add.ScatterLine(interpolation.Y, interpolation.s, color);
+                    var scatterSlope = PlotG.Plot.Add.ScatterLine(interpolation.Y, interpolation.s.Where(i => !double.IsNaN(i)).ToArray(), color); //BugFix as ScottPlot Crashes on NaNs should be fixed in future release https://github.com/ScottPlot/ScottPlot/pull/4770
                     //scatterSlope.LegendText = "Slope";
                     scatterSlope.Axes.YAxis = PlotG.Plot.Axes.Right;
                 }

TRA.Lib/TrassenElementExt.cs

@@ -52,6 +52,7 @@ public Interpolation(int num = 0)
         }
         public void Concat(Interpolation interp)
         {
+            if (interp.X ==  null || interp.Y == null || interp.S == null || interp.T == null || interp.K == null) { return; }
             X = X.Concat(interp.X).ToArray();
             Y = Y.Concat(interp.Y).ToArray();
             S = S.Concat(interp.S).ToArray();
@@ -273,7 +274,11 @@ void AddWarningCallout(string text, double X, double Y)
         public ref Interpolation Interpolate(double delta = 1.0, double allowedTolerance = 0.001)
         {
             Transform2D transform = new Transform2D(x, y, t);
-            if (TrassenGeometrie == null) { AddWarningCallout("No Gemetry for interpolation " + kz.ToString() + "set, maybe not implemented yet", Xstart,Ystart); return ref Interpolation; }
+            if (TrassenGeometrie == null) { 
+                AddWarningCallout("No Gemetry for interpolation " + kz.ToString() + "set, maybe not implemented yet", Xstart,Ystart); 
+                Interpolation = new Interpolation(0);
+                return ref Interpolation; 
+            }
             if (Double.IsNaN(l)) { AddWarningCallout("Length is NaN, no interpolation calculated", Xstart, Ystart); return ref Interpolation; }
             List<double> Xlst = new List<double>();
             List<double> Ylst = new List<double>();
@@ -320,6 +325,7 @@ public ref Interpolation Interpolate(double delta = 1.0, double allowedTolerance
 
         public double GetSAtPoint(double X, double Y, double T = double.NaN)
         {
+            if (TrassenGeometrie == null) return (0);
             Transform2D transform = new Transform2D(x, y, t);
             transform.ApplyInverse(ref X, ref Y, ref T);
             return TrassenGeometrie.sAt(X, Y, T) + s;
@@ -331,6 +337,7 @@ public double GetSAtPoint(double X, double Y, double T = double.NaN)
         /// <returns>Hochwert X,Rechtswert Y, Heading T</returns>
         public (double, double, double) GetPointAtS(double S)
         {
+            if (TrassenGeometrie == null) return (0, 0, 0);// (Double.NaN, Double.NaN, Double.NaN);
             (double X, double Y, double T, _) = TrassenGeometrie.PointAt(S - s);
             Transform2D transform = new Transform2D(x, y, t);
             transform.Apply(ref X, ref Y, ref T);

TRA.Lib/Trasseninterpolation.cs

@@ -66,7 +66,7 @@ public override (double X, double Y, double t, double k) PointAt(double s)
             double r = Math.Abs(radius);
             if (r == 0) { return (s, 0.0, 0.0, 0.0); } //Gerade
             (double X, double Y) = Math.SinCos(s / r);
-            return (X * r, sig * ((1 - Y) * r), s / r, 1 / radius);
+            return (X * r, sig * ((1 - Y) * r), s / radius, 1 / radius);
         }
 
         public override double sAt(double X, double Y, double t = double.NaN)
@@ -177,7 +177,7 @@ public override (double X, double Y, double t, double k) PointAt(double s)
 
             //Tangent at point
             double theta = gamma * s * s * 0.5 + curvature1 * s;
-            return (Cs.Real, dir * Math.Sign(gamma) * Cs.Imaginary, theta, curvature1 + gamma * s);
+            return (Cs.Real, dir * Math.Sign(gamma) * Cs.Imaginary, dir * theta, dir* (curvature1 + gamma * s));
         }
 
         // Implementing the Fresnel integrals using numerical integration

TRA.Tool/GCG2016v2023

@@ -47,4 +47,13 @@
     </EmbeddedResource>
   </ItemGroup>
 
+  <ItemGroup>
+    <None Update="bkg_license.txt">
+      <CopyToOutputDirectory>PreserveNewest</CopyToOutputDirectory>
+    </None>
+    <None Update="GCG2016v2023">
+      <CopyToOutputDirectory>PreserveNewest</CopyToOutputDirectory>
+    </None>
+  </ItemGroup>
+
 </Project>

TRA.Tool/TRA.Tool.csproj

@@ -64,7 +64,7 @@
   <data name="btn_delete.Image" type="System.Drawing.Bitmap, System.Drawing" mimetype="application/x-microsoft.net.object.bytearray.base64">
     <value>
         iVBORw0KGgoAAAANSUhEUgAAABAAAAAQCAYAAAAf8/9hAAAABGdBTUEAALGPC/xhBQAAACBjSFJNAAB6
-        JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAACXBIWXMAAA6+AAAOvgHqQrHAAAAAV0lE
+        JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAACXBIWXMAAA68AAAOvAGVvHJJAAAAV0lE
         QVQ4T2PABhQUFI4qKir+h2Eg/zBUijgA0gRlggE6HwMg20YshmpFAJCgnJycJCGMVTMIUNUAGBudhrGh
         WlDBqAHDzgB8mKYGPABJEoEfQLUAAQMDAOE0kkF4hF42AAAAAElFTkSuQmCC
 </value>

TRA.Tool/TrassenPanel.Designer.cs

@@ -0,0 +1,3 @@
+© GeoBasis-DE / BKG 2025 CC BY 4.0
+BKG: https://www.bkg.bund.de
+CC BY 4.0: https://creativecommons.org/licenses/by/4.0/