add 2nd 25th submission

Rbanism · Dec 1, 2024 · b2d665a · b2d665a
1 parent 657cd12
commit b2d665a
Show file tree

Hide file tree

Showing 17 changed files with 111 additions and 0 deletions.
diff --git a/25Nov_Heat_2/25Nov_Heat_Y.Rproj b/25Nov_Heat_2/25Nov_Heat_Y.Rproj
@@ -0,0 +1,13 @@
+Version: 1.0
+
+RestoreWorkspace: Default
+SaveWorkspace: Default
+AlwaysSaveHistory: Default
+
+EnableCodeIndexing: Yes
+UseSpacesForTab: Yes
+NumSpacesForTab: 2
+Encoding: UTF-8
+
+RnwWeave: Sweave
+LaTeX: pdfLaTeX
diff --git a/25Nov_Heat_2/PNG/Nov25_heat.png b/25Nov_Heat_2/PNG/Nov25_heat.png
diff --git a/25Nov_Heat_2/Rcode/script.R b/25Nov_Heat_2/Rcode/script.R
@@ -0,0 +1,67 @@
+#load the packages
+library(here)
+library(tidyverse)
+library(sf)
+library(magick)
+library(grid)
+
+#load logo
+rbanism_logo <- image_read('https://rbanism.org/assets/imgs/about/vi_l.jpg') # Download our logo
+
+#import shp
+points_geothermal<-st_read("../data/Vergunningen_open_bodemenergieopslagsystemenPoint.shp")
+points_geothermal
+names(points_geothermal)
+  #vergund=licensed
+  #watervoere=water feed
+  #onderkant=bottom
+  #bovenkant=top
+  #bron=source
+
+province_utrecht<-st_read("../data/province_Utrecth.shp")
+
+summary(points_geothermal)
+
+#convert charater to numeric
+points_geothermal$VERGUND_DE <- as.numeric(points_geothermal$VERGUND_DE)
+points_geothermal$ONDERKANT_ <- as.numeric(points_geothermal$ONDERKANT_)
+points_geothermal$BOVENKANT_ <- as.numeric(points_geothermal$BOVENKANT_)
+points_geothermal$BOVENKANT0 <- as.numeric(points_geothermal$BOVENKANT0)
+points_geothermal$ONDERKANT0 <- as.numeric(points_geothermal$ONDERKANT0)
+
+
+summary(points_geothermal)
+
+#plot
+ggplot()+
+  geom_sf(
+    data=province_utrecht,
+    fill="white",
+    color="black",
+  )+
+  geom_sf(
+    data=points_geothermal,
+    aes(color=ONDERKANT0),
+    size=1,
+  )+
+  scale_color_viridis_c(option="D")+
+  coord_sf(datum = st_crs(28992))+
+  labs(
+    color="temperature",
+    title="Geothermal energy",
+    subtitle = "location of the heat storage in Utrecht",
+    caption = "Data: National Georegister, Author: Youjin Lee"
+  )
+
+#add logo
+grid.raster(rbanism_logo, x = 0.9, y=0.74,  # x and y determine the position of the logo (top right)
+            width = unit(55, "points"))   # width determines the size of the logo
+
+
+#save
+ggsave(
+  filename=here("PNG","Nov_25.png"),
+  width = 8,
+  height = 8,
+  dpi=300,)
+
diff --git a/25Nov_Heat_2/Rplot.png b/25Nov_Heat_2/Rplot.png
diff --git a/25Nov_Heat_2/carbon25.png b/25Nov_Heat_2/carbon25.png
diff --git a/25Nov_Heat_2/data/Vergunningen_open_bodemenergieopslagsystemen.zip b/25Nov_Heat_2/data/Vergunningen_open_bodemenergieopslagsystemen.zip
diff --git a/25Nov_Heat_2/data/Vergunningen_open_bodemenergieopslagsystemenPoint.cst b/25Nov_Heat_2/data/Vergunningen_open_bodemenergieopslagsystemenPoint.cst
@@ -0,0 +1 @@
+ISO-8859-1
diff --git a/25Nov_Heat_2/data/Vergunningen_open_bodemenergieopslagsystemenPoint.dbf b/25Nov_Heat_2/data/Vergunningen_open_bodemenergieopslagsystemenPoint.dbf
diff --git a/25Nov_Heat_2/data/Vergunningen_open_bodemenergieopslagsystemenPoint.prj b/25Nov_Heat_2/data/Vergunningen_open_bodemenergieopslagsystemenPoint.prj
@@ -0,0 +1 @@
+PROJCS["Amersfoort / RD New", GEOGCS["Amersfoort", DATUM["Amersfoort", SPHEROID["Bessel 1841", 6377397.155, 299.1528128, AUTHORITY["EPSG","7004"]], TOWGS84[565.2369, 50.0087, 465.658, -0.40685733032239757, -0.3507326765425626, 1.8703473836067956, 4.0812], AUTHORITY["EPSG","6289"]], PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["degree", 0.017453292519943295], AXIS["Geodetic latitude", NORTH], AXIS["Geodetic longitude", EAST], AUTHORITY["EPSG","4289"]], PROJECTION["Oblique_Stereographic", AUTHORITY["EPSG","9809"]], PARAMETER["central_meridian", 5.387638888888891], PARAMETER["latitude_of_origin", 52.15616055555556], PARAMETER["scale_factor", 0.9999079], PARAMETER["false_easting", 155000.0], PARAMETER["false_northing", 463000.0], UNIT["m", 1.0], AXIS["Easting", EAST], AXIS["Northing", NORTH], AUTHORITY["EPSG","28992"]]
diff --git a/25Nov_Heat_2/data/Vergunningen_open_bodemenergieopslagsystemenPoint.shp b/25Nov_Heat_2/data/Vergunningen_open_bodemenergieopslagsystemenPoint.shp
diff --git a/25Nov_Heat_2/data/Vergunningen_open_bodemenergieopslagsystemenPoint.shx b/25Nov_Heat_2/data/Vergunningen_open_bodemenergieopslagsystemenPoint.shx
diff --git a/25Nov_Heat_2/data/province_Utrecth.cpg b/25Nov_Heat_2/data/province_Utrecth.cpg
@@ -0,0 +1 @@
+UTF-8
diff --git a/25Nov_Heat_2/data/province_Utrecth.dbf b/25Nov_Heat_2/data/province_Utrecth.dbf
diff --git a/25Nov_Heat_2/data/province_Utrecth.prj b/25Nov_Heat_2/data/province_Utrecth.prj
@@ -0,0 +1 @@
+PROJCS["RD_New",GEOGCS["GCS_Amersfoort",DATUM["D_Amersfoort",SPHEROID["Bessel_1841",6377397.155,299.1528128]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Double_Stereographic"],PARAMETER["False_Easting",155000.0],PARAMETER["False_Northing",463000.0],PARAMETER["Central_Meridian",5.38763888888889],PARAMETER["Scale_Factor",0.9999079],PARAMETER["Latitude_Of_Origin",52.1561605555556],UNIT["Meter",1.0]]
diff --git a/25Nov_Heat_2/data/province_Utrecth.qmd b/25Nov_Heat_2/data/province_Utrecth.qmd
@@ -0,0 +1,27 @@
+<!DOCTYPE qgis PUBLIC 'http://mrcc.com/qgis.dtd' 'SYSTEM'>
+<qgis version="3.34.10-Prizren">
+  <identifier></identifier>
+  <parentidentifier></parentidentifier>
+  <language></language>
+  <type></type>
+  <title></title>
+  <abstract></abstract>
+  <links/>
+  <dates/>
+  <fees></fees>
+  <encoding></encoding>
+  <crs>
+    <spatialrefsys nativeFormat="Wkt">
+      <wkt></wkt>
+      <proj4></proj4>
+      <srsid>0</srsid>
+      <srid>0</srid>
+      <authid></authid>
+      <description></description>
+      <projectionacronym></projectionacronym>
+      <ellipsoidacronym></ellipsoidacronym>
+      <geographicflag>false</geographicflag>
+    </spatialrefsys>
+  </crs>
+  <extent/>
+</qgis>
diff --git a/25Nov_Heat_2/data/province_Utrecth.shp b/25Nov_Heat_2/data/province_Utrecth.shp
diff --git a/25Nov_Heat_2/data/province_Utrecth.shx b/25Nov_Heat_2/data/province_Utrecth.shx
Original file line number	Diff line number	Diff line change
		@@ -0,0 +1 @@
		PROJCS["Amersfoort / RD New", GEOGCS["Amersfoort", DATUM["Amersfoort", SPHEROID["Bessel 1841", 6377397.155, 299.1528128, AUTHORITY["EPSG","7004"]], TOWGS84[565.2369, 50.0087, 465.658, -0.40685733032239757, -0.3507326765425626, 1.8703473836067956, 4.0812], AUTHORITY["EPSG","6289"]], PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["degree", 0.017453292519943295], AXIS["Geodetic latitude", NORTH], AXIS["Geodetic longitude", EAST], AUTHORITY["EPSG","4289"]], PROJECTION["Oblique_Stereographic", AUTHORITY["EPSG","9809"]], PARAMETER["central_meridian", 5.387638888888891], PARAMETER["latitude_of_origin", 52.15616055555556], PARAMETER["scale_factor", 0.9999079], PARAMETER["false_easting", 155000.0], PARAMETER["false_northing", 463000.0], UNIT["m", 1.0], AXIS["Easting", EAST], AXIS["Northing", NORTH], AUTHORITY["EPSG","28992"]]
Original file line number	Diff line number	Diff line change
		@@ -0,0 +1 @@
		PROJCS["RD_New",GEOGCS["GCS_Amersfoort",DATUM["D_Amersfoort",SPHEROID["Bessel_1841",6377397.155,299.1528128]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Double_Stereographic"],PARAMETER["False_Easting",155000.0],PARAMETER["False_Northing",463000.0],PARAMETER["Central_Meridian",5.38763888888889],PARAMETER["Scale_Factor",0.9999079],PARAMETER["Latitude_Of_Origin",52.1561605555556],UNIT["Meter",1.0]]