data_import.jl

Author: boriskaus

docs/src/man/Tutorial_AlpineData.md

@@ -208,10 +208,10 @@ We will only consider earthquakes with a magnitude larger than 3. The resulting
 download_data("http://www.isc.ac.uk/cgi-bin/web-db-run?request=COLLECTED&req_agcy=ISC-EHB&out_format=QuakeML&ctr_lat=&ctr_lon=&radius=&max_dist_units=deg&searchshape=RECT&top_lat=49&bot_lat=37&left_lon=4&right_lon=20&srn=&grn=&start_year=1990&start_month=1&start_day=01&start_time=00%3A00%3A00&end_year=2015&end_month=12&end_day=31&end_time=00%3A00%3A00&min_dep=&max_dep=&min_mag=3.0&max_mag=&req_mag_type=Any&req_mag_agcy=Any&min_def=&max_def=&prime_only=on&include_magnitudes=on&table_owner=iscehb","ISCData.xml")
 ```
 
-Once the data has been downloaded, we can extract `lon/lat/depth/magnitude` using the `GMG` function `GetLonLatDepthMag_QuakeML`, which returns a `GeoData` structure:
+Once the data has been downloaded, we can extract `lon/lat/depth/magnitude` using the `GMG` function `getlonlatdepthmag_QuakeML`, which returns a `GeoData` structure:
 
 ```julia
-Data_ISC = GetLonLatDepthMag_QuakeML("ISCData.xml");
+Data_ISC = getlonlatdepthmag_QuakeML("ISCData.xml");
 nothing #hide
 ```
 

docs/src/man/Tutorial_Jura.md

@@ -41,7 +41,7 @@ upperright = [8.948117154811715, 47.781282316442606, 0.0]
 We can now import the map with the `Screensho_To_GeoData` function:
 
 ```julia
-Geology  = Screenshot_To_GeoData("SchoriM_Encl_01_Jura-map_A1.png", lowerleft, upperright, fieldname=:geology_colors) # name should have "colors" in it
+Geology  = screenshotToGeoData("SchoriM_Encl_01_Jura-map_A1.png", lowerleft, upperright, fieldname=:geology_colors) # name should have "colors" in it
 ```
 
 You can "drape" this image on the topographic map with
@@ -89,7 +89,7 @@ As example, we use the cross-section
 download_data("https://zenodo.org/records/10726801/files/Schori_2020_Ornans-Miserey-v2_whiteBG.png", "Schori_2020_Ornans-Miserey-v2_whiteBG.png")
 Corner_LowerLeft = (5.92507, 47.31300, -2.0)
 Corner_UpperRight = (6.25845, 46.99550, 2.0)
-CrossSection_1 = Screenshot_To_GeoData("Schori_2020_Ornans-Miserey-v2_whiteBG.png", Corner_LowerLeft, Corner_UpperRight) # name should have "colors" in it
+CrossSection_1 = screenshotToGeoData("Schori_2020_Ornans-Miserey-v2_whiteBG.png", Corner_LowerLeft, Corner_UpperRight) # name should have "colors" in it
 ```
 
 Note that we slightly modified the image to save it with a white instead of a transparent background

docs/src/man/dataimport.md

@@ -3,9 +3,9 @@
 We have a number of ways to import data, besides using any of the additional packages in julia to read files.
 
 ```@docs
-GeophysicalModelGenerator.Screenshot_To_GeoData
-GeophysicalModelGenerator.Screenshot_To_CartData
-GeophysicalModelGenerator.Screenshot_To_UTMData
+GeophysicalModelGenerator.screenshotToGeoData
+GeophysicalModelGenerator.screenshotToCartData
+GeophysicalModelGenerator.screenshotToUTMData
 GeophysicalModelGenerator.importTopo
 GeophysicalModelGenerator.importGeoTIFF
 ```

docs/src/man/tutorial_GMT_Topography_GeologicalMap.md

