R/data0019-dragage.R
In EffetsCumulatifsNavigation/ceanav: Research compendium for project on the assessment cumulative effects of marine vessel activities in the St. Lawrence
Defines functions
get_data0019
Documented in
get_data0019
#' Data 0019 : Sites de dragage et de dépôt
#'
#' Compilation des dragages dans l’ensemble du Saint-Laurent entre 2001 et 2016.
#'
#' @keywords dragage
#' @keywords stresseurs
#'
#' @source Environnement et Changement Climatique Canada (2020) Bilan des activités de dragage dans l’Estuaire du Saint-Laurent de 2009 à 2016. Document de travail
#' @source Ministères des Transports du Québec (2016) Compilation des dragages dans l’ensemble du Saint-Laurent entre 2001 et 2016. Document de travail
#' @source Ministères des Transports du Québec (2017) Bref portrait des activités de dragage dans la voie navigable du fleuve Saint-Laurent de 2001 à 2016. 50 + iv
#'
#' @export
#'
#' @details Cette fonction formatte les données
#'
get_data0019 <- function() {
# =~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~= #
# WARNING:
# Les données brutes sont intégrées à 2 rapports et 1 tabulateur excel qui
# se prêtent mal à une intégration programmée. Le travail de formattage et
# de préparation initial des données brute a donc été effectué manuellement
# et la base de données résultante est nommée:
# `dragage_format_DavidBeauchesne_UniversiteLaval.csv`
#
# NOTE:
# Pour plusieurs sites de dépôts et dragage, les données sur l'entretien
# de la voie maritime de la garde côtière canadienne sont utilisés (id:0018)
#
# RESOURCES:
# Degrees minutes seconds to degrees decimals converter:
# https://www.rapidtables.com/convert/number/degrees-minutes-seconds-to-degrees.html
#
# Code to visualize point, draw polygon and extract coordinates
# x <- st_point(matrix(c(-67.2848, 49.4214), nrow = 1)) %>%
# st_sfc(crs = 4326)
# x <- mapview(x) %>% mapedit::editMap()
# st_coordinates(x[[1]])
# _________________________________________________________________________ #
# =~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~= #
# Import data & general parameters
# ----------------------------------------
# Output folder
output <- "data0019-dragage/"
folder <- paste0("./data/data-raw/", output)
if (!file.exists(folder)) dir.create(folder)
# BD manuellement formatée
dat <- read.csv(paste0(folder, "dragage_format_DavidBeauchesne_UniversiteLaval.csv"))
# Rayon des buffers pour les sites de depot identifiés uniquement par un point
buf <- 100
# _________________________________________________________________________ #
# =~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~= #
# Geometries dredging and deposit sites
# ----------------------------------------
# Listes pour stocker les sites de dragage et de dépôt
dragage <- depot <- list()
# ==================================================================
# Baie-Comeau
# Havre de Baie-comeau
# ==================================================================
# Dragage
coords <- rbind(c(-68.13563, 49.22962),
c(-68.13324, 49.23028),
c(-68.13186, 49.22977),
c(-68.13203, 49.22887),
c(-68.13358, 49.22840),
c(-68.13466, 49.22839),
c(-68.13563, 49.22962))
dragage[["hbc"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "hbc")
# Depot
# Site de dépôt de la Baie-des-Anglais
# Environ 1 km du quai
# Voir rapport: https://www.bape.gouv.qc.ca/fr/dossiers/programme-decennal-dragage-abords-quais-cargill-limitee-baie-comeau/documentation/
# Coordonnées du site de dépôt: 49.256618, -68.119232
coords <- matrix(c(-68.119232, 49.256618), nrow = 1)
depot[["bda"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "bda") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# ==================================================================
# Baie-Comeau
# Quai de Cargill
# ==================================================================
# Dragage
# 49°15'17.02'' N ; 68°8'0.31'' O
coords <- rbind(c(-68.13773, 49.24890),
c(-68.13546, 49.24908),
c(-68.13556, 49.24930),
c(-68.13692, 49.24920),
c(-68.13695, 49.24942),
c(-68.13527, 49.24955),
c(-68.13490, 49.24848),
c(-68.13761, 49.24834),
c(-68.13773, 49.24890))
dragage[["carg"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "carg")
# ==================================================================
# Baie-Comeau
# Quai Société Alcoa
# ==================================================================
# Dragage
coords <- rbind(c(-68.13292, 49.25245),
c(-68.13055, 49.25468),
c(-68.13348, 49.25536),
c(-68.13400, 49.25424),
c(-68.13391, 49.25348),
c(-68.13292, 49.25245))
dragage[["alco"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "alco")
# ==================================================================
# Baie-Trinité
# Havre de Baie-Trinité
# ==================================================================
# Dragage
## Coordonnées: 49º25'17.03"N; 67º17'5.27"O
coords <- rbind(c(-67.28354, 49.42130),
c(-67.28397, 49.42101),
c(-67.28451, 49.42012),
c(-67.28539, 49.42028),
c(-67.28452, 49.42157),
c(-67.28354, 49.42130))
dragage[["hbtr"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "hbtr")
# ==================================================================
# Bécancour
# Port de Bécancour
# ==================================================================
# Dragage
# https://archives.bape.gouv.qc.ca/sections/mandats/dragage-becancour/documents/PR3.2.pdf
coords <- rbind(
c(-72.37972, 46.40020),
c(-72.37594, 46.40095),
c(-72.37198, 46.40200),
c(-72.37386, 46.40660),
c(-72.37847, 46.40531),
c(-72.38006, 46.40439),
c(-72.38007, 46.40210),
c(-72.37972, 46.40020)
)
dragage[["pbec"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "pbec")
# ==================================================================
# Berthier-sur-Mer
# Havre de Berthier-sur-Mer
# ==================================================================
# Dragage
coords <- rbind(c(-70.73491, 46.93354),
c(-70.73600, 46.93315),
c(-70.73667, 46.93304),
c(-70.73715, 46.93304),
c(-70.73759, 46.93323),
c(-70.73801, 46.93391),
c(-70.73853, 46.93546),
c(-70.73797, 46.93564),
c(-70.73742, 46.93402),
c(-70.73731, 46.93382),
c(-70.73715, 46.93371),
c(-70.73575, 46.93435),
c(-70.73491, 46.93354))
dragage[["hbm"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "hbm")
# Depot
# Site de dépôt: 46˚53'01"N / 70˚54'33"O
coords <- rbind(c(-70.72694, 46.94006),
c(-70.73460, 46.93695),
c(-70.73619, 46.93857),
c(-70.72907, 46.94159),
c(-70.72694, 46.94006))
depot[["dhbm"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dhbm")
# ==================================================================
# Cacouna
# Port de Gros-Cacouna
# ==================================================================
# Dragage
# Centroide du port (manuel) : 47.93021664372833, -69.51570540967386
coords <- rbind(c(-69.50984, 47.93378),
c(-69.51883, 47.93625),
c(-69.52223, 47.93105),
c(-69.52180, 47.92788),
c(-69.51387, 47.92444),
c(-69.50984, 47.93378))
dragage[["pgc"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "pgc")
# ==================================================================
# Cap-aux-Meules
# Havre de Cap-aux-Meules
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Cap-aux-Meules
# Quai CTMA
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Chenal de navigation du fleuve Saint-Laurent
# Bécancour à Batiscan
# NOTE: Utilisation des données de l'évaluation des effets
# environnementaux (data : 0046), qui fournit des données
# entre 1985 et 2018
# ==================================================================
# ==================================================================
# Chenal de navigation du fleuve Saint-Laurent
# Cap Santé
# NOTE: Utilisation des données de l'évaluation des effets
# environnementaux (data : 0046), qui fournit des données
# entre 1985 et 2018
# ==================================================================
# ==================================================================
# Chenal de navigation du fleuve Saint-Laurent
# Champlain à Deschambault
# NOTE: Utilisation des données de l'évaluation des effets
# environnementaux (data : 0046), qui fournit des données
# entre 1985 et 2018
# ==================================================================
# ==================================================================
# Chenal de navigation du fleuve Saint-Laurent
# Lac Saint-Pierre
# NOTE: Utilisation des données de l'évaluation des effets
# environnementaux (data : 0046), qui fournit des données
# entre 1985 et 2018
# ==================================================================
# ==================================================================
# Chenal de navigation du fleuve Saint-Laurent
# Sorel
# NOTE: Utilisation des données de l'évaluation des effets
# environnementaux (data : 0046), qui fournit des données
# entre 1985 et 2018
# ==================================================================
# ==================================================================
# Chenal de navigation du fleuve Saint-Laurent
# Traverse Nord I.O.
# NOTE: Utilisation des données de l'évaluation des effets
# environnementaux (data : 0046), qui fournit des données
# entre 1985 et 2018
# ==================================================================
# ==================================================================
# Chenal de navigation du fleuve Saint-Laurent
# Trois-Rivière à Bécancour
# NOTE: Utilisation des données de l'évaluation des effets
# environnementaux (data : 0046), qui fournit des données
# entre 1985 et 2018
# ==================================================================
# ==================================================================
# Gaspé
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Godbout
# Quai de Godbout
# ==================================================================
# Dragage
# Quai de Godbout: 49° 19'21.25" N ; 67°35'30.49" O
coords <- rbind(c(-67.59252, 49.32208),
c(-67.59093, 49.32265),
c(-67.59125, 49.32301),
c(-67.59285, 49.32243),
c(-67.59252, 49.32208))
dragage[["qgod"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "qgod")
# Depot
# Coordonnées du site de dépôt: 49.2835; -67.5590
coords <- matrix(c(-67.5590, 49.2835), nrow = 1)
depot[["dqgod"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dqgod") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# ==================================================================
# Hudson/Oka
# Chenal Hudson-Oka
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Ile-aux-Coudres
# Quai de l'Île-aux-Coudres (Desserte)
# ==================================================================
# Dragage
# Quai de l’Isle-aux-Coudres : 47°25'14.75'' N, 70°23'32.48'' O
coords <- matrix(c(-70.39236, 47.42076), nrow = 1)
dragage[["qiac"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "qiac") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# Depot
coords <- rbind(c(-70.38958, 47.4391), # 47°26,346’ N, 70°23,375’ O
c(-70.3865, 47.4391), # 47°26,346’ N, 70°23,190’ O
c(-70.3865, 47.4375), # 47°26,250’ N, 70°23,190’ O
c(-70.38958, 47.4375), # 47°26,250’ N, 70°23,375’ O
c(-70.38958, 47.4391))
depot[["dqiac"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dqiac")
# ==================================================================
# Ile-aux-Grues
# Desserte
# ==================================================================
# Dragage
# Position du quai : 47° 3'18.75"N 70°31'54.43"O
coords <- matrix(c(-70.53179, 47.05521), nrow = 1)
dragage[["qiag"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "qiag") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# Depot
coords <- rbind(c(-70.49969, 47.05239), # 47°03’ 08.6’’ N / 70°29’58,9’’ O
c(-70.49767, 47.05125), # 47°03’ 04,5’’ N / 70°29’ 51,6’’ O
c(-70.49936, 47.04986), # 47°02’ 59,5’’ N / 70°29’ 57,7’’ O
c(-70.50139, 47.05103), # 47°03’ 03,7’’ N / 70°30’ 05,0’’ O
c(-70.49969, 47.05239))
depot[["dqiag"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dqiag")
# ==================================================================
# Îles-de-la-Madeleine
# Chenal de Cap-aux-Meules
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Îles-de-la-Madeleine
# Chenal de navigation de la Grande Entrée
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Les Méchins
# Quai de Les Méchins
# ==================================================================
# Dragage
# Les Méchins 49º0'17.65"N; 66º58'27.11"O
coords <- rbind(c(-66.97288, 49.00407),
c(-66.97299, 49.00442),
c(-66.97355, 49.00467),
c(-66.97383, 49.00439),
c(-66.97323, 49.00391),
c(-66.97288, 49.00407))
dragage[["qmec"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "qmec")
# ==================================================================
# Matane
# Port de Matane
# ==================================================================
# Dragage
coords <- rbind(c(-67.57563, 48.83775),
c(-67.58122, 48.84495),
c(-67.58095, 48.84660),
c(-67.57995, 48.84761),
c(-67.57752, 48.84912),
c(-67.57005, 48.84175),
c(-67.57563, 48.83775))
dragage[["pmat"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "pmat")
# Depot
# Coordonnées du site de dépôt: 48°51’ 35’’ N; 67°35’50’’ O.
coords <- matrix(c(-67.59722, 48.85972), nrow = 1)
depot[["dpmat"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dpmat") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# ==================================================================
# Matane
# Marina de Matane
# ==================================================================
# Dragage
coords <- rbind(c(-67.52713, 48.85286),
c(-67.52950, 48.85287),
c(-67.53087, 48.85256),
c(-67.53072, 48.85207),
c(-67.52710, 48.85228),
c(-67.52713, 48.85286))
dragage[["mmat"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "mmat")
# Depot
# Même que le Port de Matane: dpmat
# ==================================================================
# Montréal
# Port de Montréal
# TODO
message("Il manque les informations pour la localisation des activités de dragage recensées pour le Port de Montréal")
# ==================================================================
# Dragage
# Depot
# Bassin Vickers
# ==================================================================
# Rimouski
# Marina de Rimouski
# ==================================================================
# Dragage
coords <- rbind(c(-68.50995, 48.47923),
c(-68.51281, 48.48108),
c(-68.51503, 48.48117),
c(-68.51602, 48.48049),
c(-68.51195, 48.47784),
c(-68.50995, 48.47923))
dragage[["mrim"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "mrim")
# Depot
coords <- rbind(c(-68.55113, 48.52129), # 68° 33' 04.06" W; 48° 31' 16.65" N
c(-68.54856, 48.51892), # 68° 32' 54.80" W; 48° 31' 08.11" N
c(-68.55212, 48.51721), # 68° 33' 07.64" W; 48° 31' 01.96" N
c(-68.5547, 48.51958), # 68° 33' 16.91" W; 48° 31' 10.49" N
c(-68.55113, 48.52129))
depot[["dmrim"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dmrim")
# ==================================================================
# Rimouski
# Port de Rimouski - Bassin est
# ==================================================================
# Dragage
# Port: https://buyandsell.gc.ca/cds/public/2018/07/19/03b09a391aa4b24241247b8a1c2cb162/ABES.PROD.PW_QCM.B008.B17442.ATTA002.PDF
coords <- rbind(c(-68.51351, 48.47823),
c(-68.51531, 48.47710),
c(-68.51735, 48.47842),
c(-68.51769, 48.47820),
c(-68.52002, 48.47978),
c(-68.51932, 48.48203),
c(-68.51868, 48.48230),
c(-68.51852, 48.48703),
c(-68.51774, 48.48697),
c(-68.51801, 48.48202),
c(-68.51751, 48.48106),
c(-68.51775, 48.48094),
c(-68.51351, 48.47823))
dragage[["prim"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "prim")
# Depot
# Même que la marina: dmrim
# ==================================================================
# Rivière-du-Loup
# Desserte
# ==================================================================
# Dragage
# Quai de traversier : 47.847816, -69.570059
coords <- rbind(c(-69.56939, 47.84786),
c(-69.57083, 47.84784),
c(-69.57080, 47.84752),
c(-69.56935, 47.84756),
c(-69.56939, 47.84786))
dragage[["qrdl"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "qrdl")
# Depot 1
coords <- rbind(c(-69.57778, 47.86222), # A) 47°51’44’’ N, 69°34’40’’ O;
c(-69.59139, 47.86667), # B) 47°52’00’’ N, 69°35’29’’ O;
c(-69.57389, 47.88944), # C) 47°53’22’’ N, 69°34’26’’ O;
c(-69.56083, 47.86222), # D) 47°51’44’’ N, 69°33’39’’ O
c(-69.57778, 47.86222))
depot[["dqrdl1"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dqrdl1")
# Depot 2
coords <- rbind(c(-69.56676, 47.87908), # A) 47°52’44.6816’’N – 69°34’00.3367 O
c(-69.56442, 47.88231), # C) 47°52’56.3321’’N – 69°33’51.9119 O
c(-69.56923, 47.88389), # D) 47°53’01.9994’’ N – 69°34’09.2316 O
c(-69.57157, 47.88065), # B) 47°52’50.3485’’N – 69°34’17.6557 O
c(-69.56676, 47.87908))
depot[["dqrdl2"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dqrdl2")
# ==================================================================
# Rivière-du-Loup
# Parc maritime de Rivière-du-Loup
# ==================================================================
# Dragage
# Marina : 47.846733, -69.569109
coords <- rbind(c(-69.56964, 47.84726),
c(-69.57028, 47.84699),
c(-69.57045, 47.84628),
c(-69.57207, 47.84622),
c(-69.57215, 47.84572),
c(-69.56764, 47.84591),
c(-69.56783, 47.84738),
c(-69.56964, 47.84726))
dragage[["mrdl"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "mrdl")
# Depot
# Même site que dépôt 1 pour la desserte de Rivière-du-Loup: qrdl1
# WARNING: Vérifier s'il s'agit réellement du bon site de dépôt
# ==================================================================
# Saint-Jean-Port-Joli
# Parc nautique de Saint-Jean-Port-Joli
# ==================================================================
# Dragage
coords <- rbind(c(-70.27263, 47.21520),
c(-70.27379, 47.21420),
c(-70.27492, 47.21478),
c(-70.27549, 47.21601),
c(-70.27537, 47.21640),
c(-70.27490, 47.21701),
c(-70.27462, 47.21684),
c(-70.27448, 47.21656),
c(-70.27445, 47.21633),
c(-70.27460, 47.21614),
c(-70.27263, 47.21520))
dragage[["sjpj"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "sjpj")
# Depot
# Site de dépôt: 600m de la jetée
# Figure et coorodnnées disponible dans un rapport d'évaluation environnementale:
# http://collections.banq.qc.ca/ark:/52327/bs47822
# "Ce site de rejet des sédiments est situé à 70° 18 de longitude ouest et
# 47° 15 de latitude nord (Figure 1)"
coords <- matrix(c(-70.365030, 47.448753), nrow = 1)
depot[["dsjpj"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dsjpj") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# ==================================================================
# Saint-Joseph-de-la-Rive
# Desserte
# ==================================================================
# Dragage
# Quai de Saint-Joseph-de-la-Rive : 47.448753, -70.365030
coords <- matrix(c(-70.365030, 47.448753), nrow = 1)
dragage[["qsjr"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "qsjr") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# Depot
# Même que pour l'îles au coudres id: dqiac
# ==================================================================
# Saint-Joseph-de-Sorel
# Zone portuaire QIT-Fer et Titane inc.
# ==================================================================
# Dragage
# https://www.environnement.gouv.qc.ca/evaluations/decret/2018/799-2018-rae.pdf
coords <- rbind(
c(-73.13460, 46.04883),
c(-73.13527, 46.04986),
c(-73.14184, 46.04783),
c(-73.14122, 46.04687),
c(-73.13460, 46.04883)
)
dragage[["sjds"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "sjds")
# ==================================================================
# Saint-Laurent-de-L'Île-D'Orléans
# Port de refuge de St-Laurent-de-L'Île-D'Orléans
# ==================================================================
# Dragage
coords <- rbind(c(-71.00592195539443, 46.85883628359575),
c(-71.00401319047153, 46.85942860450257),
c(-71.00326961538109, 46.8584064537137),
c(-71.00530869758786, 46.85779315390031),
c(-71.00592195539443, 46.85883628359575))
dragage[["cnib"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "cnib")
# Depot
# Site de dépôt: 46˚53'01"N / 70˚54'33"O
coords <- matrix(c(-71.00324129280196, 46.85683960985289), nrow = 1)
depot[["dcnib"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dcnib") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# ==================================================================
# Saint-Michel-de-Bellechasse
# Halte nautique de Saint-Michel-de-Bellechasse
# ==================================================================
# Dragage
coords <- rbind(c(-70.90834133311455, 46.87879077008841),
c(-70.9100442274937, 46.87822770301899),
c(-70.90990605635466, 46.877006174053804),
c(-70.90781046074576, 46.877893989103704),
c(-70.90834133311455, 46.87879077008841))
dragage[["smdb"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "smdb")
# Depot
# Site de dépôt: 46˚53'01"N / 70˚54'33"O
coords <- matrix(c(-70.90917, 46.88361), nrow = 1)
depot[["dsmdb"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dsmdb") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# ==================================================================
# Sept-Îles
# Bassin des remorqueurs de la Mine IOC
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Sept-Îles
# Installations portuaires de la compagnie IOC
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Sept-Îles - Quais 2
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Sorel-Tracy
# Marina de Sorel
# ==================================================================
coords <- rbind(
c(-73.09273, 46.04931),
c(-73.08984, 46.05037),
c(-73.09060, 46.05133),
c(-73.09105, 46.05121),
c(-73.09160, 46.05116),
c(-73.09173, 46.05125),
c(-73.09296, 46.05085),
c(-73.09289, 46.05073),
c(-73.09355, 46.05037),
c(-73.09273, 46.04931)
)
dragage[["msor"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "msor")
# ==================================================================
# Sorel-Tracy
# Port de Sorel-Tracy
# ==================================================================
# Dragage
coords <- rbind(
c(-73.11810, 46.04653),
c(-73.11714, 46.04821),
c(-73.11568, 46.04843),
c(-73.11545, 46.04867),
c(-73.11571, 46.04937),
c(-73.11408, 46.04948),
c(-73.11360, 46.04801),
c(-73.11585, 46.04769),
c(-73.11570, 46.04715),
c(-73.11810, 46.04653)
)
dragage[["psor"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "psor")
# ==================================================================
# Tadoussac
# Desserte est
# ==================================================================
# Dragage
coords <- rbind(c(-69.72755, 48.13854),
c(-69.72688, 48.13941),
c(-69.72614, 48.13918),
c(-69.72687, 48.13833),
c(-69.72755, 48.13854))
dragage[["detd"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "detd")
# Depot
# Site d’immersion des Bergeronnes : (48°12’00’’N ; 69°34’48’’O)
coords <- matrix(c(-69.58, 48.2), nrow = 1)
depot[["drtd"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "drtd") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# ==================================================================
# Tadoussac
# Desserte ouest
# NOTE: not used
# ==================================================================
# # Dragage
# coords <- rbind(c(-69.73103, 48.12714),
# c(-69.73011, 48.12770),
# c(-69.72878, 48.12671),
# c(-69.72959, 48.12613),
# c(-69.73103, 48.12714))
#
# dragage[["dotd"]] <- list(coords) %>%
# st_polygon() %>%
# st_sfc(crs = 4326) %>%
# st_sf(id_dragage = "dotd")
#
# # Depot
# # Même que desserte est : drtd
# ==================================================================
# Tadoussac
# Marina de Tadoussac
# ==================================================================
# Dragage
# Marina de Tadoussac (48° 8'19.60"N ; 69°43'4.69"O)
coords <- rbind(c(-69.71762, 48.13832),
c(-69.71746, 48.13921),
c(-69.71452, 48.13855),
c(-69.71497, 48.13770),
c(-69.71762, 48.13832))
dragage[["mtad"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "mtad")
# Depot
# Même que desserte est : drtd
# ==================================================================
# Valleyfield
# Port de Valleyfield
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Voie Maritime
# WARNING: Not done, outside the study area
# ==================================================================
# _________________________________________________________________________ #
# =~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~= #
# Spatial dataset
# ----------------------------------------
# Dragage
dragage <- bind_rows(dragage) %>%
left_join(dat, ., by = "id_dragage") %>%
st_as_sf() %>%
mutate(type = "dragage") %>%
filter(!is.na(id_dragage))
# Depot
depot <- bind_rows(depot) %>%
left_join(dat, ., by = "id_depot") %>%
st_as_sf() %>%
mutate(type = "depot") %>%
filter(!is.na(id_depot))
# -----
data0019 <- bind_rows(dragage, depot) %>%
st_transform(crs = global_parameters()$crs)
# _________________________________________________________________________ #
# =~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~= #
# Export data
# ----------------------------------------
# Output
st_write(obj = data0019,
dsn = "./data/data-format/data0019-dragage.geojson",
delete_dsn = TRUE)
# _________________________________________________________________________ #
}
EffetsCumulatifsNavigation/ceanav documentation built on April 17, 2023, 1:02 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Defines functions get_data0019
Documented in get_data0019
#' Data 0019 : Sites de dragage et de dépôt
#'
#' Compilation des dragages dans l’ensemble du Saint-Laurent entre 2001 et 2016.
#'
#' @keywords dragage
#' @keywords stresseurs
#'
#' @source Environnement et Changement Climatique Canada (2020) Bilan des activités de dragage dans l’Estuaire du Saint-Laurent de 2009 à 2016. Document de travail
#' @source Ministères des Transports du Québec (2016) Compilation des dragages dans l’ensemble du Saint-Laurent entre 2001 et 2016. Document de travail
#' @source Ministères des Transports du Québec (2017) Bref portrait des activités de dragage dans la voie navigable du fleuve Saint-Laurent de 2001 à 2016. 50 + iv
#'
#' @export
#'
#' @details Cette fonction formatte les données
#'
get_data0019 <- function() {
# =~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~= #
# WARNING:
# Les données brutes sont intégrées à 2 rapports et 1 tabulateur excel qui
# se prêtent mal à une intégration programmée. Le travail de formattage et
# de préparation initial des données brute a donc été effectué manuellement
# et la base de données résultante est nommée:
# `dragage_format_DavidBeauchesne_UniversiteLaval.csv`
#
# NOTE:
# Pour plusieurs sites de dépôts et dragage, les données sur l'entretien
# de la voie maritime de la garde côtière canadienne sont utilisés (id:0018)
#
# RESOURCES:
# Degrees minutes seconds to degrees decimals converter:
# https://www.rapidtables.com/convert/number/degrees-minutes-seconds-to-degrees.html
#
# Code to visualize point, draw polygon and extract coordinates
# x <- st_point(matrix(c(-67.2848, 49.4214), nrow = 1)) %>%
# st_sfc(crs = 4326)
# x <- mapview(x) %>% mapedit::editMap()
# st_coordinates(x[[1]])
# _________________________________________________________________________ #
# =~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~= #
# Import data & general parameters
# ----------------------------------------
# Output folder
output <- "data0019-dragage/"
folder <- paste0("./data/data-raw/", output)
if (!file.exists(folder)) dir.create(folder)
# BD manuellement formatée
dat <- read.csv(paste0(folder, "dragage_format_DavidBeauchesne_UniversiteLaval.csv"))
# Rayon des buffers pour les sites de depot identifiés uniquement par un point
buf <- 100
# _________________________________________________________________________ #
# =~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~= #
# Geometries dredging and deposit sites
# ----------------------------------------
# Listes pour stocker les sites de dragage et de dépôt
dragage <- depot <- list()
# ==================================================================
# Baie-Comeau
# Havre de Baie-comeau
# ==================================================================
# Dragage
coords <- rbind(c(-68.13563, 49.22962),
c(-68.13324, 49.23028),
c(-68.13186, 49.22977),
c(-68.13203, 49.22887),
c(-68.13358, 49.22840),
c(-68.13466, 49.22839),
c(-68.13563, 49.22962))
dragage[["hbc"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "hbc")
# Depot
# Site de dépôt de la Baie-des-Anglais
# Environ 1 km du quai
# Voir rapport: https://www.bape.gouv.qc.ca/fr/dossiers/programme-decennal-dragage-abords-quais-cargill-limitee-baie-comeau/documentation/
# Coordonnées du site de dépôt: 49.256618, -68.119232
coords <- matrix(c(-68.119232, 49.256618), nrow = 1)
depot[["bda"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "bda") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# ==================================================================
# Baie-Comeau
# Quai de Cargill
# ==================================================================
# Dragage
# 49°15'17.02'' N ; 68°8'0.31'' O
coords <- rbind(c(-68.13773, 49.24890),
c(-68.13546, 49.24908),
c(-68.13556, 49.24930),
c(-68.13692, 49.24920),
c(-68.13695, 49.24942),
c(-68.13527, 49.24955),
c(-68.13490, 49.24848),
c(-68.13761, 49.24834),
c(-68.13773, 49.24890))
dragage[["carg"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "carg")
# ==================================================================
# Baie-Comeau
# Quai Société Alcoa
# ==================================================================
# Dragage
coords <- rbind(c(-68.13292, 49.25245),
c(-68.13055, 49.25468),
c(-68.13348, 49.25536),
c(-68.13400, 49.25424),
c(-68.13391, 49.25348),
c(-68.13292, 49.25245))
dragage[["alco"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "alco")
# ==================================================================
# Baie-Trinité
# Havre de Baie-Trinité
# ==================================================================
# Dragage
## Coordonnées: 49º25'17.03"N; 67º17'5.27"O
coords <- rbind(c(-67.28354, 49.42130),
c(-67.28397, 49.42101),
c(-67.28451, 49.42012),
c(-67.28539, 49.42028),
c(-67.28452, 49.42157),
c(-67.28354, 49.42130))
dragage[["hbtr"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "hbtr")
# ==================================================================
# Bécancour
# Port de Bécancour
# ==================================================================
# Dragage
# https://archives.bape.gouv.qc.ca/sections/mandats/dragage-becancour/documents/PR3.2.pdf
coords <- rbind(
c(-72.37972, 46.40020),
c(-72.37594, 46.40095),
c(-72.37198, 46.40200),
c(-72.37386, 46.40660),
c(-72.37847, 46.40531),
c(-72.38006, 46.40439),
c(-72.38007, 46.40210),
c(-72.37972, 46.40020)
)
dragage[["pbec"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "pbec")
# ==================================================================
# Berthier-sur-Mer
# Havre de Berthier-sur-Mer
# ==================================================================
# Dragage
coords <- rbind(c(-70.73491, 46.93354),
c(-70.73600, 46.93315),
c(-70.73667, 46.93304),
c(-70.73715, 46.93304),
c(-70.73759, 46.93323),
c(-70.73801, 46.93391),
c(-70.73853, 46.93546),
c(-70.73797, 46.93564),
c(-70.73742, 46.93402),
c(-70.73731, 46.93382),
c(-70.73715, 46.93371),
c(-70.73575, 46.93435),
c(-70.73491, 46.93354))
dragage[["hbm"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "hbm")
# Depot
# Site de dépôt: 46˚53'01"N / 70˚54'33"O
coords <- rbind(c(-70.72694, 46.94006),
c(-70.73460, 46.93695),
c(-70.73619, 46.93857),
c(-70.72907, 46.94159),
c(-70.72694, 46.94006))
depot[["dhbm"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dhbm")
# ==================================================================
# Cacouna
# Port de Gros-Cacouna
# ==================================================================
# Dragage
# Centroide du port (manuel) : 47.93021664372833, -69.51570540967386
coords <- rbind(c(-69.50984, 47.93378),
c(-69.51883, 47.93625),
c(-69.52223, 47.93105),
c(-69.52180, 47.92788),
c(-69.51387, 47.92444),
c(-69.50984, 47.93378))
dragage[["pgc"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "pgc")
# ==================================================================
# Cap-aux-Meules
# Havre de Cap-aux-Meules
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Cap-aux-Meules
# Quai CTMA
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Chenal de navigation du fleuve Saint-Laurent
# Bécancour à Batiscan
# NOTE: Utilisation des données de l'évaluation des effets
# environnementaux (data : 0046), qui fournit des données
# entre 1985 et 2018
# ==================================================================
# ==================================================================
# Chenal de navigation du fleuve Saint-Laurent
# Cap Santé
# NOTE: Utilisation des données de l'évaluation des effets
# environnementaux (data : 0046), qui fournit des données
# entre 1985 et 2018
# ==================================================================
# ==================================================================
# Chenal de navigation du fleuve Saint-Laurent
# Champlain à Deschambault
# NOTE: Utilisation des données de l'évaluation des effets
# environnementaux (data : 0046), qui fournit des données
# entre 1985 et 2018
# ==================================================================
# ==================================================================
# Chenal de navigation du fleuve Saint-Laurent
# Lac Saint-Pierre
# NOTE: Utilisation des données de l'évaluation des effets
# environnementaux (data : 0046), qui fournit des données
# entre 1985 et 2018
# ==================================================================
# ==================================================================
# Chenal de navigation du fleuve Saint-Laurent
# Sorel
# NOTE: Utilisation des données de l'évaluation des effets
# environnementaux (data : 0046), qui fournit des données
# entre 1985 et 2018
# ==================================================================
# ==================================================================
# Chenal de navigation du fleuve Saint-Laurent
# Traverse Nord I.O.
# NOTE: Utilisation des données de l'évaluation des effets
# environnementaux (data : 0046), qui fournit des données
# entre 1985 et 2018
# ==================================================================
# ==================================================================
# Chenal de navigation du fleuve Saint-Laurent
# Trois-Rivière à Bécancour
# NOTE: Utilisation des données de l'évaluation des effets
# environnementaux (data : 0046), qui fournit des données
# entre 1985 et 2018
# ==================================================================
# ==================================================================
# Gaspé
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Godbout
# Quai de Godbout
# ==================================================================
# Dragage
# Quai de Godbout: 49° 19'21.25" N ; 67°35'30.49" O
coords <- rbind(c(-67.59252, 49.32208),
c(-67.59093, 49.32265),
c(-67.59125, 49.32301),
c(-67.59285, 49.32243),
c(-67.59252, 49.32208))
dragage[["qgod"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "qgod")
# Depot
# Coordonnées du site de dépôt: 49.2835; -67.5590
coords <- matrix(c(-67.5590, 49.2835), nrow = 1)
depot[["dqgod"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dqgod") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# ==================================================================
# Hudson/Oka
# Chenal Hudson-Oka
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Ile-aux-Coudres
# Quai de l'Île-aux-Coudres (Desserte)
# ==================================================================
# Dragage
# Quai de l’Isle-aux-Coudres : 47°25'14.75'' N, 70°23'32.48'' O
coords <- matrix(c(-70.39236, 47.42076), nrow = 1)
dragage[["qiac"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "qiac") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# Depot
coords <- rbind(c(-70.38958, 47.4391), # 47°26,346’ N, 70°23,375’ O
c(-70.3865, 47.4391), # 47°26,346’ N, 70°23,190’ O
c(-70.3865, 47.4375), # 47°26,250’ N, 70°23,190’ O
c(-70.38958, 47.4375), # 47°26,250’ N, 70°23,375’ O
c(-70.38958, 47.4391))
depot[["dqiac"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dqiac")
# ==================================================================
# Ile-aux-Grues
# Desserte
# ==================================================================
# Dragage
# Position du quai : 47° 3'18.75"N 70°31'54.43"O
coords <- matrix(c(-70.53179, 47.05521), nrow = 1)
dragage[["qiag"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "qiag") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# Depot
coords <- rbind(c(-70.49969, 47.05239), # 47°03’ 08.6’’ N / 70°29’58,9’’ O
c(-70.49767, 47.05125), # 47°03’ 04,5’’ N / 70°29’ 51,6’’ O
c(-70.49936, 47.04986), # 47°02’ 59,5’’ N / 70°29’ 57,7’’ O
c(-70.50139, 47.05103), # 47°03’ 03,7’’ N / 70°30’ 05,0’’ O
c(-70.49969, 47.05239))
depot[["dqiag"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dqiag")
# ==================================================================
# Îles-de-la-Madeleine
# Chenal de Cap-aux-Meules
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Îles-de-la-Madeleine
# Chenal de navigation de la Grande Entrée
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Les Méchins
# Quai de Les Méchins
# ==================================================================
# Dragage
# Les Méchins 49º0'17.65"N; 66º58'27.11"O
coords <- rbind(c(-66.97288, 49.00407),
c(-66.97299, 49.00442),
c(-66.97355, 49.00467),
c(-66.97383, 49.00439),
c(-66.97323, 49.00391),
c(-66.97288, 49.00407))
dragage[["qmec"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "qmec")
# ==================================================================
# Matane
# Port de Matane
# ==================================================================
# Dragage
coords <- rbind(c(-67.57563, 48.83775),
c(-67.58122, 48.84495),
c(-67.58095, 48.84660),
c(-67.57995, 48.84761),
c(-67.57752, 48.84912),
c(-67.57005, 48.84175),
c(-67.57563, 48.83775))
dragage[["pmat"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "pmat")
# Depot
# Coordonnées du site de dépôt: 48°51’ 35’’ N; 67°35’50’’ O.
coords <- matrix(c(-67.59722, 48.85972), nrow = 1)
depot[["dpmat"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dpmat") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# ==================================================================
# Matane
# Marina de Matane
# ==================================================================
# Dragage
coords <- rbind(c(-67.52713, 48.85286),
c(-67.52950, 48.85287),
c(-67.53087, 48.85256),
c(-67.53072, 48.85207),
c(-67.52710, 48.85228),
c(-67.52713, 48.85286))
dragage[["mmat"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "mmat")
# Depot
# Même que le Port de Matane: dpmat
# ==================================================================
# Montréal
# Port de Montréal
# TODO
message("Il manque les informations pour la localisation des activités de dragage recensées pour le Port de Montréal")
# ==================================================================
# Dragage
# Depot
# Bassin Vickers
# ==================================================================
# Rimouski
# Marina de Rimouski
# ==================================================================
# Dragage
coords <- rbind(c(-68.50995, 48.47923),
c(-68.51281, 48.48108),
c(-68.51503, 48.48117),
c(-68.51602, 48.48049),
c(-68.51195, 48.47784),
c(-68.50995, 48.47923))
dragage[["mrim"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "mrim")
# Depot
coords <- rbind(c(-68.55113, 48.52129), # 68° 33' 04.06" W; 48° 31' 16.65" N
c(-68.54856, 48.51892), # 68° 32' 54.80" W; 48° 31' 08.11" N
c(-68.55212, 48.51721), # 68° 33' 07.64" W; 48° 31' 01.96" N
c(-68.5547, 48.51958), # 68° 33' 16.91" W; 48° 31' 10.49" N
c(-68.55113, 48.52129))
depot[["dmrim"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dmrim")
# ==================================================================
# Rimouski
# Port de Rimouski - Bassin est
# ==================================================================
# Dragage
# Port: https://buyandsell.gc.ca/cds/public/2018/07/19/03b09a391aa4b24241247b8a1c2cb162/ABES.PROD.PW_QCM.B008.B17442.ATTA002.PDF
coords <- rbind(c(-68.51351, 48.47823),
c(-68.51531, 48.47710),
c(-68.51735, 48.47842),
c(-68.51769, 48.47820),
c(-68.52002, 48.47978),
c(-68.51932, 48.48203),
c(-68.51868, 48.48230),
c(-68.51852, 48.48703),
c(-68.51774, 48.48697),
c(-68.51801, 48.48202),
c(-68.51751, 48.48106),
c(-68.51775, 48.48094),
c(-68.51351, 48.47823))
dragage[["prim"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "prim")
# Depot
# Même que la marina: dmrim
# ==================================================================
# Rivière-du-Loup
# Desserte
# ==================================================================
# Dragage
# Quai de traversier : 47.847816, -69.570059
coords <- rbind(c(-69.56939, 47.84786),
c(-69.57083, 47.84784),
c(-69.57080, 47.84752),
c(-69.56935, 47.84756),
c(-69.56939, 47.84786))
dragage[["qrdl"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "qrdl")
# Depot 1
coords <- rbind(c(-69.57778, 47.86222), # A) 47°51’44’’ N, 69°34’40’’ O;
c(-69.59139, 47.86667), # B) 47°52’00’’ N, 69°35’29’’ O;
c(-69.57389, 47.88944), # C) 47°53’22’’ N, 69°34’26’’ O;
c(-69.56083, 47.86222), # D) 47°51’44’’ N, 69°33’39’’ O
c(-69.57778, 47.86222))
depot[["dqrdl1"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dqrdl1")
# Depot 2
coords <- rbind(c(-69.56676, 47.87908), # A) 47°52’44.6816’’N – 69°34’00.3367 O
c(-69.56442, 47.88231), # C) 47°52’56.3321’’N – 69°33’51.9119 O
c(-69.56923, 47.88389), # D) 47°53’01.9994’’ N – 69°34’09.2316 O
c(-69.57157, 47.88065), # B) 47°52’50.3485’’N – 69°34’17.6557 O
c(-69.56676, 47.87908))
depot[["dqrdl2"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dqrdl2")
# ==================================================================
# Rivière-du-Loup
# Parc maritime de Rivière-du-Loup
# ==================================================================
# Dragage
# Marina : 47.846733, -69.569109
coords <- rbind(c(-69.56964, 47.84726),
c(-69.57028, 47.84699),
c(-69.57045, 47.84628),
c(-69.57207, 47.84622),
c(-69.57215, 47.84572),
c(-69.56764, 47.84591),
c(-69.56783, 47.84738),
c(-69.56964, 47.84726))
dragage[["mrdl"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "mrdl")
# Depot
# Même site que dépôt 1 pour la desserte de Rivière-du-Loup: qrdl1
# WARNING: Vérifier s'il s'agit réellement du bon site de dépôt
# ==================================================================
# Saint-Jean-Port-Joli
# Parc nautique de Saint-Jean-Port-Joli
# ==================================================================
# Dragage
coords <- rbind(c(-70.27263, 47.21520),
c(-70.27379, 47.21420),
c(-70.27492, 47.21478),
c(-70.27549, 47.21601),
c(-70.27537, 47.21640),
c(-70.27490, 47.21701),
c(-70.27462, 47.21684),
c(-70.27448, 47.21656),
c(-70.27445, 47.21633),
c(-70.27460, 47.21614),
c(-70.27263, 47.21520))
dragage[["sjpj"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "sjpj")
# Depot
# Site de dépôt: 600m de la jetée
# Figure et coorodnnées disponible dans un rapport d'évaluation environnementale:
# http://collections.banq.qc.ca/ark:/52327/bs47822
# "Ce site de rejet des sédiments est situé à 70° 18 de longitude ouest et
# 47° 15 de latitude nord (Figure 1)"
coords <- matrix(c(-70.365030, 47.448753), nrow = 1)
depot[["dsjpj"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dsjpj") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# ==================================================================
# Saint-Joseph-de-la-Rive
# Desserte
# ==================================================================
# Dragage
# Quai de Saint-Joseph-de-la-Rive : 47.448753, -70.365030
coords <- matrix(c(-70.365030, 47.448753), nrow = 1)
dragage[["qsjr"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "qsjr") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# Depot
# Même que pour l'îles au coudres id: dqiac
# ==================================================================
# Saint-Joseph-de-Sorel
# Zone portuaire QIT-Fer et Titane inc.
# ==================================================================
# Dragage
# https://www.environnement.gouv.qc.ca/evaluations/decret/2018/799-2018-rae.pdf
coords <- rbind(
c(-73.13460, 46.04883),
c(-73.13527, 46.04986),
c(-73.14184, 46.04783),
c(-73.14122, 46.04687),
c(-73.13460, 46.04883)
)
dragage[["sjds"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "sjds")
# ==================================================================
# Saint-Laurent-de-L'Île-D'Orléans
# Port de refuge de St-Laurent-de-L'Île-D'Orléans
# ==================================================================
# Dragage
coords <- rbind(c(-71.00592195539443, 46.85883628359575),
c(-71.00401319047153, 46.85942860450257),
c(-71.00326961538109, 46.8584064537137),
c(-71.00530869758786, 46.85779315390031),
c(-71.00592195539443, 46.85883628359575))
dragage[["cnib"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "cnib")
# Depot
# Site de dépôt: 46˚53'01"N / 70˚54'33"O
coords <- matrix(c(-71.00324129280196, 46.85683960985289), nrow = 1)
depot[["dcnib"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dcnib") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# ==================================================================
# Saint-Michel-de-Bellechasse
# Halte nautique de Saint-Michel-de-Bellechasse
# ==================================================================
# Dragage
coords <- rbind(c(-70.90834133311455, 46.87879077008841),
c(-70.9100442274937, 46.87822770301899),
c(-70.90990605635466, 46.877006174053804),
c(-70.90781046074576, 46.877893989103704),
c(-70.90834133311455, 46.87879077008841))
dragage[["smdb"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "smdb")
# Depot
# Site de dépôt: 46˚53'01"N / 70˚54'33"O
coords <- matrix(c(-70.90917, 46.88361), nrow = 1)
depot[["dsmdb"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "dsmdb") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# ==================================================================
# Sept-Îles
# Bassin des remorqueurs de la Mine IOC
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Sept-Îles
# Installations portuaires de la compagnie IOC
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Sept-Îles - Quais 2
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Sorel-Tracy
# Marina de Sorel
# ==================================================================
coords <- rbind(
c(-73.09273, 46.04931),
c(-73.08984, 46.05037),
c(-73.09060, 46.05133),
c(-73.09105, 46.05121),
c(-73.09160, 46.05116),
c(-73.09173, 46.05125),
c(-73.09296, 46.05085),
c(-73.09289, 46.05073),
c(-73.09355, 46.05037),
c(-73.09273, 46.04931)
)
dragage[["msor"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "msor")
# ==================================================================
# Sorel-Tracy
# Port de Sorel-Tracy
# ==================================================================
# Dragage
coords <- rbind(
c(-73.11810, 46.04653),
c(-73.11714, 46.04821),
c(-73.11568, 46.04843),
c(-73.11545, 46.04867),
c(-73.11571, 46.04937),
c(-73.11408, 46.04948),
c(-73.11360, 46.04801),
c(-73.11585, 46.04769),
c(-73.11570, 46.04715),
c(-73.11810, 46.04653)
)
dragage[["psor"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "psor")
# ==================================================================
# Tadoussac
# Desserte est
# ==================================================================
# Dragage
coords <- rbind(c(-69.72755, 48.13854),
c(-69.72688, 48.13941),
c(-69.72614, 48.13918),
c(-69.72687, 48.13833),
c(-69.72755, 48.13854))
dragage[["detd"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "detd")
# Depot
# Site d’immersion des Bergeronnes : (48°12’00’’N ; 69°34’48’’O)
coords <- matrix(c(-69.58, 48.2), nrow = 1)
depot[["drtd"]] <- st_point(coords) %>%
st_sfc(crs = 4326) %>%
st_sf(id_depot = "drtd") %>%
st_transform(32198) %>%
st_buffer(buf) %>%
st_transform(4326)
# ==================================================================
# Tadoussac
# Desserte ouest
# NOTE: not used
# ==================================================================
# # Dragage
# coords <- rbind(c(-69.73103, 48.12714),
# c(-69.73011, 48.12770),
# c(-69.72878, 48.12671),
# c(-69.72959, 48.12613),
# c(-69.73103, 48.12714))
#
# dragage[["dotd"]] <- list(coords) %>%
# st_polygon() %>%
# st_sfc(crs = 4326) %>%
# st_sf(id_dragage = "dotd")
#
# # Depot
# # Même que desserte est : drtd
# ==================================================================
# Tadoussac
# Marina de Tadoussac
# ==================================================================
# Dragage
# Marina de Tadoussac (48° 8'19.60"N ; 69°43'4.69"O)
coords <- rbind(c(-69.71762, 48.13832),
c(-69.71746, 48.13921),
c(-69.71452, 48.13855),
c(-69.71497, 48.13770),
c(-69.71762, 48.13832))
dragage[["mtad"]] <- list(coords) %>%
st_polygon() %>%
st_sfc(crs = 4326) %>%
st_sf(id_dragage = "mtad")
# Depot
# Même que desserte est : drtd
# ==================================================================
# Valleyfield
# Port de Valleyfield
# WARNING: Not done, outside the study area
# ==================================================================
# ==================================================================
# Voie Maritime
# WARNING: Not done, outside the study area
# ==================================================================
# _________________________________________________________________________ #
# =~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~= #
# Spatial dataset
# ----------------------------------------
# Dragage
dragage <- bind_rows(dragage) %>%
left_join(dat, ., by = "id_dragage") %>%
st_as_sf() %>%
mutate(type = "dragage") %>%
filter(!is.na(id_dragage))
# Depot
depot <- bind_rows(depot) %>%
left_join(dat, ., by = "id_depot") %>%
st_as_sf() %>%
mutate(type = "depot") %>%
filter(!is.na(id_depot))
# -----
data0019 <- bind_rows(dragage, depot) %>%
st_transform(crs = global_parameters()$crs)
# _________________________________________________________________________ #
# =~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~=~-~= #
# Export data
# ----------------------------------------
# Output
st_write(obj = data0019,
dsn = "./data/data-format/data0019-dragage.geojson",
delete_dsn = TRUE)
# _________________________________________________________________________ #
}
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.