Skip to content
Snippets Groups Projects
plot_map.m 1.76 KiB
Newer Older
function plot_map(hTsgGUI, hPlotAxes)
% Function to plot the earth map and ship trackline
%
% Input
% -----
% hTsgGUI ............ Handel to the main user interface
% hPlotAxes .......... Handels to the graphic axes
%
% Output
% ------
%
% Library : 'M_MAP' 
%
% PROBLEM NOT YET RESOLVED
%    Does M_MAP use 0-360 or -180-180 covention ?
%
% $Id$

% Retrieve named application data
% -------------------------------
tsg = getappdata( hTsgGUI, 'tsg_data');

% Get the Geographic limit of the TSG time series
% -----------------------------------------------
dateLim = get(hPlotAxes(1), 'Xlim');
ind = find( tsg.DAYD >= dateLim(1) & tsg.DAYD <= dateLim(2));
latMin = min( tsg.LATX(ind) );
latMax = max( tsg.LATX(ind) );
lonMin = min( tsg.LONX(ind) );
lonMax = max( tsg.LONX(ind) );

% Positionning the right axes
% ----------------------------
axes(hPlotAxes(4));

%hLines = get( hPlotAxes(4), 'Children' );
%if ~isempty( hLines )
%  delete(hLines);
%end

% Use of Mercator projection
% --------------------------
%m_proj('Mercator','lat',[-70 70],'long',[-180 180]);
m_proj('Mercator','lat',[latMin latMax],'long',[lonMin lonMax]);

% detailed coast lines
% m_gshhs_i('patch',[.7 .7 .7]);

% Less detailed coast lines
% m_gshhs_l('patch',[.7 .7 .7]);

% Lesser  detailed coast lines
% ----------------------------
m_coast('patch',[.7 .7 .7]);

%m_grid('box','fancy','tickdir','in', ...
%       'ytick', [-70:20:70],...
%       'xtick', [-180:40:180],...
%       'Fontsize', [10]);
m_grid('box','fancy','tickdir','in', 'Fontsize', [10]);
            
% Plot the ship track - Check if LAT-LON in deg-decimal
% -----------------------------------------------------
m_line(tsg.LONX,tsg.LATX,'marker','.','markersize',3,'color','b');

end