Kurt's Weblog

...
"Networked Transport of RTCM via Internet Protocol" (Ntrip) stands for
an application level proto-col streaming "Global Navigation Satellite
System (GNSS)" data over the Internet. Ntrip is a generic, stateless
protocol based on the Hypertext Transfer Protocol (HTTP/1.1) and the
RealTime Streaming Protocol (RTSP). Ntrip is a reduced subset of
HTTP/RTSP required for real-time data streaming. Ntrip communication
usually takes place over HTTP/TCP/IP or RTSP/TCP/IP and RTP/UDP/IP
connections. The default port is TCP 2101, but other ports can be
used. This does not preclude Ntrip from being implemented on top of
any other protocol on the Internet, or on other networks like modem
connections or UDP-based transport respectively connection protocols.
...

...
One reinforced ship is already in service and four are being built by
Aker Yards for delivery by mid-2009, all using new hull designs that
allow for bow- or stern-first sailing, depending on the thickness of
the ocean surface. While Norilsk and others may benefit from the
melting, Borodin says the ship orders were driven by cost calculations
rather than global warming. Most of the Siberia route is navigable
from June through September, and it may become the first trans-Arctic
channel open to reinforced vessels, say scientists.
...

elif sys.platform == 'darwin':
    lib = find_library('geos_c')
    if lib is None:
        ## try a few more locations
        lib_paths = [
            # The Framework build from Kyng Chaos:
            "/Library/Frameworks/GEOS.framework/Versions/Current/GEOS",
        ]
        for path in lib_paths:
            if os.path.exists(path):
                lib = path
                break
            else:
                raise ImportError, "Could not find geos_c library"
    lgeos = CDLL(lib)
    free = CDLL(find_library('libc')).free
else:

% fink list -i libgeos
 i   libgeos2                   2.2.3-2007          Geometry Engine - Open Source
 i   libgeos2-shlibs            2.2.3-2007          Geometry Engine - Open Source
 i   libgeos3-shlibs            3.0.0-2             Geometry Engine - Open Source
 i   libgeosc1-shlibs           3.0.0-2             Geometry Engine - Open Source
% dpkg -L libgeosc1-shlibs | grep 'dylib'
/sw/lib/libgeos_c.1.4.1.dylib
% file /sw/lib/libgeos_c.1.4.1.dylib
/sw/lib/libgeos_c.1.4.1.dylib: Mach-O dynamically linked shared library i386
% otool -L /sw/lib/libgeos_c.1.4.1.dylib
/sw/lib/libgeos_c.1.4.1.dylib:
        /sw/lib/libgeos_c.1.dylib (compatibility version 6.0.0, current version 6.1.0)
        /sw/lib/libgeos-3.0.0.dylib (compatibility version 0.0.0, current version 0.0.0)
        /usr/lib/libstdc++.6.dylib (compatibility version 7.0.0, current version 7.4.0)
        /usr/lib/libgcc_s.1.dylib (compatibility version 1.0.0, current version 1.0.0)
        /usr/lib/libSystem.B.dylib (compatibility version 1.0.0, current version 111.0.0)
% otool -hv /sw/lib/libgeos_c.1.4.1.dylib
/sw/lib/libgeos_c.1.4.1.dylib:
Mach header
      magic cputype cpusubtype  caps    filetype ncmds sizeofcmds      flags
   MH_MAGIC    I386        ALL  0x00       DYLIB    13       1648   NOUNDEFS DYLDLINK TWOLEVEL WEAK_DEFINES BINDS_TO_WEAK NO_REEXPORTED_DYLIBS

PatchScript: <<
  perl -pi -e 's|/Library/Frameworks/GEOS.framework/Versions/Current/GEOS|/sw/lib/libgeos_c.1.4.1.dylib|g' shapely/geos.py
<<

% fink install shapely-py25
% ipython
In [1]: from shapely.geometry import Point
In [2]: point = Point(-106.0, 40.0)
In [3]: p = point.buffer(10.0)
In [4]: print p
POLYGON ((-96.0000000000000000 40.0000000000000000, -96.0481527332780303 39.0198285967043930, -96.1921471959676921 38.0490967798387203, -96.4305966426779122 37.0971532274553795, -96.7612046748871251 36.1731656763491074, -97.1807873565164471 35.2860326317400279, -97.6853038769745439 34.4442976698039800, -98.2698954663726312 33.65606
...

I'm doing an experiment with my recent papers: each has a Trac
instance, so that code and the article text (a LaTeX file) can be
shared between colaborators using Subversion. Data is hosted on an
Opendap server, running my trusty pydap. But the best part is that the
whole process, from downloading data to composing the article pdf, is
done through a Makefile. The Makefile creates a virtual environment,
downloads the required Python modules, and bootstraps the
analyses. Just type make and enjoy.
.
Now imagine if all scientists embraced the open-source philosophy of
open collaboration and real meritocracy and did the same.

...
Owner Joyce Russell said she thinks business has slowed down from
construction, parking and access issues. The problems started about a
year ago when the restaurant lost some parking spaces from the River
Walk construction behind it. This past year's Exeter Street
construction decreased business because people could not get through
town as easily and the upcoming Main Street project factored into the
decision to close, she said. The "frosting on the cake" are the recent
barricades preventing access out of Joyce's parking lot through the
Rivermoor Landing parking lot. Cars can now only get in and out from
Water Street off of Main Street, Russell said.
...
Breakfast is still available at Joyce's Kitchen at 13 Water St., until
the end of March. Hours are Monday through Friday from 6 a.m. until
noon. Saturdays from 6 a.m. until 1 p.m. and Sundays from 7 a.m. until
1 p.m.

...
The NOAA-UNH Joint Hydrographic Center was founded eight years ago to
find ways to process the massive amounts of data coming from these new
sonar systems at rates commensurate with data collection; that is, to
make the data ready for chart production as rapidly as the data could
be collected. As we begin to achieve this goal, our attention is also
turning to the wonderful opportunity to utilize this massive flow of
information into products that meet a wide range of needs, such as
safer navigation, marine habitat assessments, fisheries management and
national security.  In this short period of time, we have built a
vibrant Center with over 70 employees and an international reputation
as the place, "where the cutting edge of Hydrography is now located"
(Adam Kerr, Past Director of the International Hydrographic
Organization in Hydro International).
...

var lat= 43.08472
var lon=-70.775
...
maxExtent: new OpenLayers.Bounds(-20037508.34,-20037508.34,20037508.34,20037508.34),
                      maxResolution:156543.0399, units:'meters', projection: "EPSG:900913"} );
