McIDAS Programmer's Manual

Chapter 4 - McIDAS-X Utilities

This chapter describes many of the McIDAS-X library functions that you will use when writing your applications programs. These utility functions perform a variety of common programming tasks, such as:

This chapter is organized into these sections

Overview

The McIDAS-X library functions presented in this chapter are the building blocks for your applications programs. These functions abstract the hardware, such as the display, keyboard and pointing device, and provide routines for manipulating scalar quantities, such as dates and latitudes.

When a user starts a McIDAS-X session, a display similar to the one shown in Figure 4-1 appears on the screen. The window in the foreground is the McIDAS-X Text and Command Window. The window in the background is the McIDAS-X Image Window.

Figure 4-1. At start-up, the McIDAS-X display contains the McIDAS-X Text and Command
Window and the McIDAS-X Image Window.

The visual output of your programs is displayed on these windows. Input is received from the keyboard via the command line and from the pointing device, or mouse, via the Image Window and cursor. Applications also make use of the terminal characteristics, the state of the McIDAS-X session and the application environment.

The functions described in this chapter are grouped by the tasks they perform:

User interface utilities

Users ultimately interact with McIDAS-X in one of two ways:

This section describes each part of the McIDAS-X command line, provides descriptions and examples of the command-line functions you can use in your programs, and defines the return status codes for these functions. This section also describes the standard pointing device for McIDAS-X and provides descriptions and examples of the functions available for reporting the state of the mouse buttons and changing the cursor's location, size, color and type.

 

For additional information about the functions described in this section, see the online man pages provided with the McIDAS-X software.

Command line

McIDAS-X applications are command-line driven. McIDAS-X has two basic command formats:

The multiple-letter command format is discussed in this section. The number of characters permitted in a command line is workstation dependent, although there is no practical limit.

A McIDAS-X command line may contain one or more of the following:

A sample McIDAS-X command line is shown below.

Each command line value is separated by one or more blank spaces. Surround an individual parameter requiring blank spaces with single quotes (' ') so it is treated as one parameter.

Each part of the McIDAS-X command line is discussed below, followed by a description of the functions for extracting individual values from the command line and interrogating the structure of the command line.

Command positional parameters

Command positional parameters provide input to a command and must be entered in the exact order specified. Use them in commands that have options which the user must always enter. One advantage of positional parameters is minimizing the number of keystrokes a user types. Be careful not to negate this advantage by using so many positional parameters that the user can't remember them all.

Keywords

Like command positional parameters, keywords are used for entering command input. Unlike command positional parameters, they are optional for most McIDAS-X commands and can be entered in any order as long as they follow command positional parameters and precede quoted text.

Keyword parameters are often used to clarify commands with many complicated options. Although keywords can occur in any order, their positional parameters must be entered in the order indicated. Since users don't always specify all keyword positional parameters in a command, be sure to assign them reasonable defaults.

In addition to the keywords that an application has built into it, McIDAS-X has six global keywords, which are common to all McIDAS-X commands:

You can retrieve the parameters for these global keywords in your applications, as long as you don't change their functionality.

 

For more information about the McIDAS-X global keywords, see the McIDAS User's Guide.

Quoted text

Quoted text input is most often used when strings entered by a user require blank spaces. Each application can contain only one quote string and it must be the last part of the command the user enters. Using quote fields in McIDAS-X commands is relatively easy from the McIDAS-X Text and Command Window; running McIDAS-X commands with quote strings from Unix shell prompts is more difficult. For this reason, the preferred method for entering strings with spaces in them from the command line is to use a positional or keyword parameter and surround the string with single quotes.

Functions for extracting values from a command line

The McIDAS-X library functions for extracting individual values from the command line are listed alphabetically in the table below. These functions have similar calling sequences.

C function

Fortran function

Description

Mccmddbl

mccmddbl

extracts a value as a double precision number

Mccmddhr

mccmddhr

extracts a time value as a double precision number in units of hours

Mccmddll

mccmddll

extracts a latitude or longitude value as a double precision number in units of degrees

Mccmdihr

mccmdihr

extracts a time value as an integer value in the hhmmss format

Mccmdill

mccmdill

extracts a latitude or longitude value as an integer value in the dddmmss format

Mccmdint

mccmdint

extracts a value as an integer

Mccmdiyd

mccmdiyd

extracts a day value as an integer value in the Julian day format ccyyddd

Mccmdquo

mccmdquo

extracts a character string value from the quote field

Mccmdstr

mccmdstr

extracts a character string value

You can use any of these functions, except Mccmdquo/mccmdquo , to extract information from any command or keyword positional parameter.

All of the functions in the table above, except Mccmdquo /mccmdquo , expect the following arguments:

These functions, except Mccmdstr /mccmdstr , also expect the following additional information:

All the functions above perform similar operations when retrieving data from the command line, including:

McIDAS-X uses a special syntax for identifying keyword names in these functions within an application. This syntax identifies a mandatory and an optional section, separated by a period. For example, if you use FRAME as a keyword in a command, it may be used in the above functions as FRA.ME to denote only the first three characters are mandatory. The table below shows some acceptable and unacceptable forms the user can enter.

Acceptable forms

Unacceptable forms

FRA=

FR=

FRAM=

FRAEM=

FRAMENUM=

 

Keyword names can be letters or numbers, as long as the name is unique; that is, you can't use COL.OR and COL.UMN in the same command.

Below is a description of each function for extracting values from a command line, along with examples.

Sample TEST command

