NEXRAD composites

Hello All,

I have attached a McBASI script and enhancement file to this message. 
What does the script do:
1) Composite 1 to n states of radar sites. 
2) Combines stations using clear air and precip mode (common
enhancement)
3) Adds topography image as a base map (source area must be AREA8800 -
or change the script)
4) Runs WWDISP command

The topography image has been distributed with McIDAS-XRD version 7.8. 
The file is also available via anonymous ftp:

ftp ssec.wisc.edu
cd pub
cd mug.outgoing
get AREA8800
get CLRPCP.GRX
get CLRPCP.ET
get RADUS.MCB

RADUS.MCB will combine all the sites in the US.  This takes quite a
while to run!  If you see something you'd like added to the script let
me know, I can always try to make modifications.

Rick Kohrs
address@hidden



"Craddock, Mary Ellen" wrote:
> 
> We use Rick Kohr's routine MOSAIC.MCB to do this.  It is easily modified to
> either do a composite over the entire CONUS or over a regional area.
> 
> Mary Ellen
> 
> Mary Ellen Craddock
> Logicon-TASC
> address@hidden
> 
> -----Original Message-----
> From: Russ [mailto:address@hidden]
> Sent: Thursday, June 14, 2001 8:02 PM
> To: address@hidden
> Cc: Mike Patnoe; address@hidden; address@hidden
> Subject: Re: 20010614: NEXRAD composites
> 
> Tom Yoksas wrote:
> 
> > Along the same lines, Geoff Freeman of the Australian Bureau of
> > Meteorology has been working on a routine that can be used to display
> > image data on top of any navigated display (including satellite
> > images).  The image values are written into the graphics frame, so they
> > can be easily toggled on and off.  The initial version of the routine
> > that Jeff Wilson of ABoM showed me allows one to put up multiple radar
> > displays in effect creating a mosaic.  The work on this routine is
> > ongoing and probably not yet ready for prime time.  Jeff told me that
> > he hopes to contribute it to the McIDAS Developers Forum when it is
> > further along in its development cycke.
> >
> > Tom
> 
>      I'll add that there is an old McIDAS application that writes image data
> as a
> graphic. It is called DFC and was used in much the same way as Jeff's
> program
> mentioned above.
> 
>      We also have an application that allows multiple collocated images to
> be
> viewed on a single frame as a RGB composite. COMBINE is a great way view
> multispectral colorized composites. We have used this application to
> generate true
> color images using the MODIS channels 2, 3 and 5.
> 
>      If you visit our web site at
> 
>               http://www.ssec.wisc.edu/~robo/private/radar.html
> 
>     you will see the NEXRAD composites that the McIDAS data center is
> creating.
> 
>     Russ

REM ? RADCOMP -- Creates composite image of RADAR data for specified state
REM ?      RADCOMP state lat lon res
REM ? Parameters:
REM ?      state(s)    | single quoted string of 2-letter state ids of RADAR
REM ?                  | sites to be combined
REM ?      lat         | center latitude of destination image
REM ?      lon         | center longitude of destination image
REM ?      res         | resolution in kilometers of destionation image
REM ?
REM ? Example
REM ?     RADCOMP 'CA OR NV' 39 120 1

     input l$
     input n$
     input e$
     input k$
     input q$

rem
rem
rem
     if val(k$) < 1 print "Topography map only supports resolutions > 3km" : 
stop
 
     let m$ = "Precip"
     let m = 6

rem
rem --- User can define the name of the temporary data set by changing the 
value of t$
rem 
1    let t$ = "RADCOMP/TEMP"

rem
rem --- User can define the directory name for storing temporary data required 
to create composite
rem --- Change the string r$ from temp_directory to your choice.  
rem --- The default value is $HOME/mcidas/data/RAD_TEMP
rem
     let r$ = "temp_directory"

rem
rem --- If user has changed value of r$, don't bother determining the value of 
$HOME 
rem --- and move on to creating dataset name with dsserve
rem
     if r$ <> "temp_directory" goto 5

rem
rem --- Determine value of $HOME - needed for REDIRECT commands
rem
          keyin "OS DEV=T USER_HOME R {echo $HOME}"
          open "USER_HOME",10
          input @10,a$
          h$ = mid$(a$,1,len(a$))
          r$ = h$ + "/mcidas/data/RAD_TEMP"
          close 10

rem
rem --- Create a temporary data set to store temporary images for creating 
composite
rem
5    keyin "DSSERVE ADD RADCOMP/TEMP AREA 9001 9299 {Temporary areas for 
creating RADAR COMP}"

rem
rem --- Create a temporary directory for the temporary images for creating 
composite
rem
     keyin "OS {mkdir $HOME/mcidas/data/RAD_TEMP}"

