<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Progress Toward Automation of Salmon Escapement Counts</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Mark Shortis</string-name>
          <email>mark.shortis@rmit.edu.au</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Alaska Department of Fish and Game</institution>
          ,
          <addr-line>Homer, Alaska</addr-line>
          ,
          <country country="US">USA</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>GSR_3 Geospatial Science Research 3. School of Mathematical and Geospatial Science, RMIT University</institution>
        </aff>
        <aff id="aff2">
          <label>2</label>
          <institution>School of Mathematical and Geospatial Sciences, RMIT University</institution>
        </aff>
      </contrib-group>
      <abstract>
        <p>The number of adult salmon that escape predators and fishing to successfully return to their natal streams to spawn, known as escapement counts, are vital to ensure sustainability of the salmon fisheries. Sufficient numbers of salmon must reach the spawning grounds each year to ensure the long term viability of the various species. Once a critical threshold is reached, the fishery is opened to commercial exploitation. To ensure the viability of salmon in Alaska, the Alaska Department of Fish and Game conducts escapement counts using a variety of techniques. One well established technique is to use remote video stations to record salmon passing upstream and periodically use human operators to view the video sequences to determine the escapement counts. This paper presents the initial stages of the development of a technique to automate the escapement counts and thereby relieve the demand for human operators. The technique uses a global segmentation and image morphology operations followed by a Kalman filter to identify and track the candidate salmon as they pass through the region of interest on the video recording. Candidates that complete the journey across the region of interest are counted. Finally, some of the future challenges posed by real time processing, water surface turbulence and lighting changes are discussed.</p>
      </abstract>
      <kwd-group>
        <kwd>salmon</kwd>
        <kwd>escapement</kwd>
        <kwd>counts</kwd>
        <kwd>video</kwd>
        <kwd>image processing</kwd>
        <kwd>identification</kwd>
        <kwd>tracking</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>Introduction</title>
      <p>Salmon escapement counts are vital to ensure sustainability of the fisheries in Alaska. Sufficient numbers of
salmon must reach the spawning grounds each year to ensure the long term viability of the various species. In
remote areas the salmon must escape their natural predators such as salmon sharks, seals and orca in open water,
and then bears and eagles at the entrances to streams and rivers. In more populated areas the salmon are also
subject to recreational and subsistence fishing, usually at or near stream and river mouths. Catches are limited
by daily or seasonal limits per person/household. Only when sufficient numbers have passed through the
predation zones and will reach the spawning grounds are the fisheries opened to commercial exploitation.
Commercial operators apply pressure to open the fisheries as early as possible to maximise the fishing season,
consequently escapement counts have to be both accurate and reliable.</p>
      <p>
        A number of techniques are used to monitor escapement. Traditionally, artificial weirs and manual counts by a
resident fishery officer have been the mainstay for salmon escapement
        <xref ref-type="bibr" rid="ref16 ref17">(Zimmerman and Zabkar, 2007)</xref>
        .
However, manual counts at every stream, especially in remote locations, are not practical or sustainable.
Techniques such as hydroacoustics
        <xref ref-type="bibr" rid="ref8">(Maxwell, 2007)</xref>
        and counting towers (Woody, 2007) are successful in some
applications. Aerial surveys are widely used
        <xref ref-type="bibr" rid="ref5">(Jones et al., 2007)</xref>
        , but are a relative index rather than a count of
individual fish because accurately counting thousands of salmon from the air is not feasible
        <xref ref-type="bibr" rid="ref1">(Bue et al., 1998)</xref>
        .
To bolster the accuracy and reliability of escapement counts and allow remote streams to be monitored
efficiently, the Alaska Department of Fish and Game (ADF&amp;G) has pioneered the use of remote video stations
        <xref ref-type="bibr" rid="ref10">(Otis and Dickson 2002)</xref>
        . Time lapse video is captured at a rate of 2-5 frames per second to a digital video
recorder. The video recordings are transported to the ADF&amp;G office in Homer for periodic manual counting.
Whilst this approach is accurate and reliable, it is labour intensive and staff are subject to intense pressure as
escapement counts approach the required threshold at each site. Automation of the process is a key aim to
improve the efficiency of the counts.
      </p>
      <p>This paper describes the progress made to develop an automatic video counting system that operates in real time