The sample output in the examples below is based on the following user-entered command string for the program TEST. If you see the line goto 999 in the code fragments, it means: exit the command and return an error status.

 TEST LIST 4 END=X 4:30 STA=03-APR-1996 12:30 PAR=T TD HEAD=WI 'Dane County' "Home of the Black Wolf

Extracting a character string value

Use the mccmdstr function to extract a value from the command line as a character string. The code fragment below extracts the value from the first command positional parameter and puts it in the character string option .

character*12 default
character*12 option
default = ' '
ok = mccmdstr(' ', 1, default, option)
if (ok .lt. 0)then
  goto 999
endif

When the TEST command runs, the variable option will contain the string LIST . If the user enters the string ABCDEFGHIJKLMN, an error occurs because option doesn't have room to store this string. The mccmdstr function prints the error message below and returns an error status.

TEST: first positional argument is too big --> ABCDEFGHIJKLMN
TEST: Must be a character string of no more that 12 chars.

The code below extracts the value of the second keyword positional parameter for the keyword HEAD and puts it in the character string mnhead.

character*12 default
character*24 mnhead
default = ' '
ok = mccmdstr('HEA.D', 2, default, mnhead)
if (ok .lt. 0)then
  goto 999
endif

When the TEST command runs, the variable mnhead will contain the string Dane County. This string includes whitespace because the parameter is surrounded by single quotes on the command line. Note the mandatory and optional keyword components. This call will ignore a keyword named HEAP, while it will read a keyword named HEADER.

Extracting an integer value

Use the mccmdint function to extract a value from the command line as an integer value. The following code fragment extracts the value of the second positional parameter as an integer. The call to mccmdint also performs range checking to verify that the value entered by the user is within the range 1 to 9999. If the user doesn't enter a value for the second positional parameter, the default value of 1 is used.

parameter (MINFILE = 1, MAXFILE = 9999)
integer fileno
ok = mccmdint(' ', 2, 'File Number', 1, MINFILE, MAXFILE, fileno)
if (ok .lt. 0)then
  goto 999
endif

When the TEST command runs, the variable fileno will contain the value 4. If the user enters the value 10000 for the second positional parameter, the range allowed in this mccmdint call is exceeded. Thus, mccmdint prints the following error message and returns an error status.

TEST: Invalid File Number.
TEST: 2nd positional argument is too big --> 10000
TEST: Must be valid 'File Number' integer value within range 1 thru 9999.

Extracting time information

The functions to extract time information, mccmddhr and mccmdihr , accept input in a variety of formats. To get a list of the acceptable formats, type the command ARGHELP TIME from the McIDAS-X Text and Command Window.

The following code fragment extracts the value of the second keyword positional parameter as an integer for the keyword START. Note that the call to mccmdihr also performs range checking to verify that the value entered by the user falls within the range 0 to 23:59:59. This is not done automatically, so you can also use these functions to extract an interval of time such as -18:45.

integer starttime
ok = mccmdihr('STA.RT', 2, 'Starting Time', 120000, 0, 235959, starttime)
if (ok .lt. 0)then
  goto 999
endif

When the TEST command runs, the variable starttime contains the value 123000. If the user doesn't enter a value for the second keyword positional parameter for the START keyword, the value returned in starttime is the default of 120000.

Be careful when setting ranges and default values for the mccmdihr function and using the returned values. The format expected is hhmmss . To specify 15:00:00 as the default for a command, you must enter 150000 as the fourth parameter in mccmdihr . If you specify 15 for the default, the value returned will be 00:00:15.

The following code fragment extracts the value of the second keyword positional parameter as a double precision number for the keyword END. The default value is five hours. Note that range checking is turned off in this example by setting the minimum value, 99.0, greater than the maximum value, -99.0.

double precision endtime
ok = mccmddhr('END', 2, 'Ending Time', 5.0d0, 99.0d0, -99.0d0, endtime)
if (ok .lt. 0)then
  goto 999
endif

When the TEST command runs, the variable endtime will contain the value 4.5.

Extracting date information

The function to extract date information, mccmdiyd , accepts input in a variety of formats. To get a list of the acceptable formats, type the command ARGHELP DATE from the McIDAS-X Text and Command Window.

The following code fragment extracts the value of the first keyword positional parameter for the keyword START. This call sets the default to the current day and does not perform range checking.

integer currentday
integer startday
ok = mcgetday(currentday)

ok = mccmdiyd('STA.RT', 1, 'Starting Day', currentday, 99, -99, startday)
if (ok .lt. 0)then
  goto 999
endif

When the TEST command runs, the variable startday will contain the value 1996094.

Extracting a character string value from a quote field

