Added the functionality to call the extract_ProfileData! function without an empty Named Tuple for ScreenshotData

Author: mthielma

src/ProfileProcessing.jl

@@ -17,6 +17,7 @@ Structure that holds profile data (interpolated/projected on the profile)
         VolData         ::  GeophysicalModelGenerator.GeoData
         SurfData        ::  Union{Nothing, NamedTuple}
         PointData       ::  Union{Nothing, NamedTuple}
+        ScreenshotData  ::Union{Nothing, NamedTuple}
     end
 
     Structure to store cross section data
@@ -52,7 +53,6 @@ mutable struct ProfileData
     end
 end
 
-
 function show(io::IO, g::ProfileData)
     if g.vertical
         println(io, "Vertical ProfileData")
@@ -70,7 +70,9 @@ function show(io::IO, g::ProfileData)
     if !isnothing(g.PointData)
         println(io, "        PointData: $(keys(g.PointData)) ")
     end
-
+    if !isnothing(g.ScreenshotData)
+        println(io, "        ScreenshotData: $(keys(g.ScreenshotData)) ")
+    end
     return nothing
 end
 
@@ -405,6 +407,21 @@ function extract_ProfileData!(Profile::ProfileData, VolData::Union{Nothing, GeoD
     return nothing
 end
 
+"""
+    extract_ProfileData!(Profile::ProfileData,VolData::GeoData, SurfData::NamedTuple, PointData::NamedTuple; DimsVolCross=(100,100),Depth_extent=nothing,DimsSurfCross=(100,),section_width=50)
+
+Extracts data along a vertical or horizontal profile. 
+"""
+function extract_ProfileData!(Profile::ProfileData, VolData::Union{Nothing, GeoData} = nothing, SurfData::NamedTuple = NamedTuple(), PointData::NamedTuple = NamedTuple(); DimsVolCross = (100, 100), Depth_extent = nothing, DimsSurfCross = (100,), section_width = 50km)
+
+    # call the actual function to extract the profile data
+    extract_ProfileData!(Profile, VolData, SurfData, PointData,NamedTuple(); DimsVolCross = DimsVolCross, Depth_extent = Depth_extent, DimsSurfCross = DimsSurfCross, section_width = section_width)
+    return nothing
+end
+
+
+
+
 """
 This reads the picked profiles from disk and returns a vector of ProfileData
 """

test/test_ProfileProcessing.jl

@@ -22,7 +22,6 @@ data_Vol2 = GMG_Dataset("Plomerova2022", "Volume", "https://seafile.rlp.net/f/ab
 #data_Vol2   = GMG_Dataset("Zhao2016","Volume","https://seafile.rlp.net/f/e81a6d075f6746609973/?dl=1", true)
 
 # Now load these datasets into NamedTuples
-
 SurfData = load_GMG(data_Surf)
 PointData = load_GMG(data_EQ)
 ScreenshotData = load_GMG(data_SS)
@@ -87,26 +86,26 @@ GeophysicalModelGenerator.create_profile_screenshot!(prof5, Data.Screenshot)
 
 
 # Test the main profile extraction routines:
-extract_ProfileData!(prof1, VolData_combined1, Data.Surface, Data.Point, NamedTuple())
-extract_ProfileData!(prof2, VolData_combined1, Data.Surface, Data.Point, NamedTuple())
-extract_ProfileData!(prof3, VolData_combined1, Data.Surface, Data.Point, NamedTuple())
-extract_ProfileData!(prof4, VolData_combined1, Data.Surface, Data.Point, NamedTuple())
-
-extract_ProfileData!(prof1, VolData_combined2, Data.Surface, Data.Point, NamedTuple())
-extract_ProfileData!(prof2, VolData_combined2, Data.Surface, Data.Point, NamedTuple())
-extract_ProfileData!(prof3, VolData_combined2, Data.Surface, Data.Point, NamedTuple())
-extract_ProfileData!(prof4, VolData_combined2, Data.Surface, Data.Point, NamedTuple())
-
-extract_ProfileData!(prof1, VolData_combined3, Data.Surface, Data.Point, NamedTuple())
-extract_ProfileData!(prof2, VolData_combined3, Data.Surface, Data.Point, NamedTuple())
-extract_ProfileData!(prof3, VolData_combined3, Data.Surface, Data.Point, NamedTuple())
-extract_ProfileData!(prof4, VolData_combined3, Data.Surface, Data.Point, NamedTuple())
+extract_ProfileData!(prof1, VolData_combined1, Data.Surface, Data.Point)
+extract_ProfileData!(prof2, VolData_combined1, Data.Surface, Data.Point)
+extract_ProfileData!(prof3, VolData_combined1, Data.Surface, Data.Point)
+extract_ProfileData!(prof4, VolData_combined1, Data.Surface, Data.Point)
+
+extract_ProfileData!(prof1, VolData_combined2, Data.Surface, Data.Point)
+extract_ProfileData!(prof2, VolData_combined2, Data.Surface, Data.Point)
+extract_ProfileData!(prof3, VolData_combined2, Data.Surface, Data.Point)
+extract_ProfileData!(prof4, VolData_combined2, Data.Surface, Data.Point)
+
+extract_ProfileData!(prof1, VolData_combined3, Data.Surface, Data.Point)
+extract_ProfileData!(prof2, VolData_combined3, Data.Surface, Data.Point)
+extract_ProfileData!(prof3, VolData_combined3, Data.Surface, Data.Point)
+extract_ProfileData!(prof4, VolData_combined3, Data.Surface, Data.Point)
 extract_ProfileData!(prof5, VolData_combined3, Data.Surface, Data.Point, Data.Screenshot)
 
 
 # Test that it works if only EQ's are provided:
 prof4 = ProfileData(depth = -20)
-extract_ProfileData!(prof4, nothing, NamedTuple(), Data.Point, NamedTuple())
+extract_ProfileData!(prof4, nothing, NamedTuple(), Data.Point)
 @test isnothing(prof4.VolData)
 @test isempty(prof4.SurfData)
 @test length(prof4.PointData[1].depth) == 3280