or rapidly analyses pre-recorded videos. The developed algorithm has the capability detect, track and count the
salmon passing through a defined region of interest in the video frame. The first step in the processing is to
segment the fish objects in the current video frame using a global threshold technique. Morphological image
processing (dilation, erosion and median filtering) is then used to eliminate speckle noise, eliminate small
regions which are likely not related to fish objects, and connect parts of the same fish object separated by water
surface turbulence. The third step of the system is to track the fish objects using a Kalman filter to overcome
occlusions and gaps in identification. The Kalman filter predicts the track up to a certain number of frames,
beyond which the object is eliminated as a false positive. Finally, fish objects with a full track across the region
of interest are counted.</p>
    </sec>
    <sec id="sec-2">
      <title>Remote Video Stations</title>
      <p>
        In the 1990s the ADF&amp;G at Homer began developing stand-alone remote video salmon escapement monitoring
stations to supplement aerial surveys of remote, clear water streams
        <xref ref-type="bibr" rid="ref10">(Otis and Dickson 2002)</xref>
        . A typical system
is based on an aluminium quadrapod sited next to a shallow area of the stream that affords a good view of the
salmon migration (figure 1). A video camera in a weatherproof housing is mounted at the top of the quadrapod,
viewing the stream obliquely from the south to minimise reflected surface glare from the sun. A light coloured
substrate panel is used to provide sufficient contrast to make the salmon visible in most conditions of lighting
and surface turbulence. The camera is outfitted with a circular polarizer to optimise the clarity of the view of the
substrate panel area.
      </p>
      <p>The video is recorded at 2-5 frames per second to a digital video recorder (DVR) inside a bear-proof strongbox.
The DVR receives a NTSC composite (analogue) video signal from the camera but records digital video
sequences onto a hot-swappable hard disk drive. Electrical power is provided by a bank of 12V lead-acid
batteries connected in parallel. Wind and solar power generators are used to maintain the battery charge level.
An automatic alarm system and satellite communications relay is used to alert ADF&amp;G if equipment failures
occur.</p>
      <p>The field of view of the overhead cameras is optimised to match the panel area as far as possible. Figure 2
shows the fields of view from the systems at Chenik Creek and Mikfik Creek, both on the western side of Cook
Inlet. Video sites are chosen based on a range of factors such as accessibility, stream morphology, minimisation
of surface turbulence and reflections, availability of solar, wind or hydropower, and many others. As noted
previously, the optimal direction for the camera view is to the north and this requirement can lead to a
compromise between the field of view and the resolution of the view. The orientation of the format of the
camera is physically limited by the weatherproof housing design, so in general the substrate panel is not aligned
with the highest resolution axis of the camera. At some sites other requirements, for example a desire to monitor
human traffic near the site, forces a wider than optimum field of view.</p>
      <p>
        The DVR has four input channels and some sites also have a low-light camera mounted on the quadrapod or an
underwater camera with a view across the stream bed. The underwater camera is used primarily to positively
identify fish enumerated with the above-water camera and to apportion the total escapement among co-occurring
species with overlapping run timing. The low light camera does not have a polarizing filter and is used primarily
during twilight periods to minimize missing counts caused by the 4-5 hours of darkness during June to August.
Previous research that compared video counts to those derived by a weir at Delight Creek in South-central
Alaska indicated that less than 2% of the salmon escapement passed the video counter during the hours of
darkness
        <xref ref-type="bibr" rid="ref10">(Otis and Dickson 2002)</xref>
        .
      </p>
      <p>The hard disks have sufficient capacity to record video for 2-3 weeks. Regular visits to the remote sites are
carried out to check the equipment and swap hard drives, with visits increasing in frequency as escapement
counts approach the threshold necessary to open the commercial fishery. Maintenance visits are commonly
pooled with aerial survey flights to efficiently use the float plane travel, as this is a significant component of the
budget for escapement counts.</p>
      <p>The hard disks are returned to the ADFG office for analysis. The disks are mounted in a multi-screen video
station for replay. If multiple cameras are available then all captured videos can be shown on the screens with
synchronous replay. An operator manually records the numbers of salmon passing over the substrate panel using
a hand counter. The operator can fast forward through periods of inactivity or step slowly through the individual
frames of groups of fish passing over the panel. An operator can review a full day of video in as little as 10-20
minutes during periods of low fish count. However, 1-3 hours are often needed to review peak passage days
which may involve 2000-4000 fish. Salmon runs span as many as 40 days and may require a total of 15 hours of
analysis by an operator, dependent on the frequency profile and total count. A typical profile of daily salmon
escapement counts is shown in figure 3.</p>
    </sec>
    <sec id="sec-3">
      <title>Video Imagery</title>
      <p>The quality of the video imagery recorded by the remote stations is, in most conditions, satisfactory for an
