| Re: Latitude longitude and Image Navigation [message #76744 is a reply to message #76706] |
Sun, 26 June 2011 21:00   |
AISHWARYA
Messages: 9
Registered: June 2011
|
Junior Member |
|
|
On Jun 25, 9:03 am, Brian Wolven <brian.wol...@gmail.com> wrote:
> Well, *if* I remember correctly what I did there: Subcur is the vector of the initial sub-observer point on the planet/satellite before the rotations, given a coordinate system defined with respect to the center of that body. It's just the normalized planetocentric coordinates of the vector pointing toward the observer:
>
> subcur = -1.*body.pos/length(body.pos)
>
> The -1 comes in because body.pos contains the vector from the observer location (Earth in this case) to the body. subdes is the desired sub-observer location - the one specified in theimagemetadata or by SPICE, based on the observation time. We need to rotate our body centered coordinates such that the z-axis aligns with the specified north pole position vector, then rotate the body about the pole until our line of sight is passing through the proper sub-observer point. The various 'ang' variables were only used for diagnostic purposes to check that the desired result was being obtained.
>
> To put it in more graphic terms, imagine that you have a globe of Venus sitting upright on your desk. In order to match your telescopic observations you need to tilt the globe until the pole axis is aligned properly, then rotate about the pole until you reach the correct observing geometry. We have a body centered lat-long grid instead of a globe, but once we know the rotation we can then project our rotated lat-long grid on theimageand have it properly aligned, or infer the lat-long coordinates of a pixel in theimageif we're working the other way.
>
> There may very well be an easier way to do this, but that's what I cooked up years ago as a grad student, when presumably I was smarter than I am now. ;)
I could understand the code partly. But, what is the role of z axis ?
This is a 2 dimensional image. I tried rotating the image with centre
as the pivot to 23 degree (the north pole position angle) and then
used mapset with lat long as the sub observer lat long (-5.47, -109)
with limints [-90,0,-180,180] since the image occupies southern
hemisphere. But, I'm getting an error saying latitude is out of range.
Please help.
Thank you,
Aishwarya.
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
