How can i improve marker controlled watershed segmentation algorithm for tree crown delineation using satellite image.I want neat boundaries for tree crowns.plz help me

3 次查看（过去 30 天）
madhuri kalapala
madhuri kalapala 2012-7-6
http://www.rmrkz.com/Gallery.html
clc
clear all
close all
rgb=imread('C:\Users\admin\Desktop\images\2.jpg');
imshow(rgb)
I = rgb2gray(rgb);
imshow(I)
I=adapthisteq(I);
hy = fspecial('sobel');
hx = hy';
Iy = imfilter(double(I), hy, 'replicate');
Ix = imfilter(double(I), hx, 'replicate');
gradmag = sqrt(Ix.^2 + Iy.^2);
figure, imshow(gradmag,[]), title('Gradient magnitude (gradmag)')
L = watershed(gradmag);
Lrgb = label2rgb(L);
%figure, imshow(Lrgb), title('Watershed transform of gradient magnitude (Lrgb)')
se = strel('disk',3);
Io = imopen(I, se);
figure, imshow(Io), title('Opening (Io)')
Ie = imerode(I, se);
Iobr = imreconstruct(Ie, I);
figure, imshow(Iobr), title('Opening-by-reconstruction (Iobr)')
Ioc = imclose(Io, se);
figure, imshow(Ioc), title('Opening-closing (Ioc)')
Iobrd = imdilate(Iobr, se);
Iobrcbr = imreconstruct(imcomplement(Iobrd), imcomplement(Iobr));
Iobrcbr = imcomplement(Iobrcbr);
figure, imshow(Iobrcbr), title('Opening-closing by reconstruction (Iobrcbr)')
fgm = imregionalmax(Iobrcbr);
%figure, imshow(fgm), title('Regional maxima of opening-closing by reconstruction (fgm)')
I2 = I;
I2(fgm) = 255;
figure, imshow(I2), title('Regional maxima superimposed on original image (I2)')
se2 = strel(ones(5,5));
fgm2 = imclose(fgm, se2);
fgm3 = imerode(fgm2, se2);
fgm4 = bwareaopen(fgm3, 20);
I3 = I;
I3(fgm4) = 255;
figure, imshow(I3)
title('Modified regional maxima superimposed on original image (fgm4)')
bw = im2bw(Iobrcbr, graythresh(Iobrcbr));
%figure, imshow(bw), title('Thresholded opening-closing by reconstruction (bw)')
D = bwdist(bw,'euclidean');
DL = watershed(D);
bgm = DL == 0;
%figure, imshow(bgm), title('Watershed ridge lines (bgm)')
gradmag2 = imimposemin(gradmag, bgm | fgm4);
L = watershed(gradmag2);
I4 = I;
I4(imdilate(L == 0, ones(3, 3)) | bgm | fgm4) = 255;
figure, imshow(I4)
title('Markers and object boundaries superimposed on original image (I4)')
Lrgb = label2rgb(L, 'jet', 'w', 'shuffle');
figure, imshow(Lrgb)
title('Colored watershed label matrix (Lrgb)')
figure, imshow(I), hold on
himage = imshow(Lrgb);
set(himage, 'AlphaData', 0.3);
title('Lrgb superimposed transparently on original image')
[Merged information from duplicate question]
I hav tried to implement marker controlled watershed segmentation algorithm but the trees are not delineated correctly.plz help me with code.i hav to count no.of trees in a setellite image.is there any other alg to do this
http://www.google.co.in/imgres?num=10&hl=en&biw=1366&bih=624&tbm=isch&tbnid=T-Aj9VZ5lMivpM:&imgrefurl=http://www.sciencedirect.com/science/article/pii/S0924271612000366&docid=OmnP10DcZ6jTyM&imgurl=http://ars.els-cdn.com/content/image/1-s2.0-S0924271612000366-gr13.jpg&w=579&h=481&ei=wmn6T4OMEIXRrQeK6b3cBg&zoom=1&iact=hc&vpx=877&vpy=140&dur=714&hovh=204&hovw=245&tx=123&ty=72&sig=115490177053501636257&sqi=2&page=1&tbnh=133&tbnw=160&start=0&ndsp=20&ved=1t:429,r:18,s:0,i:126
  2 个评论
Walter Roberson
Walter Roberson 2012-7-9
That google.co.in image is the Figure that the paper uses to illustrate that the paper's algorithm does not work well in some cases. Getting an algorithm that works better would be a topic for a sustained research project on your part, not something you can expect that one of the volunteers will be able to suggest in five or ten minutes thinking.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Image Analyst
Image Analyst 2012-7-7
I'd say that it is impossible for that photo. I can't even tell. It's just one big solid mass of green. You can't tell because you have tall trees and short trees and their branches overlap and penetrate neighboring trees, plus there's probably lots of underbrush in there as well. Any segmentation you did get probably wouldn't be very accurate. Why do you want the individual trees anyway?
There may be some hope if the trees are well separated or you have different species with different spectra. See the VisonBib: http://iris.usc.edu/Vision-Notes/bibliography/contentscartog.html#Cartography,%20Aerial%20Images,%20Remote%20Sensing,%20Buildings,%20Roads,%20Terrain,%20ATR
See these sections:
22.5.4 Trees, Individual Trees, Tree Crowns
22.5.4.1 Trees, Forest, Stem Volume, Biomass Measurements
22.5.4.2 Forest Extraction, Forest Analysis

类别

在 Help Center 和 File Exchange 中查找有关 Morphological Operations 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by