|[[start|Recenso PPEIC]]|[[campo|controle de campo]]|[[digita|digitalização dos dados]]|[[controle|dados PPSP-PEIC]]|

====== 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 {{solo_restinga_do_peic1_gomes_etal2007.pdf|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.

<code>
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)) 

</code>


===== Mapa Solo =====

Os dados acima produzem o seguinte mapa:

{{  mapasolopeic04nov2010.jpeg?500  |}}


==== Código Mapa ====

<code>
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")

</code>