[McIDAS #WIM-982947]: Opinion Question

Hi Rod,

I finally made some time to take a harder look at the script that you sent.
Here are some comments:

Script that you originally sent me:

REM
REM First, erase everyting in frame 1
REM

ERASE F 1 3

IMGCOPY G13/IR27APR11 LOC/AREAS.3003 TIME=20:25 LAT=30:00 86:30 SIZE=320 320 
PRO=MERC RES=6
IMGDISP LOC/AREAS.3003 1 LAT=30:00 86:30 EU=IMAGE
MAP STATE=MS AL GA COUNTY=ALL MCOL=8
MAP FILE=OUTVHPOL MCOL=1
IMGCOPY LOC/AREAS.3003 LOC/AREAS.3103 SIZE=SAME
IMGREMAP G13/ERAD27APR11 LOC/AREAS.3103 TIME=20:25
IMGLIST LOC/AREAS.3003 FORM=ALL
IMGCHA LOC/AREAS.3103 SS=180 STYPE=VISR
IMGDISP LOC/AREAS.3103 2 LAT=30:00 86:30 EU=CLW3
MAP LALO -2 INT=5 5 STATE=MS AL GA COUNTY=ALL MCOL=8
MAP FILE=OUTVHPOL MCOL=1
IMGFILT LOC/AREAS.3003 LOC/AREAS.3103 LOC/AREAS.3203 FILTER=DIS 0 255 185 255 0
IMGDISP LOC/AREAS.3103 3 LAT=30:00 86:30 EU=CLW3 SF=YES
MAP LALO -2 INT=5 5 STATE=MS AL GA COUNTY=ALL MCOL=8
MAP FILE=OUTVHPOL MCOL=1
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 OUTLINE=YES COLOR=2 MODE=N LINELE=61 121 F PLACE=ULEFT 
SIZE=11 28 EXCLUDE=0 1-184

My rework of the script: I made some changes and added some REMarks:

REM
REM First, erase everyting in frame 1
REM

SF 1
ERASE F 1 3

REM
REM Create an image subsector from an IR image and then display it 
REM with appropriate maps overlain
REM

IMGCOPY G13/IR27APR11 LOC/AREAS.3003 TIME=20:25 LAT=30:00 86:30 SIZE=320 320 
IMGDISP LOC/AREAS.3003 1 LAT=30:00 86:30 EU=IMAGE SF=YES
MAP STATE=MS AL GA COUNTY=ALL MCOL=8
MAP FILE=OUTVHPOL MCOL=1

REM 
REM Make a copy of the image just created; remap the effective radius image
REM into the copy; change the image header appropriately; and then display
REM the result with appropriate maps overlain
REM

IMGCOPY LOC/AREAS.3003 LOC/AREAS.3103 SIZE=SAME
IMGREMAP G13/ERAD27APR11 LOC/AREAS.3103 TIME=20:25
IMGLIST LOC/AREAS.3003 FORM=ALL
IMGCHA LOC/AREAS.3103 SS=180 STYPE=VISR
IMGDISP LOC/AREAS.3103 2 LAT=30:00 86:30 EU=CLW3 SF=YES
MAP LALO -2 INT=5 5 STATE=MS AL GA COUNTY=ALL MCOL=8
MAP FILE=OUTVHPOL MCOL=1

REM
REM Combine the IR sector with the effective radius sector and then display
REM it with appropriate maps overlain
REM

IMGFILT LOC/AREAS.3003 LOC/AREAS.3103 LOC/AREAS.3203 FILTER=DIS 0 255 185 255 0
IMGDISP LOC/AREAS.3203 3 LAT=30:00 86:30 EU=CLW3 SF=YES
MAP LALO -2 INT=5 5 STATE=MS AL GA COUNTY=ALL MCOL=8
MAP FILE=OUTVHPOL MCOL=1

REM
REM Do some interrogation of the combined image
REM 

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 OUTLINE=YES 6 MODE=N LINELE=61 121 F PLACE=ULEFT SIZE=11 28 
EXCLUDE=0 1-184

REM
REM Redisplay the combined image blown up by a factor of 3
REM

IMGDISP LOC/AREAS.3203 4 LAT=30:00 86:30 EU=CLW3 SF=YES MAG=3
MAP LALO -2 INT=5 5 STATE=MS AL GA COUNTY=ALL MCOL=8
MAP FILE=OUTVHPOL MCOL=1

Important changes are:

- commentary

- removal of PRO=MERC RES=6 from the first IMGCOPY as these keywords
  are ignored by IMGCOPY since it does not support them

- added SF=YES to the IMGDISP of LOC/AREAS.3103 in frame 2

- display the image created by the IMGFILT command (LOC/AREAS.3203) in
  frame 3 instead of the effective radius image (LOC/AREAS.3103)

- added a display of the filtered image blown-up by a factor of e
  in frame 4

  I did this so that the features are easier to see.

Now, on to your comments:

I said: What happens if you change PRO=MERC to PRO=RECT?

you responded:
> See attached image.

I was not paying enough attention to what you were doing.  The PRO= and RES=
keywords have no effect in an IMGCOPY invocation.


I asked: OK, this is the effective radius image that you are comparing to the 
GR2 one?
 ... Or, is it this image that you are comparing to the GR2 one?

You responded:

> Yes it's this image that is being compared to teh GR2 one.

I think you want to be comparing the image created by IMGFILT to the
GR2 image of radar reflectivity.

Give the reworked BATCH script a try and let me know if you have any
questions.

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: WIM-982947
Department: Support McIDAS
Priority: Normal
Status: Closed