layerTilesAtHome = new OpenLayers.Layer.OSM.Osmarender("Osmarender");

% sudo apt-get install update-manager-core
% sudo do-release-upgrade
% sudo apt-get install language-pack-en
% reboot

% sudo apt-get autoremove
% sudo apt-get install update-manager-core

Errors were encountered while processing:
 ubuntu-artwork

  sudo apt-get update
  sudo apt-get upgrade
  sudo apt-get upgrade -f 
  sudo apt-get -f install

Preconfiguring packages ...
(Reading database ... 141787 files and directories currently installed.)
Unpacking gutsy-wallpapers (from .../gutsy-wallpapers_0.17_all.deb) ...
dpkg: error processing /var/cache/apt/archives/gutsy-wallpapers_0.17_all.deb (--unpack):
 trying to overwrite `/usr/share/backgrounds/warty-final-ubuntu.png', which is also in package edgy-wallpapers
dpkg-deb: subprocess paste killed by signal (Broken pipe)
Preparing to replace base-files 4ubuntu2 (using .../base-files_4.0.0ubuntu5_i386.deb) ...
Unpacking replacement base-files ...
Errors were encountered while processing:
 /var/cache/apt/archives/gutsy-wallpapers_0.17_all.deb

sudo dpkg  --install --force-overwrite /var/cache/apt/archives/gutsy-wallpapers_0.17_all.deb
sudo apt-get upgrade -f 
sudo reboot

cat /etc/lsb-release 
DISTRIB_ID=Ubuntu
DISTRIB_RELEASE=7.10
DISTRIB_CODENAME=gutsy
DISTRIB_DESCRIPTION="Ubuntu 7.10"

• No longer packaging the epydoc html in the tar. Create it locally if you want it.
• Cleanup of IMO met hydro message ice field. Description name fixes for other IMOs
• database.py: more verbose
• uscg.py: add get_station, get_contents
• django: add define_ais for creating AIS binary messages. Lots of work left to do
• doc/ais.bib: lots more references
• doc/templates: now includes python template/boilerplate file
• doc/configure-ais-logging.txt - new howto
• init.d/ais-forward: cleanup messages
• init.d/logais: added -m to mark time with no messages
• ais-db-rebuild-cache-tables: cleanup magicdate usage
• ais-net-to-postgis: add magicdate usage, add periodic clean even if nothing received
• ais-port-forward: minor cleanup, fix SBNMS stations, tweak bounding box
• ais_build_sqlite3.py is not working. why?

• ActiveCaptain
• MyCruisingLog.com
• CaptainWiki
• CruiserLog.com
• Rose Point Coastal Explorer 2.0 wiki - not the actial site

...
Dr. Torkel Soma, Principal Safety Consultant in DNV Maritime, says:
"DNV's statistics shows that a ship is twice as likely to be involved
in a serious grounding, collision or contact accident today compared
to only five years ago. In addition, estimates show that also the
costs of these accidents have doubled. Since this is the general trend
for the international commercial fleet, the maritime industry needs to
act on this immediately."
...

from django.db import models
#
class Type(models.Model):
    name = models.CharField(maxlength=20)
    def __str__(self):
        return self.name
    class Admin:
        pass
#
class Unit(models.Model):
    name = models.CharField(maxlength=10)
    shortdescription = models.CharField(maxlength=30)
    description = models.TextField()
    refURL = models.URLField()
    def __str__(self):
        name = self.name
        if len(name)==0: name = 'None'
        return name + ' - ' + self.shortdescription
    class Admin:
        pass
#
class Aismsg(models.Model):
    name = models.CharField(maxlength=20)
    description = models.TextField()
    msgnum = models.PositiveIntegerField()
    dac = models.PositiveIntegerField()
    fi = models.PositiveIntegerField()
    def __str__(self):
        return self.name + ' ('+str(self.dac)+':'+str(self.fi)+')'
    class Admin:
        pass
#
class Field(models.Model):
    aismsg = models.ForeignKey(Aismsg)
    name = models.CharField(maxlength=20)
    description = models.TextField()
    numberofbits = models.PositiveIntegerField()
    type = models.ManyToManyField(Type)
    unavailable = models.CharField(maxlength=120)
    units = models.ForeignKey(Unit)
    class Admin:
        pass

% createdb aisdef
% python manage.py syncdb

% psql aisdef
# \d
 public | define_ais_aismsg                 | table    | schwehr
 public | define_ais_field                  | table    | schwehr
 public | define_ais_field_type             | table    | schwehr
 public | define_ais_type                   | table    | schwehr
 public | define_ais_unit                   | table    | schwehr
# \d define_ais_field
 id           | integer                | not null default nextval('define_ais_field_id_seq'::regclass)
 aismsg_id    | integer                | not null
 name         | character varying(20)  | not null
 description  | text                   | not null
 numberofbits | integer                | not null
 unavailable  | character varying(120) | not null
 units_id     | integer                | not null

% psql aisdef
INSERT INTO define_ais_type (name) VALUES ('bool');
INSERT INTO define_ais_type (name) VALUES ('uint');
INSERT INTO define_ais_type (name) VALUES ('int');
INSERT INTO define_ais_type (name) VALUES ('udecimal');
INSERT INTO define_ais_type (name) VALUES ('decimal');
INSERT INTO define_ais_type (name) VALUES ('aisstr6');
INSERT INTO define_ais_type (name) VALUES ('binary');
INSERT INTO define_ais_unit (name,shortdescription,description,"refURL") VALUES ('m','meters','long description','http://physics.nist.gov/cuu/Units/meter.html');
INSERT INTO define_ais_unit (name,shortdescription,description,"refURL") VALUES ('kg','kilogram','long description','http://physics.nist.gov/cuu/Units/kilogram.html');
INSERT INTO define_ais_unit (name,shortdescription,description,"refURL") VALUES ('K','kelvin','long description','http://physics.nist.gov/cuu/Units/kelvin.html');
INSERT INTO define_ais_unit (name,shortdescription,description,"refURL") VALUES ('s','second','long description','http://physics.nist.gov/cuu/Units/second.html');
INSERT INTO define_ais_unit (name,shortdescription,description,"refURL") VALUES ('','practical salinity scale','long description','http://en.wikipedia.org/wiki/Salinity');
INSERT INTO define_ais_unit (name,shortdescription,description,"refURL") VALUES ('Pa','Pascals','long description','http://en.wikipedia.org/wiki/Pressure');

% python manage.py runserver
% open http://127.0.0.1:8000/

...
Few, if any, of the people who made Wine Therapy their home kept
backup copies of their own posts and comments. I'm sure they didn't
even think of it. I don't think of it, when I post to the various
boards and blogs and forums I frequent. Of course I know better, but I
think of these forums as extensions of my own computer -- until they
disappear.
.
As we rely on others to maintain our writings and our relationships,
we lose control over their availability. Of course, we also lose
control over their security, as MySpace users learned last month when
a 17-GB file of half a million supposedly private photos was uploaded
to a BitTorrent site.
.
In the early days of the web, I remember feeling giddy over the wealth
of information out there and how easy it was to get to. "The internet
is my hard drive," I told newbies. It's even more true today; I don't
think I could write without so much information so easily
accessible. But it's a pretty damned unreliable hard drive.
...

...
Senator John Kerry has proposed speed limits for specific vessels that
pass across the migratory paths of the critically endangered North
Atlantic right whale.
.
While the rule was first proposed in February 2007, Mr Kerry has
called for immediate action to be taken to protect the whales'
population.
.
"The Bush administration needs to stop dragging their feet, follow the
best available scientific evidence, and take immediate action to
protect endangered right whales," he said.
...

...
The CASCADE's allision with the bridge at Point San Quentin (San
Rafael-Richmond Bridge) presents a different situation. The tug
PACIFIC WOLF, an 111 foot, 4100 horsepower tug and the CASCADE an OPA
90 double hull 300 foot barge are reported to be owned and operated by
K-Sea Transportation. The tug DELTA DEANNE a tractor tug,108 feet in
length with 4400 horsepower is owned by BayDelta Maritime. The grieved
appears to be the bridge i.e., the state of California. Responsibility
and fault will have to be determined. Courts have defined the duties
and responsibilities of tug masters and therefore the question of who
was in charge will be of primary concern.
.
The number of recent accidents involving petroleum product is
unquestionably going to generate political rhetoric and probably some
unnecessary rule making unless the reason these accidents occurred is
determined. There is a commonality. The existing laws, rules,
regulations and procedures are more than adequate if adhered to. The
process of determining the cause, or reasons, is complex as it
indicates that those in the pilot house are not alone in their
failures. The operational policies, procedures and the principles of
proper seamanship and management are published and essentially known,
but they are frequently over looked, or misunderstood.
...
Except for ground vehicles, one should know that maritime accidents
exceed all other forms of transportation accidents The daily number of
maritime incidents is astonishing.
...

...
In fact, Verganti finds that companies that innovate successfully
(defined by the researchers as winning a design award) use a "see and
launch" strategy, employing "some sort of research processs," usually
trend finding and scenario building rather than focus groups, and then
launching products based on that research.
.
Less successful was the "launch and see" approach, where lots of
products are put out there, relying on customers' ability to hone in
on the good stuff so market numbers can decide. This reminds me of
Google's mythical "launch early, launch often" process. The study also
found "wait and see" - imitating trends that are successful in the
market - an expensive way to innovate since getting in on the game
early enough means you can't really call the winners, so you end up
copying a lot.
.
"Apple is a company that is pushed by a vision," Verganti says. "Steve
Jobs has said that the market doesn't always know what it
wants. Companies that do radical innovation do not listen to users;
they eventually value market feedback, but first they propose things
to the users."
...

% fink install imagemagick
% convert -delay 50 foo-???.png foo.gif
% open -a /Applications/Safari.app foo.gif

convert foo-???.png foo.mpeg
sh: mpeg2encode: command not found
convert: Delegate failed `"mpeg2encode" "%i" "%o"'.

