Re: Adding geolocation to my data

[Mary Quinton <quintonmary@xxxxxxxxx>]

Hi Ben,
The topic you want to read about is "pseudo-
variables" in the Users Guide.  The definition
DEFINE AXIS/x=133.02  etc does define the
location the center of the grid box, and you
can see all the details by listing pseudo-variables
X, Y, XBOX, XBOXLO, and so on.

Also see DEFINE AXIS/EDGES if you want to use
the command to define the edges of the grid box.

Mary

--- Ben Burford <benb@restec.or.jp> wrote:
> Thank you Rob and Bill for your answers,
> 
> This worked fine.
> 
> I would like to confirm that the values in define
> axis (define axis/x=133.02
> E:135.52E:.25 xax) should be the centers of the edge
> pixels (and not edge to
>  edge).
> 
> This appears to be the case since "show axis xax"
> showed me:
> start 133.02E   end 135.52E
> Axis span (to cell edges) = 2.75
> 
> 2.75 is the full edge to edge, so things seem to be
> correct.
> 
> Thanks again,
> 
> Ben
> 
> 
> At 12:08 2004/07/13 -0500, you wrote:
> > 
> > Ben, just add the info to your axis definition:
> > 
> > DEFINE AXIS/x=10E:12.5E:0.25/UNITS=longitude xax
> > 
> > SHOW AXIS xax
> > 
> > Note that it's modulo. So for the x-axis, then
> type:
> > 
> > CANCEL AXIS/MODULO xax
> > 
> > By it won't be for the others.
> > 
> > go to 
> >
>
http://ferret.wrc.noaa.gov/Ferret/Documentation/Users_Guide/current/fer_ht
> ml.htm
> > 
> > and click on define axis to find out all the
> options under this command. 
> > It's all pretty obvious.
> > 
> > good luck,
> > Rob
> > 
> > 
> > 
> > 
> > At 09:21 PM 7/13/2004 +0900, you wrote:
> > >Hello Fellow Ferreters,
> > >
> > >I wrote the following script to read data from a
> file:
> > >
> >
>
>--------------------------------------------------------------------
> > >define axis/x=1:11:1 xax
> > >define axis/y=1:11:1 yax
> > >define axis/z=1:5:1 zax
> > >define grid/x=xax/y=yax/z=zax gax1
> >
>
>file/format=stream/type=i2/swap/skip=`2048/4`/grid=gax1/var=tb_scale
> F13S
> SMI
> > >1B_200107152131_wpo_low
> > >let tb=tb_scale/10
> > >set var/bad=9999 tb
> > >set var/title="Brightness Temperature" tb
> > >set var/units="Kelvin" tb
> > >   SAY "Var = tb"
> > >   SAY "Channels = 5"
> >
>
>-------------------------------------------------------------------------
> ---
> > >----------
> > >
> > >As you can see the data is on an 11 by 11 grid
> (which is 0.25 deg x 0.25 
> deg
> > >).  The two levels are data at different
> frequencies.
> > >
> > >I have additional information about this data:
> > >Grid size: 0.25 deg x 0.25 deg spacing (total
> area is 2.5 degree by 2.5 d
> egr
> > >ee)
> > >Lat/lon at the center of the Lower Left Pixel:
> lat 5,80  lon 133.02
> > >Observation Date Time(UT): 2001/07/15 21:31:18
> > >
> > >Can anyone show me how to insert lat/lon
> information, and possibly grid s
> pac
> > >ing (0.25 deg x 0.25 deg) information, into my
> script to geolocate the data
> > >onto a global grid?  I may also have to do
> something to specify high reso
> lut
> > >ion to look at this data.
> > >
> > >I want to be able to geolocate this data in order
> to use this data in Fer
> ret
> > >  in combination with other data, such as model
> output data.
> > >
> > >Thank you,
> > >
> > >Ben
> > 
> > Rob Scott
> > Research Associate
> > Institute for Geophysics
> > University of Texas at Austin
> > 4412 Spicewood Springs Rd.
> > Building 600
> > Austin, TX
> > 78759
> > 
> > Office: +1-512-471-0375
> >  
> 
> 



		
__________________________________
Do you Yahoo!?
Yahoo! Mail - 50x more storage than other providers!
http://promotions.yahoo.com/new_mail