eotvos(Lat,Long,H,ds)
A function to calculate the Eotvos effect using the Harlan (1968) formula
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Ref: https://en.wikipedia.org/wiki/E%C3%B6tv%C3%B6s_effect
Ref: Harlan, R.B. Eotvos Corrections for Airborne Gravimetry 1968. J. Geophysics Res.. 14,17.
Relative to the GRS80 reference ellipsoid
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
The input values required are vectors
- Latitude (in decimal degrees)
- Longitude (in decimal degrees)
- Height (in metres)
- ds (the time interval between measurements in seconds)
The Output is the Eotvos effect in mGal.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Example usage:
e.g.
Latitude=[-43;-43.001;-43.002];
Longitude=[120;120.001;120.002];
H=[100;100;100];
ds=10;
Eotvoseffect=eotvos(Latitude,Longitude,H,ds)
%% Output
Eotvoseffect =
0
89.951763118312542
0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
N.B. The first and last of the output values are zero. This is because a derivative is taken, using a central difference, to estimate the velocity.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Cite As
Jack (2024). eotvos(Lat,Long,H,ds) (https://www.mathworks.com/matlabcentral/fileexchange/57361-eotvos-lat-long-h-ds), MATLAB Central File Exchange. Retrieved .
MATLAB Release Compatibility
Platform Compatibility
Windows macOS LinuxCategories
- Mathematics and Optimization > Mapping Toolbox > Coordinate Reference Systems >
- Radar > Mapping Toolbox > Coordinate Reference Systems >
- Sciences > Earth, Ocean, and Atmospheric Sciences > Geophysics >
Tags
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!Discover Live Editor
Create scripts with code, output, and formatted text in a single executable document.
Version | Published | Release Notes | |
---|---|---|---|
1.0.0.0 |