% fink install mpeg2vidcodec

% convert -delay 50 foo-???.png foo.mpg
% open foo.mpg

% fink install ffmpeg
% ffmpeg -i foo-%d.png foo2.mpg
% open foo2.mpg

% fink apropos movie
% fink apropos mpeg
% fink apropos mpg
% fink apropos video

    Friday, February 22nd, 2:00 p.m.
      Chase Ocean Engineering Lab
           Dr. David Wiley
         Research Coordinator,
Stellwagon Bank National Marine Santuary
Endangered Whale Conservation Research in
the Stellwagon Bank National Marine Sanctuary
        Video Classroom 130
      All are welcome to attend

Platts LNG Daily [subscription required] reports that Excelerate
Energy has sent one of its LNG vessels, without an LNG cargo, to its
Northeast Gateway LNG terminal. Excelerate's COO told Platts that the
vessel would attach the submersible buoys to the vessel's onboard
regasification apparatus as part of the commissioning process for the
facility.

Now there is a report from University of Alberta which says that
Robert Boyle is the one who actually came up with the first working
model of sonar by 1917, beating out the other groups. He is being
ignored by history.

Stanford University today announced the largest increase in its
history for its financial aid program for undergraduates.
.
Under the new program, parents with incomes of less than $100,000 will
no longer pay tuition. Parents with incomes of less than $60,000 will
not be expected to pay tuition or contribute to the costs of room,
board and other expenses.
...

% ipython
from osgeo import gdal
f = gdal.Open('a.gif')
b = f.GetRasterBand(1)
c = b.GetRasterColorTable()
for i in range(c.GetCount()):
    print c.GetColorEntry(i)

(0, 0, 0, 255)
(255, 255, 255, 255)
(135, 135, 135, 255)
(220, 242, 252, 255)
(184, 229, 246, 255)
(133, 160, 165, 255)
(171, 192, 194, 255)
(89, 93, 91, 255)
(224, 227, 140, 255)
(252, 252, 148, 255)
(247, 244, 90, 255)
(53, 51, 40, 255)
(130, 113, 55, 255)
(198, 143, 80, 255)
(211, 163, 117, 255)
(255, 255, 255, 0)