operator to identify and count the salmon. The resolution of the video is limited to 720 by 480 pixels by the
NTSC broadcast signal, which was adopted in the design of the remote stations for compatibility with the low
cost surveillance cameras and security industry DVRs that are available off-the-shelf. This technology has many
advantages, particularly that it has a proven track record, is very robust and there is a vast knowledge base of
systems and technological solutions. In addition, the ability to readily inject a time code into the video image or
the audio track allows a permanent record of the date and time, and allows accurate synchronisation of multiple
cameras.</p>
      <p>Nevertheless, the NTSC resolution is a limitation that is more evident for some stations, such as the Mikfik
Creek site shown in figure 2. The field of view could be narrowed using a longer focal length lens, but in some
cases the width of the stream is too wide to permit sufficiently high resolution images.</p>
      <p>A further limitation on the quality of the video imagery is the video compression. The DVR units used in the
remote video stations use a proprietary Motion JPEG compression to maximise the recording time by the units.
The compression, combined with the relatively low number of frames per second, allows the units to be
deployed for 2-3 weeks of 24 hour recording. This is particularly important during off-peak periods when the
hard drives are not swapped out as frequently.</p>
      <p>As is typically the case for surveillance DVRs, the operator cannot disable the video compression but does have
some control over the level of compression. However the compression levels are high even within the allowable
range. For example the lowest to the highest compression ratios per frame are 30:1 to 65:1. Visually there is
little discernible difference between the extremes of quality recorded by the DVR. An apparent characteristic of
the proprietary Motion JPEG algorithm is that adjacent odd and even horizontal lines in the image are adjusted
to have similar RGB intensities, leading to an effective halving of the vertical resolution and further deterioration
of the image quality through loss of resolution. Unfortunately testing with alternative compression algorithms,
such as MPEG, produced even greater deterioration in image quality and also had a greater propensity for
dropped frames when the video was reviewed at faster replay speeds.</p>
    </sec>
    <sec id="sec-4">
      <title>Automation</title>
      <p>As discussed above, operators can review an entire 24 hour video recording in a span of 10-180 minutes
depending on the number of salmon to be counted. This level of efficiency is admirable for a human operator,
however over a full season with many remote video stations the total time for review can be measured in
hundreds of hours.</p>
      <p>
        Computer-based image analysis has been used in marine and freshwater science for many tasks associated with
detection and identification of fish from still images and video sequences. For example,
        <xref ref-type="bibr" rid="ref15">Walther et al. (2004)</xref>
        detected candidates for estimates of biomass using motion analysis from a ROV-based video camera.
        <xref ref-type="bibr" rid="ref14">Tillett et
al. (2000)</xref>
        used image analysis for fitting to a species-specific template of the body shape of fish to extract
specific metrics. There are many other examples of the use of image analysis for biological studies
        <xref ref-type="bibr" rid="ref12 ref16 ref17 ref3 ref7">(Derry and
Elliot 1997, Kocak et al. 1999, Rova et al. 2007)</xref>
        that use or build on these two fundamental processes of
candidate detection and template matching.
      </p>
      <p>
        More recently,
        <xref ref-type="bibr" rid="ref13">Spampinato et al. (2010)</xref>
        reports on a system for detection, tracking, counting and classifying fish
that is in continuous operation on a sub-tropical reef in Taiwan waters. C
        <xref ref-type="bibr" rid="ref2">haralampidis et al. (2012</xref>
        ) describe a
NOAA research project investigating recognition and tracking of fish in video sequences. Both systems broadly
use subtraction of the current image from a ‘clean’ background image; segmentation of the image using
threshold or histogram techniques; morphological operations to eliminate noise and merge adjacent regions; and
then edge detection and boundary definition to detect and identify candidate fish objects. C
        <xref ref-type="bibr" rid="ref2">haralampidis et al.
(2012</xref>
        ) then uses a Kalman filter to predict the locations of each candidate in subsequent frames.
Full automation of the salmon counts is technically feasible based on a similar process of segmentation of the
image to isolate the objects in motion, morphological operations to eliminate noise and merge regions, and then
tracking of individual fish through a region of interest (RoI) to compile the count. The RoI can be placed within
the video frame by an operator to ensure that all fish are detected and counted.
      </p>
      <p>The behaviour of the salmon facilitates the RoI approach, rather than attempting the virtually impossible task of