@@ -51,10 +51,10 @@ At this step, only the topographic data is imported. Now we have to import the t
 julia> Corner_LowerLeft  = (lon_min, lat_min , 0.0)
 julia> Corner_UpperRight = (lon_max, lat_max , 0.0)
 ```
-and import the png file with the GMG function `Screenshot_To_GeoData`:  
+and import the png file with the GMG function `screenshotToGeoData`:  
 
 ```julia
-julia> DataPNG = Screenshot_To_GeoData(filename_geo, Corner_LowerLeft, Corner_UpperRight)
+julia> DataPNG = screenshotToGeoData(filename_geo, Corner_LowerLeft, Corner_UpperRight)
 ```
 
 Next, we drape the screenshot on the topographic map with `drape_on_topo`:

docs/src/man/tutorial_Screenshot_To_Paraview.md

@@ -43,7 +43,7 @@ Once this is done, and we saved the picture under `Lippitsch_Fig13a.png`, you ca
 
 ```julia
 julia> using GeophysicalModelGenerator
-julia> data_profile1 = Screenshot_To_GeoData("Lippitsch_Fig13a.png",Corner_LowerLeft, Corner_UpperRight)
+julia> data_profile1 = screenshotToGeoData("Lippitsch_Fig13a.png",Corner_LowerLeft, Corner_UpperRight)
 Extracting GeoData from: Lippitsch_Fig13a.png
            └ Corners:         lon       lat       depth
               └ lower left  = (4.65   , 45.73  ,  -400.0 )
@@ -80,7 +80,7 @@ Corner_LowerLeft    =   ( 3.5, 43.0 , -150.0)
 Corner_UpperRight   =   (15.5, 50.0 , -150.0)
 Corner_LowerRight   =   (15.5, 43.0 , -150.0)
 Corner_UpperLeft    =   (3.5 , 50.0 , -150.0)
