LSST Applications g0f08755f38+9c285cab97,g1635faa6d4+13f3999e92,g1653933729+a8ce1bb630,g1a0ca8cf93+bf6eb00ceb,g28da252d5a+0829b12dee,g29321ee8c0+5700dc9eac,g2bbee38e9b+9634bc57db,g2bc492864f+9634bc57db,g2cdde0e794+c2c89b37c4,g3156d2b45e+41e33cbcdc,g347aa1857d+9634bc57db,g35bb328faa+a8ce1bb630,g3a166c0a6a+9634bc57db,g3e281a1b8c+9f2c4e2fc3,g414038480c+077ccc18e7,g41af890bb2+fde0dd39b6,g5fbc88fb19+17cd334064,g781aacb6e4+a8ce1bb630,g80478fca09+55a9465950,g82479be7b0+d730eedb7d,g858d7b2824+9c285cab97,g9125e01d80+a8ce1bb630,g9726552aa6+10f999ec6a,ga5288a1d22+2a84bb7594,gacf8899fa4+c69c5206e8,gae0086650b+a8ce1bb630,gb58c049af0+d64f4d3760,gc28159a63d+9634bc57db,gcf0d15dbbd+4b7d09cae4,gda3e153d99+9c285cab97,gda6a2b7d83+4b7d09cae4,gdaeeff99f8+1711a396fd,ge2409df99d+5e831397f4,ge79ae78c31+9634bc57db,gf0baf85859+147a0692ba,gf3967379c6+41c94011de,gf3fb38a9a8+8f07a9901b,gfb92a5be7c+9c285cab97,w.2024.46
LSST Data Management Base Package
Loading...
Searching...
No Matches
Functions | Variables
lsst.afw.coord._refraction Namespace Reference

Functions

 refraction (wavelength, elevation, observatory, weather=None)
 
 differentialRefraction (wavelength, wavelengthRef, elevation, observatory, weather=None)
 
 deltaN (wavelength, weather)
 
 densityFactorDry (weather)
 
 densityFactorWater (weather)
 
 humidityToPressure (weather)
 
 extractTemperature (weather, useKelvin=False)
 
 defaultWeather (altitude)
 

Variables

float deltaRefractScale = 1.0E8
 

Function Documentation

◆ defaultWeather()

lsst.afw.coord._refraction.defaultWeather ( altitude) 
Set default local weather conditions if they are missing.

Parameters
----------
weather : `lsst.afw.coord.Weather`
    Class containing the measured temperature, pressure, and humidity
    at the observatory during an observation
altitude : `astropy.units.Quantity`
    The altitude of the observatory, in meters.

Returns
-------
default : `lsst.afw.coord.Weather`
    Updated Weather class with any `nan` values replaced by defaults.

Definition at line 293 of file _refraction.py.

293def defaultWeather(altitude):
294 """Set default local weather conditions if they are missing.
295
296 Parameters
297 ----------
298 weather : `lsst.afw.coord.Weather`
299 Class containing the measured temperature, pressure, and humidity
300 at the observatory during an observation
301 altitude : `astropy.units.Quantity`
302 The altitude of the observatory, in meters.
303
304 Returns
305 -------
306 default : `lsst.afw.coord.Weather`
307 Updated Weather class with any `nan` values replaced by defaults.
308 """
309 if isinstance(altitude, units.quantity.Quantity):
310 altitude2 = altitude
311 else:
312 altitude2 = altitude*units.meter
313 p0 = 101325.*units.pascal # sea level air pressure
314 g = 9.80665*units.meter/units.second**2 # typical gravitational acceleration at sea level
315 R0 = 8.31447*units.Joule/(units.mol*units.Kelvin) # gas constant
316 T0 = 19.*units.Celsius # Typical sea-level temperature
317 lapseRate = -6.5*units.Celsius/units.km # Typical rate of change of temperature with altitude
318 M = 0.0289644*units.kg/units.mol # molar mass of dry air
319
320 temperature = T0 + lapseRate*altitude2
321 temperatureK = temperature.to(units.Kelvin, equivalencies=units.temperature())
322 pressure = p0*np.exp(-(g*M*altitude2)/(R0*temperatureK))
323 humidity = 40. # Typical humidity at many observatory sites.
324 weather = Weather((temperature/units.Celsius).value, (pressure/units.pascal).value, humidity)
325 return weather

