driftmap.m

计量工具箱

function driftmap(long,lat,variable,nl,nc,interpol,varargin)
% PURPOSE: This function adds a grid to the map and for each rectangle of the grid computes
% the mean of the spatial units included in this part.The means and Medians are also computed for each row and each column
% and plotted on the side. 
%------------------------------------------------------------------------
% USAGE: driftmap(long,lat,variable,nl,nc,interpol,theta,carte,label)
%   where : lat = n x 1 vector of coordinates on the second axis
%           long = n x 1 vector of coordinates on the first axis
%           variable = n x 1 vector of the variable to study
%           nl = number of rows of the grid.
%           mc = number of columns of the grid
%           interpol: * interpol=1 : means and medians are linearly interpolated
%                     * interpol=0 : means and medians are not interpolated
%           theta = optionnal parameter giving the angle of rotation of the map (in degree)
%           carte = n x 2 optionnal matrix giving the polygons of the edges of the spatial units
%------------------------------------------------------------------------
% uses rotation.m
%------------------------------------------------------------------------
% Christine Thomas-Agnan, Anne Ruiz-Gazen, Julien Moutel
% June 2003
% Universit