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
NOTE: All email exchanges with Unidata User Support are recorded in the Unidata inquiry tracking system and then made publicly available through the web. If you do not want to have your interactions made available in this way, you must let us know in each email you send to us.