comp.lang.idl-pvwave archive: archive » Using TV command instead of contour

Home » Public Forums » archive » Using TV command instead of contour

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

Re: Using TV command instead of contour [message #76947 is a reply to message #76865]

Wed, 13 July 2011 08:21

KRDean
Messages: 69
Registered: July 2006

Member

On Jul 12, 4:18 am, Amit <amisra.a...@gmail.com> wrote:
> Dear Friends,
>
> I have a data in form of flags, for example, dust=1,smoke=2 and so on.
> There cannot be numbers in between.
>
> The data is in the form of 2-dimensional array, with each cell
> corresponding to a latitude, longitude pair. I want to display it --
> the output would look something like a chess board. If I use the
> CONTOUR command, it interpolates at the boundaries, producing
> incorrect plots.
> This problem is solved if I use TVIMAGE program. However, I cannot
> pass the latitude,longitude data to the command (like CONTOUR). Can
> anyone suggest me how to use TVIMAGE or TV in a way similar to
> CONTOUR?
>
> Thanks for your help
> Amit

I normally use PLOTS to display index values, like landuse data were
certain numbers represent regions. For exmaple, 10:forest and
40:grass.

Below is an example how I use PLOTS.

Kelly Dean
Milliken, Colorado

;+
;
;---------------------------------------------
PRO IndexPlot

file = FILEPATH( 'worldelv.dat', SUBDIRECTORY = ['examples', 'data'])
image = READ_BINARY(file, DATA_DIMS = [360,360])

siz = SIZE( image )
LatArr = REPLICATE( 1.0, siz[1] ) # ( ( LINDGEN( siz[2] ) * 0.5 ) -
89.5 )
LngArr = ( LINDGEN( siz[1] ) * 1.0 ) # REPLICATE (1.0, siz[2] )

west = WHERE( LngArr GT 180.0, cnt )
IF ( cnt GT 0 ) THEN LngArr[west] = LngArr[west] - 360.0

; The user defined plotting symbol ( PSYM=8 ) is used by PlotS to
; add a symbol on the plot. So the IDL routine USERSYM is used to
; create a filled symbol of a circle.

A = FINDGEN(17) * (!PI*2/16.)
USERSYM, COS(A), SIN(A), /FILL

;image = ROTATE( image,7 )

DEVICE, decompose=0
WINDOW, 0, XSIZE=siz[1], YSIZE=siz[2], TITLE='Index Plot'

LOADCT, 27 ; Assign a color table

MAP_SET, 0.0, 180.0, 0.0, /Mercator

PlotS, LngArr, LatArr, COLOR = image, $
PSYM = 8, SYMsize = 0.8, /DATA

MAP_CONTINENTS

END

Report message to a moderator

[Message index]

		Using TV command instead of contour By: Amit on Tue, 12 July 2011 03:18
		Re: Using TV command instead of contour By: Amit on Thu, 14 July 2011 23:58
		Re: Using TV command instead of contour By: pgrigis on Thu, 14 July 2011 14:06
		Re: Using TV command instead of contour By: Amit on Wed, 13 July 2011 12:25
		Re: Using TV command instead of contour By: Amit on Wed, 13 July 2011 12:22
		Re: Using TV command instead of contour By: pgrigis on Wed, 13 July 2011 08:55
		Re: Using TV command instead of contour By: KRDean on Wed, 13 July 2011 08:21
		Re: Using TV command instead of contour By: Amit on Wed, 13 July 2011 05:21
		Re: Using TV command instead of contour By: pgrigis on Tue, 12 July 2011 12:01
		Re: Using TV command instead of contour By: Brian Wolven on Tue, 12 July 2011 11:08

Previous Topic:	Re: GONG SOLAR DATA ANALYSIS IN IDL
Next Topic:	Multiple-plot multiple-page landscape postscript problem

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Wed Oct 08 19:36:21 PDT 2025

Total time taken to generate the page: 0.00438 seconds