◆ deltaN()

lsst.afw.coord._refraction.deltaN ( wavelength,
weather ) 
Calculate the differential refractive index of air.

Parameters
----------
wavelength : `float`
    wavelength is in nanometers
weather : `lsst.afw.coord.Weather`
    Class containing the measured temperature, pressure, and humidity
    at the observatory during an observation

Returns
-------
deltaN : `float`
    The difference of the refractive index of air from 1.,
    calculated as (n_air - 1)*10^8

Notes
-----
The differential refractive index is the difference of
the refractive index from 1., multiplied by 1E8 to simplify
the notation and equations. Calculated as (n_air - 1)*10^8

This replicates equation 14 of [1]_

References
----------
.. [1] R. C. Stone, "An Accurate Method for Computing Atmospheric
   Refraction," Publications of the Astronomical Society of the Pacific,
   vol. 108, p. 1051, 1996.

Definition at line 123 of file _refraction.py.

123def deltaN(wavelength, weather):
124 """Calculate the differential refractive index of air.
125
126 Parameters
127 ----------
128 wavelength : `float`
129 wavelength is in nanometers
130 weather : `lsst.afw.coord.Weather`
131 Class containing the measured temperature, pressure, and humidity
132 at the observatory during an observation
133
134 Returns
135 -------
136 deltaN : `float`
137 The difference of the refractive index of air from 1.,
138 calculated as (n_air - 1)*10^8
139
140 Notes
141 -----
142 The differential refractive index is the difference of
143 the refractive index from 1., multiplied by 1E8 to simplify
144 the notation and equations. Calculated as (n_air - 1)*10^8
145
146 This replicates equation 14 of [1]_
147
148 References
149 ----------
150 .. [1] R. C. Stone, "An Accurate Method for Computing Atmospheric
151 Refraction," Publications of the Astronomical Society of the Pacific,
152 vol. 108, p. 1051, 1996.
153 """
154 waveNum = 1E3/wavelength # want wave number in units 1/micron
155 dryAirTerm = 2371.34 + (683939.7/(130. - waveNum**2.)) + (4547.3/(38.9 - waveNum**2.))
156 wetAirTerm = 6487.31 + 58.058*waveNum**2. - 0.71150*waveNum**4. + 0.08851*waveNum**6.
157 return (dryAirTerm*densityFactorDry(weather) + wetAirTerm*densityFactorWater(weather))
158
159

◆ densityFactorDry()

lsst.afw.coord._refraction.densityFactorDry ( weather) 
Calculate dry air pressure term to refractive index calculation.

Parameters
----------
weather : `lsst.afw.coord.Weather`
    Class containing the measured temperature, pressure, and humidity
    at the observatory during an observation

Returns
-------
densityFactor : `float`
    Returns the relative density of dry air
    at the given pressure and temperature.

Notes
-----
This replicates equation 15 of [1]_

References
----------
.. [1] R. C. Stone, "An Accurate Method for Computing Atmospheric
   Refraction," Publications of the Astronomical Society of the Pacific,
   vol. 108, p. 1051, 1996.

Definition at line 160 of file _refraction.py.

160def densityFactorDry(weather):
161 """Calculate dry air pressure term to refractive index calculation.
162
163 Parameters
164 ----------
165 weather : `lsst.afw.coord.Weather`
166 Class containing the measured temperature, pressure, and humidity
167 at the observatory during an observation
168
169 Returns
170 -------
171 densityFactor : `float`
172 Returns the relative density of dry air
173 at the given pressure and temperature.
174
175 Notes
176 -----
177 This replicates equation 15 of [1]_
178
179 References
180 ----------
181 .. [1] R. C. Stone, "An Accurate Method for Computing Atmospheric
182 Refraction," Publications of the Astronomical Society of the Pacific,
183 vol. 108, p. 1051, 1996.
184 """
185 temperature = extractTemperature(weather, useKelvin=True)
186 waterVaporPressure = humidityToPressure(weather)
187 airPressure = weather.getAirPressure()*units.pascal
188 dryPressure = airPressure - waterVaporPressure
189 eqn = dryPressure.to_value(cds.mbar)*(57.90E-8 - 9.3250E-4/temperature.to_value(units.Kelvin)
190 + 0.25844/temperature.to_value(units.Kelvin)**2.)
191 densityFactor = (1. + eqn)*dryPressure.to_value(cds.mbar)/temperature.to_value(units.Kelvin)
192 return densityFactor
193
194