uniquely identifying each fish, because the salmon rarely deviate from a direct path across the substrate panels.
The fish often pause at the threshold of the substrate panel because it is an unnatural feature, but then swim
directly across the panel to the natural creek bed on the far side. The salmon are driven by instinct to swim
upstream, so multiple factors militate against the fish repeatedly swimming or hesitating over the panel.
Accordingly, a feasible approach is to track the fish objects across the RoI and count each detected object that
completes the passage.</p>
    </sec>
    <sec id="sec-5">
      <title>Video Processing</title>
      <p>The video processing and tracking approach described in the previous section has been implemented using
Matlab1 as a development and testing platform. Once a mature, generalised solution is developed the intent is to
implement the video processing algorithms as a stand-alone application that can be deployed on any PC. Each
step in the process is described in more detail below.</p>
      <p>The first step in the processing is to define the RoI within the video frame. This approach is feasible because of
the fixed nature of the field of view of the remote video stations. The Matlab function roipoly allows the
operator to specify a polygonal RoI within an image from the video sequence and creates a binary image that is
used as a mask. The mask is then applied to every subsequent video frame to eliminate areas outside the RoI
from further processing. The RoI reduces the computational burden on the system, increases the effectiveness of
the segmentation performance and avoids the false detection of fish objects in the cluttered image areas outside
of the RoI.</p>
      <p>An example of the operator selected polygon and the subsequent mask are shown in figures 4 and 5. The RoI
should match the substrate panel, but should avoid the edges and shallow sections near the shore to provide
uniform background intensity. If the salmon are hesitating, or entering and then re-entering the substrate panel
area, the downstream edge of the RoI can be positioning further inside the substrate panel area to minimise
tracking errors.</p>
      <p>
        Segmentation is the next step in the video processing. This commences with converting the video frame from
RGB colour into an average greyscale intensity for each pixel using the ColorSpaceConverter function. The
medfilt2 median filter function is applied next to remove noise and smooth the image. Under the simple
assumption that the substrate panel will always be a relatively high intensity value and the fish objects will be a
relatively low intensity value, a simple average of the minimum and maximum intensities in the RoI is used to
globally segment the image and create a binary image of the candidate fish objects. This approach is based on
the principle of a bi-modal distribution of intensity used in the Otsu method
        <xref ref-type="bibr" rid="ref11">(Otsu, 1979)</xref>
        . In mathematical
terms:
and
(1)
1 http://www.mathworks.com.au/products/matlab/
where and are the maximum and minimum values of the intensities in an intensity image, dst(x,y) is the
pixel value of destination binary image and src(x,y) is the pixel value of the source greyscale image. Two
examples of original RGB images and the segmented, binary images are shown in figure 6. The binary images
include fish candidate blobs, non-fish blobs and noise.
The image morphology dilation and erosion functions Dilate and Close, using a 3 by 3 pixel square structuring
element, are then applied to remove noise and merge adjacent blobs. The next step in the process is the removal
of any remaining small blobs within the RoI. A threshold of size can be set empirically based on the apparent
size of the salmon within the video images, which is in turn based on the set up at each remote video station.
The difference in the image span of the salmon between Chenik Creek and Mikfik Creek is readily apparent in
figure 2. An example of noise and small blob removal is shown in figure 7.
      </p>
      <p>The Matlab BlobAnalysis function is then applied, with empirically determined minimum and maximum blob
sizes, to analyse the RoI and identify the fish candidates. The blob analysis identifies the connected groups of
pixels and computes characteristics such as the area, centroid and bounding box. The effectiveness of the blob
analysis algorithm relies completely on the quality of the segmentation. A poor segmentation will result in no
candidates being detected or many false positives due to lighting conditions or noise in the image. In turn this
will lead to inaccurate results from the tracking and counting processes, described in the next section.</p>
    </sec>
    <sec id="sec-6">
      <title>Tracking and Counting</title>
      <p>
        In computer vision moving objects are tracked for variety of applications such as activity recognition, traffic
monitoring and automotive safety
        <xref ref-type="bibr" rid="ref4 ref9">(Hu et al., 2004; Moeslund et al., 2006)</xref>
        . In this project, tracking is used for
