mirror of
https://github.com/TheAlgorithms/Python.git
synced 2025-02-25 02:18:39 +00:00
Compare commits
2 Commits
2cfef0913a
...
8b831cb600
Author | SHA1 | Date | |
---|---|---|---|
|
8b831cb600 | ||
|
d31750adec |
52
physics/altitude_pressure.py
Normal file
52
physics/altitude_pressure.py
Normal file
@ -0,0 +1,52 @@
|
||||
"""
|
||||
Title : Calculate altitude using Pressure
|
||||
|
||||
Description :
|
||||
The below algorithm approximates the altitude using Barometric formula
|
||||
|
||||
|
||||
"""
|
||||
|
||||
|
||||
def get_altitude_at_pressure(pressure: float) -> float:
|
||||
"""
|
||||
This method calculates the altitude from Pressure wrt to
|
||||
Sea level pressure as reference .Pressure is in Pascals
|
||||
https://en.wikipedia.org/wiki/Pressure_altitude
|
||||
https://community.bosch-sensortec.com/t5/Question-and-answers/How-to-calculate-the-altitude-from-the-pressure-sensor-data/qaq-p/5702
|
||||
|
||||
H = 44330 * [1 - (P/p0)^(1/5.255) ]
|
||||
|
||||
Where :
|
||||
H = altitude (m)
|
||||
P = measured pressure
|
||||
p0 = reference pressure at sea level 101325 Pa
|
||||
|
||||
Examples:
|
||||
>>> get_altitude_at_pressure(pressure=100_000)
|
||||
105.47836610778828
|
||||
>>> get_altitude_at_pressure(pressure=101_325)
|
||||
0.0
|
||||
>>> get_altitude_at_pressure(pressure=80_000)
|
||||
1855.873388064995
|
||||
>>> get_altitude_at_pressure(pressure=201_325)
|
||||
Traceback (most recent call last):
|
||||
...
|
||||
ValueError: Value Higher than Pressure at Sea Level !
|
||||
>>> get_altitude_at_pressure(pressure=-80_000)
|
||||
Traceback (most recent call last):
|
||||
...
|
||||
ValueError: Atmospheric Pressure can not be negative !
|
||||
"""
|
||||
|
||||
if pressure > 101325:
|
||||
raise ValueError("Value Higher than Pressure at Sea Level !")
|
||||
if pressure < 0:
|
||||
raise ValueError("Atmospheric Pressure can not be negative !")
|
||||
return 44_330 * (1 - (pressure / 101_325) ** (1 / 5.5255))
|
||||
|
||||
|
||||
if __name__ == "__main__":
|
||||
import doctest
|
||||
|
||||
doctest.testmod()
|
Loading…
x
Reference in New Issue
Block a user