◆ densityFactorWater()

lsst.afw.coord._refraction.densityFactorWater ( weather) 
Calculate water vapor pressure term to refractive index calculation.

Parameters
----------
weather : `lsst.afw.coord.Weather`
    Class containing the measured temperature, pressure, and humidity
    at the observatory during an observation

Returns
-------
densityFactor : `float`
    Returns the relative density of water vapor
    at the given pressure and temperature.

Notes
-----
This replicates equation 16 of [1]_

References
----------
.. [1] R. C. Stone, "An Accurate Method for Computing Atmospheric
   Refraction," Publications of the Astronomical Society of the Pacific,
   vol. 108, p. 1051, 1996.

Definition at line 195 of file _refraction.py.

195def densityFactorWater(weather):
196 """Calculate water vapor pressure term to refractive index calculation.
197
198 Parameters
199 ----------
200 weather : `lsst.afw.coord.Weather`
201 Class containing the measured temperature, pressure, and humidity
202 at the observatory during an observation
203
204 Returns
205 -------
206 densityFactor : `float`
207 Returns the relative density of water vapor
208 at the given pressure and temperature.
209
210 Notes
211 -----
212 This replicates equation 16 of [1]_
213
214 References
215 ----------
216 .. [1] R. C. Stone, "An Accurate Method for Computing Atmospheric
217 Refraction," Publications of the Astronomical Society of the Pacific,
218 vol. 108, p. 1051, 1996.
219 """
220 temperature = extractTemperature(weather, useKelvin=True)
221 waterVaporPressure = humidityToPressure(weather)
222 densityEqn1 = (-2.37321E-3 + 2.23366/temperature.to_value(units.Kelvin)
223 - 710.792/temperature.to_value(units.Kelvin)**2.
224 + 7.75141E-4/temperature.to_value(units.Kelvin)**3.)
225 densityEqn2 = waterVaporPressure.to_value(cds.mbar)*(1. + 3.7E-4*waterVaporPressure.to_value(cds.mbar))
226 relativeDensity = waterVaporPressure.to_value(cds.mbar)/temperature.to_value(units.Kelvin)
227 densityFactor = (1 + densityEqn2*densityEqn1)*relativeDensity
228
229 return densityFactor
230
231

◆ differentialRefraction()

lsst.afw.coord._refraction.differentialRefraction ( wavelength,
wavelengthRef,
elevation,
observatory,
weather = None ) 
Calculate the differential refraction between two wavelengths.

Parameters
----------
wavelength : `float`
    wavelength is in nm (valid for 230.2 < wavelength < 2058.6)
wavelengthRef : `float`
    Reference wavelength, typically the effective wavelength of a filter.
elevation : `lsst.geom.Angle`
    Elevation of the observation, as an Angle.
observatory : `lsst.afw.coord.Observatory`
    Class containing the longitude, latitude,
    and altitude of the observatory.
weather : `lsst.afw.coord.Weather`, optional
    Class containing the measured temperature, pressure, and humidity
    at the observatory during an observation
    If omitted, typical conditions for the observatory's elevation will be calculated.

Returns
-------
differentialRefraction : `lsst.geom.Angle`
    The refraction at `wavelength` minus the refraction at `wavelengthRef`.

Definition at line 94 of file _refraction.py.

