<p>Picture was taken somewhere near 30°22'02.0&quot;N 75°00'35.7&quot;E
Bathinda, Punjab, India.</p>
<p>Answering your question about time - how it enters into the sight reduction: RA, all by itself, isn't useful. It's the angle between the plane of the meridian of a celestial location and the plane of the meridian of the first point of Aries. That angle is more or less fixed, but the first point of Aries, and its meridian, are continually changing, going through 360 degrees in a sidereal day.</p>
<p>Assumptions: pitch is the altitude of the celestial location (its angular distance above the horizon), and it's not corrected for refraction. Yaw is the azimuth of the location, CW degrees from north.</p>
<p>General method: knowing the time, find the Greenwich hour angle (GHA) of the first point of Aries at that time, by interpolating between hourly values in an almanac. The GHA of the celestial location is the GHA of Aries minus RA. The geographical position (GP) of a celestial location, is the point on the earth's surface where the celestial location is in the zenith. The GP I found after some number crunching was 3.56 deg N, 72.655 deg E. It's in the middle of the Arabian Sea. That GP, the observer's unknown location, and the north pole form a spherical triangle, which must be solved to get the observer's location fix. Two sides and an angle are known: one of the sides is 90 deg minus the latitude of the GP, the other is 90 deg minus the observed altitude of the celestial location. The angle is the 360 deg minus the azimuth of the location. It's one of those awkward spherical triangles where the angle is not the included angle between the sides. It's solved by first using the law of sines to find the angle at the pole then using Napier's Analogies to find the third side.</p>
<p>In this particular case, the celestial location GP is almost due south of the observer, so there's a way to make a simple estimate. The observer is roughly 26.9 deg north of the GP, a little less than 30.46 deg N, and a little east of the GP 72.65 E.</p>
<p>I have a few more things to say if anyone's interested. I'd appreciate a check on my numbers. It's very easy to make a mistake, which is why sight reductions in the pre-GPS era were done by two persons independently who'd compare their results. An unintentional grounding can ruin your whole day.</p>

CC BY-SA 4.0
2020-08-04T22:37:02.687Z
@Mike G: thanks. I was waiting for someone to show some interest.
stretch
CC BY-SA 4.0
2020-08-04T17:20:50.393Z
2020-08-04T22:33:00.753Z
2020-08-04T22:33:00.753Z
stretch
<p>Given a picture of stars - assume it is possible to identify the stars and their Right Ascension (RA) and Declination (Dec).  By choosing one near the center of the image (or interpolating RA/Dec to the center) the effective "pointing angle" of the camera would be known.</p>

<p>With the RA/Dec of the camera, plus the Roll (R), Pitch (P), and Yaw (Y) of the camera, and the precise time (TAI, UTC, and/or UT1) the picture was taken... how would a location calculation of where the camera was on Earth when it was taken be possible (assume it is possible).</p>

<p>It seems like it would be a simple coordinate transformation, however understanding how time affects the calculation isn't clear.  I can't find a simple resource of algorithms to find one close enough.</p>

<p>ETA: Roll, Pitch, and Yaw could be found in a few different ways... in my current case I took a good first initial "guess" and then plotted an overlay of stars based on perturbing the initial first guess until the stars overlaid well over the image (of course, this required I knew where the picture was taken!)  However, it is also possible to know R,P,Y based on a local gravity vector and inertial mechanics (a simply INU)... which is where I would eventually like to get to.</p>

<p>The RA and Dec is currently found using <a href="http://astrometry.net/" rel="nofollow noreferrer">http://astrometry.net/</a> by uploading the picture they returned astrometric meda-data.</p>

<p>Here is an example set of numbers:
R=10.7deg, P=63.1deg, Y=-174.8deg</p>

<p>Time=2016,Oct,1,0hr,31m,28s</p>

<p>Leap Secs=36 GPS-TAI</p>

<p>UT1-UTC = -0.28</p>

<p>Center of Image (~pointing vector of camera)</p>

<p>RA = 90.7 (6h,2m,53.422s), </p>