gdalwarp -dstalpha -srcnodata 255 -dstnodata 0 \
  -s_srs EPSG:32619 -t_srs EPSG:4326 input-32619.tif output-4326.tif

• Imo_met_hydro (001:11): IMO meteorological and hydroglogical data. SN/Circ.236 Annex 2.
• imo_fairway_closed (001:13): IMO fairway closed. SN/Circ.236 Annex 2, page 4, App 3
• imo_tidal_window (001:13): IMO tidal window. SN/Circ.236 Annex 2, p5, Appl 3
• ris_waterlevel (???:24): River Information System water level report
• sls_waterlevel (366 316:1:3): St Lawrance Seaway water level message
• sls_weatherreport (366 316:1:2): St Lawrance Seaway wind information

sudo apt-get update
sudo apt-get upgrade
sudo apt-get install update-manager-core
sudo do-release-upgrade

dpkg --configure -a
apt-get update
apt-get upgrade

dpkg --force-depends -i keyspan-ubuntu-2.6.22.47.deb

mknod /dev/ttyUSB0 c 188 0
chgrp dialout /dev/ttyUSB?
chmod g+w /dev/ttyUSB?

ssh -p xxx somewhere -v
~schwehr/projects/noaadata/scripts/serial_logger.py --port=/dev/ttyUSB0

Because TiddlyWiki is a single HTML file, you've actually already
downloaded the entire software just by viewing this site. If you want
to be able to SaveChanges, you can save your own blank TiddlyWiki to
your local drive by right clicking on this link to empty.html and
selecting 'Save link as...' or 'Save target as...'. You can choose
where to save the file, and what to call it (but make sure that it's
saved in HTML format and with an HTML extension).

This graphic depicts the current status of USCG AIS projects underway
or planned, and their incorporation into the Common Operational
Picture (COP).  The Coast Guard Research and Development Center has
AIS projects in the Pacific Northwest, California, Hawaii, the Eastern
seaboard and Florida.  Vessel Traffic Service locations are
abbreviated beginning with "VTS" followed by letter abbreviations for
each location.  For instance, the Vessel Traffic Service in Prince
William Sound Alaska is identified by VTSPWS.  SCC Charleston is an
interagency Sector Command Center and SCC-J represents the Joint
Navy/Coast Guard Sector Command Centers in San Diego and Norfolk.  The
circles shown in solid green are currently incorporated into the COP.
Those shown in hatched green are operational AIS sites which do not
yet feed into the COP.  The blue circles are AIS projects not yet
operational.

On November 7, 2007, a container ship collided with a tower of the San
Francisco-Oakland Bay Bridge, spilling 58,000 gallons of the ship's
bunker fuel into the Bay. The incident has raised questions about the
role of maritime pilots and vessel traffic services (VTSs) in
preventing accidents in U.S. harbors, such as: Is there a need for
further independent oversight of pilot performance? Could VTSs operate
more like Air Traffic Control centers? Should the pilot and ship
captain be required to agree on a passage plan before transiting a
harbor? Because of the additional challenges of navigating large ships
through the constricted waters of a harbor, most U.S. ports require
shipping lines to hire a maritime pilot. Maritime pilots, through a
lengthy apprenticeship process and many years of experience, have
gained intimate knowledge of the navigational requirements of a
particular harbor. Despite the federal government's prominent role in
regulating interstate commerce, Congress has largely left it to the
various coastal states to regulate pilotage. Pilots are licensed by
the state for ships engaged in foreign trade ("registered" vessels
carrying international cargo), which accounts for the vast majority of
port calls. Typically a state pilot board oversees the hiring,
training, and performance of pilots, as well as setting pilotage rates
charged to the shipping lines. The U.S. Coast Guard has jurisdiction
only over pilots of ships engaged in domestic trade and the courts
have ruled that the Coast Guard does not have the authority to suspend
or revoke the license of pilots for violations while piloting a
foreign-trade ship. To assist pilots, the busiest U.S. ports have
established VTS stations to monitor ship traffic and provide relevant
information to pilots, such as the location of other ships. While
often compared to an air traffic control tower, a VTS is not directly
involved in the movement of vessels and is more accurately described
as an advisory service than a traffic control center. Ever-larger
ships, difficult or challenging slow-speed handling characteristics of
some of these ships, and rising port traffic that is predominantly
foreign-flagged have led to proposals concerning pilotage, VTSs, and
other safety-related navigation services in U.S. ports. Given the
federal interest in marine environmental protection and the Coast
Guard's mission to ensure the safety of shipping in U.S. waters, some
experts have advocated stronger federal oversight of pilots. They
recommend that the Coast Guard or a national commission establish
national standards for pilot training and proficiency or that the
Coast Guard be given disciplinary authority over state-licensed
pilots. State pilots resist greater federal oversight, arguing that
the unique geography and navigational requirements of each port
justifies local oversight. Whether VTSs should exert more direct
control over vessel movement is also raised as a safety measure, but
most acknowledge that an experienced on-board mariner is probably in
the best position to direct a vessel's movement. Requiring that a
pilot and ship captain first agree on a harbor passage plan,
investigating language difficulties between pilots and foreign crews,
and Coast Guard rotational staffing practices, are other issues
policymakers may examine in assessing the safety of ship navigation in
U.S. harbors.

sudo apt-get update
sudo apt-get upgrade
sudo apt-get install update-manager-core
sudo do-release-upgrade

cat /etc/lsb-release 
DISTRIB_ID=Ubuntu
DISTRIB_RELEASE=7.10
DISTRIB_CODENAME=gutsy
DISTRIB_DESCRIPTION="Ubuntu 7.10"

trac-admin /var/trac/projects/projectname upgrade

SAVANNAH, Ga. - The Coast Guard Captain of the Port for Brunswick,
Ga., has closed the port to deep-draft vessels Wednesday at 6:15
p.m. following a grounding that occurred earlier in the day.
.
A survey team from the U.S. Army Corps of Engineers is scheduled to
arrive Thursday morning to ensure the channel is safe for deep draft
vessel to transit. After the survey, the Captain of the Port will
re-evaluate the decision to keep the port closed to deep-draft
vessels.
.
At approximately 7:30 a.m. Wednesday the captain of the motor vessel
Centaurus Leader reported it had run aground on the north side of the
Brunswick River while it was transiting in-bound to the Colonels
Island terminal.
...

