misc/SpNuit_Coord.R

library(data.table)
#library(rgdal)
#library(raster)
#library(sp)
#library(ggplot2)
#library(MASS)
#library(rgeos)
#pour afficher les milisecondes
op <- options(digits.secs=3)

#récupération des données participation
Particip=fread("C:/wamp64/www/p_export.csv")
#récupération des localités
SiteLoc=fread("C:/wamp64/www/sites_localites.txt")

args="C:/wamp64/www/SpNuit2_50_DataLP_PF_exportTot"
#args="ALL"

if(args=="ALL")
{
  SelPar=Particip  
}else{
#recupération des données chiros
DataCPL3=fread(paste0(args,".csv"))

#liste des coordonnées existantes dans ce jeu de données
ListPar=levels(as.factor(DataCPL3$participation))

SelPar=subset(Particip,Particip$participation %in% ListPar)
}

SelParSL=merge(SiteLoc,SelPar,by.x=c("site","nom"),by.y=c("site","point"))
CoordCPL3=aggregate(SelParSL$participation
                ,by=c(list(SelParSL$longitude),list(SelParSL$latitude))
                ,FUN=length)

fwrite(CoordCPL3,paste0("./VigieChiro/GIS/coordWGS84_",basename(args),".csv"))
cesco-lab/Vigie-Chiro_scripts documentation built on April 4, 2024, 4:27 a.m.