<p>Dec=3.56(3deg, 33' 24.95)</p>

CC BY-SA 4.0
2020-02-27T04:28:21.230Z
How can a camera have all five of these defined?  (i.e. R.A., Dec., Pitch, Roll, and Yaw?) I think providing a specific example with numbers will be helpful. I'm also wondering if azimuth and elevation might be better,  "Roll" come from "which way is north" in the image using the stars identified which give R.A. and Dec.
uhoh
CC BY-SA 4.0
2020-02-27T08:14:05.193Z
Does this answer your question? [Is it possible to determine an observer's position on Earth from a photo of the night sky and the time of the shot?](https://astronomy.stackexchange.com/questions/33150/is-it-possible-to-determine-an-observers-position-on-earth-from-a-photo-of-the)
usernumber
CC BY-SA 4.0
2020-02-27T09:04:01.343Z
@usernumber the answer there does not answer this question so it's not a duplicate. This question asks "How to calculate..." and there's nothing like that there. *That question* asks "Is it possible..." and the answer is "Yes, it is possible" but nothing more. Now if you had voted to close as requiring clarification that would have been a different matter. But closing a new user's first question for the wrong reason isn't being very hospitable.
uhoh
CC BY-SA 4.0
2020-02-27T23:31:40.267Z
`+1` for the edit, thanks! I don't yet understand how pitch roll and yaw are defined for a camera and how to convert them to altitude and azimuth, but that's just me; perhaps someone else will.
uhoh
CC BY-SA 4.0
2020-02-28T18:50:26.420Z
@uhoh Typically RPY are Euler angles from a local navigation reference frame, North East Down (NED).  Similar to an aircraft, but for cameras the z direction is typically forward, X out the right and y down. barrycenter - I'm not sure knowing how the sky will look like at 0 RPY will help, but interested in that path!  I was anticipating a triangulation of sorts (or vector addition really), knowing 2 vectors - the center of the earth to the star closest to the center of the image and a vector from the camera to the same, and solve for the vector from the center of the earth to the camera.
Biaspoint
CC BY-SA 4.0
2020-02-28T23:48:37.037Z
@Biaspoint where does 0 RPY point, If I say "I point my camera forward" Is it North along the horizon? If so, then adding that back to the question would be great. *However* since this is asked in Astronomy SE your question might fit in better if you use azimuth, elevation, and field rotation ([1](https://themcdonalds.net/richard/wp/astrophotography-field-rotation/), [2](https://telescopemount.org/astrophotography-with-alt-az-telescope-mounts/)) rather than require answers to address one of the [six possible Euler angle conventions](https://en.wikipedia.org/wiki/Euler_angles#Conventions).
uhoh
CC BY-SA 4.0
2020-03-02T13:25:06.140Z
@uhoh - thanks, I have some mathwork to do!  I believe 0 RPY is pointing north horizontally.
Biaspoint
<photography><celestial-mechanics><navigation>
How to calculate local position, knowing the RA and Dec of the center of an image, RPY of the camera, and precise time?
Biaspoint
# How to calculate local position, knowing the RA and Dec of the center of an image, RPY of the camera, and precise time?

(Asked by: username_0 on 2020-02-27)

Given a picture of stars - assume it is possible to identify the stars and their Right Ascension (RA) and Declination (Dec). By choosing one near the center of the image (or interpolating RA/Dec to the center) the effective "pointing angle" of the camera would be known.

With the RA/Dec of the camera, plus the Roll (R), Pitch (P), and Yaw (Y) of the camera, and the precise time (TAI, UTC, and/or UT1) the picture was taken... how would a location calculation of where the camera was on Earth when it was taken be possible (assume it is possible).

It seems like it would be a simple coordinate transformation, however understanding how time affects the calculation isn't clear. I can't find a simple resource of algorithms to find one close enough.

ETA: Roll, Pitch, and Yaw could be found in a few different ways... in my current case I took a good first initial "guess" and then plotted an overlay of stars based on perturbing the initial first guess until the stars overlaid well over the image (of course, this required I knew where the picture was taken!) However, it is also possible to know R,P,Y based on a local gravity vector and inertial mechanics (a simply INU)... which is where I would eventually like to get to.

The RA and Dec is currently found using <http://astrometry.net/> by uploading the picture they returned astrometric meda-data.

Here is an example set of numbers:
R=10.7deg, P=63.1deg, Y=-174.8deg

Time=2016,Oct,1,0hr,31m,28s

Leap Secs=36 GPS-TAI

UT1-UTC = -0.28

Center of Image (~pointing vector of camera)

RA = 90.7 (6h,2m,53.422s),

Dec=3.56(3deg, 33' 24.95)

**Tags:** photography, celestial-mechanics, navigation
> **Comments:**
> **username_1**: How can a camera have all five of these defined?  (i.e. R.A., Dec., Pitch, Roll, and Yaw?) I think providing a specific example with numbers will be helpful. I'm also wondering if azimuth and elevation might be better,  "Roll" come from "which way is north" in the image using the stars identified which give R.A. and Dec.  
> **username_2**: Does this answer your question? [Is it possible to determine an observer's position on Earth from a photo of the night sky and the time of the shot?](https://astronomy.stackexchange.com/questions/33150/is-it-possible-to-determine-an-observers-position-on-earth-from-a-photo-of-the)  
> **username_1**: @usernumber the answer there does not answer this question so it's not a duplicate. This question asks "How to calculate..." and there's nothing like that there. *That question* asks "Is it possible..." and the answer is "Yes, it is possible" but nothing more. Now if you had voted to close as requiring clarification that would have been a different matter. But closing a new user's first question for the wrong reason isn't being very hospitable.  
> **username_1**: `+1` for the edit, thanks! I don't yet understand how pitch roll and yaw are defined for a camera and how to convert them to altitude and azimuth, but that's just me; perhaps someone else will.  
> **username_0**: @uhoh Typically RPY are Euler angles from a local navigation reference frame, North East Down (NED).  Similar to an aircraft, but for cameras the z direction is typically forward, X out the right and y down. barrycenter - I'm not sure knowing how the sky will look like at 0 RPY will help, but interested in that path!  I was anticipating a triangulation of sorts (or vector addition really), knowing 2 vectors - the center of the earth to the star closest to the center of the image and a vector from the camera to the same, and solve for the vector from the center of the earth to the camera.  
> **username_1**: @Biaspoint where does 0 RPY point, If I say "I point my camera forward" Is it North along the horizon? If so, then adding that back to the question would be great. *However* since this is asked in Astronomy SE your question might fit in better if you use azimuth, elevation, and field rotation ([1](https://themcdonalds.net/richard/wp/astrophotography-field-rotation/), [2](https://telescopemount.org/astrophotography-with-alt-az-telescope-mounts/)) rather than require answers to address one of the [six possible Euler angle conventions](https://en.wikipedia.org/wiki/Euler_angles#Conventions).  
> **username_0**: @uhoh - thanks, I have some mathwork to do!  I believe 0 RPY is pointing north horizontally.  

---

## 1 answers

### Answer by username_3 on 2020-08-04. Score: 2

Picture was taken somewhere near 30°22'02.0"N 75°00'35.7"E
Bathinda, Punjab, India.

Answering your question about time - how it enters into the sight reduction: RA, all by itself, isn't useful. It's the angle between the plane of the meridian of a celestial location and the plane of the meridian of the first point of Aries. That angle is more or less fixed, but the first point of Aries, and its meridian, are continually changing, going through 360 degrees in a sidereal day.

Assumptions: pitch is the altitude of the celestial location (its angular distance above the horizon), and it's not corrected for refraction. Yaw is the azimuth of the location, CW degrees from north.

General method: knowing the time, find the Greenwich hour angle (GHA) of the first point of Aries at that time, by interpolating between hourly values in an almanac. The GHA of the celestial location is the GHA of Aries minus RA. The geographical position (GP) of a celestial location, is the point on the earth's surface where the celestial location is in the zenith. The GP I found after some number crunching was 3.56 deg N, 72.655 deg E. It's in the middle of the Arabian Sea. That GP, the observer's unknown location, and the north pole form a spherical triangle, which must be solved to get the observer's location fix. Two sides and an angle are known: one of the sides is 90 deg minus the latitude of the GP, the other is 90 deg minus the observed altitude of the celestial location. The angle is the 360 deg minus the azimuth of the location. It's one of those awkward spherical triangles where the angle is not the included angle between the sides. It's solved by first using the law of sines to find the angle at the pole then using Napier's Analogies to find the third side.

In this particular case, the celestial location GP is almost due south of the observer, so there's a way to make a simple estimate. The observer is roughly 26.9 deg north of the GP, a little less than 30.46 deg N, and a little east of the GP 72.65 E.

I have a few more things to say if anyone's interested. I'd appreciate a check on my numbers. It's very easy to make a mistake, which is why sight reductions in the pre-GPS era were done by two persons independently who'd compare their results. An unintentional grounding can ruin your whole day.

> **Comments:**
> **username_3**: @Mike G: thanks. I was waiting for someone to show some interest.  

---

