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

[GEMPAK #BWT-519283]: Plotting NEXRAD Level II data in Mercator projection

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