From fa5c61184cfbb6aed2022c9db37e330718f9ded3 Mon Sep 17 00:00:00 2001
From: Hal Simpson <hasimpson@usgs.gov>
Date: Tue, 7 Jul 2015 16:18:55 -0600
Subject: [PATCH] Added get_deltaf/get_computed_f_squares routine(s) for
computing deltaf
---
geomagio/ChannelConverter.py | 38 ++++++++++++++++++++++++++++++++++++
1 file changed, 38 insertions(+)
diff --git a/geomagio/ChannelConverter.py b/geomagio/ChannelConverter.py
index 8181e42c..1058f6c9 100644
--- a/geomagio/ChannelConverter.py
+++ b/geomagio/ChannelConverter.py
@@ -378,6 +378,44 @@ def get_obs_h_from_mag(h, d, d0=0):
return numpy.multiply(h, numpy.cos(obs_d))
+def get_deltaf(fv, fs):
+ """gets the deltaf value given the scalar F and computed F values.
+
+ Parameters
+ ----------
+ fv: array_like
+ the F vector computed from the observatory instruments
+ fs: array_like
+ the measured F value
+ """
+ return numpy.subtract(fv, fs)
+
+
+def get_computed_f_using_squares(x, y, z):
+ """gets the computed f value
+
+ Parameters
+ ----------
+ x: array_like
+ the x component from the observatory
+ y: array_like
+ the y component from the observatory
+ z: array_like
+ the z component from the observatory
+ Notes
+ -----
+ This works for geographic coordinates, or observatory coordinates.
+ ie x, y, z or h, e, z
+ We're using variables x,y,z to represent generic cartisian coordinates.
+ """
+ x2 = numpy.square(x)
+ y2 = numpy.square(y)
+ z2 = numpy.square(z)
+ fv = numpy.add(x2, y2)
+ fv = numpy.add(fv, z2)
+ return numpy.sqrt(fv)
+
+
def get_radians_from_minutes(m):
"""gets the radian value given the decimal value
Parameters
--
GitLab