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