| 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/
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
