Added to geoplotFS 'Name',Value bsb

geoplotAPP now is called with 4 input arguments
geoplotAPP(Ztable,score,ShapeFile,bsb) 

biplotAPP has been modified in order to call geoplotAPP for selected units in presence of the Shape file when check box brush units is selected

Author: MarcoRianiUNIPR

toolbox/+aux/computePCA.m

@@ -210,86 +210,9 @@
     else
     end
 
+    % Show the geoplot using the Shape file
     if ~isempty(ShapeFile)
-        geoplotAPP(Ztable,score,ShapeFile)
+         geoplotAPP(Ztable,score,ShapeFile,bsb)
     end
-
-
-    % % ShapeFile
-    % if ~isempty(ShapeFile)
-    %
-    %
-    %     delete(findobj(0, 'type', 'figure','tag','pl_geoplot'));
-    %     figure('Name','geoplot','tag','pl_geoplot')
-    %     % Specify the variable for which the colors have to be set
-    %     colorIntensity = score(:,1);
-    %     Xgeotable = readgeotable(ShapeFile);
-    %
-    %     % Define a colorbar from green to red
-    %     numColors=50;
-    %     customGradient = interp1([0;1],[0 1 0; 1 0 0],linspace(0,1,numColors));
-    %
-    %
-    %     % Check if Mapping toolbox is installed
-    %     instTOO=struct2table(ver);
-    %     if any(string(instTOO{:,1})=="Mapping Toolbox")
-    %
-    %         geoaxes
-    %         % Add to the geotable the Variable which contains the colors
-    %         Xgeotable.ColorVar=colorIntensity;
-    %         geoplot(Xgeotable,'ColorVariable','ColorVar');
-    %         colormap(customGradient);
-    %         h=colorbar;
-    %         h.Label.String='Color based on first PC';
-    %         geobasemap('topographic')
-    %
-    %         % Add labels in the centroids of the shapes
-    %         % estract XY pairs coords of enclosing poligons for each row of input table
-    %         Xtab= geotable2table(Xgeotable, ["X" "Y"]);
-    %         % Obtain the averate latitute and longitude for each row
-    %         info = shapeinfo(ShapeFile);
-    %         proj = info.CoordinateReferenceSystem;
-    %         arrlonglat=zeros(n,2);
-    %         for i = 1:n
-    %             [lat, lon] = projinv(proj, Xtab.X{i},Xtab.Y{i});  % Convert to lat/lon
-    %             arrlonglat(i,1:2)=[mean(lat,'omitmissing') mean(lon,'omitmissing')];
-    %         end
-    %
-    %         % Show the labels
-    %         text(arrlonglat(:,1), arrlonglat(:,2),rownames)
-    %
-    %     else % mapping toolbox is not installed
-    %
-    %         fc=customGradient(round(rescale(colorIntensity,1,numColors)),:);
-    %         polygon=mapshow(Xgeotable, "FaceAlpha",0.9);
-    %         for i=1:n
-    %             set(polygon.Children(i), 'FaceColor', fc(i,:),"FaceAlpha",0.9);
-    %         end
-    %         colormap(customGradient);
-    %         rang=[min(colorIntensity) max(colorIntensity)];
-    %         h=colorbar;
-    %         h.Label.String='Color based on first PC';
-    %         soglie=linspace(min(rang),max(rang),6);
-    %         valori=(soglie-min(soglie))/(max(soglie)-min(soglie));
-    %         h.Ticks=valori;
-    %         h.TickLabels=string(soglie);
-    %
-    %         % Add labels in the centroids of the shapes
-    %         arrlonglat=zeros(n,2);
-    %         % estract XY pairs coords of enclosing poligons for each row of input table
-    %         Xtab= geotable2table(Xgeotable, ["X" "Y"]);
-    %
-    %         for i=1:n
-    %             % prima creiamo degli oggetti poligonali
-    %             polObj= polyshape(Xtab.X{i,1}, ...
-    %                 Xtab.Y{i,1});
-    %             [centX,centY]=centroid(polObj);
-    %             % popoliamo i vettori
-    %             arrlonglat(i,1:2)=[centX centY];
-    %         end
-    %         text(arrlonglat(:,1), arrlonglat(:,2),rownames)
-    %
-    %
-    %     end
-    % end
+  
 end

toolbox/graphics/biplotAPP.mlapp

