1 winter quarter 2003rolando v. raqueño basic idl commands simg 726

Winter Quarter 2003 Rolando V. Raqueño 11

Basic IDL Commands

SIMG 726

Objective

• Quick Tour of IDL

• Writing IDL Programs

Before we start

% source ~rvrpci/.simg726.rc

% script% set history=100

AT THE END OF THE UNIX SESSION,DON’T FORGET TO TYPE

Before we END

AT THE END OF THE UNIX SESSION,DON’T FORGET TO TYPE

You now have file called typescript which contains a log of your UNIX session

Starting IDL in UNIX

• To start command line interface IDL environment% idl

• To start graphical interface IDL environment% idlde

IDL Help

• MOST IMPORTANT IDL COMMAND OF ALL

IDL> ?

IDL Help Window

Other IDL Help Facilities

IDL> doc_library,’moment’

• Allows specially formatted code comments to be viewed without having to view the source code

General IDL Information

• N.B. Platform specific limitation – UNIX

• File name conflict between lower and upper case characters

– e.g. addnum.pro vs. AddNum.pro

• Recommend use of strictly lowercase file names

IDL Character Set• IDL is not case sensitive• Commands and arguments are delimited by

comas.– e.g. IDL> print, a

• Comments are started by a semicolon (;)– e.g. IDL> ; This is a comment

• Command options are set using a backslash (/)– e.g. IDL> tv, a, /order

Other IDL Characters

• System variables are started with an exclamation point (!)e.g. IDL> print, !PIe.g. IDL> y = sin( x*!DTOR )

• Command line continuation with $e.g. IDL> print, 180 / $

Other IDL Characters

• Multiple IDL commands can be separated by the & symbol– e.g. IDL> a=0.0 & b=1.0

• UNIX Command execution– e.g. IDL> $ls -la

Journal Facility of IDL

• The journal command allows you to record all IDL activity during a particular session

IDL> journal, ‘sequence.pro’

IDL> print, 3*5

IDL> journal

sequence.pro

; IDL Version 3.6.1a (MacOS Macintosh); Journal File for Rolando Raqueño@green; Working directory: Raqueño's Hard Drive:Applications:ENVI:; Date: Tue Sep 17 00:56:18 1996 print,3*5

Sample Commands

• Some sample commandIDL> print,3*5

IDL> print,SysTime()

IDL> num = 5*12

IDL> help,numNUM INT = 60

IDL> print,num60

Creating an Array

• Some sample commandIDL> num = fltarr(40)

IDL> help,numNUM FLOAT = Array(40)

IDL> print,num 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000

0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000

0.00000 0.00000 0.00000 0.00000 0.00000

Creating an “Indexed Array”

• Some sample commandIDL> for i=0,39 do num(i)=i

IDL> print,num 0.00000 1.00000 2.00000 3.00000 4.00000 5.00000 6.00000

7.00000 8.00000 9.00000 10.0000 11.0000 12.0000 13.0000

14.0000 15.0000 16.0000 17.0000 18.0000 19.0000 20.0000

21.0000 22.0000 23.0000 24.0000 25.0000 26.0000 27.0000

28.0000 29.0000 30.0000 31.0000 32.0000 33.0000 34.0000

35.0000 36.0000 37.0000 38.0000 39.0000

Shortcut Creating Indexed Array

• Some sample commandIDL> num = findgen(40)

IDL> print,num 0.00000 1.00000 2.00000 3.00000 4.00000 5.00000 6.00000

7.00000 8.00000 9.00000 10.0000 11.0000 12.0000 13.0000

14.0000 15.0000 16.0000 17.0000 18.0000 19.0000 20.0000

21.0000 22.0000 23.0000 24.0000 25.0000 26.0000 27.0000

28.0000 29.0000 30.0000 31.0000 32.0000 33.0000 34.0000

35.0000 36.0000 37.0000 38.0000 39.0000

Indexed Array Example

• Some sample commandIDL> num = findgen(40)*10

IDL> print,num 0.00000 10.0000 20.0000 30.0000 40.0000 50.0000 60.0000

70.0000 80.0000 90.0000 100.000 110.000 120.000 130.000

140.000 150.000 160.000 170.000 180.000 190.000 200.000

210.000 220.000 230.000 240.000 250.000 260.000 270.000

280.000 290.000 300.000 310.000 320.000 330.000 340.000

350.000 360.000 370.000 380.000 390.000

Indexed Array Example

• Some sample commandIDL> num = findgen(40)*10-195IDL> print,num