94def differentialRefraction(wavelength, wavelengthRef, elevation, observatory, weather=None):
95 """Calculate the differential refraction between two wavelengths.
96
97 Parameters
98 ----------
99 wavelength : `float`
100 wavelength is in nm (valid for 230.2 < wavelength < 2058.6)
101 wavelengthRef : `float`
102 Reference wavelength, typically the effective wavelength of a filter.
103 elevation : `lsst.geom.Angle`
104 Elevation of the observation, as an Angle.
105 observatory : `lsst.afw.coord.Observatory`
106 Class containing the longitude, latitude,
107 and altitude of the observatory.
108 weather : `lsst.afw.coord.Weather`, optional
109 Class containing the measured temperature, pressure, and humidity
110 at the observatory during an observation
111 If omitted, typical conditions for the observatory's elevation will be calculated.
112
113 Returns
114 -------
115 differentialRefraction : `lsst.geom.Angle`
116 The refraction at `wavelength` minus the refraction at `wavelengthRef`.
117 """
118 refractionStart = refraction(wavelength, elevation, observatory, weather=weather)
119 refractionEnd = refraction(wavelengthRef, elevation, observatory, weather=weather)
120 return refractionStart - refractionEnd
121
122

◆ extractTemperature()

lsst.afw.coord._refraction.extractTemperature ( weather,
useKelvin = False ) 
Thin wrapper to return the measured temperature from an observation.

Parameters
----------
weather : `lsst.afw.coord.Weather`
    Class containing the measured temperature, pressure, and humidity
    at the observatory during an observation
useKelvin : bool, optional
    Set to True to return the temperature in Kelvin instead of Celsius
    This is needed because Astropy can't easily convert
    between Kelvin and Celsius.

Returns
-------
temperature : `astropy.units.Quantity`
    The temperature in Celsius, unless `useKelvin` is set.

Definition at line 269 of file _refraction.py.

269def extractTemperature(weather, useKelvin=False):
270 """Thin wrapper to return the measured temperature from an observation.
271
272 Parameters
273 ----------
274 weather : `lsst.afw.coord.Weather`
275 Class containing the measured temperature, pressure, and humidity
276 at the observatory during an observation
277 useKelvin : bool, optional
278 Set to True to return the temperature in Kelvin instead of Celsius
279 This is needed because Astropy can't easily convert
280 between Kelvin and Celsius.
281
282 Returns
283 -------
284 temperature : `astropy.units.Quantity`
285 The temperature in Celsius, unless `useKelvin` is set.
286 """
287 temperature = weather.getAirTemperature()*units.Celsius
288 if useKelvin:
289 temperature = temperature.to(units.Kelvin, equivalencies=units.temperature())
290 return temperature
291
292

◆ humidityToPressure()

lsst.afw.coord._refraction.humidityToPressure ( weather) 
Convert humidity and temperature to water vapor pressure.

Parameters
----------
weather : `lsst.afw.coord.Weather`
    Class containing the measured temperature, pressure, and humidity
    at the observatory during an observation

Returns
-------
pressure : `astropy.units.Quantity`
    The water vapor pressure in Pascals
    calculated from the given humidity and temperature.

Notes
-----
This replicates equations 18 & 20 of [1]_

References
----------
.. [1] R. C. Stone, "An Accurate Method for Computing Atmospheric
   Refraction," Publications of the Astronomical Society of the Pacific,
   vol. 108, p. 1051, 1996.

Definition at line 232 of file _refraction.py.

232def humidityToPressure(weather):
233 """Convert humidity and temperature to water vapor pressure.
234
235 Parameters
236 ----------
237 weather : `lsst.afw.coord.Weather`
238 Class containing the measured temperature, pressure, and humidity
239 at the observatory during an observation
240
241 Returns
242 -------
243 pressure : `astropy.units.Quantity`
244 The water vapor pressure in Pascals
245 calculated from the given humidity and temperature.
246
247 Notes
248 -----
249 This replicates equations 18 & 20 of [1]_
250
251 References
252 ----------
253 .. [1] R. C. Stone, "An Accurate Method for Computing Atmospheric
254 Refraction," Publications of the Astronomical Society of the Pacific,
255 vol. 108, p. 1051, 1996.
256 """
257 humidity = weather.getHumidity()
258 x = np.log(humidity/100.0)
259 temperature = extractTemperature(weather)
260 temperatureEqn1 = (temperature + 238.3*units.Celsius)*x + 17.2694*temperature
261 temperatureEqn2 = (temperature + 238.3*units.Celsius)*(17.2694 - x) - 17.2694*temperature
262 dewPoint = 238.3*temperatureEqn1/temperatureEqn2
263 waterVaporPressure = (4.50874 + 0.341724*dewPoint + 0.0106778*dewPoint**2 + 0.184889E-3*dewPoint**3
264 + 0.238294E-5*dewPoint**4 + 0.203447E-7*dewPoint**5)*133.32239*units.pascal
265
266 return waterVaporPressure
267
268