-data_Fig13_map      =   Screenshot_To_GeoData("Fig13_mapview.png",Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight,Corner_UpperLeft=Corner_UpperLeft)
+data_Fig13_map      =   screenshotToGeoData("Fig13_mapview.png",Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight,Corner_UpperLeft=Corner_UpperLeft)
 Write_Paraview(data_Fig13_map, "Lippitsch_Fig13_mapview")
 ```
 

src/data_import.jl

@@ -6,7 +6,7 @@
 
 using LightXML
 
-export Screenshot_To_GeoData, Screenshot_To_CartData, Screenshot_To_UTMData, GetLonLatDepthMag_QuakeML
+export screenshotToGeoData, screenshotToCartData, screenshotToUTMData, getlonlatdepthmag_QuakeML
 
 # import CSV data using standard library functions
 # here we assume that the data is indeed comma separated and that comments are preceded with a "#"
@@ -160,7 +160,7 @@ end
 
 
 """
-    Screenshot_To_GeoData(filename::String, Corner_LowerLeft, Corner_UpperRight; Corner_LowerRight=nothing, Corner_UpperLeft=nothing, Cartesian=false, UTM=false, UTMzone, isnorth=true, fieldname::Symbol=:colors)
+    screenshotToGeoData(filename::String, Corner_LowerLeft, Corner_UpperRight; Corner_LowerRight=nothing, Corner_UpperLeft=nothing, Cartesian=false, UTM=false, UTMzone, isnorth=true, fieldname::Symbol=:colors)
 
 Take a screenshot of Georeferenced image either a `lat/lon`, `x,y` (if `Cartesian=true`) or in UTM coordinates (if `UTM=true`) at a given depth or along profile and converts it to a `GeoData`, `CartData` or `UTMData` struct, which can be saved to Paraview
 
@@ -170,7 +170,7 @@ The lower right and upper left corners can be specified optionally (to take non-
 
 *Note*: if your data is in `UTM` coordinates you also need to provide the `UTMzone` and whether we are on the northern hemisphere or not (`isnorth`).
 """
-function Screenshot_To_GeoData(filename::String, Corner_LowerLeft, Corner_UpperRight; Corner_LowerRight=nothing, Corner_UpperLeft=nothing, Cartesian=false, UTM=false, UTMzone=nothing, isnorth::Bool=true, fieldname::Symbol=:colors)
+function screenshotToGeoData(filename::String, Corner_LowerLeft, Corner_UpperRight; Corner_LowerRight=nothing, Corner_UpperLeft=nothing, Cartesian=false, UTM=false, UTMzone=nothing, isnorth::Bool=true, fieldname::Symbol=:colors)
 
     img         =   load(filename)      # load image
 
@@ -287,39 +287,39 @@ end
 
 
 """
-    Data = Screenshot_To_CartData(filename::String, Corner_LowerLeft, Corner_UpperRight; Corner_LowerRight=nothing, Corner_UpperLeft=nothing)
+    Data = screenshotToCartData(filename::String, Corner_LowerLeft, Corner_UpperRight; Corner_LowerRight=nothing, Corner_UpperLeft=nothing)
 
-Does the same as `Screenshot_To_GeoData`, but returns a `CartData` structure
+Does the same as `screenshotToGeoData`, but returns a `CartData` structure
 """
-function Screenshot_To_CartData(filename::String, Corner_LowerLeft, Corner_UpperRight; Corner_LowerRight=nothing, Corner_UpperLeft=nothing, fieldname::Symbol=:colors)
+function screenshotToCartData(filename::String, Corner_LowerLeft, Corner_UpperRight; Corner_LowerRight=nothing, Corner_UpperLeft=nothing, fieldname::Symbol=:colors)
 
 
     # first create a GeoData struct
-    Data_Cart = Screenshot_To_GeoData(filename, Corner_LowerLeft, Corner_UpperRight; Corner_LowerRight=Corner_LowerRight, Corner_UpperLeft=Corner_UpperLeft, Cartesian=true, fieldname=fieldname)
+    Data_Cart = screenshotToGeoData(filename, Corner_LowerLeft, Corner_UpperRight; Corner_LowerRight=Corner_LowerRight, Corner_UpperLeft=Corner_UpperLeft, Cartesian=true, fieldname=fieldname)
 
     return Data_Cart
 
 end
 
 """
-    Data = Screenshot_To_UTMData(filename::String, Corner_LowerLeft, Corner_UpperRight; Corner_LowerRight=nothing, Corner_UpperLeft=nothing, UTMzone::Int64=nothing, isnorth::Bool=true, fieldname=:colors)
+    Data = screenshotToUTMData(filename::String, Corner_LowerLeft, Corner_UpperRight; Corner_LowerRight=nothing, Corner_UpperLeft=nothing, UTMzone::Int64=nothing, isnorth::Bool=true, fieldname=:colors)
 
-Does the same as `Screenshot_To_GeoData`, but returns for UTM data
+Does the same as `screenshotToGeoData`, but returns for UTM data
 Note that you have to specify the `UTMzone` and `isnorth`
 """
-function Screenshot_To_UTMData(filename::String, Corner_LowerLeft, Corner_UpperRight; Corner_LowerRight=nothing, Corner_UpperLeft=nothing, UTMzone::Int64=nothing, isnorth::Bool=true, fieldname::Symbol=:colors)
+function screenshotToUTMData(filename::String, Corner_LowerLeft, Corner_UpperRight; Corner_LowerRight=nothing, Corner_UpperLeft=nothing, UTMzone::Int64=nothing, isnorth::Bool=true, fieldname::Symbol=:colors)
 
       # first create a GeoData struct
-      Data_UTM = Screenshot_To_GeoData(filename, Corner_LowerLeft, Corner_UpperRight; Corner_LowerRight=Corner_LowerRight, Corner_UpperLeft=Corner_UpperLeft, Cartesian=false, UTM=true, UTMzone=UTMzone, isnorth=isnorth, fieldname=fieldname)
+      Data_UTM = screenshotToGeoData(filename, Corner_LowerLeft, Corner_UpperRight; Corner_LowerRight=Corner_LowerRight, Corner_UpperLeft=Corner_UpperLeft, Cartesian=false, UTM=true, UTMzone=UTMzone, isnorth=isnorth, fieldname=fieldname)
       return Data_UTM
 end
 
 """
-    Data = GetLonLatDepthMag_QuakeML(filename::String)
+    Data = getlonlatdepthmag_QuakeML(filename::String)
 
 Extracts longitude, latitude, depth and magnitude from a QuakeML file that has been e.g. downloaded from ISC. The data is then returned in GeoData format.
 """
-function GetLonLatDepthMag_QuakeML(filename::String)
+function getlonlatdepthmag_QuakeML(filename::String)
     # The QuakeML format consists of a tree with quite a lot of branches, so we have to traverse it to quite some extent to get the desired values
     # using LightXML: extension???
     xdoc = parse_file(filename); # parse the whole file

test/test_data_import.jl

@@ -32,25 +32,25 @@ using GeophysicalModelGenerator
 filename            =   "test.png";             # fake png
 Corner_LowerLeft    =   (18.0, 51.0, -590.0)
 Corner_UpperRight   =   (9.0, 42.0,    0.0)
-data_Image          =   Screenshot_To_GeoData(filename,Corner_LowerLeft, Corner_UpperRight)
+data_Image          =   screenshotToGeoData(filename,Corner_LowerLeft, Corner_UpperRight)
 @test data_Image.lon[1000] ≈ 17.592964824120603
 @test data_Image.lat[1000] ≈ 50.59296482412061
 @test Value(data_Image.depth[1000])==-590km
 @test Write_Paraview(data_Image, "Profile_1")==nothing
 
 # test if we use a different name for the color dataset
-data_Image_newfieldname  =   Screenshot_To_GeoData(filename,Corner_LowerLeft, Corner_UpperRight, fieldname=:fake)
+data_Image_newfieldname  =   screenshotToGeoData(filename,Corner_LowerLeft, Corner_UpperRight, fieldname=:fake)
 @test  keys(data_Image_newfieldname.fields)[1] == :fake
 
 # Test in CartData
-data_Image          =   Screenshot_To_GeoData(filename,Corner_LowerLeft, Corner_UpperRight, Cartesian=true)
+data_Image          =   screenshotToGeoData(filename,Corner_LowerLeft, Corner_UpperRight, Cartesian=true)
 @test Value(data_Image.x[22]) == 18.0km
 @test Value(data_Image.y[22]) == 51.0km
 @test Value(data_Image.z[22]) ≈ -125.15151515151516km
 
 
 # Test in UTM zone [note that depth should be in m]
-data_Image          =   Screenshot_To_GeoData(filename,Corner_LowerLeft, Corner_UpperRight, UTM=true, UTMzone=33, isnorth=true)
+data_Image          =   screenshotToGeoData(filename,Corner_LowerLeft, Corner_UpperRight, UTM=true, UTMzone=33, isnorth=true)
 @test data_Image.EW.val[22] == 18.0
 @test data_Image.NS.val[22] == 51.0
 @test Value(data_Image.depth[22]) ≈ -125.15151515151516m
@@ -61,27 +61,27 @@ Corner_LowerLeft    =   (2.0,  40.0, -15.0)
 Corner_UpperRight   =   (22.0, 51.0, -15.0)
 Corner_LowerRight   =   (20.0, 40.0, -15.0)
 Corner_UpperLeft    =   (0.0,  51.0, -15.0)
-data_Image          =   Screenshot_To_GeoData(filename,Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight, Corner_UpperLeft=Corner_UpperLeft)
+data_Image          =   screenshotToGeoData(filename,Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight, Corner_UpperLeft=Corner_UpperLeft)
 @test data_Image.lon[1000]  ≈  2.814070351758794
 @test data_Image.lat[1000]  ≈ 40.00000000000001
 @test Value(data_Image.depth[1000])==-15km
 @test Write_Paraview(data_Image, "MapView_1") == nothing
 
 # MapView in CartData
-data_Image          =   Screenshot_To_CartData(filename,Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight, Corner_UpperLeft=Corner_UpperLeft)
+data_Image          =   screenshotToCartData(filename,Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight, Corner_UpperLeft=Corner_UpperLeft)
 @test Value(data_Image.x[22]) ≈ 0.42424242424242425km
 @test Value(data_Image.y[22]) ≈ 48.666666666666664km
 @test Value(data_Image.z[22]) ≈ -15km
 
 # MapView in UTMData
-data_Image          =   Screenshot_To_UTMData(filename,Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight, Corner_UpperLeft=Corner_UpperLeft, UTMzone=33, isnorth=true)
+data_Image          =   screenshotToUTMData(filename,Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight, Corner_UpperLeft=Corner_UpperLeft, UTMzone=33, isnorth=true)
 @test data_Image.EW.val[22] ≈ 0.42424242424242425
 @test data_Image.NS.val[22] ≈ 48.666666666666664
 @test Value(data_Image.depth[22]) ≈ -15.0m
 
 # test the import of xml files from ISC
 # the search criteria are set in a way that only one event should be found
 download_data("http://www.isc.ac.uk/cgi-bin/web-db-run?request=COLLECTED&req_agcy=ISC-EHB&out_format=QuakeML&ctr_lat=&ctr_lon=&radius=&max_dist_units=deg&searchshape=RECT&top_lat=49&bot_lat=37&left_lon=4&right_lon=20&srn=&grn=&start_year=2000&start_month=1&start_day=01&start_time=00%3A00%3A00&end_year=2005&end_month=12&end_day=31&end_time=00%3A00%3A00&min_dep=&max_dep=&min_mag=5.8&max_mag=&req_mag_type=Any&req_mag_agcy=Any&min_def=&max_def=&include_magnitudes=on&include_links=on&include_headers=on&include_comments=on&table_owner=iscehb","ISCTest.xml")
-Data_ISC = GetLonLatDepthMag_QuakeML("ISCTest.xml");
+Data_ISC = getlonlatdepthmag_QuakeML("ISCTest.xml");
 @test Value(Data_ISC.depth[1])==-13.0km
 @test Data_ISC.fields.Magnitude[1]==5.8

tutorials/Lippitsch_Screenshots.jl

@@ -12,53 +12,53 @@ using GeophysicalModelGenerator
 
 # Process cross-sections of Figure 13. Note that we estimated some of the lon/lat locations
 
-data_Fig13a         =   Screenshot_To_GeoData("Lippitsch_Fig13a.png",( 4.65,45.73, -400.0), (17.23, 43.80, 0.0))
+data_Fig13a         =   screenshotToGeoData("Lippitsch_Fig13a.png",( 4.65,45.73, -400.0), (17.23, 43.80, 0.0))
 Write_Paraview(data_Fig13a, "Lippitsch_Fig13a") 
 
-data_Fig13b         =   Screenshot_To_GeoData("Lippitsch_Fig13b.png",( 5.51,51.53, -400.0), (12.04, 43.68 , 0.0))
+data_Fig13b         =   screenshotToGeoData("Lippitsch_Fig13b.png",( 5.51,51.53, -400.0), (12.04, 43.68 , 0.0))
 Write_Paraview(data_Fig13b, "Lippitsch_Fig13b") 
 
-data_Fig13c         =   Screenshot_To_GeoData("Lippitsch_Fig13c.png",(17.78,50.95, -400.0), (11.66, 43.68, 0.0))
+data_Fig13c         =   screenshotToGeoData("Lippitsch_Fig13c.png",(17.78,50.95, -400.0), (11.66, 43.68, 0.0))
 Write_Paraview(data_Fig13c, "Lippitsch_Fig13c") 
 
 # Mapview images
 Corner_LowerLeft    =   ( 3.5, 43.0 , -150.0)
 Corner_UpperRight   =   (15.5, 50.0 , -150.0)
 Corner_LowerRight   =   (15.5, 43.0 , -150.0)
 Corner_UpperLeft    =   (3.5 , 50.0 , -150.0)
-data_Fig13_map      =   Screenshot_To_GeoData("Fig13_mapview.png",Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight,Corner_UpperLeft=Corner_UpperLeft)
+data_Fig13_map      =   screenshotToGeoData("Fig13_mapview.png",Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight,Corner_UpperLeft=Corner_UpperLeft)
 Write_Paraview(data_Fig13_map, "Lippitsch_Fig13_mapview") 
 
 Depth                   =   -90.0;
 Corner_LowerLeft        =   (Corner_LowerLeft[1],   Corner_LowerLeft[2],    Depth)
 Corner_UpperRight       =   (Corner_UpperRight[1],  Corner_UpperRight[2],   Depth)
 Corner_LowerRight       =   (Corner_LowerRight[1],  Corner_LowerRight[2],   Depth)
 Corner_UpperLeft        =   (Corner_UpperLeft[1],   Corner_UpperLeft[2],    Depth)
-data_Fig12_90km    =   Screenshot_To_GeoData("Fig12_90km.png",Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight,Corner_UpperLeft=Corner_UpperLeft)
+data_Fig12_90km    =   screenshotToGeoData("Fig12_90km.png",Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight,Corner_UpperLeft=Corner_UpperLeft)
 Write_Paraview(data_Fig12_90km, "Lippitsch_Fig12_90km") 
 
 Depth                   = -180.0;
 Corner_LowerLeft        =   (Corner_LowerLeft[1],   Corner_LowerLeft[2],    Depth)
 Corner_UpperRight       =   (Corner_UpperRight[1],  Corner_UpperRight[2],   Depth)
 Corner_LowerRight       =   (Corner_LowerRight[1],  Corner_LowerRight[2],   Depth)
 Corner_UpperLeft        =   (Corner_UpperLeft[1],   Corner_UpperLeft[2],    Depth)
-data_Fig12_180km   =   Screenshot_To_GeoData("Fig12_180km.png",Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight,Corner_UpperLeft=Corner_UpperLeft)
+data_Fig12_180km   =   screenshotToGeoData("Fig12_180km.png",Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight,Corner_UpperLeft=Corner_UpperLeft)
 Write_Paraview(data_Fig12_180km, "Lippitsch_Fig12_180km") 
 
 Depth                   = -300.0;
 Corner_LowerLeft        =   (Corner_LowerLeft[1],   Corner_LowerLeft[2],    Depth)
 Corner_UpperRight       =   (Corner_UpperRight[1],  Corner_UpperRight[2],   Depth)
 Corner_LowerRight       =   (Corner_LowerRight[1],  Corner_LowerRight[2],   Depth)
 Corner_UpperLeft        =   (Corner_UpperLeft[1],   Corner_UpperLeft[2],    Depth)
-data_Fig12_300km   =   Screenshot_To_GeoData("Fig12_300km.png",Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight,Corner_UpperLeft=Corner_UpperLeft)
+data_Fig12_300km   =   screenshotToGeoData("Fig12_300km.png",Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight,Corner_UpperLeft=Corner_UpperLeft)
 Write_Paraview(data_Fig12_300km, "Lippitsch_Fig12_300km") 
 
 Depth                   = -400.0;
 Corner_LowerLeft        =   (Corner_LowerLeft[1],   Corner_LowerLeft[2],    Depth)
 Corner_UpperRight       =   (Corner_UpperRight[1],  Corner_UpperRight[2],   Depth)
 Corner_LowerRight       =   (Corner_LowerRight[1],  Corner_LowerRight[2],   Depth)
 Corner_UpperLeft        =   (Corner_UpperLeft[1],   Corner_UpperLeft[2],    Depth)
-data_Fig12_400km   =   Screenshot_To_GeoData("Fig12_400km.png",Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight,Corner_UpperLeft=Corner_UpperLeft)
+data_Fig12_400km   =   screenshotToGeoData("Fig12_400km.png",Corner_LowerLeft, Corner_UpperRight, Corner_LowerRight=Corner_LowerRight,Corner_UpperLeft=Corner_UpperLeft)
 Write_Paraview(data_Fig12_400km, "Lippitsch_Fig12_400km")