-195.000 -185.000 -175.000 -165.000 -155.000 -145.000 -135.000 -125.000 -115.000 -105.000 -95.0000 -85.0000 -75.0000 -65.0000 -55.0000 -45.0000 -35.0000 -25.0000 -15.0000 -5.00000 5.00000 15.0000 25.0000 35.0000 45.0000 55.0000 65.0000 75.0000 85.0000 95.0000 105.000 115.000 125.000 135.000 145.000 155.000 165.000 175.000 185.000 195.000

• In general you will be using a linear formy=m*x+b

Extracting Array Subset

• Some sample commandIDL> num1 = num(0:36)

IDL> print,num1 -195.000 -185.000 -175.000 -165.000 -155.000 -145.000 -135.000

-125.000 -115.000 -105.000 -95.0000 -85.0000 -75.0000 -65.0000

-55.0000 -45.0000 -35.0000 -25.0000 -15.0000 -5.00000 5.00000

15.0000 25.0000 35.0000 45.0000 55.0000 65.0000 75.0000

85.0000 95.0000 105.000 115.000 125.000 135.000 145.000

155.000 165.000

IDL> help,num1 NUM1 FLOAT = Array(37)

Scalar Operation on Arrays

• Some sample commandsIDL> num = findgen(37)*10

IDL> print,num0.00000 10.0000 20.0000 30.0000 40.0000 50.0000 60.0000

70.0000 80.0000 90.0000 100.000 110.000 120.000 130.000

140.000 150.000 160.000 170.000 180.000 190.000 200.000

210.000 220.000 230.000 240.000 250.000 260.000 270.000

280.000 290.000 300.000 310.000 320.000 330.000 340.000

350.000 360.000

IDL> help,num NUM FLOAT = Array(37)

Two-Dimensional Plotting

• Displaying Two-Dimensional GraphicsIDL> line = Sin(num * !DtoR)

IDL> plot,line

Two-Dimensional Plotting

• Displaying Two-Dimensional GraphicsIDL> plot,line IDL> plot,num,line

Two-Dimensional Plotting Options

IDL> plot,num,line,/Psym

Two-Dimensional Plotting Options

IDL> plot,num,line,Psym=1

Other Plotting OptionsIDL> plot,num,line,linestyle=1

Labeling IDL PlotsIDL> plot,num,line,xrange=[0,360],$

ytitle='Sin(x)',Title='Sine Plot'

Other Plot FacilitiesIDL> plotplot,num,line,xrange=[0,360],$

ytitle='Sin(x)',Title='Sine Plot'

IDL> plots,[0,360],[0,0]

Polygon FillingIDL> loadct, 5IDL> polyfill,num(0:18),line(0:18),color=120IDL> polyfill,num(18:36),line(18:36),color=50

Data File Inputs

• Copy the following files from ~rvrpci/pub/IDL/training

surface.dat

ctscan.dat

Three-Dimensional Plotting• Displaying Three-Dimensional Graphics

– To get the data into IDL, you will need to open a file and read it into an array using the following commands (N.B. this is specific for this particular data set)

IDL> file=‘surface.dat’IDL> openr,lun,file,/get_lunIDL> peak=fltarr(40,40)IDL> readf,lun,peakIDL> free_lun, lun

Surface Plots• Displaying Three-Dimensional Graphics

IDL> print, max(peak), min(peak)

IDL> surface,peak

Surface PlotsIDL> x = findgen(40)*2.5 & y=findgen(40)+50

IDL> surface,peak,x,y

Color Look-up Tables– The following command lets you view the

current colors loaded into the color look-up tables

IDL> cindex

Background/Foreground Colors– The following surface command lets you

manipulate the background and foreground colors

IDL>surface,peak,background=200,color=45,title='Surface Projection',charsize=2.0

Shaded Surface PlottingIDL> shade_surf, peak

Draping Data Over a Surface– We can drape a given surface over another one by the

following command (we read in the data snow the same manner as peak)

IDL>shade_surf,peak,$

shade=snow

Surface Plot Colored by Altitude– We can generate a simple elevation color map.

IDL> shade_surf, peak, shade=bytscl(peak)

Simple Contour Plots– We can generate a simple contour map.

IDL> contour,peak

Contour Plots Spacing– We can generate a simple contour map and adjust spacing.

IDL> contour, peak, nlevels=10

Contour Plots Labelling– We can generate a simple contour map and label them

IDL> contour, peak, nlevels=10, /follow

Contour Plots Spacing Adjustment– We can generate a simple contour map and adjust spacing

along with line styles.IDL> contour, peak, nlevels=10, c_linestyle = indgen(5)

