Description Usage Arguments Details Value Author(s) References See Also Examples
La fonction extrait la liste, par type d'objet des fichiers IFN de donnees brutes, presents dans le dossier 'doss'
1 2 3 4 5 6 7 8 | getFich_IFN(
obj = c("Pla", "Arb", "Eco", "Flore", "Couv"),
Peup = FALSE,
Morts = FALSE,
Doc = TRUE,
ans = 2005:2013,
doss
)
|
obj |
vecteur des principaux 'objets' IFN dont on cherche a recuperer les donnees :
|
Peup |
faut-il charger les donnees Peuplier ? Vaut |
Morts |
faut-il charger les donnees des arbres morts ? Vaut |
Doc |
faut-il charger les tables de documentation ? Vaut |
ans |
vecteur (numerique) des annees charger. |
doss |
le nom du dossier contenant les fichiers zip de l'IFN (non modifies, avec leurs noms en 20xx-v.zip, avec '20xx' l'annee et 'v' l'eventuelle version du fichier) |
Attention, le temps d'execution de la fonction peut etre un peu long. Il est lie au temps d'acces au dossier (typiquement sous T:) et au temps de lecture, ouverture et chargement des donnees IFN. Par ailleurs les donnees d'un meme type d'objet restent separees dans le resultat produit. C'est a l'utilisateur de veiller au fait qu'il peu, ou non, les rassembler dans un seul et unique data.frame.
Attention, ce n'est pas parce qu'une variable existe plusieurs annees de suite qu'elle peut etre regroupee : les modalites peuvent avoir change ou (pire) le sens de ces modalite !!
Une seule regle : lire et relire la doc !!
Voir aussi les exemples.
une liste contenant les types d'objets, elle même subdivisee en data.frame par annee
François Morneau, Pascal Obstetar
Voir le portail internet IFN et toute la documentation associee au donnees. Ainsi que, par exemple : IGN Les donnees de l'inventaire forestier : etat des lieux et evolution, l'iF, La feuille de l'Inventaire forestier, 34, 17 p., nov. 2014
Other IFN:
getListFich_IFN()
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 | ## Not run:
# Recuperation des seules donnees 'placette' des points foret et Peupleraie,
# sans la doc et pour qq annees
Plac_IFN <- getFich_IFN(obj = 'Pla', Peup = TRUE, Doc = FALSE, ans = c(2006, 2005, 2013, 2008))
# Sous selection des seuls identifiants de points et des coordonnees
PlaPeup <- lapply(Plac_IFN$PlaPeup, function(x) subset(x, select = c(idp, xl93, yl93)))
PlaFor <- lapply(Plac_IFN$Pla, function(x) subset(x, select = c(idp, xl93, yl93)))
rm(Plac_IFN)
# Creation d'un data.frame unique
PlaPeup <- do.call('rbind', PlaPeup)
PlaFor <- do.call('rbind', PlaFor)
Pla <- rbind(data.frame(type_pla = 'Foret', PlaFor), data.frame(type_pla = 'Peup', PlaPeup))
Pla$an <- as.factor(floor(Pla$idp/ 100000) + 2005) # annee d'inventaire
xtabs(~ an , Pla)
rm(PlaPeup, PlaFor)
# Changement de projection
require(sp)
require(rgdal)
coordinates(Pla) <- ~ xl93 + yl93
# Definition de la projection utilisee
proj4string(Pla) <- lambert93
# Transformation en lambert2 etendu
Pla <- spTransform(Pla, CRS(lambert2e))
Pla <- data.frame(Pla)
names(Pla)[names(Pla) %in% c('xl93', 'yl93')] <- c('x_l2e', 'y_l2e')
## Carte des points IFN
require(ggplot2)
# Creation du graphique par ajout des couches 'france.df' et 'Pla'
p <- ggplot(Pla) +
\t\taes(x = x_l2e, y = y_l2e, colour = type_pla) +
facet_wrap(~ an) +
geom_point() +
scale_color_brewer('Type de placette', palette = 'Set1') +
coord_equal() +
\t\ttheme_bw()
p
### Exemple 3
# Recuperation des donnees 'Eco'. Ces dernieres n'ont change qu'en 2013
# avec l'ajout de 4 variables : 'obsriv' 'obsriv2' 'distriv' 'denivriv'
Eco_IFN <- getFich_IFN(obj = 'Eco', Doc = FALSE, ans = 2010:2013)$Eco
str(Eco_IFN, 1)
# tp <- do.call('rbind', Eco_IFN) # Conduit a une erreur
tp2 <- do.call('rbind', Eco_IFN[1:3]) # Fonctionne !
str(tp2)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.