comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: create an UTM grid
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: create an UTM grid [message #78411 is a reply to message #78410] Mon, 28 November 2011 12:02 Go to previous messageGo to previous message
natha is currently offline  natha
Messages: 482
Registered: October 2007
Senior Member
I think the following code is better :

x_grid=REBIN(x_grid,xdim,ydim,/SAMPLE)
y_grid=TRANSPOSE(REBIN(y_grid,ydim,xdim,/SAMPLE))

result=MAP_PROJ_INVERSE(x_grid, y_grid, MAP_STRUCTURE=map_utm)

res_lon=REFORM(result[0,*],xdim,ydim)
res_lat=REFORM(result[1,*],xdim,ydim)

By the way, I am missing something here...

UTM is a cartesian coordinate system. Theoretically, the distance
between two points on the grid could be performed using the
Pythagorean theorem (because all distances should be equal). What I
realize is that the distances are different and never equal to 1km. I
checked that using the routine MAP_2POINTS. Do you know why ?
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Re: Writing arrays to text file - format code trickery?
Next Topic: Search single column of array - removing nasty loop

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Fri Oct 10 21:35:24 PDT 2025

Total time taken to generate the page: 0.48004 seconds