[
Date Prev][
Date Next][
Thread Prev][
Thread Next][
Date Index][
Thread Index]
[McIDAS #WIM-982947]: Opinion Question
- Subject: [McIDAS #WIM-982947]: Opinion Question
- Date: Thu, 06 Dec 2012 14:06:14 -0700
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