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

[Support #FWQ-357482]: Re: conversion to image format

> Dr. Chiswell,
> 
> I have a question of gd2img.  The satellite observations that I'm using to
> compare with my model-simulated satellite are all in celsius.  I have both
> celsius and kelvin versions of cloud top temperatures in my gempak files
> before they're converted to gini with gd2img.  Can I also use gd2img to
> convert the celsius cloud top temperatures?

Yes, you can use any variable for GFUNC as long as you create a nex2gini.tbl
entry such as:
CTMPC     15    2 K,255,176,-110,-31;175,0,-31.5,15         TICZ99 TMPK

The band number of n=15 is what is used in imgtyp.tbl column 6
as 2**n-1 (so that would be 2**14), where I had previously
sent you 2**15 for band=16 (or you can just redefine CTMPK to CTMPC
using the same band number).

Steve Chiswell
Unidata User Support

> I attempted to put a table entry for a celsius variable (CTMPC) in
> nex2gini.tbl, by subtracting 273 off all the values for the kelvin cloud
> top temps. Is that correct?  I'm also not sure what values to put into
> imgtyp.tbl for CTMPC.
> Thanks!
> Megan
> 
> > Megan,
> >
> > I put together a completely standalone tarfile (attached to this message)
> > that should run without any upgrade of your previous versions.
> >
> > First thing is to see if your Redhat can run the binary I built here under
> > fedora core 5.
> >
> > Unpack the tarfile (it will create its own "ncsu" subdirectory where you
> > unpack it.
> > In the directory are the nex2gini.tbl file that defines the mapping from
> > temperature to pixel value:
> > CTMPK     16    2 K,255,176,163,242;175,0,242.5,330         TICZ99 TMPK
> >
> > The 2 part IR mapping is pixel values 255 -> 176 convert to 163K -> 242K,
> >                                       175 ->   0 convert to 242.5K -> 330K
> > These match the standard IR curve used for10.7 micron imagery in
> > $GEMPAK/source/gemlib/im/imbtot.f
> > and GEMPAK/source/gemlib/im/imttob.f.
> >
> > The CTMPK in this entry must match the name of your grid function! Band
> > used will be "16" which is used in the $GEMTBL/imgtyp.tbl entry I sent
> > you previously as 2**(band-1) (aka 2**15).
> >
> > You will be able to run the program within the ncsu directory since all
> > tables are searched
> > in the current working directory. Otherwise, the nex2gini.tbl would be
> > found in
> > $GEMTBL/unidata/nex2gini.tbl, the gd2img.pdf would be found in
> > $GEMPDF/gd2img.pdf.
> >
> > The program can be run as:
> >> gd2img
> >  GDATTIM   Grid date/time                    f024
> >  GDFILE    Grid file                         satir_f24.gem
> >  GLEVEL    level                             0
> >  GVCORD    Vertical coordinate               none
> >  GFUNC     Gridded parameter                 CTMPK
> >  SCALE     Output scaling                    0
> >  PROJ      Map projection/angles/margins|dr
> >  GRDAREA   Area covered by grid
> >  KXKY      Number of grid points in x;y
> >  CPYFIL    Grid file whose navigation is to  satir_f24.gem
> >  SATFIL    output file name                  testimg.gini
> >
> > Since your grid is a 4x4km LCC grid, I used that in cpyfil for
> > the output image projection. Otherwise, a GRDAREA, PROJ, and KXKY could
> > be specified but it has to be a lambert tangent cone with dx=dy, but
> > you probably don't want to have to figure that out unless you change the
> > model resoltion!
> >
> > The GFUNC name CTMPK must match the entry in nex2gini.tbl
> >
> > Go ahead and give this a try. As I mentioned, its under development so it
> > spits out some debug information
> > and the help file is basically non-existant, but improving 100% all the
> > time ;-)
> >
> > Let me know if you can get the program to run, and I'll keep working on
> > cleaning up
> > the interface and input. Of course, you will eventually want to run a
> > production version in an
> > up to date GEMPAK release in the future.
> >
> > Steve Chiswell
> > Unidata User Support
> >
> >
> >
> >
> >> Dr. Chiswell,
> >> We are running Red Hat Enterprise Linux, 4 or 5.  I am currently using
> >> gempak v5.6.  If that's too old, I think 5.8.3 is also compiled on the
> >> cluster I'm using, so I could use it instead.
> >> Thanks again for your help,
> >> Megan
> >>
> >> >> Mr. Chiswell,
> >> >> I took a look at that test image you sent, and that's exactly what I
> >> >> wanted.
> >> >> I'd like to generate a lot more of these images from my model output.
> >>  I
> >> >> was wondering if you could send me the code you used.  Could I
> >> install
> >> >> gd2img by itself or would I need a whole new gempak version?
> >> >> Thanks,
> >> >> Megan
> >> >
> >> > Megan,
> >> >
> >> > I have the gd2img directory in the source distribution of 5.10.1, but
> >> > didn't
> >> > configure it to be built in the binary distributions since it is still
> >> > under development.
> >> > I have been testing it out to create topo images from dem data, but
> >> > created the
> >> > TMPK to pixel value entry in the table for your data. I will make it
> >> part
> >> > of the standard
> >> > build for 5.10.2; I just need to add some user friendliness to the
> >> > program.
> >> >
> >> > I can send you the program to run without upgrading your distribution
> >> if
> >> > you can tell me what OS etc you are running. Do you know what version
> >> of
> >> > GEMPAK you are
> >> > currently running?
> >> >
> >> > Dr. Steve Chiswell
> >> > NC State 87, 90, 94
> >> > Unidata User Support
> >> >
> >> >
> >> >
> >> >>
> >> >>
> >> >> > Megan,
> >> >> >
> >> >> > I created the attached GINI format satellite image from your data.
> >> >> > To display the image, you will need to add the following line to
> >> your
> >> >> > $GEMTBL/sat/imgtyp.tbl file (or imgtyp.tbl file in your working
> >> >> > directory):
> >> >> >
> >> >> > CTMPK                IR_TMPK       0    255     11  2**15      1
> >> >> > ir_upc3.tbl
> >> >> >
> >> >> > Note use only spaces in that file, not tabs.
> >> >> >
> >> >> > The conversion process is somewhat manual for this process since I
> >> >> have to
> >> >> > generate a
> >> >> > mapping of pixel values to gridded data values. This is something I
> >> >> need
> >> >> > to add for
> >> >> > users in a working program, but doing this with your data helps me
> >> get
> >> >> to
> >> >> > that
> >> >> > point - I did discuss this type of project with Gary Lackmann here
> >> >> last
> >> >> > summer,
> >> >> > so its a work in progress.
> >> >> >
> >> >> > I need to find out from you if this mapping of brightness
> >> temperature
> >> >> to
> >> >> > pixel
> >> >> > is what you are looking for (its the standard 2 part temperature
> >> curve
> >> >> > used for
> >> >> > IR satellites), and if not, I'll need to get this worked out to
> >> >> provide
> >> >> > you with the
> >> >> > program.
> >> >> >
> >> >> > I displayed this with your grid points and contours in gdplot2 and
> >> the
> >> >> > image lines up
> >> >> > with the points, so my first pass seems to be in the right
> >> direction.
> >> >> >
> >> >> > Let me know how this works for you.....
> >> >> >
> >> >> > Steve Chiswell
> >> >> > Unidata User Support
> >> >> >
> >> >> >
> >> >> >> Here's the grid navigation from gdinfo...
> >> >> >> GDFILE    Grid file                         satir_f24.gem
> >> >> >> LSTALL    Full list flag                    YES
> >> >> >> OUTPUT    Output device/filename            T
> >> >> >> GDATTIM   Grid date/time                    last
> >> >> >> GLEVEL    Grid level                        0
> >> >> >> GVCORD    Grid vertical coordinate          none
> >> >> >> GFUNC     Scalar grid                       CTMPK
> >> >> >> Parameters requested:
> >> >> GDFILE,LSTALL,OUTPUT,GDATTIM,GLEVEL,GVCORD,GFUNC.
> >> >> >> GEMPAK-GDINFO>r
> >> >> >>
> >> >> >> GRID FILE: satir_f24.gem
> >> >> >>
> >> >> >> GRID NAVIGATION:
> >> >> >> PROJECTION:          LCC
> >> >> >> ANGLES:                27.0   -85.3    27.0
> >> >> >> GRID SIZE:          159 159
> >> >> >> LL CORNER:              15.83    -82.64
> >> >> >> UR CORNER:              21.20    -76.45
> >> >> >>
> >> >> >> GRID ANALYSIS BLOCK:
> >> >> >> ANALYSIS TYPE:        BARNES
> >> >> >> DELTAN:               0.073
> >> >> >> DELTAX:           -9999.000
> >> >> >> DELTAY:           -9999.000
> >> >> >> GRID AREA:            15.00  -84.00   22.00  -75.00
> >> >> >> EXTEND AREA:          15.00  -84.00   23.00  -75.00
> >> >> >> DATA AREA:            15.00  -84.00   23.00  -75.00
> >> >> >>
> >> >> >> Also, the file (satir_f24.gem) is posted at..
> >> >> >> http://tempest.meas.ncsu.edu/megan
> >> >> >> if you'd like to have the actual file.
> >> >> >>
> >> >> >> Basically, I'm working on generating a model-simulated satellite
> >> IR
> >> >> >> product, kind of like model-simulated radar.  But I'd like to use
> >> the
> >> >> >> satellite color tables to make an observational comparison.  I'd
> >> >> >> eventually like to try and use the model-simulated IR on grids of
> >> >> >> different navigation, but I started out with this grid because
> >> it's
> >> >> 4-km
> >> >> >> grid spacing, similar to satellite resolution.
> >> >> >>
> >> >> >> Thanks for your help! If there's any other information you'd like,
> >> >> just
> >> >> >> let me know.
> >> >> >>
> >> >> >> Megan
> >> >> >>
> >> >> >>
> >> >> >> > Megan,
> >> >> >> > I have been working on a gd2img program in the source
> >> distribution
> >> >> >> which
> >> >> >> > can convert a grid to an image using a specified satellite id
> >> and
> >> >> >> > channel number with a projection remapped into one that the
> >> image
> >> >> >> format
> >> >> >> > can handle.
> >> >> >> >
> >> >> >> > Can you send me some information about what grid
> >> >> projection/navigation
> >> >> >> > you have so I can see if this will be sufficient for your
> >> >> application?
> >> >> >> > The GINI format is somewhat limited specifying a projection and
> >> >> pixel
> >> >> >> > shape, so might have to use an AREA format instead. Please send
> >> >> this
> >> >> >> to
> >> >> >> > address@hidden.
> >> >> >> >
> >> >> >> > Steve Chiswell
> >> >> >> > Unidata User Support
> >> >> >> >
> >> >> >> > On Fri, 2007-03-30 at 10:25 -0400, Megan Gentry wrote:
> >> >> >> >> I have a GEMPAK grid file that I'd like to convert to image
> >> format
> >> >> >> and
> >> >> >> >> view using the satellite IR color tables.  I was wondering what
> >> >> the
> >> >> >> best
> >> >> >> >> way is to convert a grid file to an image file that can be
> >> viewed
> >> >> >> using
> >> >> >> >> the IR color tables?
> >> >> >> >> Thanks,
> >> >> >> >> Megan
> >> >> >> >>
> >> >> >> >> ----------------------------------------------------
> >> >> >> >> Megan S. Gentry
> >> >> >> >> Graduate Research Assistant, Forecasting Lab
> >> >> >> >> Department of Marine, Earth and Atmospheric Sciences
> >> >> >> >> North Carolina State University
> >> >> >> >>
> >> >> >> >> ===============================================================================
> >> >> >> >> To unsubscribe gembud, visit:
> >> >> >> >> http://www.unidata.ucar.edu/mailing-list-delete-form.html
> >> >> >> >> ===============================================================================
> >> >> >> > --
> >> >> >> > Steve Chiswell <address@hidden>
> >> >> >> > Unidata
> >> >> >> >
> >> >> >>
> >> >> >>
> >> >> >> ----------------------------------------------------
> >> >> >> Megan S. Gentry
> >> >> >> Graduate Research Assistant, Forecasting Lab
> >> >> >> Department of Marine, Earth and Atmospheric Sciences
> >> >> >> North Carolina State University
> >> >> >>
> >> >> >>
> >> >> >
> >> >> >
> >> >> > Ticket Details
> >> >> > ===================
> >> >> > Ticket ID: FWQ-357482
> >> >> > Department: Support
> >> >> > Priority: Normal
> >> >> > Status: Closed
> >> >>
> >> >>
> >> >> ----------------------------------------------------
> >> >> Megan S. Gentry
> >> >> Graduate Research Assistant, Forecasting Lab
> >> >> Department of Marine, Earth and Atmospheric Sciences
> >> >> North Carolina State University
> >> >>
> >> >>
> >> >
> >> >
> >> > Ticket Details
> >> > ===================
> >> > Ticket ID: FWQ-357482
> >> > Department: Support GEMPAK
> >> > Priority: Normal
> >> > Status: Closed
> >> >
> >> >
> >>
> >>
> >> ----------------------------------------------------
> >> Megan S. Gentry
> >> Graduate Research Assistant, Forecasting Lab
> >> Department of Marine, Earth and Atmospheric Sciences
> >> North Carolina State University
> >>
> >>
> >
> >
> > Ticket Details
> > ===================
> > Ticket ID: FWQ-357482
> > Department: Support GEMPAK
> > Priority: Normal
> > Status: Closed
> 
> 
> ----------------------------------------------------
> Megan S. Gentry
> Graduate Research Assistant, Forecasting Lab
> Department of Marine, Earth and Atmospheric Sciences
> North Carolina State University
> 
> 


Ticket Details
===================
Ticket ID: FWQ-357482
Department: Support GEMPAK
Priority: Normal
Status: Closed