...
Bathymetric soundings taken last year showed the foot of Alaska's
continental slope extending more than 100 nautical miles farther from
the U.S. coast than previously believed, according to the National
Oceanic and Atmospheric Administration.
.
"We found evidence that the foot of the slope was much farther out
than we thought," said Larry Mayer, the chief scientist for the
expedition last year. "That was the big discovery."
...

BOSTON - The Unified Command at Coast Guard Sector Boston is
continuing to monitor and assist the Liquefied Natural Gas tanker
Catalunya Spirit, which lost propulsion and became disabled east of
Cape Cod Monday morning.
.
The Catalunya Spirit is about 33 miles east of Chatham, Mass., and is
being towed by tugs Freedom and Iona McAllister, traveling about four
miles-per-hour in a northwesterly direction, parallel to
shore. Additionally, the Coast Guard Cutter Escanaba, the Maine
Responder and Atlantic Salvor are also on scene.
.
The Catalunya Spirit has restored some electrical power to the
vessel's switchboards; however, full power has not yet been restored
to the vessel.
...

...
Teekay Corporation, which operates the Catalunya Spirit, contracted
four tugboats and three technical representatives to assist the
tanker. The tugboats Freedom and Liberty arrived on scene with the
tanker about 10 p.m., Monday. The Freedom placed the Catalunya Spirit
under tow, stopping the tanker's southeasterly movement, and is slowly
heading west to a temporary offshore location, while conducting
repairs.
.
The tug Iona McAllister is also on scene standing by. The salvage tug,
Atlantic Salvor, is expected to arrive at 1 p.m., Tuesday to assist as
necessary.
...

% xdpyinfo
vendor string:    The X.Org Foundation
vendor release number:    10300000
.Org version: 1.3.0
...

% uname -a
Darwin laptop 9.2.0 Darwin Kernel Version 9.2.0: Tue Feb  5 16:13:22 PST 2008; root:xnu-1228.3.13~1/RELEASE_I386 i386 i386 MacBookPro3,1 Darwin

BOSTON - The Coast Guard is monitoring and assisting a Liquefied
Natural Gas tanker, which lost propulsion and is disabled and adrift
with 29 people aboard approximately 35 miles east of Chatham, Mass.,
today.
.
The vessel is not in immediate danger and is drifting away from Cape
Cod.
.
The 933-foot tanker Catalunya Spirit, carrying a full load of LNG from
Trinidad and Tobago to Boston, became disabled about 3 a.m., this
morning.
.
The tanker crew is making hourly situation reports to Sector
Boston. Additionally, Air Station Cape Cod is planning to deliver a
Coast Guard marine inspector and a technical representative to the
Catalunya Spirit today. The Coast Guard Cutter Escanaba is on scene
with the tanker to assist with communication.
.
Teekay Corporation, the company that operates Catalunya Spirit, has
contracted two tug boats to assist. The tugs are expected to arrive on
scene with the Catalunya Spirit at 11:30 p.m. A second Coast Guard
marine inspector is aboard one of the tugs.
.
Current weather is 30-35 knot winds with 12 foot seas and is forecast
to diminish overnight.
.
"As part of our response plan, we've notified our National Strike
Team, and we're coordinating salvage and pollution response assets,"
said Capt. Gail Kulisch, Captain of the port of Boston. "A very
comprehensive safety system has been developed by the Coast Guard in
conjunction with port partners and the shipping company to minimize
the risk to the marine environment and public safety."

New Arctic sea floor data released today by the University of New
Hampshire and the National Oceanic and Atmospheric Administration
suggests that the foot of the continental slope off Alaska is more
than 100 nautical miles farther from the U.S. coast than previously
assumed.
.
The data, gathered during a recent mapping expedition to the Chukchi
Cap some 600 nautical miles north of Alaska, could support U.S. rights
to natural resources of the sea floor beyond 200 nautical miles from
the coast.
.
"We found evidence that the foot of the slope was much farther out
than we thought," said Larry Mayer, expedition chief scientist and
co-director of the Joint Hydrographic Center at UNH. "That was the big
discovery."
...

...
The $4.1 billion request represents a 7.7 percent increase over the
President's FY 2008 request and a 5.2 percent increase over the FY
2008 budget passed by Congress.
...
Specific budget highlights include:
    * $31.8 million to implement new requirements under the
      Magnuson-Stevens Fisheries Conservation and Management
      Reauthorization Act of 2006 designed to end overfishing
    * $7 million to improve our understanding of the oceans through
      the Integrated Ocean Observing System
    * $4 million to clean up dangerous marine debris
    * $5.4 million to restore streams and rivers by removing obsolete dams
    * $5.3 million to improve hurricane forecasts

opendiff template.py template2.py

SYNOPSIS
     opendiff file1 file2 [-ancestor ancestorFile] [-merge mergeFile]
     opendiff dir1 dir2 [-ancestor ancestorDirectory] [-merge mergeDirectory]
...
     If the -ancestor flag is given, FileMerge will compare the two files or
     directories to a common ancestor.  This is useful if two people indepen-
     dently modify copies of a single original file or directory.

open /Developer/Applications/Utilities/FileMerge.app

    def __str__(self):
        return '%s last position' % self.name

    class Admin:
        list_display = ('userid','name','cog','cg_timestamp')

Use Pythonesque coding pattern
This is very related to the previous item, obviously. Python has some
well-known constructs to handle some situations. Get to know and
understand them.  An example: as you might know, Python has no
switch/case construct. There's a very neat way to implement this
though by simply using a dict and function objects (or lambda
functions). An example:



def handle_one():
    return 'one'
def handle_two():
    return 'two'
def handle_default():
    return 'unknown'
cases = {
    'one': handle_one,
    'two': handle_two,
    'three': lambda: 'three',
}
for i in ('one', 'two', 'three', 'four', ):
    handler = cases.get(i, handle_default)
    print handler()
We're using the dict.get method here, which can take an optional
'default' argument. Pretty neat, huh?

