angleplotmap.m

计量工具箱

function [out]=angleplotmap(long,lat,variable,opt,varargin)
% PURPOSE: This function links a map and an angle plot (only the angle plot is active)
% angle plot: on the x-axis, the angle in radians between the vector joining a couple of sites
% and the x-axis; on the y-axis,the absolute difference of variable between these two sites
%------------------------------------------------------------------------
% USAGE: out=angleplotmap(long,lat,variable,opt,quantiles,carte,label,symbol)
%   where : 
%           long = n x 1 vector of coordinates on the first axis
%           lat = n x 1 vector of coordinates on the second axis
%           variable = n x 1 vector of the variable to study
%           opt = parameter that tells how to draw the scatterplot
%                   * opt=1 : the scatterplot alone is drawn
%                   * opt=2 : a kernel mean estimator is added
%                   * opt=3 : quantiles are added ( see following optionnal parameter )
%           quantiles = optionnal vector of quantile orders to draw
%           carte = n x 2 optionnal matrix giving the polygons of the edges of the spatial units
%           label = n x 1 optionnal variable used to label selected observations
%           symbol = * symbol=1 : selected spatial units are marked with a different symbol
%                    * symbol=0 : selected spatial units are marked with a different color only (default)
%------------------------------------------------------------------------
% OUTPUTS: out = (n x 1) 0-1 variable: selected spatial units are marked with a 1  
%------------------------------------------------------------------------
% NOTES: This fonction uses the function fastbinsmooth.m to calculate the mean estimator
%------------------------------------------------------------------------
% MANUAL: Select points on the angle plot by clicking with the left mouse button
%         You can select points inside a polygon: - right click to set the first point of the polygon
%                                                 - left click to set the other points
%                                                 - right click to close the polygon
%         Selection is lost when you click on the map
%         To quit, click the middle button or press 'q'
%------------------------------------------------------------------------
% uses setdens.m, fastbinsmooth.m, quant.m, selectstat.m
%------------------------------------------------------------------------
% Christine Thomas-Agnan, Anne Ruiz-Gazen, Julien Moutel
% June 2003
% Universit