[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[McIDAS #RDQ-637778]: Re: New Questions Regarding My Research



Hi Rod,

re:
> Yes I wrote this:
> 
> REM
> REM First, erase everyting in frame 1
> REM
> 
> SF 1
> ERASE F 1 4
> 
> 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.3103 LOC/AREAS.3003 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 RICE OUTLINE=YES 6 MODE=N LINELE=61 121 F PLACE=ULEFT 
> SIZE=11 28 EXCLUDE=-109.57 -6.77

OK, you are trying to probe the image created by the IMGFILT command
invocation and then displayed by IMGDISP.  Since you are getting the
error you sent in the previous email, it is obvious that the calibration
was carried along during the various manipulations.  This means that
you need to calibrate the image referenced by LOC/AREAS.3203; you
do this by adding the appropriate lines immediately after the IMGFILT
line (and before the IMGDISP line):

IMGCHA LOC/AREAS.3203 OTYPE=PRD
PRDUTIL ADD LOC/AREAS.3203 1 RICE -1095 515 0 255 ERAD 10 -2139062144

> REM
> REM Redisplay the combined image blown up by a factor of 3
> REM
> 
> IMGDISP LOC/AREAS.3203 4 LAT=34:00 88:20 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
> 
> So the I displayed the image before using IMGPROBE stat so I guess
> this means its not calibrated properly, weirdly.

Not so weird. Remember the need to run IOMGCHA after the IMGREMAP
step.

re:
> So I downloaded the script on my desktop, transfer it to
> /home/rod/mcidas/bin then I chmod +x CalibrateEffrad.sh followed by cd
> $MCDATA
> CalibrateEffrad.sh $MCDATA/archive/27apr2011
> 
> And Got this:
> 
> Current Calibration Block Segments In File:
> Output         Input
> Product  Product  Brit  Brit
> Num Type low      high     low   high  Units Scale Missing Val
> --- ---- -------- -------- ----- ----- ----- ----- -----------
> 1: RICE    -1095      515     0   255  ERAD    10 -2139062144
> 
> So it should have been calibrated. Dont see the problem there.

The original image(s) are calibrated, but the one created by
the sequence in your BATCH file is not.  All you need to do
is calibrate it.

By the way, after the IMGFILT invocation, your new image will
not have any values below 185 BRITness.   This means that you
should change the calibration for the combination to be two
part:

PRDUTIL ADD LOC/AREAS.3203 1 RICE     0   0   0 184 ERAD 10 -2139062144
PRDUTIL ADD LOC/AREAS.3203 2 RICE -1095 515 185 255 ERAD 10 -2139062144

Since I haven't tried this, I am not 100% sure that it will work the
way that is intended.  If it doesn't work, you will need to play around
with the PRDUTIL invocations to get what you want.

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