Willy Tadema FrieseWoudloper

## tile.R
setwd("D:/temp")

# Functie om op basis van latitude, longitude en zoomniveau de bijbehorende OpenStreetMap tile te berekenen
# Zie http://wiki.openstreetmap.org/wiki/Slippy_map_tilenames

deg2num <- function(lat_deg, lon_deg, zoom) {
    lat_rad <- lat_deg * pi /180
    n <- 2.0 ^ zoom
    xtile <- floor((lon_deg + 180.0) / 360.0 * n)
    ytile = floor((1.0 - log(tan(lat_rad) + (1 / cos(lat_rad))) / pi) / 2.0 * n)

## Geschiktheid panden provincie Groningen voor opwekken zonne-energie.RMD
---
title: "Welke provinciale panden zijn geschikt voor het opwekken van zonne-energie?"
author: "Willy Bakker"
date: "5 januari 2016"
output: html_document
---

## Vraagstelling en uitgangspunten
Een collega bij de Provincie Groningen vroeg zich af: Welke panden van de provincie Groningen zijn geschikt voor het opwekken van zonne-energie? Een interessante vraag! Bovendien een uitgelezen kans om me verder te verdiepen in geografische analyses met R. In de kerstvakantie ben ik er mee aan de slag gegaan. Dit document bevat de code, documentatie en resultaten van mijn analyse.

## test.sld
<?xml version="1.0" encoding="UTF-8" ?>
 <sld:StyledLayerDescriptor version="1.0.0"
xmlns:sld="http://www.opengis.net/sld" xmlns:ogc="http://www.opengis.net/ogc"
xmlns:gml="http://www.opengis.net/gml"
xmlns:xlink="http://www.w3.org/1999/xlink"
xmlns:deegreeogc="http://www.deegree.org/ogc"
xmlns:app="http://www.deegree.org/app">
 <sld:NamedLayer>
 <sld:Name>BP:Bestemmingsplangebied</sld:Name>
 <sld:LayerFeatureConstraints>

## index.html
<!-- Voorbeeld waarin dynamische gegevens uit OpenStreetMap met behulp van Overpass API en                   -->
<!-- de osmtogeojson JavaScript library worden weergeggeven in Leaflet.                                      -->
<!-- Afgeleid van: Alex Morega, https://github.com/mgax/workshop-geocj2015-overpass/blob/master/cafes.html   -->

<!doctype html>
<meta charset="utf-8">
<link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/leaflet/0.7.3/leaflet.css">
<style>
html, body, #map { height: 100%; margin: 0; }
</style>

## test.ipynb

      
              1 file
            
          
              0 forks
            
          
              0 comments
            
          
              0 stars
            
          
                FrieseWoudloper
                / test.ipynb
            
            
              Created
              September 27, 2016 09:12
            
          
      Sorry, something went wrong. Reload?
      Sorry, we cannot display this file.
      Sorry, this file is invalid so it cannot be displayed.
      
          Viewer requires iframe.
      
    
## NGR_CSW.MD

      
              1 file
            
          
              0 forks
            
          
              0 comments
            
          
              0 stars
            
          
                FrieseWoudloper
                / NGR_CSW.MD
            
            
              Last active
              July 5, 2018 06:24
            
              
                CSW Nationaal Georegister (NGR)
              
          
    Op deze plek bewaar ik mijn experimentjes met de Catalog Service for the Web (CSW) van het Nationaal Georegister (NGR).
De documentatie van PDOK vind je hier.

Ook handig: documentatie van GeoNetwork en deze site.
Het endpoint van de service is http://nationaalgeoregister.nl/geonetwork/srv/dut/csw?

Let op: er is ook een endpoint voor INSPIRE-services. Dat is http://nationaalgeoregister.nl/geonetwork/srv/dut/inspire?
De requests in dit document zijn URL encoded. Bij sommige requests is dit namelijk een vereiste, dus heb ik het maar voor alle gedaan. Met deze [online tool](http://mey

  
## metadata_test.xml
<?xml version="1.0" encoding="UTF-8"?>
<metadata xml:lang="en">
	<dataIdInfo>
		<idAbs>Testing if I can update the metadata of an item in AOL</idAbs>
	</dataIdInfo>
</metadata>

## Parse_WFS_response.R
install.packages("XML")
library(XML)

doc <- xmlParse("http://geodata.nationaalgeoregister.nl/cbsgebiedsindelingen/wfs?request=GetCapabilities", isURL = TRUE)
output_formats <- unlist(xpathApply(doc, "//ows:Operation[@name='GetFeature']/ows:Parameter[@name='outputFormat']/ows:AllowedValues/ows:Value/text()", xmlValue))
print(output_formats)

feature_types  <- unlist(xpathApply(doc,"//wfs:FeatureType/wfs:Name/text()", xmlValue))
head(feature_types)

## verblijfsobjecten_in_zuidhorn.R
library(sf)
library(tmap)
library(magrittr)
library(httr)
library(xml2)

url <- parse_url("geodata.nationaalgeoregister.nl/bag/wfs?")
url$scheme <- "https"

# Request server limit

## AGOL.R
# Living Atlas of the world
# Rijksmonumenten
# https://www.arcgis.com/home/item.html?id=da0ee4c59ec44cff9683016c58a67470
# https://services.arcgis.com/nSZVuSZjHpEZZbRo/arcgis/rest/services/Rijksmonumenten/FeatureServer

library(sf)
library(tmap)
library(httr)
library(dplyr)
library(jsonlite)
	setwd("D:/temp")

	# Functie om op basis van latitude, longitude en zoomniveau de bijbehorende OpenStreetMap tile te berekenen
	# Zie http://wiki.openstreetmap.org/wiki/Slippy_map_tilenames

	deg2num <- function(lat_deg, lon_deg, zoom) {
	lat_rad <- lat_deg * pi /180
	n <- 2.0 ^ zoom
	xtile <- floor((lon_deg + 180.0) / 360.0 * n)
	ytile = floor((1.0 - log(tan(lat_rad) + (1 / cos(lat_rad))) / pi) / 2.0 * n)
	---
	title: "Welke provinciale panden zijn geschikt voor het opwekken van zonne-energie?"
	author: "Willy Bakker"
	date: "5 januari 2016"
	output: html_document
	---

	## Vraagstelling en uitgangspunten
	Een collega bij de Provincie Groningen vroeg zich af: Welke panden van de provincie Groningen zijn geschikt voor het opwekken van zonne-energie? Een interessante vraag! Bovendien een uitgelezen kans om me verder te verdiepen in geografische analyses met R. In de kerstvakantie ben ik er mee aan de slag gegaan. Dit document bevat de code, documentatie en resultaten van mijn analyse.
	<?xml version="1.0" encoding="UTF-8" ?>
	<sld:StyledLayerDescriptor version="1.0.0"
	xmlns:sld="http://www.opengis.net/sld" xmlns:ogc="http://www.opengis.net/ogc"
	xmlns:gml="http://www.opengis.net/gml"
	xmlns:xlink="http://www.w3.org/1999/xlink"
	xmlns:deegreeogc="http://www.deegree.org/ogc"
	xmlns:app="http://www.deegree.org/app">
	<sld:NamedLayer>
	<sld:Name>BP:Bestemmingsplangebied</sld:Name>
	<sld:LayerFeatureConstraints>
	<!-- Voorbeeld waarin dynamische gegevens uit OpenStreetMap met behulp van Overpass API en -->
	<!-- de osmtogeojson JavaScript library worden weergeggeven in Leaflet. -->
	<!-- Afgeleid van: Alex Morega, https://github.com/mgax/workshop-geocj2015-overpass/blob/master/cafes.html -->

	<!doctype html>
	<meta charset="utf-8">
	<link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/leaflet/0.7.3/leaflet.css">
	<style>
	html, body, #map { height: 100%; margin: 0; }
	</style>
	<?xml version="1.0" encoding="UTF-8"?>
	<metadata xml:lang="en">
	<dataIdInfo>
	<idAbs>Testing if I can update the metadata of an item in AOL</idAbs>
	</dataIdInfo>
	</metadata>
	install.packages("XML")
	library(XML)

	doc <- xmlParse("http://geodata.nationaalgeoregister.nl/cbsgebiedsindelingen/wfs?request=GetCapabilities", isURL = TRUE)
	output_formats <- unlist(xpathApply(doc, "//ows:Operation[@name='GetFeature']/ows:Parameter[@name='outputFormat']/ows:AllowedValues/ows:Value/text()", xmlValue))
	print(output_formats)

	feature_types <- unlist(xpathApply(doc,"//wfs:FeatureType/wfs:Name/text()", xmlValue))
	head(feature_types)
	library(sf)
	library(tmap)
	library(magrittr)
	library(httr)
	library(xml2)

	url <- parse_url("geodata.nationaalgeoregister.nl/bag/wfs?")
	url$scheme <- "https"

	# Request server limit
	# Living Atlas of the world
	# Rijksmonumenten
	# https://www.arcgis.com/home/item.html?id=da0ee4c59ec44cff9683016c58a67470
	# https://services.arcgis.com/nSZVuSZjHpEZZbRo/arcgis/rest/services/Rijksmonumenten/FeatureServer

	library(sf)
	library(tmap)
	library(httr)
	library(dplyr)
	library(jsonlite)