[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
20010614: NEXRAD composites
- Subject: 20010614: NEXRAD composites
- Date: Thu, 14 Jun 2001 14:41:54 -0600
>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