variocloudmap.m

计量工具箱

function [out]=variocloudmap(long,lat,variable,robust,opt,varargin)
% PURPOSE: This function links a map and a variogram cloud (only the variogram cloud is active)
%------------------------------------------------------------------------
% USAGE: out=variocloudmap(long,lat,variable,robust,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
%           robust = parameter that tells how to calculate the variogram
%                   * robust=0: usual version
%                   * robust=1: robust version
%           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
%------------------------------------------------------------------------
% 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 variogram cloud 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, selectmap.m, selectstat.m
%------------------------------------------------------------------------
%------------------------------------------------------------------------
% Christine Thomas-Agnan, Anne Ruiz-Gazen, Julien Moutel
% June 2003
% Universit