◆ refraction()

lsst.afw.coord._refraction.refraction ( wavelength,
elevation,
observatory,
weather = None ) 
Calculate overall refraction under atmospheric and observing conditions.

Parameters
----------
wavelength : `float`
    wavelength is in nm (valid for 230.2 < wavelength < 2058.6)
elevation : `lsst.geom.Angle`
    Elevation of the observation, as an Angle.
observatory : `lsst.afw.coord.Observatory`
    Class containing the longitude, latitude,
    and altitude of the observatory.
weather : `lsst.afw.coord.Weather`, optional
    Class containing the measured temperature, pressure, and humidity
    at the observatory during an observation
    If omitted, typical conditions for the observatory's elevation will be calculated.

Returns
-------
refraction : `lsst.geom.Angle`
    The angular refraction for light of the given wavelength,
    under the given observing conditions.

Notes
-----
The calculation is taken from [1]_.

References
----------
.. [1] R. C. Stone, "An Accurate Method for Computing Atmospheric
   Refraction," Publications of the Astronomical Society of the Pacific,
   vol. 108, p. 1051, 1996.

Definition at line 36 of file _refraction.py.

36def refraction(wavelength, elevation, observatory, weather=None):
37 """Calculate overall refraction under atmospheric and observing conditions.
38
39 Parameters
40 ----------
41 wavelength : `float`
42 wavelength is in nm (valid for 230.2 < wavelength < 2058.6)
43 elevation : `lsst.geom.Angle`
44 Elevation of the observation, as an Angle.
45 observatory : `lsst.afw.coord.Observatory`
46 Class containing the longitude, latitude,
47 and altitude of the observatory.
48 weather : `lsst.afw.coord.Weather`, optional
49 Class containing the measured temperature, pressure, and humidity
50 at the observatory during an observation
51 If omitted, typical conditions for the observatory's elevation will be calculated.
52
53 Returns
54 -------
55 refraction : `lsst.geom.Angle`
56 The angular refraction for light of the given wavelength,
57 under the given observing conditions.
58
59 Notes
60 -----
61 The calculation is taken from [1]_.
62
63 References
64 ----------
65 .. [1] R. C. Stone, "An Accurate Method for Computing Atmospheric
66 Refraction," Publications of the Astronomical Society of the Pacific,
67 vol. 108, p. 1051, 1996.
68 """
69 if wavelength < 230.2:
70 raise ValueError("Refraction calculation is valid for wavelengths between 230.2 and 2058.6 nm.")
71 if wavelength > 2058.6:
72 raise ValueError("Refraction calculation is valid for wavelengths between 230.2 and 2058.6 nm.")
73 latitude = observatory.getLatitude()
74 altitude = observatory.getElevation()
75 if weather is None:
76 weather = defaultWeather(altitude*units.meter)
77 reducedN = deltaN(wavelength, weather)/deltaRefractScale
78 temperature = extractTemperature(weather, useKelvin=True)
79 atmosScaleheightRatio = 4.5908E-6*temperature.to_value(units.Kelvin)
80
81 # Account for oblate Earth
82 # This replicates equation 10 of Stone 1996
83 relativeGravity = (1. + 0.005302*np.sin(latitude.asRadians())**2.
84 - 0.00000583*np.sin(2.*latitude.asRadians())**2. - 0.000000315*altitude)
85
86 # Calculate the tangent of the zenith angle.
87 tanZ = np.tan(np.pi/2. - elevation.asRadians())
88 atmosTerm1 = reducedN*relativeGravity*(1. - atmosScaleheightRatio)
89 atmosTerm2 = reducedN*relativeGravity*(atmosScaleheightRatio - reducedN/2.)
90 result = float(atmosTerm1*tanZ + atmosTerm2*tanZ**3.)*lsst.geom.radians
91 return result
92
93

Variable Documentation

◆ deltaRefractScale

float lsst.afw.coord._refraction.deltaRefractScale = 1.0E8

Definition at line 33 of file _refraction.py.

Generated on Wed Nov 20 2024 09:02:25 for LSST Applications by doxygen 1.10.0