Colored Contour Plots– We can generate a simple contour map and adjust spacing

along with line color.IDL> contour, peak, nlevels=10, c_colors = indgen(5)*20+80

Images Data• Displaying Image Data

– You can read in image data in a manner similar to other forms of two-dimensional data. You will need to ascertain whether the data is an ASCII file or a binary file

IDL> openr, lun,’ctscan.dat’,/get_lunIDL> scan=bytarr(256,256)IDL> readu, lun, scanIDL> free_lun, lun

Image Data• Displaying Image Data

IDL> window, /free, xsize=256, ysize=256

IDL> tvscl, scan

Image Data• Displaying Image Data

IDL> tvscl, scan IDL> tvscl,hist_equal(scan)

Image Processing• Displaying Image Data

IDL> tvscl, scan IDL> tvscl, sobel(scan)

XLOADCT Routine• Displaying Image Data

IDL> xloadctIDL> tvscl, scan

Writing IDL Programs

General Form of an IDL Program• Parts of an IDL program

pro program_name [, arg1, arg2, ..., argn]

IDL statement #1IDL statement #2

.IDL statement #n

Compiling an IDL Program

• Create a program to take celsius arguments and print out fahrenheit valuespro celsius_to_fahrenheit, celsius

print, 9/5*celsius+32

• To compile:IDL> .compile celsius_to_fahrenheit.pro

Executing an IDL Program

• To use the program in scalar mode:

IDL> c = 100.0

IDL> celsius_to_fahrenheit, c

• To use the program in array mode:

IDL> c = findgen(201)/2

IDL> celsius_to_fahrenheit, c

General Form of an IDL Function• Parts of an IDL function

function program_name [, arg1, arg2, ..., argn]

IDL statement #1IDL statement #2

.IDL statement #nreturn, answer

Compiling an IDL Function• Create a function to take celsius

arguments and print out fahrenheit valuesfunction celsius_to_fahrenheit, celsius

answer=9/5*celsius+32

return, answer

• To compile:IDL> .compile celsius_to_fahrenheit.pro

Executing an IDL Function

• To use the program:IDL> c = findgen(201)/2IDL> f = celsius_to_fahrenheit(c)IDL> print, fIDL> plot, c, f

• Equivalently:IDL> print, celsius_to_fahrenheit(c)IDL> plot, c, celsius_to_fahrenheit(c)

Program Designprogram file vs. function file

• Use functions in case of mathematical functions to allow composition of functions

IDL> k=celsius_to_kelvin(fahrenheit_to_celsius(f))

• Use program file in case of process not requiring the return of data

IDL> write_to_file,’temperature.dat’,k

Parameter Passing in IDL

• IDL passes parameters into functions and programs by reference

• This means parameters can be modified through “side effect”

pro celsius_to_fahrenheit, celsius

celsius =9/5*celsius+32

Problems with “Side Effect” Parameter Passing

IDL> a=32.0

IDL> celsius_to_fahrenheit,a

% Compiled module: CELSIUS_TO_FAHRENHEIT.

IDL> print,a

64.0000

A Debugging Exercise

• The IDL function and program examples celsius_to_fahrenheit all compile, but generate incorrect answers in some cases.

• This is left as a debugging exercise

Hints for Statistics Problem

• Study the IDL functionsn_elements()

size()

• Try to reuse as many of the functions as possible.

• Do not be concerned with efficiency at this point.

1 winter quarter 2003rolando v. raqueño basic idl commands simg 726

Documents

specialist international medical graduate (simg

“anzca can take pride in the contributions we make to...

simg. società italiana di medicina generale e delle...

combined specialist international medical graduate (simg...

rolando v. raqueñothursday, july 16, 2015 1 computing for...

j. kerekes, c. raqueño, salvaggio, j. van aardt, t. bauch,...

digital watermarking simg 786 advanced digital image...

1 winter quarter 2003rolando v. raqueño final project...

«technical using of measurement...

1 winter quarter 2005rolando v. raqueño quick latex...

preliminary programme - simg. società italiana di ... ·...

digital imaging and remote sensing laboratory hyperspectral...

04-simg-217-20043-history of astronomical imaging … stars...

1 winter quarter 2005rolando v. raqueño twas the lecture...

idl gui for digital halftoning final project for simg-726...

advances in cardiac electrophysiology -...

specialist international medical graduate (simg ...anzca cv...

simg-782 introduction to digital image processing …...

simg assessment guide - ranzco

1 winter quarter 2004rolando v. raqueño quick latex...