infnanmean added to noise plots in SkyTEM

GeoscienceAustralia · Sep 15, 2024 · f3063ff · f3063ff
1 parent 165b94a
commit f3063ff
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/src/SkyTEM1DInversion.jl b/src/SkyTEM1DInversion.jl
@@ -329,8 +329,8 @@ function plotsoundingdata(nsoundings, LM_times, HM_times, d_LM, d_HM,
     axLM = ax[1].twiny()
     # axLM.semilogy(LM_noise, LM_times)
     # axLM.set_xlabel("high alt noise")
-    axLM.semilogy(mean(LM_noise./abs.(d_LM), dims=1)[:], LM_times, "r")
-    axLM.semilogy(mean(LM_noise./abs.(d_LM), dims=1)[:], LM_times, "--w")
+    axLM.semilogy(infnanmean(LM_noise./abs.(d_LM), 1)[:], LM_times, "r")
+    axLM.semilogy(infnanmean(LM_noise./abs.(d_LM), 1)[:], LM_times, "--w")
     axLM.set_xlabel("avg LM noise fraction")
     ax[2] = subplot(2,2,3,sharex=ax[1], sharey=ax[1])
     plot_dHM = permutedims(d_HM)
@@ -344,8 +344,8 @@ function plotsoundingdata(nsoundings, LM_times, HM_times, d_LM, d_HM,
     axHM = ax[2].twiny()
     # axHM.semilogy(HM_noise, HM_times)
     # axHM.set_xlabel("high alt noise")
-    axHM.semilogy(mean(HM_noise./abs.(d_HM), dims=1)[:], HM_times, "r")
-    axHM.semilogy(mean(HM_noise./abs.(d_HM), dims=1)[:], HM_times, "--w")
+    axHM.semilogy(infnanmean(HM_noise./abs.(d_HM), 1)[:], HM_times, "r")
+    axHM.semilogy(infnanmean(HM_noise./abs.(d_HM), 1)[:], HM_times, "--w")
     axHM.set_xlabel("avg HM noise fraction")
     ax[3] = subplot(2,2,2, sharex=ax[1])
     ax[3].plot(1:nsoundings, zRx, label="Rx")