the purpose of counting the fish candidates. Each object is tracked from the point it enters the RoI to the point
where it exits the RoI. The foundation of this approach is that migrating salmon are the only likely cause of an
object transiting the RoI from the downstream edge to the upstream edge. If an object completes the journey
across the RoI it is classified as a salmon and added to the count. Static objects and transitory candidates caused
by surface ripples or reflections are automatically eliminated.
      </p>
      <p>
        The multiple objects in the RoI are tracked with a Kalman filter
        <xref ref-type="bibr" rid="ref6">(Kalman, 1960)</xref>
        . The Kalman filter has been
used for many applications in trajectory modelling, both within the 2D image space for computer vision and
within 3D world space for positioning. The Kalman filter is a recursive solution to discrete-data linear filtering,
used to estimate the state of a linear system assumed to be distributed as a Gaussian function. The filter supports
estimation of past, present and future states, and can compute the estimates when the precise nature of the
modelled system is unknown
        <xref ref-type="bibr" rid="ref16">(Welch and Bishop, 2004)</xref>
        .
      </p>
      <p>In this project a Kalman filter has been used for object tracking by predicting the position from the previous
information and verifying the existence of the object at the predicted location. In Matlab, the estimates are
computed using the vision.KalmanFilter system object. The association of candidate detection to the same
object is based solely on motion. The motion of each track is estimated by the Kalman filter which determines
the likelihood of each candidate detection being assigned to each track.</p>
      <p>Track maintenance becomes an important aspect of counting the salmon. In any given frame, some candidate
detections may be assigned to tracks, while other candidates and tracks may remain unassigned. The assigned
tracks are updated using the corresponding detections. An unassigned candidate detection begins a new track.
Candidates must be tracked for a certain number of consecutive frames to be assigned to a track. Tracks that are
unassigned for a certain number of consecutive frames are deleted. False positives caused by noise or surface
turbulence tend to result in short-lived tracks. The limit is set empirically and in most circumstances a value of
five consecutive frames is effective.</p>
      <p>Assigning object detections in the current frame to existing tracks is done by minimizing a cost function. The
cost is defined as the negative log-likelihood of a detection corresponding to a track. The algorithm involves
two steps:</p>
      <sec id="sec-6-1">
        <title>Step 1:</title>
      </sec>
      <sec id="sec-6-2">
        <title>Step 2:</title>
        <p>Compute the cost of assigning every detection to each track using the distance method of the
vision.KalmanFilter system object. The cost takes into account the Euclidean distance between the
predicted centroid of the track and the centroid of the detection. It also includes the confidence of
the prediction, which is maintained by the Kalman filter. The results are stored in an MxN matrix,
where M is the number of tracks and N is the number of detections.</p>
        <p>Solve the assignment problem represented by the cost matrix using the assignDetectionsToTracks
function. The function takes into account the cost matrix and the cost of not assigning any detections
to a track.</p>
        <p>The value for the cost of not assigning detection to a track depends on the range of values returned by the
distance method of the vision.KalmanFilter. This value must be tuned experimentally. Setting it too low
increases the likelihood of creating a new track, and may result in track fragmentation. Setting it too high may
result in a single track corresponding to a series of separate moving candidates.</p>
        <p>Figure 8 shows an example of fish candidates mapped onto the original frame, so that the results of segmentation
and counting can be shown within the RoI. The top left annotation of the bounding box shows the fish counter
and top right annotation shows the life of the candidate in number of frames. When the candidate is occluded or
fragmented by surface turbulence, or is touching an edge of the RoI, the top right annotation indicates that the
location is ‘predicted’. The green cross marks show the centroids of the candidate blobs.</p>
        <p>Candidates that complete an unbroken track across the RoI are added to an accumulation count of salmon.</p>
        <p>In the limited testing of good quality video sequences that has been conducted to date, the counter has been
100% reliable. However this perfect reliability has to be strongly qualified by the quality of the video and the
heuristic nature of the approach. To establish a baseline of effectiveness, the initial testing has concentrated on
video sequences in which the salmon are clearly visible and counts can be easily verified by a human operator.
Moreover, factors such as the limits on blob size, the upper limit for the elimination of small blobs and the
selected number of frames used to eliminate false positives must all be chosen very carefully.</p>
      </sec>
    </sec>
    <sec id="sec-7">
      <title>Conclusions</title>
      <p>The approach to salmon escapement counts described in this paper has the potential to relieve the demand for
