Difference between revisions of "HOWTO Drifter Data Processing"

From Locl
Jump to navigationJump to search
Line 17: Line 17:
  
 
<h2>Folder:  .../kirk/user/0_Info</h2>
 
<h2>Folder:  .../kirk/user/0_Info</h2>
 +
<b>Purpose: Track drifter usage history</b>
 
<li><b>On deployment day</b>: be sure to fill out .../kirk/user/0_Info/<b>drifter_hist.xls</b> (drifter history) with the
 
<li><b>On deployment day</b>: be sure to fill out .../kirk/user/0_Info/<b>drifter_hist.xls</b> (drifter history) with the
 
performance of the drifters to track if there are any consistent problems (<span style="border:1px solid red;background-color:#FFEEEE;padding:1px;">Note that this file is now out of date</span>).<br>
 
performance of the drifters to track if there are any consistent problems (<span style="border:1px solid red;background-color:#FFEEEE;padding:1px;">Note that this file is now out of date</span>).<br>
Line 24: Line 25:
  
 
<h2>Folder:  .../kirk/user/1_get_data_from_web</h2>
 
<h2>Folder:  .../kirk/user/1_get_data_from_web</h2>
 +
<b>Purpose: Fetch the appropriate data from pacificgyre.com</b>
 
<li>Edit <b>data_dates.m</b> with the start/stop date/times to get.<br>
 
<li>Edit <b>data_dates.m</b> with the start/stop date/times to get.<br>
Notes:<br>
+
Notes:<br>
 
<ul>
 
