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

Home » Public Forums » archive » Re: Map point
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Switch to threaded view of this topic Create a new topic Submit Reply
Re: Map point [message #25425] Mon, 18 June 2001 05:37
btt is currently offline  btt
Messages: 345
Registered: December 2000
Senior Member
Hello,

>
> and now I want to have a point or a cross in the place of some
> important spot like (42.938,-0.6)
>

The Map_Set procedure has set up the coordinate conversions for you.

Now simply plot a point with the PLOTS command:

PlotS, -0.06, 42.938, PSYM = 1, ....

>
> Also I have a 2d array with value of concentration that I want to plot
> on the map and the array have the same (lat-lon) dimension of the map,
> how can I do it?
>

You have a couple of choices, you could contour the array (CONTOUR
command) or use Map_Patch (or Map_Image to warp the image properly to
your projection then use TV (or TVImage or ImDisp, etc.) to display it
on you 'map'.

Ben

--
Ben Tupper
Bigelow Laboratory for Ocean Sciences
180 McKown Point Rd.
W. Boothbay Harbor, ME 04575
btupper@bigelow.org
Re: Map point [message #25434 is a reply to message #25425] Fri, 15 June 2001 19:52 Go to previous message
david[2] is currently offline  david[2]
Messages: 100
Registered: June 2001
Senior Member
Francisco writes:

> I have produced a map in this way
>
> latmin=39.74 & latmax=45.63
> lonmin=-5.32 & lonmax=4.79
>
> latcenter=42.8 & loncenter=-0.517
>
> map_set,latcenter,loncenter,/lambert,/hires,
> color=0,limit=[latmin,lonmin,latmax,lonmax]
>
> map_continents,/countries
>
> and now I want to have a point or a cross in the place of some
> important spot like (42.938,-0.6)
>
> How can I do it?

PLOTS, 42.939, -0.6, PSYM=1, SymSize=3.0

> Also I have a 2d array with value of concentration that I want to plot
> on the map and the array have the same (lat-lon) dimension of the map,
> how can I do it?

You could try warping the 2D array on to the image
with either Map_Image or Map_Patch. Keywords allow
you to specify the latitude and longitude coordinates
of the edges of the image.

warped = Map_Image(array, xx, yy, xs, ys, LonMin=-5.32, $
LonMax = 4.79, LatMin=39.74, LatMax=45.63)
TV, warped, xx, yy, XSize=xs, YSize=ys
Map_Continents, /Countries

Or, you could put the 2D array over the image with
the CONTOUR command. Just be sure to specify latitude
and longitude vectors along with the array, and use
the OVERPLOT keyword.

Best Regards,

David

--
David Fanning, Ph.D.
Fanning Software Consulting
Phone: 970-221-0438 E-Mail: davidf@dfanning.com
Coyote's Guide to IDL Programming: http://www.dfanning.com/
Toll-Free IDL Book Orders: 1-888-461-0155
  Switch to threaded view of this topic Create a new topic Submit Reply
Previous Topic: Re: regression with poisson models
Next Topic: Re: IDL and UNIX setenv

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

Current Time: Wed Oct 08 13:26:24 PDT 2025

Total time taken to generate the page: 0.00447 seconds