toplight

datasets/DEM/ETOPO1_Ice_g_geotiff_v3.R

+# Packages
+
+library(sf)
+library(raster)
+library(sp)
+library(tiff)
+
+prj <- "+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"
+fact <- 100
+dem <- raster("ETOPO1_Ice_g_geotiff.tif")
+dem <- aggregate(dem, fact=fact, fun=max)
+dots <-  as(dem, 'SpatialPointsDataFrame')
+dots <- st_as_sf(dots)
+dots <- dots[dots$ETOPO1_Ice_g_geotiff > 0,]
+coords <- as.data.frame(st_coordinates(dots))
+dots$X <- coords$X
+dots$Y <- coords$Y
+plot(st_geometry(dots))
+
+write.csv(dots, 'ETOPO1_light.csv')