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

Re: Radar Composites

As Daryl noted, CED isn't one of the projections that the GINI image
format supports (see the GINI specification for its description).
For LCC, a tangent cone is supported (eg Lat1 = Lat2).

Here's a quick script snippet to use a Central lat/lon input to create a
region for nex2gini (CLAT is input command line arg #2, CLON
is command line arg #2). KXKY of 755;605 chosen because of a model
domain ( eg 5km * 151x121 ).


#!/bin/csh -f

set LAT=$1
set LON=$2

set DLAT="3.025"
set DLON=`echo "scale=3; 3.775/c(${LAT}*4*a(1)/180)" | bc -l`

exec nex2gini << EOF
   restore $RADNTS/nex2gini.nts
   KXKY = 755;605
   grdarea = #${LAT};${LON};${DLAT};${DLON}
   PROJ = LCC/${LAT};${LON};${LAT}
   r

   e
EOF


The output can be seen under:
http://my.unidata.ucar.edu/content/software/gempak/wseta/index.html

Steve Chiswell
Unidata User Support



On Mon, 2005-02-28 at 14:41, Daryl Herzmann wrote:
> Good afternooon,
> 
> On Mon, 28 Feb 2005, Ryan Jewell wrote:
> 
> > Yes I know what the MAN says, but apparently I am too stupid to know how 
> > many pixels are going to be in the image produced when using bounds of 
> > 35;-105;45;-90 and wanting 1km resolution.  Please tell me again how 
> > many pixels per degree lat/lon there are for 1km res?
> 
> Well, a hand waving estimate is 0.01 degree of lat/lon is 1 km.  So if you 
> want to go from 35 to 45 in the y direction, you need 10 x 100 or 1000 
> pixels.  and -105 to -90, that is 15 x 100 or 1500 pixels.
> 
> 
> daryl