Re: Polar-projection multiple plot with image+contours [message #46810] |
Wed, 04 January 2006 12:45  |
David Fanning
Messages: 11724 Registered: August 2001
|
Senior Member |
|
|
Emre Ishik writes:
> I'm in trouble with creating a multiple (even single) plot which is
> warped onto stereographic polar projection.
>
> ; Initially I have a rectangular array
> A FLOAT = Array[360, 180]
> ; and I set the mapping as
> map_set, /stereographic, /isotropic, 90, latdel=20, londel=20,
> /horizon, /grid
> ; Then I map the image A
> bScala = BYTSCL(a,minScala,maxScala) ; rescale the intensity
> mapped_imagea = MAP_IMAGE(bscala, compress=1)
>
> What I want to do is to plot some (say two) images A, B into a
> multiple- plot PS-output by setting !P.MULTI = [0,2,1] in this case. I
> use the routine TVSCL. Then I want to overplot contours onto each
> pole-on image such that the levels are the constants minScala and
> maxScala. Even with a single plot, the contour doesn't coincide with
> the patterns on the image. For multiple plots with !P.MULTI, things are
> even worse. Isn't there any easy way? I tried as much ways as I can..
Humm. Yes, quite a few things wrong here. :-)
I'd start by reading these articles:
http://www.dfanning.com/tips/map_pmulti.html
http://www.dfanning.com/tips/map_image.html
http://www.dfanning.com/color_tips/fill_colors.html
Digest the information there, then let's see where you are.
> Another minor problem is that the grids that I set in map_set are not
> seen on the TVSCL outputs of the mapped image.
Redraw the grids, with MAP_GRID, after you have done everything else.
Cheers,
David
--
David Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming: http://www.dfanning.com/
|
|
|
Re: Polar-projection multiple plot with image+contours [message #46861 is a reply to message #46810] |
Mon, 09 January 2006 11:08  |
Emre Ishik
Messages: 1 Registered: January 2006
|
Junior Member |
|
|
Thanks a lot for your help.. I've just removed my post because I
received many spams due to my e-mail address explicit on top of the
message (nickname); now at least I changed my nickname. Maybe it was
useless, but anyway, it's still quoted in your post. (oh, I've just
read the remark on top of this page: "The group you are posting to is a
Usenet group. Messages posted to this group will make your email
visible to anyone on the Internet." So, no way of escaping spams..)
I've read the your articles. Those have given me insight, thank you..
Finally I've accomplished it by forgetting about displaying the image,
and considering contours with /cell_fill option. What I did is just
"map_set" and then "contour" directly, with choosing the levels and the
corresponding colors from the color table which were loaded by
"loadct". Everything works well with multiple plots, too.
Greetings,
Emre Ishik
MPl for Solar System Research, Germany
|
|
|