Use the mccmdquo function sparingly when extracting a quote string that contains whitespace. The preferred method for retrieving parameters containing whitespace is to use the mccmdstr function and have the user enter the parameter with single quotes (') as shown in the mccmdstr section above.

The following code fragment extracts the value of the quote string.

character*256 quote
ok = mccmdquo(quote)
if (ok .lt. 0)then
  goto 999
endif

When the TEST command runs, the variable quote will contain the string Home of the Black Wolf . If this field is missing from the user command, the variable is set to blank characters.

Functions for interrogating the command line structure

Sometimes applications need to get information from the command line beyond the actual values entered by the user. The functions described in the table below provide these services.

C function

Fortran function

Description

Mccmd

mccmd

returns the entire user command in one string

Mccmdkey

mccmdkey

verifies that the keywords entered by the user are valid

Mccmdnam

mccmdnam

returns a list of the keywords entered by the user

Mccmdnum

mccmdnum

returns the number of parameters entered for a particular keyword

Each function is described below, along with an example. The sample output is based on the user-entered command string for the program TEST, as shown in the section titled Sample TEST command above.

Verifying the validity of a keyword

Because of the complexity of the McIDAS-X command line, users will sometimes enter an incorrect keyword that is ignored and then not be able to determine why their command doesn't run properly. The mccmdkey function will verify that the keywords entered by the user are valid for the command. It compares the keywords the user enters with a list of acceptable keywords for a command and prints an error message if a discrepancy occurs. Additionally, mccmdkey checks your application to verify that it contains no ambiguous keywords; for example, COL.OR and COL.UMN.

The global keywords described earlier in this section are implicitly included in mccmdkey , so you don't need to include them in your list of keywords to be tested. Your cannot use the mccmdkey function if the keyword names in your command are not fixed; for example, cases where keyword names are related to data structures. Your program cannot predetermine the names the user may have entered; therefore, it cannot use mccmdkey .

The code fragment below shows how to validate the sample TEST command for the keywords END, FILE, PARAMETER, START and TITLE.

parameter (MAXKW = 5)
character*12 validkw(MAXKW)
data validkw/'END','FIL.E','PAR.AMETER','STA.RT','TIT.LE'/

ok = mccmdkey(MAXKW, validkw)
if (ok .lt. 0)then
  goto 999
endif

If the user enters the additional keyword argument COUNTRY=UK , the TEST command will output the following line to the McIDAS-X Text and Command Window and exit.

TEST: Invalid command keywords: COUNTRY=UK

Note that it is the programmer's responsibility to keep the keyword list used in mccmdkey consistent with the keywords actually used by the application. This checking can be turned off by the user by setting syskey word 15 to a non-zero value.

Retrieving a list of user-entered keywords

Sometimes an application needs to know exactly which keywords the user entered. The mccmdnam function will return this information. The code fragment below shows how to retrieve the list of user-entered keywords.

parameter (MAXKW = 5)
character*12 entered_kw(MAXKW)

n_ent = mccmdnam (MAXKW, entered_kw)
if (n_ent .lt. 0)then
  goto 999
endif

When the TEST command runs, the variables in this code fragment will contain the following values:


n_ent           = 4
entered_kw(1)       = 'END'
entered_kw(2)       = 'STA'
entered_kw(3)       = 'PAR'
entered_kw(4)       = 'TITLE'

Note that mccmdkey returns abbreviated values for the START and PARAMETERS keywords, since that is how they were entered on the command line.

Verifying the number of keyword parameters

When applications allow for multiple keyword positional parameters, you must be able to tell how many parameters the user actually included in the command. The mccmdnum function returns this information. The following code fragment extracts the number of keyword positional parameters that the user entered for the PARAMETER keyword.

integer n_parms

n_parms = mccmdnum('PAR.AMETER')

When the sample TEST command runs, n_parms will contain the value 2.

To find out how many command positional parameters the user entered, the code fragment would look like this:

n_parms = mccmdnum(' ')

Status codes

The command-line functions all return status codes. A non-negative value indicates a successful return. A value greater than zero indicates a successful return, but with additional information such as the parameter returned is the default value. A value less than zero indicates a problem, for example in the formatting or range checking. The argument fetching status codes are defined in the table below.

The status codes use the format abcd. The value in the a position indicates the parameter's location; for example, the command line or string table. The values in the b and c positions indicate the form of the value; for example, time or angle. The value in the d position indicates the status.

Status code

Definition

[-]0bcd

argument comes from the default

[-]1bcd

argument comes from the command line

[-]2bcd

argument comes from the McIDAS-X string table

[-]a00d

character string argument

[-]a01d

quote field string argument

[-]a10d

integer argument

[-]a11d

integer hexadecimal argument

[-]a20d

double decimal argument

[-]a21d

double hexadecimal argument

[-]a30d

date argument

[-]a31d

current date argument

- a32d

year within date argument is invalid

- a33d

mon month within date argument is invalid

- a34d

mm month within date argument is invalid

- a35d

day of month (dd) within date argument is invalid

- a36d

day of year (ddd) within date argument is invalid

[-]a40d

integer time argument

[-]a41d

current integer time argument

- a42d

hours within integer time argument are invalid

- a43d

minutes within integer time argument are invalid

- a44d

seconds within integer time argument are invalid

[-]a45d

double time argument

[-]a46d

current double time argument

- a47d

hours within double time argument are invalid

- a48d

minutes within double time argument are invalid

- a49d

seconds within double time argument are invalid

[-]a50d

integer lat/lon argument

- a52d

degrees within integer lat/lon argument are invalid

- a53d

minutes within integer lat/lon argument are invalid

- a54d

seconds within integer lat/lon argument are invalid

[-] a55d

double lat/lon argument

- a57d

degrees within double lat/lon argument are invalid

- a58d

minutes within double lat/lon argument are invalid

- a59d

seconds within double lat/lon argument are invalid

[-] 90d

keyword status

abc0

argument is ok

- abc1

argument has an invalid character

- abc2

integer argument can't contain a fraction

- abc3

argument exceeds system limits for desired format

- abc4

out-of-range argument < given min

- abc5

out-of-range argument > given max

Pointing device

The standard pointing device is usually a three-button mouse. The leftmost button is used by the window manager and the middle and right buttons are used by the McIDAS-X mouse interface.

If you have a two-button mouse, the left button is used by the window manager and the right button performs the tasks normally performed by the right button on a three-button mouse. Press both buttons to perform the tasks normally performed by the middle mouse button. Check the mouse settings in our X-server configuration if you see unexpected behavior.

Users can interact with the McIDAS-X Image Window and McIDAS-X commands via the mouse-driven cursor. McIDAS-X defines the cursor's size, color and shape as long as the cursor resides on the McIDAS-X Image Window. If the cursor is moved outside the window, the cursor's size, color and shape revert to the workstation's settings.

The McIDAS-X function for interfacing with the mouse is mcmoubtn , which performs the following tasks:

You can use the McIDAS-X functions putcur , mcgetcur , sizcur , colcur and typcur in your applications to define the cursor's location, size, color and type. The function mcgetcur will return the cursor position of the cursor when your command starts.

Reporting the state of the mouse buttons

The mcmoubtn function reports the state of the mouse buttons as soon as a mouse button event occurs. Button events include the following:

Applications specify which of the above events triggers a response via the first argument of the mcmoubtn function. Once the event occurs, the state of each mouse button and the cursor's line and element position in the McIDAS-X Image Window are returned to the calling application.

Calling mcmoubtn inherently causes the application to sleep or idle. Thus, you can call it within tight, intensive loops without fear of saturating the processor. If, at any time, the desired mouse event does not occur within a two-minute time span, mcmoubtn returns a timeout status to the calling application. The application is free to call mcmoubtn again if that two-minute time span is not considered unusual.

The code fragment below, from the McIDAS-X ZLM command, monitors the mouse event status.

     parameter button_release=3
     ....
     
c---wait for mouse release
10  status=mcmoubtn(button_release,button(1),button(2),tvlin,tvele)

c---end mouse button sampling
   if((button(1).ne.0.and.button(2).ne.0).or.status.ne.0) then
       call beep(50,10)
       goto 100

c---on left mouse button press

   else if(button(1).ne.0) then

   ......

Note that the mcmoubtn function neither tests for nor protects against two applications using it at the same time. If several applications that call mcmoubtn are started concurrently, the results may be unpredictable.

Reporting the status of Alt G and Alt Q

The mcmoubtn function also indicates the status of the keyboard toggles Alt G and Alt Q, which are monitored by McIDAS-X similar to the mouse buttons. If you hold down the Alt key on the keyboard and press either the G or Q key, mcmoubtn will report that Alt G or Alt Q was pressed. These keyboard toggles are useful when you need to precisely position the cursor on the McIDAS-X Image Window and the act of clicking a mouse button may cause the mouse to move.

Defining cursor location, size, color and type

Use the functions below to change the appearance of the McIDAS-X cursor.

Fortran function

Description

putcur

moves the cursor to a specified location

sizcur

changes the cursor size

colcur

changes the cursor color

typcur

changes the cursor type; for example, crosshair or box

The values that these functions place into User Common are used by the mcimage program to define the appearance of the McIDAS-X cursor. The example below illustrates the use of these cursor functions.

    integer line, element
    integer height, width
    integer color
    integer type

    ..
C Define the new location of the cursor
    line = 200
    element = 351

C Define the new height and width; note: these values must be ODD
    height = 23
    width = 13

C Reassign the cursor to be displayed in color level 3
    color = 3

C Define the new type; 1=box, 2=cross-hair, 3=both, 4=solid,
C 5=bullseye
    type = 3

C Now make the change:

    call putcur ( line, element )
    call sizcur ( height, width )
    call colcur ( color )
    call typcur ( type )

...

Display utilities

The McIDAS-X library provides many functions for displaying text on the McIDAS-X Text and Command Window and displaying images and graphics on the McIDAS-X Image Window.

This section provides descriptions and examples of the functions that you can use in your programs to display text, images and graphics on the McIDAS-X windows, and to obtain information about the McIDAS-X display.

 

For additional information about the functions described in this section, see the online man pages provided with the McIDAS-X software.

Text messages

McIDAS-X has three types of text messages that you can use in your applications:

You should insert debug messages to help trace progress in applications, avoiding extraneous or misleading messages. Also avoid placing debug messages in code where thousands of lines of output could be generated. You will probably not have sufficient room on the McIDAS-X Text and Command Window to view that much output and even if the debug messages are suppressed, they will increase CPU usage unnecessarily.

Figure 4-2. The McIDAS-X Text and Command Window displays standard, error and debug text messages

Routing text

Text can be sent, or routed, to several destinations, as shown in the table below. Users enter the McIDAS-X global keyword, DEV, to specify the destination device of the text output generated by a command.

This DEV= option

Sends the output to

DEV=C

McIDAS-X Text and Command Window (default for standard and error messages)

DEV=P

printer

DEV=F

disk file; 80 characters per line

DEV=N

all output is discarded (default for debug messages)

DEV=T

text file

The user enters the disposition of the text messages as a series of three, contiguous disposition types in this order: standard, error, debug . For example, if the user enters DEV=CCC, all standard, error and debug messages are written on the McIDAS-X Text and Command Window. Entering DEV=CNN displays only the standard messages on the window; the error and debug messages are suppressed. The default is DEV=CCN, meaning the standard and error messages are written on the McIDAS-X Text and Command Window and the debug messages are suppressed.

 

For more information about McIDAS-X global keywords, see the McIDAS User's Guide.

Generating text

To display standard, error and debug messages in C or Fortran applications, use the library functions shown in the table below.

C function

Fortran function

Message type

Mcprintf

sdest

standard text

Mceprintf

edest

error

Mcdprintf

ddest

debug

The sample code below illustrates the use of each message type in a Fortran application.

      subroutine main0
C $ ABCD - illustrates message disposition

      implicit none

c --- external functions
      integer mccmdint   ! integer arg fetch

c --- internal variables
      integer parameter  ! 1st parameter

c --- get the first positional parameter from the command line
      if( mccmdint(' ',1,'1st Parameter',1,1,0,parameter).lt.0)
    & return

      call sdest('STANDARD -- First parameter =',parameter)
      call edest('ERROR -- First parameter =',parameter)
      call ddest('DEBUG -- First parameter =',parameter)

      return
      end

The following table shows where text would be routed, based on the application, ABCD, presented in the code above.

Entering this command

Results in

Output goes to

ABCD 100

STANDARD -- First parameter =100
ABCD: ERROR -- First parameter =100

both messages go to the McIDAS-X Text and Command Window

ABCD 200 DEV=CCC

STANDARD -- First parameter =200
ABCD: ERROR -- First parameter =200
ABCD*DEBUG -- First parameter =200

all messages go to the McIDAS-X Text and Command Window

ABCD 300 DEV=PPF YY

STANDARD -- First parameter =300
ABCD: ERROR -- First parameter =300
ABCD*DEBUG -- First parameter =300

printer
printer
disk file named YY

ABCD 400 DEV=NNN

 

no output

If you don't have an integer value to be printed at the end of sdest, edest or ddest calls, specify zero as the second argument. These functions don't print an integer value of zero. To print a zero, format it into the string. The lines below show the use of a Fortran write to format a text message.

   write(cline,fmt='(a18,1x,f9.2)')'the temperature is',temp
   call sdest(cline,0)

You can create the same formatted output with the following C call.

(void) Mcprintf ("the temperature is %9.2f\n ",temp);

Although sdest, edest and ddest automatically insert a new line when they are called, you must explicitly include the new line character, \n, in applications calling Mcprintf, Mceprintf and Mcdprintf.

Images

An image is information that is usually represented as shades of gray in a two-dimensional matrix, such as a satellite image, a radar image or an image derived from grids. Images are displayed on the McIDAS-X Image Window, as shown in Figure 4-3.

Figure 4-3. The McIDAS-X Image Window displays McIDAS-X-generated images .

You will use the mcline function to write a line of image data into a frame object for display on the McIDAS-X Image Window. A frame object is a memory-based collection of information that completely describes the contents and appearance of a frame to the mcimage process, which creates a visible picture. A frame object contains the actual image data, along with navigation and color enhancement information. Frame objects are stored in McIDAS-X shared memory. The number and size of the frame objects allowed per McIDAS-X session are user-configurable.

The sample code below illustrates how to use mcline to write a frame object. Note that the mcline function writes an entire line. Use mclineseg to write only a partial line.

      integer array(1000)
      .
      .
      .
C --- get the image frame number
      frame = luc(51)

C --- get the size of the frame
      status = mcfsize(frame, num_lines, num_elems)
      if( status.lt.0 ) then
         call edest('Invalid image frame number=',frame)
         return
      endif

C --- fill the output array with a bounded gray scale
      do 10 elem = 1,num_elems
10    array(elem) = MOD( elem-1, 255 )

C --- write array to the frame object
      do 100 line = 1,num_lines
100   call mcline( frame, line, array )

The output from this code displays a series of black-to-white gray shades. It may not map to the McIDAS-X Image Window quite as you expect, however. McIDAS-X applications envision the McIDAS-X Image Window as being n lines by m elements, and capable of displaying 256 colors. Each element in the array holds a 1-byte value (0 to 255) representing a brightness value.

McIDAS-X applications use the mcline function to write image data, or brightness levels, into a frame object. The mcimage program then takes this 8-bit image data and maps it into the McIDAS-X Image Window. Because hardware configurations impose limitations on the number of displayable brightness levels, the image data in the frame object must pass through a filter that maps the brightness levels into the range of display levels. Consequently, the image brightness levels in the frame object rarely correspond one-to-one with the display levels in the McIDAS-X Image Window.

 

For more information about McIDAS-X image display characteristics, see the McIDAS-X TERM command with the SCALE option in the McIDAS User's Guide. Also see the section titled Displaying Images and Graphics in the Introduction chapter of the McIDAS User's Guide.

For more information about McIDAS-X shared memory and frame objects, see the section titled Shared memory in Chapter 2, Learning the Basics.

Graphics

McIDAS-X graphics appear to the user as text, symbols and line segments drawn in color on the McIDAS-X Image Window. Though displayed as part of the same frame object as images, graphics are manipulated through a separate, vector-based API. Vector-based means that graphics are produced by an imaginary drawing pen that can be moved about on the frame by function calls. The pen traces a line segment of the specified color as it moves, and the graphics subsystem converts this into the appropriate pixels and places it in the frame object. McIDAS-X also provides a simple interface for writing text and numbers as vector graphics.

This section describes how to use the graphics subsystem, providing both basic and advanced McIDAS-X graphics techniques. An extended example, GRAF.PGM, and the output it produces, accompanies much of the discussion. The listing for GRAF.PGM appears later in this section.

The table below alphabetically lists the McIDAS-X graphics functions used on the following pages. Additional graphics functions for displaying wind barbs, wind vectors, and weather symbols on a frame are provided later in this section.

Fortran function

Description

dshoff

turns the dash mode off

dshon

turns the dash mode on

endplt

closes or binds off the graphics frame

enpt

flushes the graphics buffer

initpl

initializes the graphics package to write to a frame object

newplt

erases the current graphics frame

page

resizes the viewport or world

plot

draws a line

pltdig

writes an integer number

qgdash

returns the current dash mode

qscale

returns the current scaling mode

sclhgt

converts text height from world to frame coordinates; use only if scaling is turned on; does not reduce the height of the text

scloff

turns scaling off; subsequent graphics calls use frame coordinates

sclon

turns scaling on; subsequent graphics calls use world coordinates

sclpnt

converts a point from world to frame coordinates

wrtext

writes a text string

Basic McIDAS-X graphics techniques

This section describes the techniques that you will most often use when working with the McIDAS-X graphics subsystem. Figure 4-4 below shows
a graphic produced by the sample GRAF application when run from a McIDAS-X session with no arguments and a frame size of 480 lines by
640 elements.

Figure 4-4. This GRAF.PGM output is the result of a McIDAS-X session run with
no arguments and a frame size of 480 lines by 640 elements.

To produce such a graphic, an application must do the following:

These procedures are described below.

Initializing the graphics subsystem

Before you can use the McIDAS-X graphics subsystem, you must initialize it with a call to initpl as shown in the following code fragment from GRAF.PGM, lines 88 to 93; frame is the number of the frame where the graphics will appear and width is the line width, in pixels. Specify zero to use the current frame and width.

C     // Initialize the plot and fetch the color

      frame = luc(-1)
      call initpl( frame, width )

      if( mccmdint('COL.OR', 1, 'Graphic Color', 3, 1, 3,
     &  color ) .lt. 0 ) return

Drawing lines

All graphics objects are composed of line segments built using the plot function. The graphics subsystem lets you specify the following line segment attributes:

Line color

The diagonal line in Figure 4-4 is generated by these lines of code (91 to 132) in GRAF.PGM:

      call initpl( frame, width )
      ...
      call plot( w_ullin, w_ulele, PEN_UP ) 
      call plot( w_lrlin, w_lrele, color  ) 
      call enpt

The second plot call draws the line in color level color from the previous cursor position to line w_lrlin and element w_lrele , here the lower-right corner of the frame. The previous position is always the one specified by the most recent plot call; the PEN_UP defined constant, with a value of zero, moves the pen from its previous, undefined location to the start of the line at (w_ullin,w_ulele ) without drawing.

To use a particular graphics form repeatedly, you can write a routine to generate it at the desired location with desired attributes. A simple example from GRAF.PGM, lines 263 to 286, is shown below.

      subroutine box( ullin, ulele, lrlin, lrele, color )

      <...variable declaractions deleted...>

      call plot( ullin, ulele, PEN_UP )
      call plot( ullin, lrele, color  )
      call plot( lrlin, lrele, color  )
      call plot( lrlin, ulele, color  )
      call plot( ullin, ulele, color  ) 
      call enpt
      ...

This example draws rectangles on the display with a single call and makes applications easier to write and understand.

Dash mode

The dash mode is set by the dshon and dshoff functions. To determine the current mode at any time, call the qgdash function.

When a McIDAS-X session starts, the dash pattern is preset to 10 pixels on and 10 pixels off. You can change this pattern using the McIDAS-X GD command, either from the command line before running any graphics applications, or within an application using the keyin or Mckeyin functions. These values are stored in User Common words 47 to 49.

Lines 146 to 151 in GRAF.PGM draw a dashed box to outline the viewport, which is described later, using the calls below.

      if( border .gt. 0 ) then
          call qgdash( old_dash )
         call dshon 
          call box( v_ullin, v_ulele, v_lrlin, v_lrele, color )
          if( old_dash .eq. 0 ) call dshoff 
      end if

Though not strictly necessary in this context, this sequence also illustrates a helpful idiom: capturing an aspect of graphics status (here the dashing state old_dash ), changing it to a known state (dshon) and then restoring the original state when done. Such behavior is especially important when writing library functions, since they should be careful to leave things the way they found them.

Line width

The line width is established when the graphics subsystem is initialized with initpl . To change the line width, you must close the current graphic and restart another, specifying a new width. Previously drawn graphics remain on the frame with the former line width and new graphics are drawn with the new line width.

Displaying text and numbers

The McIDAS-X library provides two graphics functions, wrtext and pltdig , for writing text and numbers on a graphics frame. When using these functions, you must designate the following:

Lines 361 to 386 in GRAF.PGM use wrtext to produce the alphanumeric characters shown in Figure 4-4 through its own plttext subroutine, as shown below.

C --- Determine the center point in frame coordinates. Because
C --- sclpnt() always works whether the world is active or not,
C --- make this call only if scaling is actually on.

      call qscale( sclstat )
      if( sclstat .eq. 1 ) then
         call sclpnt( clin, cele, txtclin, txtcele )
      else
          txtclin = clin
          txtcele = cele
      end if
C --- Determine the length of the text message in characters.
C --- This information is then used, together with the actual
C --- heights, to compute the upper left corner in frame coordinates
      nchar = len_trim( text )
      ullin = txtclin - txthgt/2
      ulele = txtcele - txthgt*nchar/2 + txthgt/10

C --- Display the text. Turn scaling off first (we have figured
C --- the text height and location in frame coordinates ourselves,
C --- right?) but be sure to restore it to its original state
C --- according to the result 'sclstat' of the qscale() call above
      call scloff
      call wrtext( ullin, ulele, txthgt, text, nchar, color )
      if( sclstat .eq. 1 ) call sclon

For titles and labels, it is often easier to decide where text should be centered rather than where its left margin should be. The above code uses the number of characters nchar and height in pixels txthgt to compute the offsets in line and element from the center to the upper-left. The internal plttext subroutine encapsulates the logic and other computations to account for graphics scaling so GRAF.PGM generates a centered line of text with a single call, just as box draws a graphics object in a clear, abstract way.

The pltdig function is similar to wrtext except that it also allows the user to specify the number of spaces. If the value to be plotted contains fewer digits, the result is left-padded with zeros.

Flushing the buffers

Because the graphics subsystem is buffered, graphics may not appear as they are being drawn. Interactive applications may need to force graphics to appear at various stages in the process. To do this, you must enter the call below to flush the buffers without ending the application's ability to draw additional graphics.

      call enpt

In the listing for GRAF.PGM, this is done at the end of the box and plttext subroutines (lines 286 and 388). All graphics drawn to that point become visible and the graphics package, including all graphics options the application may have set, remains active.

Binding off graphics

When you're done creating graphics, call the endplt function (line 224 in GRAF.PGM) to bind off the graphics as shown below.

      call endplt

Calling endplt displays all graphics on the McIDAS-X Image Window and frees all associated resources. To create more graphics on this or other frames, you will have to call initpl again.

Advanced McIDAS-X graphics techniques

This section describes some advanced techniques that you can use with the McIDAS-X graphics subsystem. It builds upon the example and the McIDAS-X session started in the previous section, Basic McIDAS-X graphics techniques.

Using world coordinates

The example presented in the previous section positioned and drew all graphics in terms of pixel locations on the screen, assuming the typical frame size of 480 lines by 640 elements. However, McIDAS-X frames don't have a fixed size. If you restart the previous McIDAS-X session with a frame size of 350 by 480 and rerun GRAF.PGM, the display will look like Figure 4-5 below.

Figure 5-5. This GRAF.PGM output is the result of a McIDAS-X session run with no arguments and a frame size of 350 lines by 480 elements.

Although the text is still centered on (240,320), this is no longer the center of the frame. The text message describing the frame size gets its values from this line of code (line 195) in GRAF.PGM:

      call mcfsize( frame, f_nlins, f_neles )

The upper-left corner is always (1,1) and the lower-right corner (f_lrlin, f_lrele ) can be readily computed from the size.

The problem of positioning graphics in frames of various sizes remains, however. An obvious solution is for the application to create its graphics objects in some normal coordinate system in which the frame always occupies the same range, whatever its physical size in pixels. McIDAS-X refers to such systems as scaled or world coordinates because they let graphics applications run in their own little world, which is always the size that the application finds convenient. McIDAS-X supports the use of world coordinates with the page, sclon, scloff, sclpnt and sclhgt functions.

To use world coordinates properly, you must understand the concepts of defining and activating worlds. The sclon function activates world coordinates (scaling) and scloff deactivates them. You can use only one coordinate system, either world or frame, at a time. All functions that generate graphics using a line or element pixel coordinate (plot, pltdig, wrtext) use the current system, whether frame or world. Frame coordinates are determined entirely by the hardware and/or the McIDAS-X session; world coordinates, on the other hand, may be defined to be convenient for the application. The initpl function defines a default world of 480 lines by 640 pixels; this world may be redefined using the page function. Keep in mind that page also activates the world it defines, whereas initpl does not.

In Figure 4-5, the world is defined to be 480 by 640 by initpl but not activated, as noted on the display. If you instead enter this command:

GRAF WORLD=1 1 480 640

GRAF.PGM looks at the four WORLD parameters and defines and activates a world using the code below (lines 99 to 117).

C   // If the user has specified the necessary four arguments to
C   // WORLD= read them and define the world; page()
C   // automatically activates it.
C   // The defaults are 1 1 480 640 but there is no significance
C   // to this here; the argfetchers will trigger a return
C   // if the user doesn't actually enter a legal value for
C   // each of the four positions

    if( mccmdnum('WORLD') .eq. 4 ) then
        if( mccmdint('WORLD',1, 'World upper left line',
   &      1, 1, 0, w_ullin ) .lt. 1100 ) return

    ...three more argfetches deleted...
    
    call page( w_ullin, w_ulele, w_lrlin, w_lrele, SET_WORLD )
    end if

The resulting output is shown in Figure 4-6 below.

Figure 4-6. This GRAF.PGM output uses world coordinates to position the text.

The diagonal now goes from corner to corner as in Figure 4-4 and the text is again centered even though the frame sizes are different. Making applications that can generate attractive, properly positioned output regardless of the physical size of the frame is the primary use of world coordinates and the page function. If a 480 by 640 world is suitable, you can just enter the calls below to activate initpl 's default world.

      call initpl( frame, width )
      call sclon

In general, you should use a world that is as large as or larger than your largest frame size to minimize truncation errors resulting from a single world pixel occupying multiple pixels of the frame.

One complication of world coordinates is the positioning and sizing of text. The wrtext function can be given positions and a height in either frame or world coordinates. If the world is smaller than the frame, the text is enlarged; if the world is larger than the frame, the text is not reduced. This prevents the text from becoming too small to read. The difficulty is that the application must be able to account for both cases when centering text. One solution is to do the positioning and output directly in frame coordinates even if scaling is on; the sclpnt and sclhgt functions allow explicit conversion of points and text heights from world to frame coordinates. The plttext subroutine in GRAF.PGM, lines 354 to 386, illustrates this technique as shown in the code samples below.

The actual text size in frame coordinates is determined first:

C --- Determine the actual (frame coordinates) height of the
C --- text. (The use of local variable 'txthgt' is needed to keep
C --- plttext() from modifying its inputs.)

      txthgt = hgt
      call sclhgt( txthgt )

The center point is then transformed to frame coordinates:

C --- Determine the center point in frame coordinates. Because
C --- sclpnt() always works whether the world is active or not,
C --- make this call only if scaling is actually on.
      call qscale( sclstat )
      if( sclstat .eq. 1 ) then
          call sclpnt( clin, cele, txtclin, txtcele ) 
      else
          txtclin = clin
          txtcele = cele
      end if

Then the upper-left corner is computed, again in frame coordinates:

C --- Determine the length of the text message in characters.
C --- This information is then used, together with the actual
C --- heights, to compute the upper-left corner in frame
C --- coordinates.
      nchar = len_trim( text )
      ullin = txtclin - txthgt/2
      ulele = txtcele - txthgt*nchar/2 + txthgt/10

Given that the position is in frame coordinates already, you must verify that scaling is off before calling wrtext and then be sure to restore the scaling to its original setting, whether off or on, before returning control to the caller:

C --- Display the text. Turn scaling off first (we have figured
C --- the text height and location in frame coordinates
C --- ourselves, right?) but be sure to restore it to its
C --- original state according to the result 'sclstat' of the
C --- qscale() call above.
      call scloff 
      call wrtext( ullin, ulele, txthgt, text, nchar, color )
      if( sclstat .eq. 1 ) call sclon

In the above examples, you used world coordinates to make graphics applications independent of frame size. You can also change the orientation and origin of world coordinates as well as their extent.

Some graphics such as scatter plots are easier to generate in an (x,y ) coordinate system in which x increases leftward and y increases upward. The page function allows this, but there are two caveats. First, you should make your world domains as large as, or larger than, the frame to avoid truncation. Second, inverting the world by making either coordinate decrease, rather than increase, from the upper-left corner of the frame to the lower-right may cause confusion if you want to display text in world coordinates or use viewports (see below). The difficulty with text is that a specific, negative scale factor for line may result in text becoming very small or invisible.

 

For more information about frame coordinates, see the section titled Coordinate systems in Chapter 2, Learning the Basics.

Clipping and viewports

Sometimes it is desirable to limit graphics to a particular region of the frame. Consider Figure 4-7 below, which is identical to Figure 4-6 except the line segment is only drawn when it is more than 130 pixels, in world coordinates, from the edge of the world.

Figure 4-7. This GRAF.PGM output shows a viewport, or clipping region.

One solution is to determine the intersection of the diagonal with the 130-pixel border to get two new endpoints for the line segment. This algorithm is cumbersome because the diagonal may intersect with any two of the borders. The simple solution, supported by the McIDAS-X graphics subsystem, is to use a viewport . Also called a clipping region, a viewport defines a region of the frame outside of which graphics will not appear even if drawn. Like worlds, viewports are defined using page , but with a fifth argument of zero instead of one. Lines 139 to 151 in GRAF.PGM define the above viewport and draws its extent, as shown below.

      if( mccmdint('BOR.DER', 1, 'Border Width', 0, 0,
     &  (w_lrlin-w_ullin)/3, border ).lt. 0 ) return
      v_ullin = w_ullin + border
      v_ulele = w_ulele + border
      v_lrlin = w_lrlin - border
      v_lrele = w_lrele - border
      call page( v_ullin, v_ulele, v_lrlin, v_lrele, SET_VIEW ) 
      if( border .gt. 0 ) then
          call qgdash( old_dash )
          call dshon
          call box( v_ullin, v_ulele, v_lrlin, v_lrele, color )
          if( old_dash .eq. 0 ) call dshoff
      end if

Note that the viewport also clips text. The main use of viewports in McIDAS-X programming is to generate graphical output in regions set by the LINE and ELEM keywords, as is done by the McIDAS-X commands PTLIST and GRDDISP. Some users have written simple, frame-based graphical interfaces that set aside a portion of the frame as a toolbar and use a viewport to prevent meteorological data from plotting there when the rest of the frame is generated.

Listing for GRAF.PGM

  1 C   THIS IS SSEC PROPRIETARY SOFTWARE - ITS USE IS RESTRICTED.
  2 
  3 C *** McIDAS Revision History ***
  4 C *** McIDAS Revision History ***
  5 
  6 C ? GRAF -- Demo of basic McIDAS graphics
  7 C ?   GRAF <keywords>
  8 C ? Parameters:
  9 C ?   (none)
 10 C ? Keywords:
 11 C ?   BORDER   = | border width, in pixels
 12 C ?   COLOR    = | graphics color level
 13 C ?   WORLD    = | mn_lin mn_ele mx_lin mx_ele  to define and
 14 C ?              | activate a world
 15 C ? Remarks
 16 C ?       This demo nominally draws a diagonal line from 
 17 C ?   the upper left to the lower right of the display.
 18 C ?
 19 C ?       If WORLD= arguments are not specified, the initpl()
 20 C ?   default world of 480 lines by 640 elements is defined
 21 C ?   but not activated. The diagonal line will be correct
 22 C ?   only if the display size is 480 by 640. If WORLD= arguments
 23 C ?   are specified the world will be defined and activated and
 24 C ?   the diagonal correct regardless of world size.
 25 C ? 
 26 C ?   BORDER= defines a viewport (clipping region) the indicated
 27 C ?   number of pixels (world coordinates) from the edge of the
 28 C ?   world and outlines it with a dashed line.
 29 C ? ----------
 30       subroutine main0
 31
 32       implicit      NONE
 33
 34 C