depts france

datasets/DEM/.Rhistory

-colnames(output) <- c("lat","lon","elevation")
-View(output)
-library(elevatr)
-write.csv(output, 'dem.csv')
-# 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 <- 20
-dem <- raster("alwdgg.tif")
-dem <- aggregate(dem, fact=fact, fun=mean)
-dots <-  as(dem, 'SpatialPointsDataFrame')
-dots <- st_as_sf(dots)
-dots$alwdgg[dots$alwdgg <0] <- 0
-coords <- as.data.frame(st_coordinates(dots))
-dots$X <- coords$X
-dots$Y <- coords$Y
-plot(st_geometry(dots))
-colnames(output)
-output <- dots %>% st_drop_geometry()
-output <- output[,c("Y","X","alwdgg")]
-colnames(output) <- c("lat","lon","elevation")
-View(output)
-write.csv(output, 'dem.csv')
-# 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 <- 10
-dem <- raster("alwdgg.tif")
-dem <- aggregate(dem, fact=fact, fun=mean)
-dots <-  as(dem, 'SpatialPointsDataFrame')
-dots <- st_as_sf(dots)
-dots$alwdgg[dots$alwdgg <0] <- 0
-coords <- as.data.frame(st_coordinates(dots))
-dots$X <- coords$X
-dots$Y <- coords$Y
-plot(st_geometry(dots))
-colnames(output)
-output <- dots %>% st_drop_geometry()
-output <- output[,c("Y","X","alwdgg")]
-colnames(output) <- c("lat","lon","elevation")
-View(output)
-write.csv(output, 'dem.csv')
-# 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 <- 10
-dem <- raster("alwdgg.tif")
-dem <- aggregate(dem, fact=fact, fun=mean)
 dots <-  as(dem, 'SpatialPointsDataFrame')
 dots <- st_as_sf(dots)
 dots$alwdgg[dots$alwdgg <0] <- 0
@@ -510,3 +453,60 @@ dots$Y <- coords$Y
 plot(st_geometry(dots))
 dots <- dots %>% st_drop_geometry()
 write.csv(dots, "ETOPO1_light2.csv")
+# Packages
+library(sf)
+library(raster)
+library(sp)
+library(tiff)
+fact <- 500
+dem <- raster("ETOPO1_Ice_g_geotiff.tif")
+dem1 <- dem
+dem1[dem1$ETOPO1_Ice_g_geotiff < 0] <- 0
+dem1 <- aggregate(dem1, fact=fact, fun=max)
+dots1 <-  as(dem1, 'SpatialPointsDataFrame')
+dots1 <- st_as_sf(dots1)
+dots1 <- dots1[dots1$layer >= 0,]
+dem2 <- dem
+dem2[dem2$ETOPO1_Ice_g_geotiff >= 0] <- 0
+dem2 <- aggregate(dem2, fact=fact, fun=min)
+dots2 <-  as(dem2, 'SpatialPointsDataFrame')
+dots2 <- st_as_sf(dots2)
+dots2 <- dots2[dots2$layer < 0,]
+dots <- rbind(dots1, dots2)
+plot(dots)
+# dots$alwdgg[dots$alwdgg <0] <- 0
+coords <- as.data.frame(st_coordinates(dots))
+dots$X <- coords$X
+dots$Y <- coords$Y
+output <- dots %>% st_drop_geometry()
+output <- output[,c("Y","X","layer")]
+colnames(output) <- c("lat","lon","elevation")
+write.csv(output, 'ETOPO1_minmax_fact500.csv')
+# Packages
+library(sf)
+library(raster)
+library(sp)
+library(tiff)
+fact <- 50
+dem <- raster("ETOPO1_Ice_g_geotiff.tif")
+dem1 <- dem
+dem1[dem1$ETOPO1_Ice_g_geotiff < 0] <- 0
+dem1 <- aggregate(dem1, fact=fact, fun=max)
+dots1 <-  as(dem1, 'SpatialPointsDataFrame')
+dots1 <- st_as_sf(dots1)
+dots1 <- dots1[dots1$layer >= 0,]
+dem2 <- dem
+dem2[dem2$ETOPO1_Ice_g_geotiff >= 0] <- 0
+dem2 <- aggregate(dem2, fact=fact, fun=min)
+dots2 <-  as(dem2, 'SpatialPointsDataFrame')
+dots2 <- st_as_sf(dots2)
+dots2 <- dots2[dots2$layer < 0,]
+dots <- rbind(dots1, dots2)
+# dots$alwdgg[dots$alwdgg <0] <- 0
+coords <- as.data.frame(st_coordinates(dots))
+dots$X <- coords$X
+dots$Y <- coords$Y
+output <- dots %>% st_drop_geometry()
+output <- output[,c("Y","X","layer")]
+colnames(output) <- c("lat","lon","elevation")
+write.csv(output, 'ETOPO1_minmax_fact50.csv')