Just trying it out, using pg_plotimage
http://hea-www.harvard.edu/~pgrigis/idl_stuff/pg_plotimage.p ro
;linear test array, easier to deal with than dist
a=[[0.,1,2,1,0],[1,2,3,2,1],[2,3,4,3,2],[1,2,3,2,1],[0,1,2,1 ,0]]
loadct,5
pg_plotimage,a,findgen(5),findgen(5)
contour,a,/over,levels=[0.5,1.5,2.5,3.5],color=0
(The contours should have a vertex exactly at the pixel
boundaries by construction of the image)
The half pixel at the boundary is missing however.
Paolo
David Fanning wrote:
> David Fanning writes:
>
>> Wayne Landsman writes:
>>
>>> Does anyone have a routine to overlay a contour plot on a small image
>>> that is both accurate and aesthetically pleasing?
>>
>> Isn't this what you want:
>>
>> im = dist(9)
>> thispostion = [0.1, 0.1, 0.9, 0.9]
>> TVIMAGE, bytscl(im), POSITION=thisPosition, /KEEP_ASPECT_RATIO, $
>> /NOINTERP, /ERASE
>> contour,im, Position=thisPosition,/xsty,/ysty,/noerase,$
>> scale_vector(findgen(9), 0, 9), scale_vector(findgen(9), 0, 9)
>
> I don't know. What do you think of this:
>
> im = dist(9)
> im[2,2] = 15
> im[4:5, 4:5] = 10
> thispostion = [0.1, 0.1, 0.9, 0.9]
> TVIMAGE, bytscl(im), POSITION=thisPosition, /KEEP_ASPECT_RATIO, $
> /NOINTERP, /ERASE
> contour,im, Position=thisPosition,/xsty,/ysty,/noerase,$
> scale_vector(findgen(9), 0, 9), scale_vector(findgen(9), 0,9), $
> level = 6
>
> isocontour, im, v, c, c_value = 6
> Plots, v[0,0:6]+0.5, v[1,0:6]+0.5, color=fsc_color('dodger blue')
> Plots, v[0,7:*]+0.5, v[1,7:*]+0.5, color=fsc_color('dodger blue')
>
>
> Note how I have to add a half pixel to the ISOCONTOUR to get
> it anywhere near where I want it. Weird.
>
> Cheers,
>
> David
> --
> David Fanning, Ph.D.
> Fanning Software Consulting, Inc.
> Coyote's Guide to IDL Programming: http://www.dfanning.com/
> Sepore ma de ni thui. ("Perhaps thou speakest truth.")
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
