WIP:dfn file and proj

GeoscienceAustralia · Mar 7, 2024 · 469edf3 · 469edf3
1 parent c760049
commit 469edf3
Show file tree

Hide file tree

Showing 2 changed files with 21 additions and 5 deletions.
diff --git a/zz_portalcurtains/RDP.jl b/zz_portalcurtains/RDP.jl
@@ -312,22 +312,28 @@ function writevtkfromxyzrhodir(nlayers::Int; src_dir="", dst_dir="", src_epsg=0,
 end
 
 function writeaseggdffromxyzrho(nlayers::Int; src_dir="", dst_dir="", 
-         prefix="", epsg=0)
+         fname="", src_epsg=0)
     isdir(dst_dir) || mkpath(dst_dir)
-    sfmt = ["%15.3f", "%15.3f", "%15.3f", "%15.3f", "%10.5f", "%10.5f", "%10.5f", "%10.5f", "%10.5f", "%10.5f"], 
+    sfmt = ["%15i", "%15.3f", "%15.3f", "%15.3f", "%15.3f", "%10.5f", "%10.5f", "%10.5f", "%10.5f", "%10.5f", "%10.5f"] 
     channel_names = [["Line", "X", "Y", "Z", "zcenter", "log10_cond_low", "log10_cond_mid", "log10_cond_high", "log10_cond_avg", 
                       "phid_mean", "phid_sdev"], 
                      ["", "m", "m", "m", "m", "Log10_Siemens_per_m", "Log10_Siemens_per_m", "Log10_Siemens_per_m", "Log10_Siemens_per_m",
                       "", ""],
                      ["Line", "X", "Y", "Z", "zcenter", "log10_cond_low", "log10_cond_mid", "log10_cond_high", "log10_cond_avg", 
                       "phid_mean", "phid_sdev"]
                     ]
-    outfile =                 
+    outfile = joinpath(dst_dir, fname*"_EPSG_$src_epsg")
+    lines = transD_GP.getprobabilisticlinesfromdirectory(src_dir)
     map(lines) do ln
+        @info "Doing Line $ln"
         X, Y, Z, zall, ρlow, ρmid, ρhigh, ρavg, ϕmean, ϕsdev = transD_GP.readxyzrhoϕ(ln, nlayers; pathname=src_dir)
         for i in 1:length(X)
-            vonerow = [ln, X[i], Y[i], Z[i], zall, -ρhigh[:,i], -ρmid[:,i], -ρlow[:,i], ρavg[:,i], ϕmean[i], ϕsdev[i]]
-            transD_GP.CommonToAll.writeasegdat(vonerow, sfmt, outfile, mode::String)
+            vonerow = [ln, X[i], Y[i], Z[i], zall, -ρhigh[:,i], -ρmid[:,i], -ρlow[:,i], -ρavg[:,i], ϕmean[i], ϕsdev[i]]
+            mode = i == 1 ? "w" : "a"
+            transD_GP.CommonToAll.writeasegdat(vonerow, sfmt, outfile, mode)
+            transD_GP.CommonToAll.writeasegdfn(vonerow, channel_names, sfmt, outfile)
         end    
     end
+end
+
 end # module
diff --git a/zz_portalcurtains/makeaseggdf_from_xyzrho.jl b/zz_portalcurtains/makeaseggdf_from_xyzrho.jl
@@ -0,0 +1,10 @@
+using HiQGA, PyPlot
+includet("RDP.jl")
+src_epsg = 28354
+nlayers = 52
+src_dir = "/Users/anray/Documents/work/projects/largeaem/final_01/summaries_all/AusAEM_03_ERC/2021_ERC_01/summary"
+dst_dir = "/Users/anray/Documents/work/projects/curtainstuff/ERC_01_dat"
+parts = split(src_dir,"/")
+fname = parts[end-2]*"_"*parts[end-1]
+RDP.writeaseggdffromxyzrho(nlayers; src_dir, dst_dir, 
+         fname, src_epsg)