toplight

datasets/DEM/ETOPO1_Ice_g_geotiff_v3.R

@@ -6,7 +6,7 @@ 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
+fact <- 200
 dem <- raster("ETOPO1_Ice_g_geotiff.tif")
 dem <- aggregate(dem, fact=fact, fun=max)
 dots <-  as(dem, 'SpatialPointsDataFrame')
@@ -16,5 +16,5 @@ 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')
+dots <- dots %>% st_drop_geometry()
+write.csv(dots, "ETOPO1_light2.csv")