#!/usr/bin/env python
# Consider this template file public domain.  The __license__ is just an example.
# Replace values as needed
__author__    = 'Kurt Schwehr'
__version__   = '$Revision: 4799 $'.split()[1]
__revision__  = __version__ # For pylint
__date__ = '$Date: 2006-09-25 11:09:02 -0400 (Mon, 25 Sep 2006) $'.split()[1]
__copyright__ = '2008'
__license__   = 'GPL v3'
__contact__   = 'kurt at ccom.unh.edu'
#__deprecated__
__doc__ ='''
Example python file that is all tricked out.  Designed for epydoc,
unittest and doctest.  Please keep updating to make this the best
possible template file. Decimate these requirements to meet what you need.
@requires: U{Python<http://python.org/>} >= 2.5
@requires: U{epydoc<http://epydoc.sourceforge.net/>} >= 3.0.1
@requires: U{psycopg2<http://http://initd.org/projects/psycopg2/>} >= 2.0.6
@undocumented: __doc__ parser
@since: 2008-Feb-09
@status: under development
@organization: U{CCOM<http://ccom.unh.edu/>} 
@bug: epydoc --check complains about unittest additions
@bug: Bug two
@todo: logging, fuzzing, etc
@todo: To make your code executable after check-in (to allow for doc and unittest), do svn propset svn:executable ON yourfile.py
@todo: "svn propedit svn:keywords yourfile.py", then put in "Date Revision"
@todo: What to make of "pylint --with-line-length=150" ... does this really lead to better code?
@todo: Run pychecker on your files.
@see: U{python idioms<http://jaynes.colorado.edu/PythonIdioms.html>} - Read this before modifying any code here.
'''

Previously charted as a small sea knoll, Mount Healy abruptly rises
more than 3,000 meters from the floor of the Arctic Ocean.

• Real-time AIS works, but needs a little bit more tweaking to get time thresholding
• Dropped the .py for an initial batch of scripts including serial-logger, xy2kml, xy_pos_filter, and xymt2kml. Need to do the rest
• Moved from scripts to data subdir: gsc.dat, merg5.dat, sbnms-*.dat, stellwagen*.dat
• Created a contrib directory. Added doc from Maitland Bottoms
• Changed around the MANIFEST.in. Can anyone suggest a better way to build the MANIFEST.in?
• Working on switching from GPL v2 to v3
• Minor updates to the waterlevel2 message based on feedback from ACOE - additional datums
• Starting on updates for epydoc 3.0.1 - better front matter tagging in python files. e.g. see database.py
• database.py: added rebuild_track_lines and rebuild_last_position functions
• uscg.py: added __eq__ method and a return BitVector; can now rebuild the NMEA
• Renaming of the ais-db* and ais-port-forward command. Moved from scraps to scripts
• Moved aisxml2wxpy and aisxmlbinmsg2py from ais to scripts. They should not a part of the ais module!
• Added initial django-ais subdir that doesn't yet do anything
• Much improved fink info file. Still needs a lot of work

0020000001000010E6C051D30925D1DA0B4044A79AE924F228

#!/usr/bin/env python
import psycopg2
import psycopg2.extensions
import GeoTypes
#
factory = GeoTypes.OGGeoTypeFactory()
parser = GeoTypes.HEXEWKBParser(factory)
#
c = psycopg2.connect("dbname='ais' user='schwehr'")
cu = c.cursor()
#
cu.execute('SELECT position FROM position LIMIT 1;')
row = cu.fetchone()
print row[0]
#
parser.parseGeometry(row[0])
geom = factory.getGeometry()
#
print type(geom)
print geom
print 'x:',geom.getX()
print 'y:',geom.getY()
#
print 'Should match...'
cu.execute('SELECT AsText(position) FROM position LIMIT 1;')
row = cu.fetchone()
print row[0]

0020000001000010E6C051D30925D1DA0B4044A79AE924F228
<class 'GeoTypes._OGPoint.OGPoint'>
GeometryFromText('POINT(-71.297433 41.309415)',4326)
x: -71.2974333333
y: 41.309415
Should match...
POINT(-71.2974333333333 41.309415)

% python manage.py sqlall msgs

ais=# \d track_lines
                                   Table "public.track_lines"
 Column  |         Type          |                           Modifiers                           
=========+=====================--+===============================================================
 ogc_fid | integer               | not null default nextval('track_lines_ogc_fid_seq'::regclass)
 userid  | integer               | 
 name    | character varying(20) | 
 track   | geometry              | 
Indexes:
    "track_lines_pkey" PRIMARY KEY, btree (ogc_fid)
Check constraints:
    "enforce_dims_track" CHECK (ndims(track) = 2)
    "enforce_geotype_track" CHECK (geometrytype(track) = 'LINESTRING'::text OR track IS NULL)
    "enforce_srid_track" CHECK (srid(track) = 4326)

% django-admin.py startproject mysite
% cd mysite
% python manage.py runserver

% grep DATA settings.py
DATABASE_ENGINE = 'postgresql_psycopg2'
DATABASE_NAME = 'ais'
DATABASE_USER = ''
DATABASE_PASSWORD = ''
DATABASE_HOST = ''
DATABASE_PORT = ''

% python manage.py inspectdb > foo.py

from django.db import models
#
class Position(models.Model):
    key = models.IntegerField(primary_key=True)
    messageid = models.IntegerField()
    repeatindicator = models.IntegerField()
    userid = models.IntegerField()
    navigationstatus = models.IntegerField()
    rot = models.IntegerField()
    sog = models.FloatField(max_digits=4, decimal_places=1)
    positionaccuracy = models.IntegerField()
    cog = models.FloatField(max_digits=4, decimal_places=1)
    trueheading = models.IntegerField()
    timestamp = models.IntegerField()
    regionalreserved = models.IntegerField()
    spare = models.IntegerField()
    raim = models.BooleanField()
    state_syncstate = models.IntegerField()
    state_slottimeout = models.IntegerField()
    state_slotoffset = models.IntegerField()
    cg_r = models.CharField(maxlength=15)
    cg_sec = models.IntegerField()
    cg_timestamp = models.DateTimeField()
    position = models.TextField() # This field type is a guess.
    class Meta:
        db_table = 'position'