manual counts by operators. The approach has been successfully applied to pre-recorded video sequences.
However it is yet to be tested on live video and has been applied in a limited range of water surface conditions.
Further, testing is needed to validate the results against reference data provided by experienced human observers.
The most difficult circumstances for salmon counts are twilight periods, wind-induced surface turbulence,
suninduced surface glare and discoloration of the substrate panel. Figure 9 illustrates the range of different
conditions experienced over the course of a single day at Chenik Creek. It is evident that, as conditions
deteriorate, identification and counts of salmon become commensurately more difficult. Despite the use of a
polarizing filter, surface turbulence from the wind and upstream obstacles can exacerbate sun-induced surface
glare to the point where the ability to see and count fish is severely compromised. The substrate panels can
become discoloured by deposition of silt or algae growth, or may be interfered with by bears or floods.
Fortunately, interference with the panels is infrequent and discoloration is rectified during regular maintenance
visits.</p>
      <p>Further development and testing is necessary to extend the effectiveness of the approach to the full range of
conditions that will be experienced at the remote video stations and decrease the reliance on empirically chosen
factors to control the algorithms. The critical step in the video processing is the segmentation within the RoI, as
this determines the success of the following steps. Different approaches, such as segmentation within the colour
space and frame differencing from an averaged background image will be the subject of future research.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          <string-name>
            <surname>Bue</surname>
            ,
            <given-names>B.G.</given-names>
          </string-name>
          ,
          <string-name>
            <given-names>S.M.</given-names>
            <surname>Fried</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S.</given-names>
            <surname>Sharr</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D.G.</given-names>
            <surname>Sharp</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.A.</given-names>
            <surname>Wilcock</surname>
          </string-name>
          , and
          <string-name>
            <given-names>H.J.</given-names>
            <surname>Geiger</surname>
          </string-name>
          .
          <year>1998</year>
          .
          <article-title>Estimating salmon escapement using area-under-the-curve, aerial observer efficiency, and stream-life estimates: the Prince William Sound pink salmon example</article-title>
          .
          <source>North Pacific Anadromous Fish Commission Bulletin</source>
          <volume>1</volume>
          :
          <fpage>240</fpage>
          -
          <lpage>250</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          <string-name>
            <surname>H.</surname>
          </string-name>
          ,
          <year>2012</year>
          .
          <article-title>Feature analysis for classification of fish in underwater video</article-title>
          .
          <source>Final Report, LA Board of Regents Contract NASA</source>
          (
          <year>2011</year>
          )
          <article-title>-</article-title>
          STENNIS-02,
          <year>December 2012</year>
          .
          <volume>19</volume>
          pages.
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          <string-name>
            <surname>Derry</surname>
            ,
            <given-names>J. F.</given-names>
          </string-name>
          and
          <string-name>
            <surname>Elliott</surname>
            ,
            <given-names>C. J. H.</given-names>
          </string-name>
          (
          <year>1997</year>
          ). Automated 3
          <article-title>-D tracking of video-captured movement using the example of an aquatic mollusk</article-title>
          .
          <source>Behavior Research Methods, Instruments and Computers</source>
          <volume>29</volume>
          (
          <issue>3</issue>
          ):
          <fpage>353</fpage>
          -
          <lpage>357</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          <string-name>
            <surname>Hu</surname>
            ,
            <given-names>W.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Tan</surname>
            ,
            <given-names>T.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Wang</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          and
          <string-name>
            <surname>Maybank</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <year>2004</year>
          .
          <article-title>A survey on visual surveillance of object motion and behaviors</article-title>
          ,
          <source>IEEE Transactions on Systems, Man, and Cybernetics</source>
          ,
          <volume>34</volume>
          (
          <issue>3</issue>
          ):
          <fpage>334</fpage>
          -
          <lpage>352</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          <string-name>
            <surname>Jones</surname>
            <given-names>III</given-names>
          </string-name>
          , E.L.,
          <string-name>
            <given-names>S.</given-names>
            <surname>Heinl</surname>
          </string-name>
          , and
          <string-name>
            <given-names>K.</given-names>
            <surname>Pahlke</surname>
          </string-name>
          .
          <year>2007</year>
          .
          <string-name>
            <given-names>Aerial</given-names>
            <surname>Counts</surname>
          </string-name>
          .
          <source>Pages</source>
          <volume>399</volume>
          -410 in D.H.
          <string-name>
            <surname>Johnson</surname>
            ,
            <given-names>B.M.</given-names>
          </string-name>
          <string-name>
            <surname>Shrier</surname>
            ,
            <given-names>J.S.</given-names>
          </string-name>
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          <string-name>
            <surname>Kalman</surname>
            ,
            <given-names>R. E.</given-names>
          </string-name>
          ,
          <year>1960</year>
          .
          <article-title>A new approach to linear filtering and prediction problems</article-title>
          .
          <source>Transaction of the ASME Journal of Basic Engineering</source>
          ,
          <fpage>35</fpage>
          -
          <lpage>45</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          <string-name>
            <surname>Kocak</surname>
            ,
            <given-names>D.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Vitoria-Lobo</surname>
            ,
            <given-names>N.</given-names>
          </string-name>
          and
          <string-name>
            <surname>Widder</surname>
            ,
            <given-names>E.</given-names>
          </string-name>
          ,
          <year>1999</year>
          .
          <article-title>Computer Vision Techniques for Quantifying, Tracking, and Identifying Bioluminescent Plankton</article-title>
          .
          <source>IEEE Journal of Oceanic Engineering</source>
          <volume>24</volume>
          (
          <issue>1</issue>
          ):
          <fpage>81</fpage>
          -
          <lpage>95</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          <string-name>
            <surname>Maxwell</surname>
            ,
            <given-names>S.L.</given-names>
          </string-name>
          <year>2007</year>
          . Hydroacoustics: Rivers. Pages
          <fpage>133</fpage>
          -152 in D.H.
          <string-name>
            <surname>Johnson</surname>
            ,
            <given-names>B.M.</given-names>
          </string-name>
          <string-name>
            <surname>Shrier</surname>
            ,
            <given-names>J.S. O</given-names>
          </string-name>
          <string-name>
            <surname>'Neil</surname>
            ,
            <given-names>J.A.</given-names>
          </string-name>
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          <string-name>
            <surname>Moeslund</surname>
            ,
            <given-names>T. B.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Hilton</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          and
          <string-name>
            <surname>Krüger</surname>
            ,
            <given-names>V.</given-names>
          </string-name>
          ,
          <year>2006</year>
          .
          <article-title>A survey of advances in vision-based human motion capture and analysis</article-title>
          .
          <source>Computer Vision</source>
          and Image Understanding,
          <volume>104</volume>
          (
          <issue>2-3</issue>
          ):
          <fpage>90</fpage>
          -
          <lpage>126</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          <string-name>
            <surname>Otis</surname>
            ,
            <given-names>E.O.</given-names>
          </string-name>
          and
          <string-name>
            <given-names>M.</given-names>
            <surname>Dickson</surname>
          </string-name>
          .
          <year>2002</year>
          .
          <article-title>Improved salmon escapement enumeration using remote video and timelapse recording technology</article-title>
          ,
          <source>Exxon Valdez Restoration Final Report (Restoration Project 00366)</source>
          , Alaska Department of Fish and Game, Division of Commercial Fisheries, Homer, Alaska.
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          <string-name>
            <surname>Otsu</surname>
            ,
            <given-names>N.</given-names>
          </string-name>
          ,
          <year>1979</year>
          .
          <article-title>A threshold selection method from gray-level histograms</article-title>
          .
          <source>IEEE Trans. Systems, Man, and Cybernetics</source>
          <volume>9</volume>
          (
          <issue>1</issue>
          ):
          <fpage>62</fpage>
          -
          <lpage>66</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          <string-name>
            <surname>Rova</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Mori</surname>
            ,
            <given-names>G.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Dill</surname>
            ,
            <given-names>L. M.</given-names>
          </string-name>
          ,
          <year>2007</year>
          .
          <article-title>One fish, two fish, butterfish, trumpeter: Recognizing fish in underwater video</article-title>
          .
          <source>Proceedings, IAPR Conference on Machine Vision Applications</source>
          , May 16-18,
          <year>2007</year>
          , Tokyo, Japan.
          <source>ISBN 978-4-901122-07-8</source>
          pp
          <fpage>404</fpage>
          -
          <lpage>407</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          <string-name>
            <surname>Spampinato</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Giordano</surname>
            ,
            <given-names>D.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Di</surname>
            <given-names>Salvo</given-names>
          </string-name>
          ,
          <string-name>
            <given-names>R.</given-names>
            ,
            <surname>Chen-Burger</surname>
          </string-name>
          , J., Fisher, R.,
          <string-name>
            <surname>Nadarajan</surname>
            ,
            <given-names>G.</given-names>
          </string-name>
          ,
          <year>2010</year>
          .
          <article-title>Automatic fish classification for underwater species behaviour understanding</article-title>
          .
          <source>Proceedings of First ACM International Workshop on Analysis and Retrieval of Tracked Events and Motion in Imagery Streams. Firenze</source>
          , Italy, ACM:
          <fpage>45</fpage>
          -
          <lpage>50</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          <string-name>
            <surname>Tillett</surname>
            ,
            <given-names>R.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>McFarlane</surname>
            ,
            <given-names>N.</given-names>
          </string-name>
          and
          <string-name>
            <surname>Lines</surname>
            ,
            <given-names>J.</given-names>
          </string-name>
          <year>2000</year>
          .
          <article-title>Estimating dimensions of free-swimming fish using 3D point distribution models</article-title>
          .
          <source>Computer Vision</source>
          and Image Understanding,
          <volume>79</volume>
          ,
          <fpage>123</fpage>
          -
          <lpage>141</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref15">
        <mixed-citation>
          <string-name>
            <surname>Walther</surname>
            ,
            <given-names>D.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Edgington</surname>
            ,
            <given-names>D. R.</given-names>
          </string-name>
          and
          <string-name>
            <surname>Koch</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          <year>2004</year>
          .
          <article-title>Detection and tracking of objects in underwater video</article-title>
          .
          <source>In IEEE International Conference on Computer Vision</source>
          and Pattern
          <string-name>
            <surname>Recognition</surname>
            ,
            <given-names>F.M.</given-names>
          </string-name>
          Titsworth (ed.). New York: Institute of Electrical and
          <source>Electronics Engineers</source>
          <volume>1</volume>
          ,
          <fpage>544</fpage>
          -
          <lpage>549</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref16">
        <mixed-citation>
          <string-name>
            <surname>Welch</surname>
            ,
            <given-names>G.</given-names>
          </string-name>
          and
          <string-name>
            <surname>Bishop</surname>
            ,
            <given-names>G.</given-names>
          </string-name>
          ,
          <year>2004</year>
          .
          <article-title>An introduction to the Kalman Filter</article-title>
          , In University of North Carolina at Chapel Hill, Department of Computer Science.
          <source>Tech. Rep</source>
          .
          <volume>95</volume>
          -
          <fpage>041</fpage>
          , available at http://www.cs.unc.edu/~welch/media/pdf/kalman_intro.pdf Woody,
          <string-name>
            <surname>C.A.</surname>
          </string-name>
          <year>2007</year>
          .
          <article-title>Tower counts</article-title>
          .
          <source>Pages</source>
          <volume>363</volume>
          -384 in D.H.
          <string-name>
            <surname>Johnson</surname>
            ,
            <given-names>B.M.</given-names>
          </string-name>
          <string-name>
            <surname>Shrier</surname>
            ,
            <given-names>J.S. O</given-names>
          </string-name>
          <string-name>
            <surname>'Neil</surname>
            ,
            <given-names>J.A.</given-names>
          </string-name>
          <string-name>
            <surname>Knutzen</surname>
            ,
            <given-names>X.</given-names>
          </string-name>
        </mixed-citation>
      </ref>
      <ref id="ref17">
        <mixed-citation>
          <string-name>
            <surname>Zimmerman</surname>
            ,
            <given-names>C.E.</given-names>
          </string-name>
          and
          <string-name>
            <given-names>L.M.</given-names>
            <surname>Zabkar</surname>
          </string-name>
          .
          <year>2007</year>
          . Weirs. Pages
          <fpage>385</fpage>
          -398 in D.H.
          <string-name>
            <surname>Johnson</surname>
            ,
            <given-names>B.M.</given-names>
          </string-name>
          <string-name>
            <surname>Shrier</surname>
            ,
            <given-names>J.S. O</given-names>
          </string-name>
          <string-name>
            <surname>'Neil</surname>
            ,
            <given-names>J.A.</given-names>
          </string-name>
          <string-name>
            <surname>Knutzen</surname>
            ,
            <given-names>X.</given-names>
          </string-name>
          <string-name>
            <surname>Augerot</surname>
          </string-name>
          , T.A.
          <string-name>
            <surname>O'Neil</surname>
          </string-name>
          , and T.N. Pearsons, editors.
          <source>Salmonid Field</source>
          Protocols Handbook:
          <article-title>Techniques for assessing status and trends in salmon and trout populations</article-title>
          .
          <source>American Fisheries Society</source>
          , Bethesda, Maryland.
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>