rem
rem --- Set up a file redirection for temporary images
rem
     keyin "REDIRECT ADD AREA90* {"; r$ ;"}"
     keyin "REDIRECT ADD AREA91* {"; r$ ;"}"
     keyin "REDIRECT ADD AREA92* {"; r$ ;"}"
     keyin "IMGDEL RADCOMP/TEMP.1 189"

rem
rem --- Loop through list of states
rem
     let c = 1
     let p = 10
     x = 1
     l = len(l$)
15   if x > l goto 100
     s$ = mid$(l$,x,2)
     
rem
rem --- Determine the stations for the state requested
rem
     keyin "STNLIST ST="; s$ ;" TYPE=N DEV=T RAD_SITE.DOC R"

rem
rem --- Loop through all the stations
rem
     open "RAD_SITE.DOC",10
     input @10,a$
10   input @10,a$
     print a$
     if mid$(a$,1,4)=EOF$ close 10 : x = x + 3 : goto 15
     if mid$(a$,1,4)="Numb" close 10 : x = x + 3 : goto 15
     p$ = mid$(a$,7,3)

rem 
rem --- Get a listing for the ID and check mode to see if we process the data
rem
     keyin "IMGLIST NEXRAD/BREF1 FORM=EXP ID="; p$ ;" DEV=T RAD_PROD.DOC R"
     open "RAD_PROD.DOC",20
     input @20,b$
     input @20,b$
     input @20,b$
     input @20,b$
     input @20,b$
     input @20,b$
     input @20,b$
     print m$
     print mid$(b$,69,m)
     if mid$(b$,69,m) <> m$ close 20 : goto 10

rem
rem --- Copy the data and process
rem
     keyin "IMGCOPY NEXRAD/BREF1 "; t$ ;".1 SIZE=SAME STYPE=VISR ID="; p$
     keyin "IMGREMAP "; t$ ;".1 "; t$ ;"."; p ;" PRO=MERC LATLON="; n$,e$ ;" 
RES="; k$
     let z=status
     if z=0 if c=1 keyin "IMGCOPY "; t$ ;"."; p ;" "; t$ ;".199 SIZE=ALL" : let 
c=10 : close 20 : let p = p + 1 : goto 10
     if z=0 if c<>1 keyin "IMGOPER "; t$ ;"."; p ;" "; t$ ;".199 "; t$ ;".200 
ZERO=DATA FORM=CMAX SIZE=ALL"
     if z=0 if c<>1 keyin "IMGCOPY "; t$ ;".200 "; t$ ;".199 SIZE=ALL"
     let p = p + 1
     close 20
     goto 10

100  keyin "IMGOPER "; t$ ;".200 "; t$ ;"."; p ;" "; t$ ;".199 ZERO=DATA 
FORM=CMAX SIZE=ALL"
   
rem
rem --- Scale the precip/clear area to -28 to 75 dbz then scale to britness 
values 0 128
rem
     if q$ = "NO" goto 101
     if m$ = "Precip" let m$ = "Clear" : let m = 5 : keyin "IMGOPER "; t$ 
;".199 "; t$ ;".190 SCALE=16 240 82 255 SIZE=ALL"  : goto 1
101  keyin "IMGOPER "; t$ ;".199 "; t$ ;".191 SCALE=16 240 0 139 SIZE=ALL"
     keyin "IMGOPER "; t$ ;".190 "; t$ ;".191 "; t$ ;".192 SCALE=0 255 0 128 
SIZE=ALL FORM=CMAX ZERO=DATA"

rem 
rem --- Set up navigation for topo map
rem
     keyin "IMGCOPY "; t$ ;".190 "; t$ ;".193  SIZE=ALL"
     keyin "IMGREMAP A/A.8800 "; t$ ;".193"
     keyin "IMGOPER "; t$ ;".193 "; t$ ;".194 SCALE=0 255 0 255 LLMT=149"
     keyin "IMGOPER "; t$ ;".194 "; t$ ;".192 "; t$ ;".195 SCALE=0 255 0 255 
SIZE=ALL FORM=CMIN ZERO=DATA"

rem 
rem --- Display the final image, add map and weather boxes

     keyin "IMGDISP "; t$ ;".195 EU=CLRPCPT"
     keyin "PRDUTIL ADD "; t$ ;".195 2 ECHO -28 75 0 128 DBZ 0 0"
     keyin "LWU POKE AREA9195 0 20"
     keyin "GU REST CLRPCP"
     keyin "ERASE G"
     keyin "MAP"
     keyin "WWDISP FILL=TRAN 2 COLOR=2 4 3 5 TYPE=SEVWARN SEVWATCH TORWARN 
TORWATCH PLOT=BOX"
     keyin "GU MAKE 2 FIREBRICK"
     keyin "GU MAKE 3 RED"
     keyin "GU MAKE 4 FIREBRICK"
     keyin "GU MAKE 5 ORANGE"

     stop

Attachment: Clrpcp.et
Description: Binary data