segysqlshiptrack (version 1.6, 2005/10/13)
index
/Users/schwehr/projects/xcore/src/segy-py/segysqlshiptrack.py

Pull the bathymetry from a GMT grid file for the ship track of a seismic line.
 
FIX: add error checking
FIX: allow inserting the z values into the database
FIX: get my PyGMT changes into the official PyGMT release and into fink
 
     Copyright (C) 2005  Kurt Schwehr
 
    This program is free software; you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation; either version 2 of the License, or
    (at your option) any later version.
 
    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.
 
    You should have received a copy of the GNU General Public License
    along with this program; if not, write to the Free Software
    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 
$Id: segysqlshiptrack.py,v 1.6 2005/10/13 20:22:11 schwehr Exp $

 
Modules
       
PyGMT
os
sqlite
sys
verbosity

 
Functions
       
interpolateXY(points)
Cheesy linear interpolation between the first and last x,y pairs
 
e.g. interpolateXY ([[0,1],[1,1],[2,1],[3,],[4,2]]) will only look
at [0,1] and [4,2].
lonlat2xy(grd, lon, lat)
Get the x,y cell offsets for a lon lat position in a gmt grid file
 
grd - PyGMT grid file
lon - decimal longitude (-180 to +180)
lat - decimal latitude  (-90 to +90)
 
FIX: this needs some error checking.  What if we are off the grid?
shipTrack(options, cx)
Pull the ship track out of a database and get the z values from a
gmt grid file.  Optionally interpolate between navigation (GPS) updates.
 
options - command line options plus a projection instance
cx - database connection

 
Data
        ALWAYS = 0
BOMBASTIC = 4
TERSE = 1
TRACE = 2
VERBOSE = 3
VERSION = '0.31'
__author__ = 'Kurt Schwehr'
__credits__ = 'Paul Henkart for the idea of using the grid file.\nMagnus Hagdorn for writing PyGMT.'
__date__ = '2005/10/13'
__version__ = '1.6'
myparser = <optparse.OptionParser instance>

 
Author
        Kurt Schwehr

 
Credits
        Paul Henkart for the idea of using the grid file.
Magnus Hagdorn for writing PyGMT.