comp.lang.idl-pvwave archive: archive » map

Home » Public Forums » archive » map_set stereographic projection

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: map_set stereographic projection [message #48582 is a reply to message #48520]

Fri, 28 April 2006 12:41

mattie
Messages: 4
Registered: January 2006

Junior Member

kuyper@wizard.net writes:

> mattie wrote:
>> kuyper@wizard.net writes:
>>
>>> dvila wrote:
> ...
>>> ; These are the u-v values corresponding to pixel centers along each
>>> edge.
>>> topv = (uv[1,1]+uv[1,2])*0.5
>>> botv = (uv[1,0]+uv[1,3])*0.5
>>> leftu = (uv[0,0]+uv[0,1])*0.5
>>> rightu = (uv[0,2]+uv[0,3])*0.5
>>
>>
>> Aren't these the u-v values corresponding to the _center_ of the
>> gridcell along each edge? (assuming the initial data gave the
>> centerpoints of the gridcell.)
>
> Yes, that's what my comment line says.

Ah, so it does. I read it as the center point of each edge.

>>> ; U-V coordinates of midpoints of outer edges
>>> u = [leftu-0.5*xscale, 0.5*(leftu+rightu), rightu+0.5*xscale, $
>>> 0.5*(leftu+rightu)]
>>> v = [0.5*(botv+topv), topv+0.5*yscale, 0.5*(botv+topv),
>>> botv-0.5*yscale]
>>> lonlat = MAP_PROJ_INVERSE(u, v, MAP_STRUCTURE=stereo)
>>> limit = [lonlat[1,*],lonlat[0,*]]
>>
>> Again, I'm not sure, but don't you have to add half a gridcell to each
>> direction to get the outter limit of each grid cell?
>
> Yes, that is precisely why I wrote such things as "leftu-0.5*xscale",

Yup again. When I looked at it, I kept thinking u was upper, rather
than the
uv direction, and it's completely obvious now. I was getting thrown by
the
u[1] factor as (leftu + rightu) * .5.

So I don't have any reason to say "But" when I said:

"But this was a very informative tutorial on how to register images
when
corner points are known."

Thanks for the lesson.

Matt

--
Matthew Savoie - Scientific Programmer
National Snow and Ice Data Center
(303) 735-0785 http://nsidc.org

Report message to a moderator

[Message index]

		map_set stereographic projection By: dvila on Thu, 27 April 2006 05:50
		Re: map_set stereographic projection By: dvila on Mon, 01 May 2006 08:12
		Re: map_set stereographic projection By: James Kuyper on Fri, 28 April 2006 21:18
		Re: map_set stereographic projection By: David Fanning on Fri, 28 April 2006 16:06
		Re: map_set stereographic projection By: James Kuyper on Fri, 28 April 2006 15:59
		Re: map_set stereographic projection By: David Fanning on Fri, 28 April 2006 12:52
		Re: map_set stereographic projection By: mattie on Fri, 28 April 2006 12:41
		Re: map_set stereographic projection By: James Kuyper on Fri, 28 April 2006 10:24

Previous Topic:	Re: How to rotate a figure in PLOT?
Next Topic:	Upgrade to IDL 6.3 breaks IDLtoAVI

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sat Nov 29 11:00:46 PST 2025

Total time taken to generate the page: 0.00215 seconds