#
class LastPosition(models.Model):
    key = models.IntegerField(primary_key=True)
    userid = models.IntegerField()
    name = models.CharField(maxlength=20)
    cog = models.IntegerField()
    sog = models.TextField() # This field type is a guess.
    cg_timestamp = models.DateTimeField()
    position = models.TextField() # This field type is a guess.
    class Meta:
        db_table = 'last_position'

% psql ais
ais=# \d last_position
                                       Table "public.last_position"
    Column    |            Type             |                          Modifiers                          
============--+===========================--+============================================================
 key          | integer                     | not null default nextval('last_position_key_seq'::regclass)
 userid       | integer                     | 
 name         | character varying(20)       | 
 cog          | integer                     | 
 sog          | real                        | 
 cg_timestamp | timestamp without time zone | 
 position     | geometry                    | 
Indexes:
    "last_position_pkey" PRIMARY KEY, btree ("key")
Check constraints:
    "enforce_dims_position" CHECK (ndims("position") = 2)
    "enforce_geotype_position" CHECK (geometrytype("position") = 'POINT'::text OR "position" IS NULL)
    "enforce_srid_position" CHECK (srid("position") = 4326)

Traceback (most recent call last):
  File "manage.py", line 11, in <module>
    execute_manager(settings)
  File "/sw/lib/python2.5/site-packages/django/core/management.py", line 1672, in execute_manager
    execute_from_command_line(action_mapping, argv)
  File "/sw/lib/python2.5/site-packages/django/core/management.py", line 1563, in execute_from_command_line
    from django.utils import translation
  File "/sw/lib/python2.5/site-packages/django/utils/translation/__init__.py", line 3, in 
    if settings.USE_I18N:
  File "/sw/lib/python2.5/site-packages/django/conf/__init__.py", line 28, in __getattr__
    self._import_settings()
  File "/sw/lib/python2.5/site-packages/django/conf/__init__.py", line 55, in _import_settings
    self._target = Settings(settings_module)
  File "/sw/lib/python2.5/site-packages/django/conf/__init__.py", line 83, in __init__
    raise EnvironmentError, "Could not import settings '%s' (Is it on sys.path? Does it have syntax errors?): %s" % (self.SETTINGS_MODULE, e)
EnvironmentError: Could not import settings 'ais.settings' (Is it on sys.path? Does it have syntax errors?): No module named settings

I've just uploaded a new version to the Cheeseshop. The syntax now
allows the combined use of
  x weeks ago
  x days ago
  x hours ago
  x minutes ago
  x seconds ago
So it works with your examples but also "1 week and 2 days ago", "2
days, 1 hour and 34 seconds ago", etc.

ais-db-rebuild-cache-tables --all --track-start-time='2008-01-22' \
--last-position-start-time='1 hour ago'

  logging.basicConfig(filename = options.log_file
                      , level  = options.log_level
                      )

...
But NSF rocked the ocean sciences community by deciding to pull its
$331 million Ocean Observatories Initiative (OOI) off the table
pending a final review later this year.
...

Traceback (most recent call last):
  File "./make_waterlevel_ais.py", line 38, in 
    stations = Stations.ActiveStations(forceCache=True)
  File "/Users/schwehr/Desktop/production/src/noaadata/noaadata/stations.py", line 233, in __init__
  File "lxml.etree.pyx", line 2342, in lxml.etree.parse
  File "parser.pxi", line 1700, in lxml.etree._parseDocument
  File "parser.pxi", line 1732, in lxml.etree._parseMemoryDocument
  File "parser.pxi", line 1633, in lxml.etree._parseDoc
  File "parser.pxi", line 785, in lxml.etree._BaseParser._parseDoc
  File "parser.pxi", line 450, in lxml.etree._ParserContext._handleParseResultDoc
  File "parser.pxi", line 534, in lxml.etree._handleParseResult
  File "parser.pxi", line 476, in lxml.etree._raiseParseError
lxml.etree.XMLSyntaxError: Namespace prefix soapenv on Envelope is not defined, line 1, column 54
Traceback (most recent call last):

/ssh:username@myhost.unh.edu:/home/schwehr/projects/noaadata/scripts/Makefile

pylint serial_logger.py  --max-line-length=200
No config file found, using default configuration
************* Module serial_logger
C:  1: Missing docstring
W:  1: Redefining built-in '__doc__'
C: 34:createDaemon: Invalid name "createDaemon" (should match [a-z_][a-z0-9_]{2,30}$)
C: 39:createDaemon: Invalid name "e" (should match [a-z_][a-z0-9_]{2,30}$)
C: 48:createDaemon: Invalid name "e" (should match [a-z_][a-z0-9_]{2,30}$)
W: 52:createDaemon: Access to a protected member _exit of a client class
W: 55:createDaemon: Access to a protected member _exit of a client class
C: 64:createDaemon: Invalid name "fd" (should match [a-z_][a-z0-9_]{2,30}$)
W: 67:createDaemon: Except doesn't do anything
C: 78:dateStr: Invalid name "dateStr" (should match [a-z_][a-z0-9_]{2,30}$)
C: 84:dateStr: Operator not preceded by a space
    t=time.gmtime()
     ^
C: 84:dateStr: Invalid name "t" (should match [a-z_][a-z0-9_]{2,30}$)
C: 85:dateStr: Operator not preceded by a space
    d='%04d-%02d-%02d' % t[:3]
     ^
C: 85:dateStr: Invalid name "d" (should match [a-z_][a-z0-9_]{2,30}$)
W: 89:start: Redefining name 'options' from outer scope (line 217)
C: 95:start: Comma not followed by a space
    ser = serial.Serial(options.port,options.baud,timeout=options.timeout)

+===========+=======+
|type       |number |
+===========+=======+
|convention |49     |
+===========+=======+
|refactor   |2      |
+===========+=======+
|warning    |5      |
+===========+=======+
|error      |0      |
+===========+=======+

Your code has been rated at 0.17/10

Your code has been rated at 8.39/10 (previous run: 8.05/10)

