pprmap.m

计量工具箱

function [out]=pprmap(long,lat,data,varargin)
% PURPOSE: This function links a map and two scatterplots of projections found by ppr
%------------------------------------------------------------------------
% USAGE: [out]=pprmap(long,lat,data,direct,C,half,N,M,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
%           data = n x 1 vector of the variable to study on the first axis
%           direct = optional 1 x 2 vector  of the projections to be plotted. Default is [1 2]
%           C = optional size of the starting neighborhood for each search. Default is C=tan(80*pi/180)
%           half = optional number of steps without an increase in the index, at which time the neighborhood is halved. Default is 30
%           N = optional number of structure removal, N >=2. Default is N=2
%           M = optional number of random starts. Default is M=4
%           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  
%------------------------------------------------------------------------
% MANUAL: Select points on either the map or the ppr plots 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
%         Change the active figure by pressing the 1-2 buttons
%         Selection is lost when you switch graph
%         To quit, click the middle button or press 'q'
%------------------------------------------------------------------------
% Christine Thomas-Agnan, Anne Ruiz-Gazen, Julien Moutel
% June 2003
% Universit