[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[GEMPAK #BWT-519283]: Plotting NEXRAD Level II data in Mercator projection
- Subject: [GEMPAK #BWT-519283]: Plotting NEXRAD Level II data in Mercator projection
- Date: Wed, 16 Oct 2013 10:43:29 -0600
Hi Greg,
I played with this last night to confirm what I had found out before, the GP
programs are unable to re-project the site radar projection written into the
file. I saw Daryl's suggestions on gembud and I recommend either: use his
re-projected level 3 or create your own with nex2img.
I have used nex2img/nex2gini before on single stations by controlling what is
in a local (pwd) copy of the table nexrad.tbl for STNFIL.
Michael
> Michael,
> I'm running a script using GPNEXR2_GF to generate individual NEXRAD radar
> plots. I'd like to use these plots in a GIS application so I need them
> plotted with the Mercator projection rather the the RAD one. Is that
> possible? I tried setting proj=mer but the output is a plot with the map and
> lat lon lines on it but no radar data. GEMPAK doesn't generate any error
> message but I get output like:
>
> radar_plot.tmp*
> MAP Map color/dash/width/filter flag 4/1/1+4/1/2+4/1/1
> GAREA Graphics area 38.25;-106.75;41.625;-102.25
> PROJ Map projection/angles/margins|dr mer
> RADFIL Radar image filename(s)
> /export/ldm/data/gempak/nexrad/craft/KFTG/KFTG_20131015_1930
> LATLON Line color/dash/width/freq/inc/l 3/1/2/1/0.5;0.5
> PANEL Panel loc/color/dash/width/regn 0
> TITLE Title color/line/title 2/-1/KFTG NEXRAD Level II
> Reflectivity dz 2013/10/15 1930 UTC
> TEXT Size/fnt/wdth/brdr/N-rot/just/hw 1
> CLEAR Clear screen flag yes
> DEVICE Device|name|x size;y size|color
> gif|ops.NEXRAD.201310151930.l2_KFTG_Reflectivity.gif|900;900
> LUTFIL Enhancement lookup table filenam upc_radimg_white.tbl
> IMCBAR Color/ornt/anch/x;y/ln;wd/freq|t 4/V/ /.05;.05/1
> TILT Radar beam elevation/tilt number 0.5
> RADPARM Radar parameter (dz, vr, sw) dz
> RADTIM Radar composite current/dattim last
> Parameters requested: MAP,GAREA,PROJ,RADFIL,LATLON,PANEL,TITLE,TEXT,CLEAR,
> DEVICE,LUTFIL,IMCBAR,TILT,RADPARM,RADTIM.
> GEMPAK-GPNEXR2> GEMPAK-GPNEXR2> GEMPAK-GPNEXR2> GEMPAK-GPNEXR2>
> GEMPAK-GPNEXR2> GEMPAK-GPNEXR2> GEMPAK-GPNEXR2> GEMPAK-GPNEXR2>
> GEMPAK-GPNEXR2> GEMPAK-GPNEXR2> GEMPAK-GPNEXR2> GEMPAK-GPNEXR2>
> GEMPAK-GPNEXR2> GEMPAK-GPNEXR2> GEMPAK-GPNEXR2> GEMPAK-GPNEXR2>
> GEMPAK-GPNEXR2> PROJ should be set to RAD
> GPMAP PARAMETERS:
>
> Device:
> gif|ops.NEXRAD.201310151930.l2_KFTG_Reflectivity.gif|900;900
> Projection: mer
> Graphics area name: 38.25;-106.75;41.625;-102.25
> Map: 4/1/1+4/1/2+4/1/1
> Title: 2/-1/KFTG NEXRAD Level II Reflectivity dz 2013/10/15
> 1930 UTC
> Panel: 0
> Latlon: 3/1/2/1/0.5;0.5
> Clear: YES
> Image file: ÿÿÿÿÿÿÿÿ°cüüÿâvßV: V:ÿÿÿÿ
> Enter <cr> to accept parameters or type EXIT: Parameters requested:
> MAP,GAREA,PROJ,RADFIL,LATLON,PANEL,TITLE,TEXT,CLEAR,
> DEVICE,LUTFIL,IMCBAR,TILT,RADPARM,RADTIM.
>
> So something goes screwy with the Image file as I get all those weird
> characters. If proj=rad this script works file and the Image file item lists
> an actual file.
>
> Can you help me with this? I need level II data plotted in a mercator
> projection for use with Open Street Maps.
>
> Thanks,
> Greg
> -----------------------------------------------------------
> Greg Stossmeister
> Head - Collaborative Tools and Metadata Group
> Computing Data and Software Facility
> NCAR/Earth Observing Laboratory (EOL)
> Telephone: (303)497-8692 Facsimilie:(303)497-2044
> e-mail: address@hidden
> Mailing Address: Shipping Address:
> NCAR/EOL UCAR Shipping and Receiving
> P.O. Box 3000 3090 Center Green Dr.
> Boulder, CO 80307 Boulder, CO 80301
> ------------------------------------------------------------
>
>
Ticket Details
===================
Ticket ID: BWT-519283
Department: Support GEMPAK
Priority: Normal
Status: Open