comp.lang.idl-pvwave archive: archive » Re: create an UTM grid

Home » Public Forums » archive » Re: create an UTM grid

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

Re: create an UTM grid [message #78413 is a reply to message #78411]

Mon, 28 November 2011 11:25

natha
Messages: 482
Registered: October 2007

Senior Member

Thank you David. Just to be sure, if I want a grid of 1000x1000, 1km
resolution and centered on 45.42419lat -73.9372lon, is the following
code right ?

center_lat=45.42419
center_lon=-73.9372
xdim=1000
ydim=1000
resolution=1000.

map_utm=MAP_PROJ_INIT('UTM', /GCTP, CENTER_LON=center_lon,
CENTER_LAT=center_lat, ELLIPSOID=24)

xycenter=MAP_PROJ_FORWARD(center_lon, center_lat,
MAP_STRUCTURE=map_utm)

x_grid=xycenter[0] + FINDGEN(xdim)*resolution - xdim/2.*resolution
y_grid=xycenter[1] + FINDGEN(ydim)*resolution - ydim/2.*resolution

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

lon=REFORM(result[0,*])
lat=REFORM(result[1,*])

Thank you again ! I really appreciate your support !
nata

Report message to a moderator

[Message index]

		Re: create an UTM grid By: David Fanning on Mon, 28 November 2011 14:35
		Re: create an UTM grid By: David Fanning on Mon, 28 November 2011 14:33
		Re: create an UTM grid By: natha on Mon, 28 November 2011 12:02
		Re: create an UTM grid By: natha on Mon, 28 November 2011 11:25
		Re: create an UTM grid By: David Fanning on Mon, 28 November 2011 10:12
		Re: create an UTM grid By: David Fanning on Mon, 28 November 2011 14:38

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) ] [

] [

]

Current Time: Sat Oct 11 05:47:39 PDT 2025

Total time taken to generate the page: 0.86718 seconds