from django.http import HttpResponse
import datetime
#
def last_position(request):
    import psycopg2 as psycopg
    connectStr = "dbname='ais' user='schwehr' host='localhost'"
    cx = psycopg.connect(connectStr)
    cu = cx.cursor()
    q = 'SELECT userid,name,cog,cg_timestamp,AsText(position) FROM last_position ORDER by userid;'
    cu.execute(q)
    htmlLines=['<html><body>',]
    htmlLines.append('<table border="1">')
    htmlLines.append('<tr bgcolor="orange"><td>userid</td><td>name</td><td>cog (&deg;)</td><td>timestamp (UTC)</td><td>position</td></tr>')
    rownum=0
    for row in cu.fetchall():
        if rownum%2==0: htmlLines.append('<tr bgcolor="gray">')
        else: htmlLines.append('<tr>')
        rownum+=1
        userid,name,cog,cg_timestamp,pos = row
        htmlLines.append('<td>'+str(userid)+'</td>')
        htmlLines.append('<td>'+str(name)+'</td>')
        htmlLines.append('<td>'+str(cog)+'</td>')
        htmlLines.append('<td>'+str(cg_timestamp)+'</td>')
        htmlLines.append('<td>'+str(pos)+'</td>')
        htmlLines.append('</tr>')
    htmlLines.append('</table>')
    htmlLines.append('</body></html>')
    html = '\n'.join(htmlLines)
    print html
    return HttpResponse(html)

from django.conf.urls.defaults import *
import ablog.views 
urlpatterns = patterns(
    ''
    ,(r'^vessels/latest/$',ablog.views.last_position)
)

fink install coreutils-default
isightcapture `date +'%Y%m%d-%H%M'`.jpg
isightcapture -t png `date +'%Y%m%d-%H%M'`.png

0.25000
0
0
-0.25000
360246.83000
4771048.51000

# WGS 84 / UTM zone 19N
<32619> +proj=utm +zone=19 +ellps=WGS84 +datum=WGS84 +units=m +no_defs  <>

gdalwarp -s_srs EPSG:32619 -t_srs EPSG:4326 backscatter-utmz19.tif backscatter-4326.tif

In a combined effort of the National Oceanic and Atmospheric
Administration and private industry (Applanix and Riegl), two
different but complimentary technologies were tested in the marine
survey field to find a more efficient, higher-accuracy and safer
alternative for traditional shoreline verification methods. This
article describes the tests and the encouraging good results with this
videogrammetry and laser scanning technology.
...
Results
Results from this initial test in an urban port setting were very
encouraging. While sea conditions for these tests were ideal, there
were no noticeable motion artefacts found in the data. Riegl has found
over and over again that 10-centimetre point cloud accuracy is quite
typical for their scanner data integrated with an IMU. In this Norfolk
test, we have accomplished better results: less than 5 centimetres (1
sigma) with respect to reference objects, after post-processing. This
is typically not range-dependent under 1,000 metres.
.
It was found that in this environment it was rarely necessary to
position things at distances of greater than 250-300 metres. The
resolution of these systems was impressive, providing almost
photographic quality point clouds. The resolution was such that
cleats, bits and other pier-side hardware could easily be identified
from the point cloud.
.
A data logger on the laser system recorded raw data. The Riegl
processing includes an extraction step that produces a file with
points referenced to the laser's coordinate system. This export file
is then merged with the POS MV's file of the smooth best estimated
trajectory using the Applanix utility POS point cloud. The merge is
based on GPS time to the microsecond.
.
Summary and Conclusion
Two different but complementary technologies were tested aboard the
NOAA vesselBay Hydrographerin an effort to find a more efficient,
higher accuracy and safer alternative for traditional shoreline
verification. Both systems performed well, were relatively easy to
integrate into the existing hydrographic suite of equipment and
provided a means of positioning various items from a boat at a safe
distance offshore. The OCR capabilities of the videogrammetry system
have great productivity potential for surveying port facility
infrastructure. The long-range laser system holds promise for those
areas where the shoreline structure is entirely natural and items such
as rocks and islets dominate the survey. Each provides benefits where
the other is weakest (i.e. using the video to classify features in the
point cloud for which there may be insufficient data).  The final
result of this experiment is the production of what will inevitably
become the nautical chart of the future, a merge of previously
acquired Norfolk Harbor airborne laser data with the vessel-based
laser data. The 2D rendering in Figure 4 does not do it justice;
however, the picture is still quite stunning. Visualisation in three
dimensions, on a high-definition colour monitor is truly impressive.
...

Digipilot is an official Digital Nautical Publication, developed and
issued by the Hydrographer of the Royal Netherlands Navy.  Digipilot
provides the mariner with the necessary information to plan his voyage
in an interactive and versatile manner.  It contains Sailing
Directions and supplementary Port Information in various interactive
geographic and thematic layers, but also tidal information, Deep
Draught Planning Guide (Greenwich Buoy to Europoort) and the Mariners
Routing Guide for the Southern North Sea.  The product covers the
Netherlands and Belgian coasts including the main inland waterways up
to and including the ports of Antwerpen (Belgium) and Emden (Germany).
The content of this publication has been derived from traditional
existing paper nautical publications issued by the Hydrographer of the
Royal Netherlands Navy and contains data derived from additional
sources.  Digipilot is an official equivalent of paper publications
HP1 (Netherlands Coast Pilot) and HP33 (tide tables and stream
atlases) in accordance with SOLAS V2.2 and SOLAS V19.2.1.5.  It may
replace these traditional paper publications provided that appropriate
back-up arrangements are available (e.g. print facility or 2nd
installed program).

#!/bin/bash
xterm -r -sl 1000 -j -sb $* &

#!/bin/bash
/sw/bin/emacs -fn 6x13 -bg black -fg white $* &

echo "xterm*geometry: 80x40" | xrdb -merge
echo "xterm*foreground: white" | xrdb -merge
echo "xterm*background: black" | xrdb -merge
echo "emacs*foreground: white" | xrdb -merge
echo "emacs*background: black" | xrdb -merge
echo "emacs*font: 6x13" | xrdb -merge

! Is the comment character
! XTerm settings that are easier on the eyes
xterm*background: black
xterm*foreground: wheat
xterm*cursorColor: orchid
xterm*scrollBar: true
xterm*reverseWrap: true
xterm*eightBitInput: false
!xterm*titleBar: false