20010614: NEXRAD composites

>From:  Mike Patnoe <address@hidden>
>Organization:  Boeing
>Keywords:  200106142026.f5EKQ9p10894 McIDAS NEXRAD composites

Mike,

>Has anyone done NEXRAD composites on McIDAS?

Yes, Rick Kohrs of SSEC put together a script to composite the NEXRAD stuff
some time ago.  I will let Rick fill folks in on the details.

>That is, composites of two or
>more different sites. I don't need a lengthly reply explaining how it is
>done, I am just trying to verify that it is possible at this point, but
>maybe a short clue would be useful. For example, "use IMGREMAP, then
>IMGOPER."

Yup, that is the idea.

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

>From address@hidden Thu Jun 14 18:00:58 2001
>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

>From address@hidden Fri Jun 15 05:25:12 2001
>Subject: RE: 20010614: NEXRAD composites

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

From: Rick Kohrs <address@hidden>
Date: Fri, 15 Jun 2001 09:16:19 -0500
To: "Craddock, Mary Ellen" <address@hidden>
CC: "'address@hidden'" <address@hidden>,
   Mike Patnoe <address@hidden>, address@hidden,
   address@hidden
Subject: 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

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