This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.
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