https://wiki-pnb.eri.ucsb.edu/index.php?action=history&feed=atom&title=Make_meta_hs6.m 2026-04-04T18:03:45Z Revision history for this page on the wiki MediaWiki 1.34.0 https://wiki-pnb.eri.ucsb.edu/index.php?title=Make_meta_hs6.m&diff=182&oldid=prev 2010-08-02T23:56:46Z

<p>Created page with &#039;&lt;pre&gt; function make_meta_hs6(cruise) %Function to automatyically construct the metafile.txt %needed for HS-6 data processing %Shpuld work on PC and UNIX %Author: Tihomir Kostadi…&#039;</p> <p><b>New page</b></p><div>&lt;pre&gt;<br /> <br /> function make_meta_hs6(cruise)<br /> %Function to automatyically construct the metafile.txt<br /> %needed for HS-6 data processing<br /> %Shpuld work on PC and UNIX<br /> %Author: Tihomir Kostadinov<br /> %Date: November 20, 2003<br /> %modified:<br /> <br /> %if isunix<br /> % eval(strcat(['cd /home/data65/pb/HYDROSCAT/', cruise])); <br /> % ctd_path = strcat(['/home/data65/pb/CTD/', cruise,'/hdr/']); <br /> % <br /> % else %PC<br /> % eval(strcat([cd ,'k:\links\PB\HYDROSCAT\', cruise]));<br /> % ctd_path = strcat(['k:\links\PB\CTD', cruise,'\hdr\']);<br /> % end<br /> %Open metafile for wrting and print header; <br /> <br /> <br /> %ctd_path = strcat(['/home/eriks/',cruise,'/CTD/',cruise,'/hdr/'])<br /> % temp ctd_path for pb215<br /> ctd_path = strcat(['/home/eriks/pb218/215_test/',cruise,'/CTD/',cruise,'/hdr/']);<br /> <br /> <br /> fid = fopen('metafile.txt', 'w'); <br /> fprintf(fid,'%s', 'File Cruise_ID Date Station Latitude Longitude'); %print header; <br /> fprintf(fid, '\n');<br /> %Extract current directory info. <br /> %k = dir ----------------------------------------original<br /> %altered above if isunix look at original when ready to insert<br /> <br /> pattern = fullfile(pwd,'*.dat')<br /> k = dir(pattern)<br /> filenames = {};<br /> for i = 1:size(k,1)<br /> <br /> filenames(i) = cellstr(k(i).name)<br /> <br /> end<br /> <br /> <br /> <br /> <br /> %Filter for .asc files only;<br /> asc_files = {}; <br /> for i = 1:length(filenames)<br /> if ~isempty(findstr(char(filenames(i)),'dat'))<br /> asc_files = [asc_files, filenames(i)];<br /> end<br /> end<br /> <br /> <br /> %Loop thropugh each h*.asc file to create a line for it in the metafile<br /> %Extract LAT/LON from the CTD .asc files (NMEA) <br /> %cd ../..<br /> pwd<br /> for i = 1:length(asc_files)<br /> line = '';<br /> file = char(asc_files(i));<br /> cdate = file(2:7);<br /> line = strcat([line,file(1:end-4),' ',cruise,' ',cdate, ' PnB', num2str(pb_let2num(file(8)))]); %used to have 's' before file<br /> <br /> %Extract lat/lon and append to line<br /> ctd_file = strcat(['X',cdate,upper(file(8)),'.hdr']);<br /> ctd_fullfile = strcat(ctd_path, ctd_file);<br /> <br /> disp(['Reading NMEA lat/lon from CTD file ', ctd_fullfile])<br /> <br /> <br /> ctdfid = fopen(ctd_fullfile,'r')<br /> ctdline = '';<br /> <br /> <br /> <br /> while isempty(findstr(ctdline,'NMEA Latitude'))<br /> ctdline = fgetl(ctdfid);<br /> end<br /> <br /> <br /> %Note that the line '* NMEA Latitude = 34 04.95 N' ansd the follwoing<br /> %line have to have equal sign '=' brigh tbeore lat/lon in degrees,<br /> %space, decimal minutes format. <br /> equal_sign = findstr(ctdline, '=');<br /> lat_raw = ctdline(equal_sign+1:end);<br /> filter = isletter(lat_raw)<br /> <br /> <br /> lat = '';<br /> for i = 1:length(lat_raw)<br /> if filter(i)==0<br /> lat = [lat, lat_raw(i)];<br /> end<br /> end<br /> <br /> <br /> lat = str2num(lat);<br /> latitude = lat(1) + lat(2)/60;<br /> <br /> <br /> <br /> %Longitude processing<br /> ctdline = fgetl(ctdfid);<br /> equal_sign = findstr(ctdline, '=');<br /> lon_raw = ctdline(equal_sign+1:end);<br /> filter = isletter(lon_raw); <br /> lon = '';<br /> for i = 1:length(lon_raw)<br /> if filter(i)==0<br /> lon = [lon, lon_raw(i)];<br /> end<br /> end<br /> lon = str2num(lon);<br /> longitude = lon(1) + lon(2)/60;<br /> fclose(ctdfid); <br /> line = strcat([line, ' ',num2str(latitude),' ',num2str(longitude)]);<br /> fprintf(fid,'%s',char(line)); <br /> fprintf(fid,'\n');<br /> atend = pwd<br /> end<br /> fclose(fid); <br /> %end<br /> &lt;/pre&gt;<br /> <br /> <br /> &lt;big&gt;[[make_meta_hs6]]&lt;big&gt;</div>

128.111.101.185