<ul>
 
   <li>Use the [http://pacificgyre.com http://pacificgyre.com] website (log: <i>ohlmann</i>, pass: <i>drifterguy</i>) if you need to get an estimation of
 
   <li>Use the [http://pacificgyre.com http://pacificgyre.com] website (log: <i>ohlmann</i>, pass: <i>drifterguy</i>) if you need to get an estimation of
Line 36: Line 38:
 
<li>Run <b>get_data.m</b>
 
<li>Run <b>get_data.m</b>
 
<ul>
 
<ul>
<li>Be sure the <b>IRIDIUM</b> value in it is set correctly
+
  <li>Be sure the <b>IRIDIUM</b> value in it is set correctly
<li>Data will be saved in '..../kirk/data/Drifter_Anal/2_data/<i>yyyymmdd</i>' folder as a <i>yyyymmdd_drifter_dat.log</i> file.
+
  <li>Data will be saved in '..../kirk/data/Drifter_Anal/2_data/<i>yyyymmdd</i>' folder as a <i>yyyymmdd_drifter_dat.log</i> file.
<li>Now, move that folder into a 'deployment type' folder (e.g. .../kirk/data/Drifter_Anal/2_data/<b>SBCLS</b>)<br>
+
  <li>Now, move that folder into a 'deployment type' folder (e.g. .../kirk/data/Drifter_Anal/2_data/<b>SBCLS</b>)<br>
 
</ul>
 
</ul>
 
<br>
 
<br>
 
Notes:<br>
 
Notes:<br>
 
<ul>
 
<ul>
<li>If IRIDIUM is set, you may want to check the settings of <b>REMOVE_GPS_01</b> and <b>CONVERT_GPS_2_TO_2</b> in the <b>convert_Iridium_reports_to_Microstar_format.m</b> file.
+
  <li>If IRIDIUM is set, you may want to check the settings of <b>REMOVE_GPS_01</b> and <b>CONVERT_GPS_2_TO_2</b> in the <b>convert_Iridium_reports_to_Microstar_format.m</b> file.
<li>Data format (column numbers):
+
  <li>Data format (column numbers):
 
</ul>
 
</ul>
 
<pre>
 
<pre>
Line 67: Line 69:
  
 
<h2>Folder:  .../kirk/user/2_clean_up</h2>
 
<h2>Folder:  .../kirk/user/2_clean_up</h2>
 +
<b>Purpose: Two things: the drifters sometimes get redeployed and so tracks need to be split, also the in-water flags are sometimes wrong (e.g. due to faulty sensor or water in bucket). This step is to fix both of these problems.</b>
 +
 
<li>Edit <b>clean_up_drifter_data.m</b> with the correct <b>area</b> and <b>deploy_date</b> at the top.<br>
 
<li>Edit <b>clean_up_drifter_data.m</b> with the correct <b>area</b> and <b>deploy_date</b> at the top.<br>
 +
 +
 +
 +
 +
  
  
Line 79: Line 88:
  
 
<h2>folder: 3_process_drifter_data</h2>
 
<h2>folder: 3_process_drifter_data</h2>
<li><input type="checkbox">Edit top portion of 'process_drifter_data.m'.<br>
+
<b>Purpose: This is the step which takes the data and calculates all the derived products (e.g. dz/dt, u & v velocities, etc.) and outputs a standard formatted file of good data (optionally in interpolated round 10 minute steps).</b>
 +
 
 +
<li>Edit top portion of 'process_drifter_data.m'.<br>
 
&nbsp;Mostly just need to set one of these:<br>
 
&nbsp;Mostly just need to set one of these:<br>
 
<pre>
 
<pre>
<span class="mygeneral1-identifier">SBCIS_SBCSS </span><span class="mygeneral1-symbol">=</span><span class="mygeneral1-space"> </span><span class="mygeneral1-number">0</span><span class="mygeneral1-symbol">;</span><span class="mygeneral1-comment">%UCSB inner shelf or small-scale CODAR area work</span>
+
SBCIS_SBCSS = 0;%UCSB inner shelf or small-scale CODAR area work
<span class="mygeneral1-identifier">SBC_FULL </span><span class="mygeneral1-symbol">=</span><span class="mygeneral1-space"> </span><span class="mygeneral1-number">0</span><span class="mygeneral1-symbol">;</span><span class="mygeneral1-space"> </span><span class="mygeneral1-comment">%all of Santa Barbara Channel for Large Scale Deployments or SCR</span>
+
SBC_FULL = 0; %all of Santa Barbara Channel for Large Scale Deployments or SCR
<span class="mygeneral1-identifier">UCSD </span><span class="mygeneral1-symbol">=</span><span class="mygeneral1-space"> </span><span class="mygeneral1-number">1</span><span class="mygeneral1-symbol">;</span><span class="mygeneral1-space"> </span><span class="mygeneral1-comment">%San Diego Imperial Beach</span>
+
UCSD = 1; %San Diego Imperial Beach
<span class="mygeneral1-identifier">HBDP </span><span class="mygeneral1-symbol">=</span><span class="mygeneral1-space"> </span><span class="mygeneral1-number">0</span><span class="mygeneral1-symbol">;</span><span class="mygeneral1-space"> </span><span class="mygeneral1-comment">%Huntington Beach or Dana Point (DPIO) area</span>
+
HBDP = 0; %Huntington Beach or Dana Point (DPIO) area
<span class="mygeneral1-identifier">SAR </span><span class="mygeneral1-symbol">=</span><span class="mygeneral1-space"> </span><span class="mygeneral1-number">0</span><span class="mygeneral1-symbol">;</span><span class="mygeneral1-space"> </span><span class="mygeneral1-comment">%Santa Ana River</span>
+
SAR = 0; %Santa Ana River
  
<span class="mygeneral1-identifier">GEODAS </span><span class="mygeneral1-symbol">=</span><span class="mygeneral1-space"> </span><span class="mygeneral1-number">1</span><span class="mygeneral1-symbol">;</span><span class="mygeneral1-space"> </span><span class="mygeneral1-comment">%1=use bathy from GEODAS CD (manditory if UCSD or SAR is set above), 0=use original brian chinn sbc_bathy</span>
+
GEODAS = 1; %1=use bathy from GEODAS CD (manditory if UCSD or SAR is set above), 0=use original brian chinn sbc_bathy
 
</pre>
 
</pre>
<li><input type="checkbox">Run 'process_drifter_data.m' to process all drifters from a single deployment<br>
+
<li>Run 'process_drifter_data.m' to process all drifters from a single deployment<br>
 
Notes:<br><br>
 
Notes:<br><br>
 
<ul>
 
<ul>
<li>'process_all.m' will go thru and run 'process_drifter_data.m' on all subfolders.
+
  <li>'process_all.m' will go thru and run 'process_drifter_data.m' on all subfolders.
<li>Will generate '<i>yyyymmdd_drifter#</i>_summary.csv' files with this format:
+
  <li>Will generate '<i>yyyymmdd_drifter#</i>_summary.csv' files with this format:
 
</ul>
 
</ul>
 
<pre>
 
<pre>
Line 126: Line 137:
  
 
<h2>folder: 3b_plot_drifter_data</h2>
 
<h2>folder: 3b_plot_drifter_data</h2>
<li><input type="checkbox">Edit 'deployment_settings.m' to set max depth, dz/dt, lat/lon and velocity limits on the associated plots.
+
<b>Purpose: Generate the 4 different types of plots:  All tracks, all tracks animated, individual tracks, and summary plots.</b>
<li><input type="checkbox">Edit top portion of 'do_plots.m'.<br>
+
 
&nbsp;&nbsp;&nbsp;Mostly, to set any of these:<br>
+
<li>Edit 'deployment_settings.m' to set max depth, dz/dt, lat/lon and velocity limits on the associated plots.
<pre>
+
<li>Edit top portion of 'do_plots.m', mostly, to set any of these:<br>
<span class="mygeneral1-identifier">do_all </span><span class="mygeneral1-symbol">=</span><span class="mygeneral1-space"> </span><span class="mygeneral1-number">0</span><span class="mygeneral1-symbol">;</span><span class="mygeneral1-space">   </span><span class="mygeneral1-comment">%0 = do plot, 1 = skip</span>
+
 
<span class="mygeneral1-identifier">do_all_anim </span><span class="mygeneral1-symbol">=</span><span class="mygeneral1-space"> </span><span class="mygeneral1-number">1</span><span class="mygeneral1-symbol">;</span>
+
<pre>
<span class="mygeneral1-space">  SHOW_CODAR_TRACKS </span><span class="mygeneral1-symbol">=</span><span class="mygeneral1-space"> </span><span class="mygeneral1-number">0</span><span class="mygeneral1-symbol">;</span><span class="mygeneral1-space"> </span><span class="mygeneral1-comment">%show interpolated 'simulated drifter' tracks based on CODAR data (uses '*_compare.mat' files) for anim plot</span>
+
do_all = 0;  %0 = do plot, 1 = skip
<span class="mygeneral1-identifier">do_individ_plots </span><span class="mygeneral1-symbol">=</span><span class="mygeneral1-space"> </span><span class="mygeneral1-number">0</span><span class="mygeneral1-symbol">;</span>
+
do_all_anim = 1;
<span class="mygeneral1-identifier">do_summaries </span><span class="mygeneral1-symbol">=</span><span class="mygeneral1-space"> </span><span class="mygeneral1-number">0</span><span class="mygeneral1-symbol">;</span>
+
  SHOW_CODAR_TRACKS = 0; %show interpolated 'simulated drifter' tracks based on CODAR data (uses '*_compare.mat' files) for anim plot
<span class="mygeneral1-identifier">SHOW_CODAR </span><span class="mygeneral1-symbol">=</span><span class="mygeneral1-space"> </span><span class="mygeneral1-number">0</span><span class="mygeneral1-symbol">;</span><span class="mygeneral1-space"> </span><span class="mygeneral1-comment">%draw the hourly codar vectors on the anim plot</span>
+
do_individ_plots = 0;
 +
do_summaries = 0;
 +
SHOW_CODAR = 0; %draw the hourly codar vectors on the anim plot
 +
 
 
</pre>
 
</pre>
  
Line 142: Line 156:
  
 
<h2>folder: 4_make_web</h2>
 
<h2>folder: 4_make_web</h2>
<ol>
+
<b>Purpose: Generate a "stats" page (e.g. start/stop times, average speeds, average depth, number of drifters, etc.)</b>
 +
 
 
<li>Run 'make_stats_page.m' (generates all 'stats.php' pages for a deployment area)
 
<li>Run 'make_stats_page.m' (generates all 'stats.php' pages for a deployment area)
 
<li>Copy this deployment's folder (e.g. 20050425) to root folder of deployment type (e.g. SDSS) on the website, but
 
<li>Copy this deployment's folder (e.g. 20050425) to root folder of deployment type (e.g. SDSS) on the website, but
Line 153: Line 168:
 
<li>Syncronize E:\User\Carter with /home/data318/kirk
 
<li>Syncronize E:\User\Carter with /home/data318/kirk
 
<li>Have a beer.
 
<li>Have a beer.
</ol>
+
 
  
  
 
</ol><!-- from way near top -->
 
</ol><!-- from way near top -->

Revision as of 08:38, 17 July 2009

Kirk's HOWTO on processing Drifter data



Warning.jpg WARNING: THIS PAGE IS NOT DONE



Follow these steps to process 'raw' drifter data from pacificgyre to final posting on our website:


Notes:

  • Once proficient, it takes about 20 min to process a typical day's deployment.
  • Folders/files referenced here can be found on the ICESS network at /home/locl/users/kirk/ (hereafter referred to as '..../kirk/') in the data or user subfolders.
  • You should copy the kirk folder contents to your local computer and run the scripts there instead of on the server.


    Folder: .../kirk/user/0_Info

    Purpose: Track drifter usage history

  1. On deployment day: be sure to fill out .../kirk/user/0_Info/drifter_hist.xls (drifter history) with the performance of the drifters to track if there are any consistent problems (Note that this file is now out of date).

    Folder: .../kirk/user/1_get_data_from_web

    Purpose: Fetch the appropriate data from pacificgyre.com

  2. Edit data_dates.m with the start/stop date/times to get.
    Notes:
    • Use the http://pacificgyre.com website (log: ohlmann, pass: drifterguy) if you need to get an estimation of when the drifters went in and came out of the water (look at in-water flag values), or go to http://www.icess.ucsb.edu/drifter/realtime/BETA.php and enter custom times.
    • Keep adding date/times to the bottom of this file (or in the appropriate "deployment type" area (e.g. SBCLS), and commented out), then copy these uncommented deployment times to the top of it to fetch the data. Note, when in the "deployment type" groups, I number the 'dates' sequentially, but when copying one or two of them to the top to fetch, you need to make sure those 'dates' numbering start at 1.
    • You don't have to be exact on the times (cleanup happens in following steps), just be sure you have a time at least before or at start and a time after or at ending deployment.
  3. Run get_data.m
    • Be sure the IRIDIUM value in it is set correctly
    • Data will be saved in '..../kirk/data/Drifter_Anal/2_data/yyyymmdd' folder as a yyyymmdd_drifter_dat.log file.
    • Now, move that folder into a 'deployment type' folder (e.g. .../kirk/data/Drifter_Anal/2_data/SBCLS)


    Notes:

    • If IRIDIUM is set, you may want to check the settings of REMOVE_GPS_01 and CONVERT_GPS_2_TO_2 in the convert_Iridium_reports_to_Microstar_format.m file.
    • Data format (column numbers):
          #1 drifter ID
          #2 year
          #3 month
          #4 day
          #5 hour
          #6 minute
          #7 second
          #8 lat
          #9 long
         #10 satelites availiable
         #11 satelites used
         #12 com flag
         #13 Secs to fix
         #14 Mobitex RSSI
         #15 in-water flag
    


    Folder: .../kirk/user/2_clean_up

    Purpose: Two things: the drifters sometimes get redeployed and so tracks need to be split, also the in-water flags are sometimes wrong (e.g. due to faulty sensor or water in bucket). This step is to fix both of these problems.

  4. Edit clean_up_drifter_data.m with the correct area and deploy_date at the top.








    LEFT OFF HERE!



    folder: 3_process_drifter_data

    Purpose: This is the step which takes the data and calculates all the derived products (e.g. dz/dt, u & v velocities, etc.) and outputs a standard formatted file of good data (optionally in interpolated round 10 minute steps).

  5. Edit top portion of 'process_drifter_data.m'.
     Mostly just need to set one of these:
    SBCIS_SBCSS = 0;%UCSB inner shelf or small-scale CODAR area work
    SBC_FULL = 0; %all of Santa Barbara Channel for Large Scale Deployments or SCR
    UCSD = 1; %San Diego Imperial Beach
    HBDP = 0; %Huntington Beach or Dana Point (DPIO) area
    SAR = 0; %Santa Ana River
    
    GEODAS = 1; %1=use bathy from GEODAS CD (manditory if UCSD or SAR is set above), 0=use original brian chinn sbc_bathy
    
  6. Run 'process_drifter_data.m' to process all drifters from a single deployment
    Notes:

    • 'process_all.m' will go thru and run 'process_drifter_data.m' on all subfolders.
    • Will generate 'yyyymmdd_drifter#_summary.csv' files with this format:
             #1 drifter ID
             #2 year
             #3 month
             #4 day
             #5 hour
             #6 minute
             #7 second
             #8 lat
             #9 lon
        *   #10 u vel (cm/s)
        *   #11 v vel (cm/s)
        *   #12 along velocity (cm/s)
        *   #13 across (offshore) velocity (cm/s)
        *   #14 depth (m)
        *   #15 dz/dt (cm/s)
        *   #16 serial time
            #17 satellites available
            #18 satellites used
            #19 com flag/Receiving status
            #20 Secs to fix
            #21 Mobitex RSSI
    
        * denotes a derived product
    


    folder: 3b_plot_drifter_data

    Purpose: Generate the 4 different types of plots: All tracks, all tracks animated, individual tracks, and summary plots.

  7. Edit 'deployment_settings.m' to set max depth, dz/dt, lat/lon and velocity limits on the associated plots.
  8. Edit top portion of 'do_plots.m', mostly, to set any of these:
    do_all = 0;   %0 = do plot, 1 = skip
    do_all_anim = 1;
      SHOW_CODAR_TRACKS = 0; %show interpolated 'simulated drifter' tracks based on CODAR data (uses '*_compare.mat' files) for anim plot
    do_individ_plots = 0;
    do_summaries = 0;
    SHOW_CODAR = 0; %draw the hourly codar vectors on the anim plot
    
    


    folder: 4_make_web

    Purpose: Generate a "stats" page (e.g. start/stop times, average speeds, average depth, number of drifters, etc.)

  9. Run 'make_stats_page.m' (generates all 'stats.php' pages for a deployment area)
  10. Copy this deployment's folder (e.g. 20050425) to root folder of deployment type (e.g. SDSS) on the website, but only the *.png, stats.php & *.gif files (and 'anim' subfolder).
  11. Copy the '*_summary.log' files to the appropriate deployment type subfolder at E:\User\Carter\data\Drifter_Anal\2_data\AllData
  12. Right-click on that folder you just copied those files to, and chose 'compress to *.zip'.
  13. Move that generated *.zip folder to U:\public_html\drifter\data\data.
  14. On web, update 'mod_date.txt' with correct (today's) modified date.
  15. Check <a href="http://drifterdata.com">website</a> looks good.
  16. Syncronize E:\User\Carter with /home/data318/kirk
  17. Have a beer.