dbledensitymap.m

计量工具箱

function [out]=dbledensitymap(long,lat,var1,var2,a1,a2,varargin)
% PURPOSE: This function links a map and two density estimators
%------------------------------------------------------------------------
% USAGE: out=dbledensitymap(long,lat,var1,var2,alpha1,alpha2,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
%           var1 = n x 1 vector of the first variable to study
%           var2 = n x 1 vector of the second variable to study
%           alpha1 = used to calculate the bandwidth of the first density. band1=(alpha1/100)*(max(var1)-min(var1))/2
%           alpha2 = used to calcule the bandwidth of the second density. band2=(alpha2/100)*(max(var2)-min(var2))/2
%           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 function uses the function kern_den to calculate the densities
%------------------------------------------------------------------------
% MANUAL: Select points on the map by clicking with the left mouse button
%         Select intervals on a density graph by clicking with the left mouse button
%         You can select points inside a polygon on the map: - 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 switch graph
%         To quit, click the middle button or press 'q'
% -----------------------------------------------------------------------
% uses kern_den.m, setdens3.m, selectmap.m, selectstat.m
%------------------------------------------------------------------------
% Christine Thomas-Agnan, Anne Ruiz-Gazen, Julien Moutel
% June 2003
% Universit