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

[GEMPAK #IBI-740956]: LCC projection

This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.

  • Subject: [GEMPAK #IBI-740956]: LCC projection
  • Date: Mon, 29 Jan 2007 09:52:19 -0700
Jamie,

The PROJ variable must be "sat" to display a SATFIL, and
"rad" to display a RADFIL. The projection will be that of the image itself.

GAREA can be either "dset" to use the entire image region, or a subset region.

Steve Chiswell
Unidata User Support




> I am trying to run a script that creates a png image of the current satellite
> data. The problem that I am having is that when I set the PROJ to LCC I get a
> blank image. Not sure why?
> 
> =========Script===========
> setenv DISPLAY :1.0
> #
> set timenow=`/bin/date`
> echo "sat2png called at $timenow"  >>/home/jlahowet/cronsat/gpcron.log
> cd /home/jlahowet/cronsat/png #put PNGs in a folder
> source /home/jlahowet/Gemenviron
> #
> #   set gempak parameters
> #
> 
> set x = 1200 ; set y = 1200
> set png = "vis_sat.png"
> rm -f $png
> #  enhancement tables, list of which are located in $GEMTBL
> #  located at: /apps/nawips/NAWIPS/gempak/tables/luts
> 
> set lutfil = $GEMTBL/visgray.tbl
> 
> #  gemexe parameter set to blank but could be used to
> #  overlay data on satellite
> set sfparm = "blnk"
> set clr =  "2;15;25;24;24;23;23;11;32;2"
> set txt = "0.9/13/1/hw"
> set map = "31/10"
> set filt = .5
> # edited latlon
> # this is used for the grid lines, which I dont want
> # set latlon = "6/4/1/2/10;10"
> #
> #  map boundaries
> #  I need to edit this is that it zeros in on colorado
> set maparea = "35.10;-109.10;42.10;-101.90"
> #  map projection which needs to be set to sat
> set map_proj = "lcc/40;-105;40"
> #set map_proj = "sat"
> 
> #
> # set data files and time for current DTG
> #
> #  current surface map
> #
> set cursao = `ls -t $SAO/*sao.gem | head -1`
> #   current satellite image
> #
> set cursat = `ls -t /data/images/GINI/WEST-CONUS/1km/VIS/* |head -1`
> #  get the hour/minute that the sat file is valid
> set curtim = `ls -t /data/images/GINI/WEST-CONUS/1km/VIS/* |head -1 | cut -d_ 
> -f3`
> #  get the current hour from the filename
> set curhour = `ls -t /data/images/GINI/WEST-CONUS/1km/VIS/* |head -1 | cut -d_
> -f3 | cut -c1-2`
> 
> echo "$cursao"
> echo "$cursat"
> echo "$curtim"
> 
> #  get from system current year,month, and day
> set curyear   = `date -u +%Y`
> set curmon    = `date -u +%h`
> set curday    = `date -u +%d`
> 
> sfmap << eof
> AREA     = GAREA
> GAREA    = $maparea
> SATFIL   = $cursat
> RADFIL   =
> SFPARM   = $sfparm
> DATTIM   = $curday/$curhour
> SFFILE   = $cursao
> COLORS   = $clr
> MAP      = $map
> LATLON   = 0
> MARKER   = 0
> TITLE    = 2/-1/WEST-CONUS VIS $curday$curmon$curyear $curtim UTC
> CLEAR    = yes
> PANEL    = 0
> DEVICE   = png | $png | $x;$y
> PROJ     = $map_proj
> FILTER   = $filt
> TEXT     = $txt
> LUTFIL   = $lutfil
> SKPMIS   = NO
> STNPLT   =
> 
> run
> 
> exit
> 
> eof
> 
> gpend
> 
> ===================================
> Jamie Lahowetz
> Undergraduate of Meteorology
> University of Nebraska - Lincoln
> P: 402-304-0766
> M:address@hidden
> 
> 


Ticket Details
===================
Ticket ID: IBI-740956
Department: Support GEMPAK
Priority: Normal
Status: Closed