{"id":753,"date":"2024-01-06T13:23:13","date_gmt":"2024-01-06T20:23:13","guid":{"rendered":"https:\/\/blogs.ubc.ca\/yiwang28\/?p=753"},"modified":"2024-01-06T13:23:13","modified_gmt":"2024-01-06T20:23:13","slug":"my-r-learning-notes-3-how-to-calculate-the-shortest-distance-to-coast","status":"publish","type":"post","link":"https:\/\/blogs.ubc.ca\/yiwang28\/2024\/01\/06\/my-r-learning-notes-3-how-to-calculate-the-shortest-distance-to-coast\/","title":{"rendered":"My R learning notes 3: how to calculate the shortest distance to coast?"},"content":{"rendered":"

Recently, I have been working on processing data from multiple sites, which involved calculating the distance to the nearest coastline for the study sites. After exploring various methods in R and ArcGIS to compute the shortest distances, I found that many were complex and tedious. However, by integrating ideas from different sources<\/a>, I’ve developed a simplest and fastest approach for calculating the shortest distance for each site by using R. This method is capable of processing multiple sites at the same time. I’d happy to share my codes below. Please feel free to let me know if you have any questions \ud83d\ude42<\/p>\n

# Load library<\/span>
\nlibrary(ggOceanMaps) # in order to use its function dis2land, the output are in kilometers.<\/span>
\nlibrary(rnaturalearth)<\/span>
\nlibrary(sf)<\/span><\/p>\n

# Load coastline data<\/span>
\ncoast <- ne_coastline(scale=10, returnclass = “sf”)<\/span>
\nYour_sites <- read.csv(“Your_sites.csv”)<\/span><\/p>\n

# For example, I have a csv file called \"Global sites\", with three columns: sites_codes, Longitude (LON), and latitude (LAT):\r\n\"\"<\/pre>\n
# Calculate the shortest distance<\/span>
\nggOceanMaps::dist2land(Your_sites,shapefile =coast)<\/span><\/p>\n
# The results will look like:\"\"\r\nldist is the shortest distance (kilometers).<\/pre>\n
 <\/p>\n
# Additional explanations:<\/p>\n
The usage of ggOceanMaps library: We will use its function called “dis2land” to calculate the shortest distance. The function was originally designed for calculating the distance from a point in the ocean to the land (reference: R\/dist2land.R<\/a> ). But we\u00a0 can replace the shapefile of the ocean by the shapefile of the coastline, so that we can calculate distance to the closest shoreline. The output is “ldist” and it is in kilometers.<\/p>\n
The usage of rnaturalearth library: we can get the global coastline shapefile with the available finest resolution (scale =10) directly from this library.<\/p>\n
The method presents here is inspired by Mikko’s answer at https:\/\/stackoverflow.com\/questions\/51837454\/r-measuring-distance-from-a-coastline.<\/p>\n","protected":false},"excerpt":{"rendered":"
Recently, I have been working on processing data from multiple sites, which involved calculating the distance to the nearest coastline for the study sites. After exploring various methods in R and ArcGIS to compute the shortest distances, I found that many were complex and tedious. However, by integrating ideas from different sources, I’ve developed a … Continue reading My R learning notes 3: how to calculate the shortest distance to coast?<\/span> →<\/span><\/a><\/p>\n","protected":false},"author":40393,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[67],"tags":[],"class_list":["post-753","post","type-post","status-publish","format-standard","hentry","category-blogs"],"_links":{"self":[{"href":"https:\/\/blogs.ubc.ca\/yiwang28\/wp-json\/wp\/v2\/posts\/753","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/blogs.ubc.ca\/yiwang28\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/blogs.ubc.ca\/yiwang28\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/blogs.ubc.ca\/yiwang28\/wp-json\/wp\/v2\/users\/40393"}],"replies":[{"embeddable":true,"href":"https:\/\/blogs.ubc.ca\/yiwang28\/wp-json\/wp\/v2\/comments?post=753"}],"version-history":[{"count":1,"href":"https:\/\/blogs.ubc.ca\/yiwang28\/wp-json\/wp\/v2\/posts\/753\/revisions"}],"predecessor-version":[{"id":756,"href":"https:\/\/blogs.ubc.ca\/yiwang28\/wp-json\/wp\/v2\/posts\/753\/revisions\/756"}],"wp:attachment":[{"href":"https:\/\/blogs.ubc.ca\/yiwang28\/wp-json\/wp\/v2\/media?parent=753"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/blogs.ubc.ca\/yiwang28\/wp-json\/wp\/v2\/categories?post=753"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/blogs.ubc.ca\/yiwang28\/wp-json\/wp\/v2\/tags?post=753"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}