@@ -1,10 +1,10 @@
-function geoplotFS(Y,score,ShapeFile)
+function geoplotFS(Y,score,ShapeFile, varargin)
 %geoplotFS calls geoplotAPP.mlapp to show an interactive geoplot with polygons
 %
 %<a href="matlab: docsearchFS('geoplotFS')">Link to the help function</a>
-% 
 %
-% In the app it is also possible choose the variabile of Y to show, the
+%
+% In the app it is also possible choose the variable of Y to show, the
 % type of plot (i.e. in latitude and longitude coordinates or in the
 % original coordinates of ShapeFile), the direction of the colorbar.
 % Note that this function requires that the "mapping toolbox" is installed.
@@ -21,7 +21,7 @@ function geoplotFS(Y,score,ShapeFile)
 %               computations.
 %                Data Types - single|double
 %
-%      score  : variable to show or nxp matrix containing variables 
+%      score  : variable to show or nxp matrix containing variables
 %               to show.
 %               Numeric scalar, of string or character identifying the
 %               initial variable to use in the geoplot. Alternatively,
@@ -30,7 +30,7 @@ function geoplotFS(Y,score,ShapeFile)
 %                   Example - 'ShapeFile','shapefileName'
 %                    Data Types - char or string numeric scalar or nxp array.
 %
-%  ShapeFile  : name of ShapeFile or geotable containg shapes. 
+%  ShapeFile  : name of ShapeFile or geotable containing shapes.
 %               Character or string or geotable.
 %               Name of the ShapeFile containing the containing the
 %               geometric details of the rows. The ShapeFile, which is
@@ -40,11 +40,16 @@ function geoplotFS(Y,score,ShapeFile)
 %               "Mapping toolbox" is installed.
 %                   Example - 'ShapeFile','shapefileName'
 %                    Data Types - char or string or geotable.
-% 
+%
 %  Optional input arguments:
 %
-% Output:  
-%         
+%  bsb     :    list of units forming the initial subset. Vector. If bsb is
+%               not specified all unis are considered
+%                   Example - 'bsb',1:20
+%                    Data Types - single|double
+%
+% Output:
+%
 %
 %
 % See also: biplotFS
@@ -105,7 +110,59 @@ function geoplotFS(Y,score,ShapeFile)
     geoplotFS(X,out.score,ShapeFile)
 %}
 
+%{
+    %% Example of use of option bsb.
+    load citiesItaly2024.mat
+    X=citiesItaly2024;
+    out=pcaFS(X,'biplot',0,'dispresults',0,'plots',0);
+    ShapeFile=X.Properties.UserData{1};
+    % Color based on PCs and use of option bsb
+    % Select the provinces of Emilia-Romagna
+    bsb=[33:40 99];
+    geoplotFS(X,out.score,ShapeFile,'bsb',bsb)
+%}
+
+
 %% Beginning of code
-  geoplotAPP(Y,score,ShapeFile)
+bsb=1:size(Y,1);
+
+if nargin>3
+    options=struct('bsb',bsb);
+
+
+    [varargin{:}] = convertStringsToChars(varargin{:});
+    UserOptions=varargin(1:2:length(varargin));
+    if ~isempty(UserOptions)
+
+
+        % Check if number of supplied options is valid
+        if length(varargin) ~= 2*length(UserOptions)
+            error('FSDA:geoplotFS:WrongInputOpt','Number of supplied options is invalid. Probably values for some parameters are missing.');
+        end
+
+        % Check if all the specified optional arguments were present
+        % in structure options
+        % Remark: the nocheck option has already been dealt by routine
+        % chkinputR
+        inpchk=isfield(options,UserOptions);
+        WrongOptions=UserOptions(inpchk==0);
+        if ~isempty(WrongOptions)
+            disp(strcat('Non existent user option found->', char(WrongOptions{:})))
+            error('FSDA:geoplotFS:NonExistInputOpt','In total %d non-existent user options found.', length(WrongOptions));
+        end
+    end
+
+
+    % Write in structure 'options' the options chosen by the user
+    for i=1:2:length(varargin)
+        options.(varargin{i})=varargin{i+1};
+    end
+
+    bsb=options.bsb;
+    
+end
+
+geoplotAPP(Y,score,ShapeFile,bsb)
+
 end
 %FScategory:VIS-Mult