MSER极值区域提取

MSER代码，实现最大稳定极值区域提取，经测试可用

function regionList = chooseRegions（ storage num_regions im ）
%
% regionList = chooseRegions（ storage num_regions ）
%
% PURPOSE:
%
% Cull out regions that are not maximally stable.  Return regions that are
% most likely to give us good matches.
%
% INPUT:
%
% im          - Grayscale image.
% num_regions - Number of regions over the entire image.
% storage     - Storage structure with an entry for each intensity level of image.
%               Fields associated with each intensity level are:
%                  storage（intensityLevel）.pixels 
%                  storage（intensityLevel）.region_labels 
%                  storage（intensityLevel）.areas 
%               “pixels“ stores a representative pixel for each region at the
%                        given intensity level. （Just the first one in that
%                        region.）
%               “region_labels“ stores the number label for each region at the given
%                        intensity level.
%               “areas“  stores the number of pixels associated with each region
%                        at the given intensity level.
%
% OUTPUT: 
%
% regionList  - N x 2 list of regions chosen from the storage structure.
%               The first column holds the region label number of a chosen
%               region and the second column holds the threshold for the region.
%

imageHeight = size（im1）;
imageWidth = size（im2）;

minLifespan   = round（0.10 * length（storage））;       % must exist over at least 10 percent of intensities
minSize       = 15;                                  % minimum area in pixels
maxSize       = round（（imageHeight*imageWidth）/40）;  % maximum area in pixels
maxDeriv      = 7;                                   % maximum change in size we‘ll consider to be flat
minStableTime = round（0.05 * length（storage））;       % must be flat over at least 5 percent of intensities
maxAxisRatio  = 10;

% Create a list of all regions over all intensity levels.
% Store the size of each region at each intensity level.

regionList = zeros（num_regions1+length（storage））;
for k=1:num_regions
    regionList（k1） = -k;  % store group number at head （storing as negative allows later ops to ignore it）
end

for k=1:length（storage） % 1 to maxIntensity+1
    fprintf（‘Processing regions for intensity: %d \n‘ k-1）;
    
    region_labels = storage（k）.region_labels;
    region_labels = region_labels（ find（region_labels>0） ）;        % ignore background group...
    for g = 1:length（region_labels） 
        regionList（region_labels（g）1+k） = storage（k）.areas（g+1）;  % ...and offset for background  
    end        
end

disp（ ‘ ‘ ）;
fprintf（‘Master region list has %d entries.\n‘ num_regions）;

% image（regionList（:2:end））;
% pause;

% Cull out regions that existed at fewer than N intensities.

lifespans = [];
for k=1:num_regions
    lifespans（k） = length（find（regionList（k:） > 0））;
end
regionList = regionList（find（lifespans > minLifespan） :）;
num_regions = size（regionList1）;

disp（ ‘ ‘ ）;
fprintf

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      97261  2012-11-19 11:38  MSER\11.jpg

     文件     119371  2012-11-19 11:38  MSER\12.jpg

     文件       5976  2007-03-05 10:44  MSER\chooseRegions.m

     文件       4207  2007-03-05 10:44  MSER\computeRegions.m

     文件       3429  2007-03-05 10:29  MSER\corner_detector.m

     文件        521  2007-03-05 10:29  MSER\displayMatchedRegions.m

     文件        380  2007-03-05 10:29  MSER\displayRegions.m

     文件       2829  2007-03-05 10:44  MSER\drawEllipse.m

     文件        763  2007-03-05 10:44  MSER\drawMatchedCorners.m

     文件       1947  2007-03-05 10:44  MSER\drawMatchedRegions.m

     文件       1113  2007-03-05 10:46  MSER\drawRegion.m

     文件      10289  2007-03-05 10:46  MSER\fit_ellipse.m

     文件        305  2007-03-05 10:20  MSER\getRegions.m

     文件        164  2007-03-05 10:46  MSER\loadimages.m

     文件        191  2007-03-05 10:46  MSER\makeStorageStructure.m

     文件       2750  2007-03-05 10:46  MSER\matchCorners.m

     文件       3811  2007-03-05 10:46  MSER\matchRegions.m

     文件       2300  2007-03-05 10:36  MSER\rectifyRegions.m

     文件        484  2007-03-05 10:34  MSER\regionMaxSize.m

     文件        255  2007-03-05 10:31  MSER\regionSlideshow.m

     文件       2678  2007-03-05 10:25  MSER\show.m

     文件        517  2007-03-05 10:25  MSER\warppatch.m

     目录          0  2012-11-19 13:46  MSER

----------- ---------  ---------- -----  ----

               261541                    23