getFich_IFN: Liste les fichiers IFN par type d'objets

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/IFN.R

Description

La fonction extrait la liste, par type d'objet des fichiers IFN de donnees brutes, presents dans le dossier 'doss'

Usage

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
)

Arguments

obj

vecteur des principaux 'objets' IFN dont on cherche a recuperer les donnees : 'Pla' (placettes), 'Arb' (arbres), 'Eco' (ecologie), 'Flore' (flore) et 'Couv' (couverts)

Peup

faut-il charger les donnees Peuplier ? Vaut FALSE par defaut

Morts

faut-il charger les donnees des arbres morts ? Vaut FALSE par defaut

Doc

faut-il charger les tables de documentation ? Vaut TRUE par defaut

ans

vecteur (numerique) des annees charger. 2005:2013 par defaut

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)

Details

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.

Value

une liste contenant les types d'objets, elle même subdivisee en data.frame par annee

Author(s)

François Morneau, Pascal Obstetar

References

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

See Also

Other IFN: getListFich_IFN()

Examples

 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)

pobsteta/gftools documentation built on March 28, 2020, 8:25 p.m.