<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META content="text/html; charset=iso-8859-1" http-equiv=Content-Type>
<META name=GENERATOR content="MSHTML 8.00.6001.23520">
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT color=#0000ff><FONT face=Arial>Tentei baixar esta biblioteca
library(EBImage) mas diz que não está disponível para versão R 3.0.1 do
Windows.</FONT> <BR></FONT><FONT size=2 face=Arial>.</FONT><BR></DIV>
<BLOCKQUOTE
style="BORDER-LEFT: #000000 2px solid; PADDING-LEFT: 5px; PADDING-RIGHT: 0px; MARGIN-LEFT: 5px; MARGIN-RIGHT: 0px">
<DIV dir=ltr>
<DIV>
<DIV>
<DIV>Caros amigos, a dúvida é a seguinte: Estou tentando por meio do pacote
EBImage contar um número de áreas de pólen na colméia. Para isso estou
utilizando uma imagem retirada da internet apenas para treino. Nesta imagem os
pontos de tons amarelos (onde há o pólen), alaranjados e brancos são de meu
interesse. No entanto acho que os locais com pólen são mais difíces por
estarem separados na área, por isso estou tentando contá-los e no caso de
sucesso, prosseguirei com a contagem das demais áreas.<BR><BR></DIV>Pois bem,
o código abaixo, conseguiu isolar razoavelmente os pontos de polén (amarelos),
no entanto ficou alguns "resíduos" o qual não consegui eliminá-los para contar
apenas o pólen.<BR><BR></DIV>Alguém sabe como fazer isso? Alguma idéia melhor
de como fazer?<BR><BR></DIV>OBS: A imagem utilizada encontra-se em anexo. Peço
desculpar por não fornecer um Link direto, pois não consegui fazer isso
utilizando nenhuma das hospedagens (box, ubuntu one ou 4 shared), a função
read.Image() sempre retorna erro: "Erro em readJPEG(x, ...) : <BR> JPEG
decompression error: Not a JPEG file: starts with 0x0a 0x0a" e não sei o que
isso significa. grato pela atenção.<BR>
<DIV>
<DIV>
<DIV>
<DIV><BR><BR>library(EBImage)<BR>imagem <- readImage("<A
href="http://www.4shared.com/download/GYeHK30-/colmeia.jpg">http://www.4shared.com/download/GYeHK30-/colmeia.jpg</A>")<BR>display(imagem)<BR>str(imagem)<BR>imagem<-imagem[1:604,60:343,1:3]<BR>display(imagem)<BR>hist(imagem)<BR>imagem2
<- imageData(channel(imagem, mode="red")) # canal vermelho destacou melhor
os pontos desejados<BR>imagem2 <-
1-imagem2
# inverte as tonalidades<BR>display(imagem2)<BR>plot(density(imagem2),
xlim=c(0,1))<BR>filled.contour(imagem2, asp=1)<BR>imagem2[imagem2>=0.2]
<- 1<BR>imagem2[imagem2<0.2] <- 0<BR>display(imagem2) # os pontos de
interesse destacaram, no entanto há muitas interferências<BR># a duúvida é
como removê-los?<BR></DIV></DIV></DIV></DIV></DIV>
<P><FONT size=2 face=Arial></FONT></P></BLOCKQUOTE></BODY></HTML>