Skip to Content

profil avec les couleurs des points sur la carte géologique

  1. # lecture de la table "profil géologie"
  2. > m2 <- read.table('profilgeol.pts', header=F, sep=" ", col.names=c('distance',
  3. + 'cat','r','g','b'))
  4.  
  5. # fusion des 2 tables sur base de la variable distance qui est la même:
  6. # le data frame résultant aura donc les paramètres distance et élévation de la table
  7. # "profil MNT" et les couleurs r,g,b de la table "profil géologie"
  8. > res <- merge(m1[,c('distance','elev')],m2[,c('distance','r','g','b')],
  9. + by="distance", all=T)
  10.  
  11. #plot
  12. > plot(res$distance,res$elev,pch=21,col = rgb(res$r/255,res$g/255,res$b/255),
  13. + bg = rgb(m2$r/255,m2$g/255,m2$b/255),xlab='Distance (m)', ylab='Elevation (m)',
  14. + main='Profil géologie')