Classificação de Solos

Esses dados foram repassados pelo Renato Lima ao Ale para as análises da Mari. Foram reajustados segundo mapa disponível em Gomes et. al 2007. As análises da Márcia foram feitas com a versão antiga. Depois de pequenos ajustes a versão atual é:

• solo.class04nov2010.txt

Essa versão está bem parecida com o mapa de Gomes. Entretanto, nossos códigos tem algumas diferenças que devem ser ajustadas para publicar (caso queiram adotar a classificação!):

• esp.hist = Organossolo + Espodossolo organossólico
• esp.org = Espodossolo organossólico com horizonte (Cgi)

No artigo ele comenta que são muito parecidos: os Organossolos como os espdossolos organossólicos (que para ele inclui o hístico). Todos estão em terrenos mais baixos (“abaciados” no organossolo) e tem forte influência da flutuação do lençol. Parece que são os solos encharcados e alagados da parcela. Por essa razão estou classificando como sendo o mesmo solo.

solo=read.table("solo.class04nov2010.txt",header=T, sep="\t", as.is=T)
str(solo)
unique(solo$soil)
## ajustes segundo mapa de Gomes (2007): na classificação de gomes: esp.hist = Organossolo + Esp organossolico e esp.org = esp.org com horizonte (Cgi)
## portanto vou juntar esp.hist + esp.org como um único tipo
## minha classificação:
solo$tipo[solo$soil=="neo"]=1
solo$tipo[solo$soil=="trans"]=2
solo$tipo[solo$soil=="esp.aren"]=3
solo$tipo[solo$soil=="esp.hist"]=4
solo$tipo[solo$soil=="esp.org"]=4
solo$tipo
# mapa
matrix(solo$tipo,16,16,byrow=F)-> matsolo
image(t(matsolo)) 

Os dados acima produzem o seguinte mapa:

matsolo=matrix(solo$tipo,16, 16, byrow=F)

image(0:16,0:16 ,t(matsolo), xaxp=c(0,16,16), yaxt="n", xaxt="n", xlab="", ylab="")

axis(1, 0.5:15.5, labels=LETTERS[1:16], tick=FALSE)

axis(1, 0:16, labels=FALSE, tick=TRUE)

axis(2, 0:16, labels=FALSE, tick=TRUE)

axis(2, seq(0.5,15.5, by=1), labels=seq(1,16,by=1), las=1, tick=FALSE)

text(locator(1), label="neossolo")

text(locator(1), label="espodosol espessarênico + diúrico")

text(locator(1), label="transição")

text(locator(1), label="espodosol organossólico + Organossolo")