[McIDAS #PIR-584729]: McIDAS Research Questions

Hi Rod,

re: I need to see the full BATCH file
> Ok that should be no problem:

The gist of what I see in the BATCH file that you sent appears to be:

1) create an IR image centered on a region on interest and use the new image
   as a remap target for an effective radius image:

IMGCOPY G13/IR27APR11 LOC/AREAS.3003 TIME=20:25 LAT=30:00 86:30 SIZE=320 320

IMGCOPY LOC/AREAS.3003 LOC/AREAS.3103 SIZE=SAME
IMGREMAP G13/ERAD27APR11 LOC/AREAS.3103 TIME=20:25
IMGCHA LOC/AREAS.3103 SS=180 STYPE=VISR

2) create a new image that is the result of filtering the effective radius
   with the IR image in your region of interest:

IMGFILT LOC/AREAS.3103 LOC/AREAS.3003 LOC/AREAS.3203 FILTER=DIS 0 255 185 255 0
IMGCHA LOC/AREAS.3203 OTYPE=PRD
PRDUTIL ADD LOC/AREAS.3203 1 RICE -1095 515 0 255 ERAD 10 -2139062144

3) display the resultant image and do some probing:

IMGDISP LOC/AREAS.3203 3 LAT=30:00 86:30 EU=CLW3 SF=YES

IMGPROBE LIST POINT MODE=N LATLON=34.689243 88.234695
IMGPROBE LIST POINT MODE=N LATLON=34.685452 87.354070
IMGPROBE LIST POINT MODE=N LATLON=34.177437 88.247765
IMGPROBE LIST POINT MODE=N LATLON=34.179134 87.397865

IMGPROBE STAT BOX RICE OUTLINE=YES 6 MODE=N LINELE=61 121 F PLACE=ULEFT SIZE=11 
28 EXCLUDE=-109.57 -6.77

4) create another image that does not seem to have anything to do with
   things later on (i.e., I don't know why you are doing the following):

IMGFILT LOC/AREAS.3103 LOC/AREAS.3003 LOC/AREAS.3301 FILTER=MERGE 185 255 0

IMGDISP LOC/AREAS.3301 4 LAT=34:00 88:20 GRAY=X X X 389 EU=CLW4 SF=YES MAG=3

5) redisplay the filtered image of the region of interest in another frame and
   then use IMGPROBE to overlay it with contours of the calibrated unit of
   the image, RICE:

IMGDISP LOC/AREAS.3203 5 LAT=34:00 88:20 GRAY=X X X 910 EU=CLW4 SF=YES MAG=3

IMGPROBE CONT BOX RICE SIZE=500 700 CINT=3 CLIM=30 50 COL=7 DASH=YES LAT=34:45 
87:49 MODE=N

6) remap the effective radius image into the result of step 4)

   I don't understand why you are doing this as it does not seem to have any 
use in
   the rest of the BATCH file:

IMGREMAP G13/ERAD27APR11 LOC/AREAS.3103 TIME=20:25

7) create a grid of the BRIT values of the filtered effective radius image
   in the region of interest:

LWU DEL GRID1000
IMGGRD AREA=3203 GRIDF=1000

8) list statistics for the grid just created from the filtered effective radius
   image in the region of interest:

GRDINFO LOC/GRIDS.1000 STAT GRID=1

Comments:

- I don't understand why you are creating LOC/AREAS.3301

- it strikes me that you really want to create a grid of the RICE values
  in the filtered effective radius image in the region of interest, not
  a grid of BRIT values

  This is done by specifying the UNIT= in the IMGGRD invocation:

LWU DEL GRID1000
IMGGRD AREA=3203 GRIDF=1000 UNIT=RICE

- Once you have a grid of effective radius values, you can display the
  values on top of a radar reflectivity display that is close to the same
  time as the effective radius image.  By doing this in McIDAS, you will
  be able to really see the juxtaposition of effective radius values
  and radar reflectivities.

  I grabbed the NEXRAD Level 3 images for GWX for the times that seem to be
  of interest to you some time ago and made them available via the ADDE dataset
  LEVEL3/N0R in your account on iset2:

<as 'rod' on iset2>
cd $MCDATA
imglist.k LEVEL3/N0R.ALL ID=GWX
Image file directory listing for:LEVEL3/N0R
 Pos Satellite/         Date       Time      Center   Band(s)
     sensor                                 Lat  Lon
 --- -------------  ------------  --------  ---- ---- ------------
   1  RADAR         27 APR 11117  19:00:12     GWX    1
   2  RADAR         27 APR 11117  19:10:27     GWX    1
   3  RADAR         27 APR 11117  19:15:05     GWX    1
   4  RADAR         27 APR 11117  19:19:42     GWX    1
   5  RADAR         27 APR 11117  19:28:55     GWX    1
   6  RADAR         27 APR 11117  19:33:32     GWX    1
   7  RADAR         27 APR 11117  19:38:08     GWX    1
   8  RADAR         27 APR 11117  19:42:45     GWX    1
   9  RADAR         27 APR 11117  19:51:58     GWX    1
  10  RADAR         27 APR 11117  19:56:34     GWX    1
  11  RADAR         27 APR 11117  20:01:11     GWX    1
  12  RADAR         27 APR 11117  20:05:47     GWX    1
  13  RADAR         27 APR 11117  20:10:24     GWX    1
  14  RADAR         27 APR 11117  20:15:01     GWX    1
  15  RADAR         27 APR 11117  20:19:39     GWX    1
  16  RADAR         27 APR 11117  20:24:17     GWX    1
  17  RADAR         27 APR 11117  20:33:30     GWX    1
  18  RADAR         27 APR 11117  20:38:07     GWX    1
  19  RADAR         27 APR 11117  20:42:43     GWX    1
  20  RADAR         27 APR 11117  20:47:20     GWX    1
  21  RADAR         27 APR 11117  20:51:57     GWX    1
  22  RADAR         27 APR 11117  20:56:33     GWX    1
  23  RADAR         27 APR 11117  21:05:50     GWX    1
  24  RADAR         27 APR 11117  21:10:26     GWX    1
  25  RADAR         27 APR 11117  21:15:03     GWX    1
  26  RADAR         27 APR 11117  21:24:17     GWX    1
  27  RADAR         27 APR 11117  21:28:53     GWX    1
  28  RADAR         27 APR 11117  21:33:30     GWX    1
  29  RADAR         27 APR 11117  21:38:06     GWX    1
  30  RADAR         27 APR 11117  21:42:43     GWX    1
  31  RADAR         27 APR 11117  21:47:20     GWX    1
  32  RADAR         27 APR 11117  21:51:57     GWX    1
  33  RADAR         27 APR 11117  21:56:33     GWX    1
imglist.k: done

Cheers,

Tom
--
****************************************************************************
Unidata User Support                                    UCAR Unidata Program
(303) 497-8642                                                 P.O. Box 3000
address@hidden                                   Boulder, CO 80307
----------------------------------------------------------------------------
Unidata HomePage                       http://www.unidata.ucar.edu
****************************************************************************


Ticket Details
===================
Ticket ID: PIR-584729
Department: Support McIDAS
Priority: Normal
Status: Closed