32:T41d,Die Daten für diese Analyse stammen aus verschiedenen Quellen: Die Schläge wurden vom [Thünen-Institut für Fernerkundung](https://www.thuenen.de/de/thuenen-institut/verbundstrukturen/thuenen-fernerkundung) erfasst und mit den Daten aus dem Integrierten Verwaltungs- und Kontrollsystem (Invekos) für Nordrhein-Westfalen, Niedersachsen und Brandenburg abgeglichen; die Autobahnen wurden aus [OpenStreetMap](https://www.openstreetmap.org/) extrahiert; das Schienennetz wurde von der [Deutschen Bahn](https://data.deutschebahn.com/dataset/geo-strecke.html) bereitgestellt; die Naturschutzgebiete wurden aus dem European Environment Agency Dataset "[Nationally designated areas](https://www.eea.europa.eu/data-and-maps/data/nationally-designated-areas-national-cdda-17)" übernommen; die Strahlungsintensitäten wurden dem [Global Solar Atlas](https://globalsolaratlas.info/) entnommen; und die Gemeinden wurden vom [Bundesamt für Kartographie und Geodäsie](https://gdz.bkg.bund.de/index.php/default/digitale-geodaten/verwaltungsgebiete.html) bezogen.33:T30c7,freiflaechen_autobahnen_bahnschienen_karte_vTI = { const container = yield html`

`; // add the PMTiles plugin to the maplibregl global. let protocol = new pmtiles.Protocol(); maplibregl.addProtocol("pmtiles", protocol.tile); const map = (container.value = new maplibregl.Map({ //interactive: false, boxZoom: true, pitch: 0, bearing: 0, maplibreLogo: false, customAttribution: "Pahmeyer et. al (2023), Thünen-Institut", container, center: [10.3853, 51.5431], zoom: 4.5, style: await FileAttachment("positron@1.json").json(), cooperativeGestures: { windowsHelpText: 'Halte die Taste "strg" beim Scrollen gedrückt, um die Karte zu vergrößern', macHelpText: 'Halte die Taste "⌘" beim Scrollen gedrückt, um die Karte zu vergrößern', mobileHelpText: "Use two fingers to move the map" } // scrollZoom: false })); var geocoder_api = { forwardGeocode: async (config) => { const features = []; try { let request = "https://nominatim.openstreetmap.org/search?q=" + config.query + "&format=geojson&polygon_geojson=1&addressdetails=1&countrycodes=de"; const response = await fetch(request); const geojson = await response.json(); for (let feature of geojson.features) { let center = [ feature.bbox[0] + (feature.bbox[2] - feature.bbox[0]) / 2, feature.bbox[1] + (feature.bbox[3] - feature.bbox[1]) / 2 ]; let point = { type: "Feature", geometry: { type: "Point", coordinates: center }, place_name: feature.properties.display_name, properties: feature.properties, text: feature.properties.display_name, place_type: ["place"], center: center }; features.push(point); } } catch (e) { console.error(`Fehler beim Suchen nach Adressdaten: ${e}`); } return { features: features }; } }; map.addControl( new MaplibreGeocoder(geocoder_api, { maplibregl: maplibregl, showResultsWhileTyping: true, debounceSearch: 300 }) ); map.addControl(new maplibregl.NavigationControl()); /* const legend = new LegendControl({ // Show all properties in selected layers layers: [ "Naturschutz", "Autobahnen", "Bahnstrecken", "Landwirtschaftliche_Flächen" ] }); map.addControl(legend, "top-left"); */ map.on("load", async function () { const layers = map.getStyle().layers; let firstSymbolId; for (let i = 0; i < layers.length; i++) { if (layers[i].type === "symbol") { firstSymbolId = layers[i].id; break; } } map.addSource("wms-test-source", { type: "raster", // use the tiles option to specify a WMS tile source URL // https://maplibre.org/maplibre-gl-js-docs/style-spec/sources/ tiles: [ "https://sgx.geodatenzentrum.de/wms_sentinel2_de?bbox={bbox-epsg-3857}&format=image/png&service=WMS&version=1.1.1&request=GetMap&srs=EPSG:3857&transparent=true&width=256&height=256&layers=rgb" ], tileSize: 256 }); viewof hintergrund.addEventListener("input", () => { evalHintergrund(); }); function evalHintergrund() { const selectedHintergrund = viewof hintergrund.value; if (selectedHintergrund === "Satellit") { if (map.getLayer("wms-satellit-layer")) return; map.addLayer( { id: "wms-satellit-layer", type: "raster", source: "wms-test-source", paint: {} }, "Bahnstrecken" ); } else { if (map?.getLayer("wms-satellit-layer")) { map.removeLayer("wms-satellit-layer"); } } } map.addSource("data", { type: "vector", maxzoom: 10, tiles: [ `pmtiles://https://humusklimanetz-couch.thuenen.de/datasets/freiflächen-pv-autobahnen-bahnstrecken_500.pmtiles/{z}/{x}/{y}` ] }); map.addSource("nuts3", { type: "geojson", data: geojson_nuts3 }); map.addLayer( { id: "Autobahnen", type: "fill", source: "data", "source-layer": "autobahnenbahnstrecken_buffered_clipped_500dissolved", layout: {}, paint: { "fill-color": "#37464A", "fill-opacity": 1 }, filter: ["==", "type", "Autobahn"] }, firstSymbolId ); map.addLayer( { id: "Bahnstrecken", type: "fill", source: "data", "source-layer": "autobahnenbahnstrecken_buffered_clipped_500dissolved", layout: {}, paint: { "fill-color": "#DF7C28", "fill-opacity": 1 }, filter: ["!=", "type", "Autobahn"] }, firstSymbolId ); // Satelliten layer anzeigen evalHintergrund(); // load CanvasPattern const image = await asyncLoadImage(naturschutz_pattern_png); map.addImage("natura2000_pattern", image); map.addLayer( { id: "Naturschutz", type: "fill", source: "data", "source-layer": "NaturaSite_polygon_renamedndjson", layout: {}, paint: { "fill-pattern": "natura2000_pattern", //"fill-color": "#7DBB39", "fill-opacity": [ "interpolate", ["linear"], ["zoom"], 0, 0, 5, 0, 6, 0.5 ] } }, firstSymbolId ); map.addLayer( { id: "naturschutz-label", type: "symbol", source: "data", "source-layer": "NaturaSite_polygon_renamedndjson", layout: { "text-field": [ "format", ["get", "natura2000"], { "text-font": ["literal", ["DIN Offc Pro Italic"]], "font-scale": 0.8 } ], "text-font": ["DIN Offc Pro Medium", "Arial Unicode MS Bold"], "text-size": 12 }, paint: { "text-opacity": [ "interpolate", ["linear"], ["zoom"], 0, 0, 9.5, 0, 10.5, 1 ] } }, firstSymbolId ); map.addLayer( { id: "Landwirtschaftliche_Flächen", type: "fill", minzoom: 10, source: "data", "source-layer": "lf_an_autobahnen_bahnstrecken_inkl_daten_500ndjson", layout: {}, paint: { "fill-color": "#0B95C7", "fill-opacity": [ "interpolate", ["linear"], ["zoom"], 0, 0, 10, 0, 10.5, 1 ] }, ...(nat_ausschliessen.length ? { filter: [ "all", ["==", ["typeof", ["get", "natura2000"]], "null"], [">", ["get", "streifen_flaeche_ha"], min_size] ] } : { filter: [">", ["get", "streifen_flaeche_ha"], min_size] }) }, firstSymbolId ); map.addLayer( { id: "flaechen-line", type: "line", source: "data", "source-layer": "lf_an_autobahnen_bahnstrecken_inkl_daten_500ndjson", minzoom: 10, paint: { "line-color": "white", "line-opacity": [ "interpolate", ["linear"], ["zoom"], 0, 0, 10, 0, 10.5, 1 ], "line-width": 2 }, ...(nat_ausschliessen.length ? { filter: [ "all", ["==", ["typeof", ["get", "natura2000"]], "null"], [">", ["get", "streifen_flaeche_ha"], min_size] ] } : { filter: [">", ["get", "streifen_flaeche_ha"], min_size] }) }, firstSymbolId ); map.addLayer( { id: "flaeche-label", type: "symbol", source: "data", "source-layer": "lf_an_autobahnen_bahnstrecken_inkl_daten_500ndjson", minzoom: 10, layout: { "text-field": [ "format", [ "number-format", ["get", "streifen_flaeche_ha"], { "max-fraction-digits": 2 } ], " ha", {}, "\n", {}, "", { "text-font": ["literal", ["DIN Offc Pro Italic"]], "font-scale": 0.8 }, ["get", "crop_type"], { "text-font": ["literal", ["DIN Offc Pro Italic"]], "font-scale": 0.8 } ], "text-font": ["DIN Offc Pro Medium", "Arial Unicode MS Bold"], "text-size": 12 }, paint: { "text-opacity": [ "interpolate", ["linear"], ["zoom"], 0, 0, 11, 0, 12, 1 ] }, ...(nat_ausschliessen.length ? { filter: [ "all", ["==", ["typeof", ["get", "natura2000"]], "null"], [">", ["get", "streifen_flaeche_ha"], min_size] ] } : { filter: [">", ["get", "streifen_flaeche_ha"], min_size] }) }, firstSymbolId ); viewof nuts3_agg.addEventListener("input", () => evalNUTS3Aggregate()); function evalNUTS3Aggregate() { const curSelection = viewof nuts3_agg.value; if (map.getLayer("NUTS3")) { map.removeLayer("NUTS3"); } let steps, prop; if (curSelection == "Installierbare Leistung (MWp)") { steps = stepsMWp; prop = "Installierbare Leistung (MWp)"; } else if (curSelection == "Stromerzeugungspotenzial (MWh)") { steps = stepsMWh; prop = "Stromerzeugungspotenzial (MWh)"; } if (steps) { map.addLayer( { id: "NUTS3", type: "fill", source: "nuts3", layout: {}, paint: { "fill-color": [ "interpolate", ["linear"], ["get", prop], ...steps ], "fill-opacity": [ "interpolate", ["linear"], ["zoom"], 0.75, 0.75, 7, 0.75, 8, 0.0 ] } }, firstSymbolId ); } } evalNUTS3Aggregate(); const popup = new maplibregl.Popup({ closeButton: false, closeOnClick: false }); const maxzoom_popup = 8; map.on("mousemove", "NUTS3", function (e) { map.getCanvas().style.cursor = "pointer"; const props = e.features[0].properties; if (map.getZoom() < maxzoom_popup) { popup .setLngLat(e.lngLat) .setHTML( ` ${props.NUTS_NAME}
Installierbare Leistung: ${props[ "Installierbare Leistung (MWp)" ].toLocaleString()} MWp
Stromerzeugungspotenzial: ${props[ "Stromerzeugungspotenzial (MWh)" ].toLocaleString()} MWh
Gesamtfläche${min_size ? " >" + min_size + " ha" : ""}: ${props[ "Geeignete Flächen (ha)" ].toLocaleString()} ha ` ) .addTo(map); } }); map.on("mouseleave", "NUTS3", function () { map.getCanvas().style.cursor = ""; popup.remove(); }); map.on("zoom", () => { if (map.getZoom() > 8 && popup) popup.remove(); }); }); function asyncLoadImage(url) { return new Promise((resolve, reject) => { map.loadImage(url, (err, img) => { if (err) return reject(err); return resolve(img); }); }); } invalidation.then(() => map.remove()); }34:[{"id":"e795e3403a54716a85eb8c951e46cdb723478ea46945bcbc871219b97d46d641b23d918d34b37d9860fead427abee3e7a5c15621a55958ad854052e3a4d33339","href":"https://static.observableusercontent.com/files/e795e3403a54716a85eb8c951e46cdb723478ea46945bcbc871219b97d46d641b23d918d34b37d9860fead427abee3e7a5c15621a55958ad854052e3a4d33339","name":"lf_an_autobahnen_bahnstrecken_inkl_daten_500.parquet","type":"application/octet-stream","lastModified":1679524357062.03,"size":12050153},{"id":"041610c95fe0043fe196737c23443aa327388dc0e19c0d527943affae28e1b95a8a308132831869b1e1c753931a752070e917c4f5d1062ab33778699b315d7a8","href":"https://static.observableusercontent.com/files/041610c95fe0043fe196737c23443aa327388dc0e19c0d527943affae28e1b95a8a308132831869b1e1c753931a752070e917c4f5d1062ab33778699b315d7a8","name":"gemeinden_nuts3@3.csv","type":"text/csv","lastModified":1679524109728.714,"size":466683},{"id":"05ab906de92bf33f2719681df3bcd6be3c994cbf3a4423991f04d53090628fbee0b76143ea3ef83c2c5029a0339d8dc3afae674081ec85a467504de70c65e06e","href":"https://static.observableusercontent.com/files/05ab906de92bf33f2719681df3bcd6be3c994cbf3a4423991f04d53090628fbee0b76143ea3ef83c2c5029a0339d8dc3afae674081ec85a467504de70c65e06e","name":"ID_3697_Sonnenkollektoren_IMG_1539_large.jpeg","type":"image/jpeg","lastModified":1679524109728.714,"size":621630},{"id":"e2fc9eb7c257a5164296e37bbffb39dd35d35ec88313ac61c366d2b71c04c63cf63b44135de4a4c2c4b957444743cbcf61fcbea112e7d58703ac4ad42643d41c","href":"https://static.observableusercontent.com/files/e2fc9eb7c257a5164296e37bbffb39dd35d35ec88313ac61c366d2b71c04c63cf63b44135de4a4c2c4b957444743cbcf61fcbea112e7d58703ac4ad42643d41c","name":"NUTS250_N3 (1).json","type":"application/json","lastModified":1679524109728.714,"size":3547497},{"id":"f0579f93d66f30361e9c60238ed133b28d6900fd07e4622d55fcef9ce7d70443cdc1f711e85e7ef47b9893453fbfcfc0f89b49103fc4687cdf1f3b6b70df079c","href":"https://static.observableusercontent.com/files/f0579f93d66f30361e9c60238ed133b28d6900fd07e4622d55fcef9ce7d70443cdc1f711e85e7ef47b9893453fbfcfc0f89b49103fc4687cdf1f3b6b70df079c","name":"positron@1.json","type":"application/json","lastModified":1679524109728.714,"size":186531}] 30:{"title":"Freiflächen-PV an Autobahnen & Schienen in Deutschland (500m - EEG 2023)","cells":[{"id":0,"value":"# Freiflächen-PV an Autobahnen & Schienen in Deutschland (500m - EEG 2023)","pinned":false,"hidden":false,"mode":"md"},{"id":882,"value":"[➔ Hier geht's zur Analyse des 200m Puffers nach BGBl 2023](https://observablehq.com/@thuenen-institute/freiflachen-pv-an-autobahnen-schienen-in-deutschland-200m?collection=@thuenen-institute/freiflachen-pv)","pinned":false,"hidden":false,"mode":"md"},{"id":5,"value":"md`von *Christoph Pahmeyer, Jonas Böhm, Stefan Erasmi, Gideon Tetteh, Alexander Gocht*`","pinned":false,"hidden":false,"mode":"ojs"},{"id":171,"value":"md`Ackerflächen, die in einem Abstand von bis zu 500 Metern längs von Autobahnen und Schienenwegen verlaufen, haben einen Anspruch auf Einspeisevergütung nach dem Erneuerbaren Energien Gesetz (2023).`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1026,"value":"htl.html` $\"Sonnenkollektoren$ `","pinned":false,"hidden":false,"mode":"ojs"},{"id":1046,"value":"md`##### Freiflächen-PV Anlage entlang einer Bahnstrecke. © Michael Welling`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1132,"value":"Um mögliche Flächenkokurrenzen besser abschätzen zu können, die sich aus der Gesetzesänderung ergeben, hat das Thünen-Institut für Betriebswirtschaft eine erste Analyse der betroffenen Flächen anhand von Geodaten vorgenommen. ","pinned":false,"hidden":false,"mode":"md"},{"id":1134,"value":"Die Grundfläche im Abstand von 500 Metern längs von Autobahnen und Schienenwegen beträgt ${(\n (flaeche_autobahn_ha + flaeche_schienen_ha) /\n 100\n).toLocaleString()} km² (${perc_autobahn}% Autobahn, ${perc_schiene}% Schienenwege). Dabei wurde bei den Autobahnen ein Puffer von 40 Metern abgezogen, um das generelle Bauverbot von 40 Meter neben der Autobahn zu berücksichtigen.","pinned":false,"hidden":false,"mode":"md"},{"id":1136,"value":"In dieser Flächenkulisse wurde eine landwirtschaftlich genutzte Fläche (LF) von **${total_area.toLocaleString()} ha** auf ${tb\n .totalRows()\n .toLocaleString()} Schlägen identifiziert. Das entspricht ca. ${perc_of_total}% der Grundfläche im Abstand von 500 Metern längs von Autobahnen und Schienenwegen und ca. ${(\n (total_area / 16600000) *\n 100\n).toFixed(1)}% der gesamten LF Deutschlands.","pinned":false,"hidden":false,"mode":"md"},{"id":1139,"value":"Die restlichen ${\n 100 - perc_of_total\n}% der Grundfläche sind andere Nutzungsarten wie Wald, Siedlung oder Wasser.","pinned":false,"hidden":false,"mode":"md"},{"id":1141,"value":"Von der LF entfallen ca. ${\n areas_per_land_use_perc.Ackerland\n}% der Fläche auf Ackerland (${areas_per_land_use.Ackerland.toLocaleString()} ha), ${\n areas_per_land_use_perc.Dauergrünland\n}% auf Dauergrünland (${areas_per_land_use.Dauergrünland.toLocaleString()} ha) und ${\n areas_per_land_use_perc.Andere\n}% auf andere Flächen (${areas_per_land_use.Andere.toLocaleString()} ha, nicht identifizierte Flächen und Dauerkulturen). Etwa ${(\n Math.round((natura2000_area / total_area) * 1000) / 10\n).toLocaleString()}% der LF (${natura2000_area.toLocaleString()} ha) liegen innerhalb von Naturschutzgebieten.","pinned":false,"hidden":false,"mode":"md"},{"id":76,"value":"$32","pinned":false,"hidden":false,"mode":"md"},{"id":367,"value":"md`### Berechnung der Flächenkulisse und des theoretischen Ausbaupotentials`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1082,"value":"hinweis = md`**Hinweis:** Die in der Analyse betrachteten Schläge sind zusammenhängende Flächen, die einheitlich mit einer Kultur bebaut und per Fernerkundung identifiziert wurden. Eigentumsverhältnisse sowie mögliche Nutzungskooperationen zwischen Eigentümer:innen sind nicht bekannt und können daher nicht berücksichtigt werden.\n\nIn der praktischen Umsetzung könnten einige der identifizierten Flächen zusammengelegt werden, wodurch sich größere (und damit wirtschaftlichere) Flächengrößen erreichen ließen.`","pinned":false,"hidden":false,"mode":"ojs"},{"id":290,"value":"viewof min_size = Inputs.range([0, 10], {\n label:\n \"Annahme über min. Flächengröße für Freiflächen-PV (ha) ➔ Hinweis beachten!\",\n value: 0,\n step: 1,\n submit: true\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":298,"value":"viewof nat_ausschliessen = Inputs.checkbox([\" \"], {\n label: \"Naturschutzflächen ausschließen?\",\n value: [\" \"]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":461,"value":"md`Werden Naturschutzflächen ${nat_ausschliessen.length ? \"\" : \"nicht \"} ${\n min_size ? `und Schläge mit einer Größe von < ${min_size} ha` : \"\"\n} kategorisch ausgeschlossen, kommen ca. ${total_area_filtered.toLocaleString()} ha (${_.round(\n (total_area_filtered / total_area) * 100,\n 2\n).toLocaleString()}% der Gesamtfläche) auf ${tb_gt_treshold\n .totalRows()\n .toLocaleString()} Schlägen in Betracht für die PV-Kulisse. Die durchschnittliche Schlaggröße beträgt dabei ${mean_area.mean_area.toLocaleString()} ha (Median: ${mean_area.median_area.toLocaleString()} ha).`","pinned":false,"hidden":false,"mode":"ojs"},{"id":864,"value":"viewof fia_simple = Inputs.range([0, 5], {\n label: \"Flächenbedarf je MWp (ha/MWp)\",\n value: 1.1,\n step: 0.1\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":502,"value":"md`Unter der vereinfachten Annahme, dass je Hektar ${_.round(\n 1 / fia_simple,\n 2\n)} MWp Freiflächen-PV installiert werden können, ergibt sich ein theoretisches Potenzial von bis zu ${_.round(\n energy_sums.totalInstallierbarMWp / 1000,\n 2\n).toLocaleString()} GWp installierter Leistung auf den ausgewählten Flächen. Dies entspricht ca. ${_.round(\n (energy_sums.totalInstallierbarMWp / 1000 / 400) * 100,\n 1\n)}% des im EEG 2023 definierten Ausbauziels von 400 GW installierter Leistung von Solaranlagen bis 2040 ([EEG, 2023](https://www.gesetze-im-internet.de/eeg_2014/__4.html)).`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1150,"value":"md`Auf Basis der räumlich differenzierten Strahlungsintensitäten ergibt sich daraus ein Stromerzeugungspotenzial von bis zu **${total_TWh.toLocaleString()} TWh pro Jahr**. Dies entspricht ca. ${perc_energie_verbrauch_d.toLocaleString()}% der Bruttostromerzeugung in Deutschland im Jahr 2021 ([BMWK, 2023](https://www.bmwk.de/Redaktion/DE/Dossier/strommarkt-der-zukunft.html#:~:text=In%20Deutschland%20wurden%20im%20Jahr,das%20sind%20588%20Milliarden%20kWh%20.)).`","pinned":false,"hidden":false,"mode":"ojs"},{"id":545,"value":"md`### Karte`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1020,"value":"viewof hintergrund = Inputs.radio([\"Basiskarte\", \"Satellit\"], {\n label: \"Hintergrundkarte\",\n value: \"Basiskarte\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1319,"value":"viewof nuts3_agg = Inputs.radio(\n [\"Keine\", \"Installierbare Leistung (MWp)\", \"Stromerzeugungspotenzial (MWh)\"],\n {\n label: \"Aggregation auf Kreisebene anzeigen?\",\n value: \"Installierbare Leistung (MWp)\"\n }\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":659,"value":"Swatches(\n d3.scaleOrdinal(\n [\n \"Autobahn (40-200m Streifen)\",\n \"Bahnstrecken (200m Streifen)\",\n `Landwirtschaftliche Flächen ${min_size ? \"> \" + min_size + \" ha\" : \"\"}`,\n \"Naturschutzgebiete\"\n ],\n [\"#37464A\", \"#DF7C28\", \"#0B95C7\", `url(${naturschutz_pattern_png})`]\n )\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1294,"value":"{\n if (nuts3_agg === \"Installierbare Leistung (MWp)\") {\n return Legend(\n d3.scaleSequential(\n [0, maxMWp],\n d3.interpolateRgbBasis([\"#FFFFFF\", \"#479EDB\", \"#204F75\"])\n ),\n {\n title: \"Installierbare Leistung nach Landkreisen (MWp)\"\n }\n );\n } else if (nuts3_agg === \"Stromerzeugungspotenzial (MWh)\") {\n return Legend(\n d3.scaleSequential(\n [0, maxMWh],\n d3.interpolateRgbBasis([\"#FFFFFF\", \"#E9C54B\", \"#47332A\"])\n ),\n {\n title: \"Stromerzeugungspotenzial nach Landkreisen (MWh)\"\n }\n );\n }\n return md``;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":548,"value":"$33","pinned":false,"hidden":false,"mode":"ojs"},{"id":1357,"value":"md`**Hinweis**: Zoomen Sie in der Karte auf eine Autobahn oder Bahnstrecke, um einzelne Flächen anzusehen.`","pinned":false,"hidden":false,"mode":"ojs"},{"id":552,"value":"md`### Theoretisches Ausbaupotential nach Gemeinden`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1191,"value":"data","pinned":true,"hidden":false,"mode":"ojs"},{"id":1155,"value":"### Einordnung der Ergebnisse ###","pinned":false,"hidden":false,"mode":"md"},{"id":1166,"value":"Die im Notebook angegebenen Flächennutzungspotentiale sind eine maximale Obergrenze und stellen kein realistisches Ausbaupotential dar. Weitere Einschränkungen müssten beachtet werden, wie die Verschattung, Struktur der Eigentumsrechte der Flächen, Netzanschlussmöglichkeiten, Mindestabstand zu Siedlungsflächen und vieles mehr. Zudem muss eine Mindestflächengröße angenommen werden, um ökonomisch sinnvoll die Anlage betreiben zu könne. Für eine detaillierte Diskussion solcher Aspekte siehe auch [Kelm et al. (2019)](https://www.erneuerbare-energien.de/EE/Redaktion/DE/Downloads/bmwi_de/zsv-boschundpartner-vorbereitung-begleitung-eeg.pdf?__blob=publicationFile&v=7). Zudem weisen die Autoren darauf hin, dass mit der Berechnung und grafischen Darstellung der Flächen keine Bewertung zu der generellen Vorzüglichkeit vorgenommen wurde. Dazu sind weitere Schritte notwendig, wie beispielweise der Verschneidung von Ertragspotentialen der Landwirtschaft oder der Analyse der Eigentümerrechte.","pinned":false,"hidden":false,"mode":"md"},{"id":1055,"value":"### Warum Observable? ###\n","pinned":false,"hidden":false,"mode":"md"},{"id":1059,"value":"Ein Observable notebook ist ein interaktives Dokument, das es Ihnen ermöglicht, Daten zu erforschen, zu analysieren und zu visualisieren. Sie können Code ändern und selber, Daten importieren, Grafiken erstellen und Ihre Ergebnisse mit anderen teilen. Sie können auch Fragen stellen, Feedback geben und so neue Erkenntnisse gewinnen.","pinned":false,"hidden":false,"mode":"md"},{"id":1173,"value":"### Weiterführende Links bzw. Literatur ###","pinned":false,"hidden":false,"mode":"md"},{"id":1176,"value":"- Feuerbacher, A., Herrmann, T., Neuenfeldt, S., Laub, M., & Gocht, A. (2022). Estimating the economics and adoption potential of agrivoltaics in Germany using a farm-level bottom-up approach. Renewable and Sustainable Energy Reviews, 168, 112784., [DOI:10.1016/j.rser.2022.112784](https://www.sciencedirect.com/science/article/abs/pii/S1364032122006682)\n- [Photovoltaik auf Agrarflächen – für eine schnelle Energiewende, 2022](https://www.thuenen.de/de/themenfelder/langfristige-politikkonzepte/pv-auf-agrarflaechen)","pinned":false,"hidden":false,"mode":"md"},{"id":1117,"value":"md`### Anhang`","pinned":false,"hidden":false,"mode":"ojs"},{"id":195,"value":"md`#### Tabelle A1: Identifizierte Gesamtflächen nach Kulturen innerhalb einer Entfernung von 200m von Autobahnen und zweigleisigen Schienenwegen in Deutschland`","pinned":false,"hidden":false,"mode":"ojs"},{"id":83,"value":"Inputs.table(_.sortBy(areas_per_crop_arr, \"Fläche (ha)\").reverse(), {\n format: {\n \"Fläche (ha)\": (d) => Math.round(d).toLocaleString()\n }\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1126,"value":"md`#### Abbildung A1: Verteilung der Schlaggrößen innerhalb einer Entfernung von 200m von Autobahnen und zweigleisigen Schienenwegen in Deutschland. `","pinned":false,"hidden":false,"mode":"ojs"},{"id":898,"value":"Plot.plot({\n x: {\n label: \"Schlaggröße (ha)\"\n },\n marks: [\n Plot.barY(hist_data, {\n x: \"bin\",\n y: \"count\",\n sort: { x: \"y\", reverse: true }\n }),\n Plot.ruleY([0])\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":527,"value":"md`### Technischer Anhang - Formeln und Daten`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1261,"value":"minMWp = 0","pinned":false,"hidden":false,"mode":"ojs"},{"id":1267,"value":"maxMWp = d3.max(\n geojson_nuts3.features.map(\n (f) => f.properties[\"Installierbare Leistung (MWp)\"]\n )\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1277,"value":"minMWh = 0","pinned":true,"hidden":false,"mode":"ojs"},{"id":1279,"value":"maxMWh = d3.max(\n geojson_nuts3.features.map(\n (f) => f.properties[\"Stromerzeugungspotenzial (MWh)\"]\n )\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1270,"value":"stepsMWp = d3.range(9).reduce((bla, val, i, arr) => {\n bla.push(\n Math.round(((maxMWp - minMWp) / (arr.length + 1)) * (val + 1)) + minMWp\n );\n bla.push(\n d3\n .color(\n d3.interpolateRgbBasis([\"#FFFFFF\", \"#479EDB\", \"#204F75\"])(\n (val + 1) / (arr.length + 1)\n )\n )\n .formatHex()\n // d3.color(d3.interpolateYlOrBr((val + 1) / (arr.length + 1))).formatHex()\n );\n return bla;\n}, [])","pinned":true,"hidden":false,"mode":"ojs"},{"id":1325,"value":"stepsMWh = d3.range(9).reduce((bla, val, i, arr) => {\n bla.push(\n Math.round(((maxMWh - minMWh) / (arr.length + 1)) * (val + 1)) + minMWh\n );\n bla.push(\n d3\n .color(\n d3.interpolateRgbBasis([\"#FFFFFF\", \"#E9C54B\", \"#47332A\"])(\n (val + 1) / (arr.length + 1)\n )\n )\n .formatHex()\n // d3.color(d3.interpolateYlOrBr((val + 1) / (arr.length + 1))).formatHex()\n );\n return bla;\n}, [])","pinned":true,"hidden":false,"mode":"ojs"},{"id":555,"value":"data = energy_tb\n .groupby(\"GEN\", \"AGS\")\n .rollup({\n \"Geeignete Flächen (ha)\": (d) => Math.round(op.sum(d.streifen_flaeche_ha)),\n \"Installierbare Leistung (MWp)\": (d) => Math.round(op.sum(d.MWp)),\n \"Stromerzeugungspotenzial (MWh)\": (d) => Math.round(op.sum(d.MWh))\n })\n .objects()","pinned":false,"hidden":false,"mode":"ojs"},{"id":1211,"value":"nuts3_gemeinden = d3.group(\n await FileAttachment(\"gemeinden_nuts3@3.csv\").csv({ typed: true }),\n (g) => g.NUTS_CODE\n)","pinned":true,"hidden":false,"mode":"ojs"},{"id":1224,"value":"nuts3_topo = FileAttachment(\"NUTS250_N3 (1).json\").json()","pinned":false,"hidden":false,"mode":"ojs"},{"id":1218,"value":"geojson_nuts3 = {\n const fc = topojson.feature(nuts3_topo, nuts3_topo.objects.NUTS250_N3);\n fc.features.forEach((region, i) => {\n // Gemeinden für NUTS 3 Region abrufen\n const gemeinden = nuts3_gemeinden.get(region.properties.NUTS_CODE);\n // Potentiale aggregieren\n const gemeinde_daten = gemeinden.flatMap((g) =>\n data.filter((r) => Number(r.AGS) == g.AGS)\n );\n fc.features[i].properties[\"Installierbare Leistung (MWp)\"] = d3.sum(\n gemeinde_daten.map((r) => r[\"Installierbare Leistung (MWp)\"])\n );\n fc.features[i].properties[\"Stromerzeugungspotenzial (MWh)\"] = d3.sum(\n gemeinde_daten.map((r) => r[\"Stromerzeugungspotenzial (MWh)\"])\n );\n fc.features[i].properties[\"Geeignete Flächen (ha)\"] = d3.sum(\n gemeinde_daten.map((r) => r[\"Geeignete Flächen (ha)\"])\n );\n });\n return fc;\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":379,"value":"total_TWh = _.round(energy_sums.totalEnergieMWh / (1000 * 1000), 2) //","pinned":true,"hidden":false,"mode":"ojs"},{"id":394,"value":"perc_energie_verbrauch_d = _.round((total_TWh / 588.8) * 100, 2) // Gesamtproduktion Deutschland https://www.bmwk.de/Redaktion/DE/Dossier/strommarkt-der-zukunft.html#:~:text=In%20Deutschland%20wurden%20im%20Jahr,das%20sind%20588%20Milliarden%20kWh%20.","pinned":true,"hidden":false,"mode":"ojs"},{"id":269,"value":"hist_data = tb_gt_treshold\n .filter((d) => d.streifen_flaeche_ha < 40)\n .groupby({ bin: aq.bin(\"streifen_flaeche_ha\") })\n .count()\n .impute(\n { count: () => 0 }, // set imputed counts to zero\n { expand: { bin: (d) => op.sequence(...op.bins(d.bin)) } } // include rows for all bin values\n // Note: op.bins(d.bin) -> [start, stop, step]\n // so op.sequence(...op.bins(d.bin)) -> op.sequence(start, stop, step)\n )\n .orderby(\"bin\")\n .objects()","pinned":false,"hidden":false,"mode":"ojs"},{"id":420,"value":"function findClosestIndex(arr, goal) {\n const closestVal = arr.reduce((prev, curr) =>\n Math.abs(curr - goal) < Math.abs(prev - goal) ? curr : prev\n );\n return arr.findIndex((v) => v === closestVal);\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":480,"value":"mean_area = tb_gt_treshold\n .rollup({\n mean_area: (d) => Math.round(op.mean(d.streifen_flaeche_ha) * 100) / 100,\n median_area: (d) => Math.round(op.median(d.streifen_flaeche_ha) * 100) / 100\n })\n .object()","pinned":false,"hidden":false,"mode":"ojs"},{"id":208,"value":"energy_tb = tb_gt_treshold.derive({\n MWp: aq.escape((d) =>\n calc_installierbare_leistung_MWp(d.streifen_flaeche_ha)\n ),\n pvout: aq.escape((d) => Number(d.PVOUT_centroid)),\n MWh: aq.escape((d) =>\n calc_energie_pro_zeiteinheit_MWh(d.streifen_flaeche_ha, d.PVOUT_centroid)\n ),\n ...monate.reduce((arr, monat) => {\n arr[monat.eng_key] = aq.escape((d) =>\n calc_energie_pro_zeiteinheit_MWh(d.streifen_flaeche_ha, d[monat.eng_key])\n );\n return arr;\n }, {})\n})","pinned":true,"hidden":false,"mode":"ojs"},{"id":252,"value":"tb_gt_treshold = tb\n .filter(aq.escape((d) => d.streifen_flaeche_ha > min_size))\n .filter(aq.escape((d) => (nat_ausschliessen.length ? !d.natura2000 : d)))","pinned":true,"hidden":false,"mode":"ojs"},{"id":237,"value":"energy_sums = energy_tb\n .rollup({\n totalInstallierbarMWp: (d) => Math.round(op.sum(d.MWp)),\n avgPVOUT: (d) => Math.round(op.average(d.pvout)),\n totalEnergieMWh: (d) => Math.round(op.sum(d.MWh)),\n PVOUT_01_centroid: (d) => Math.round(op.sum(d.PVOUT_01_centroid)),\n PVOUT_02_centroid: (d) => Math.round(op.sum(d.PVOUT_02_centroid)),\n PVOUT_03_centroid: (d) => Math.round(op.sum(d.PVOUT_03_centroid)),\n PVOUT_04_centroid: (d) => Math.round(op.sum(d.PVOUT_04_centroid)),\n PVOUT_05_centroid: (d) => Math.round(op.sum(d.PVOUT_05_centroid)),\n PVOUT_06_centroid: (d) => Math.round(op.sum(d.PVOUT_06_centroid)),\n PVOUT_07_centroid: (d) => Math.round(op.sum(d.PVOUT_07_centroid)),\n PVOUT_08_centroid: (d) => Math.round(op.sum(d.PVOUT_08_centroid)),\n PVOUT_09_centroid: (d) => Math.round(op.sum(d.PVOUT_09_centroid)),\n PVOUT_10_centroid: (d) => Math.round(op.sum(d.PVOUT_10_centroid)),\n PVOUT_11_centroid: (d) => Math.round(op.sum(d.PVOUT_11_centroid))\n })\n .object()","pinned":true,"hidden":false,"mode":"ojs"},{"id":206,"value":"function calc_installierbare_leistung_MWp(ha) {\n // Flächeninanspruchnahme berechnen\n // const values = [\n // (750 / 1000) * fia.fia_750,\n // (2000 / 1000) * fia.fia_2000,\n // (5000 / 1000) * fia.fia_5000,\n // (20000 / 1000) * fia.fia_20000\n // ];\n // // Flächeninanspruchnahme auswählen, die der Schlaggröße am nächsten kommt\n // const fias = [fia.fia_750, fia.fia_2000, fia.fia_5000, fia.fia_20000];\n // const index = findClosestIndex(values, ha);\n // const MWp = ha / fias[index];\n const MWp = ha / fia_simple;\n // wenn die installierbare Leistung negativ wird (zu kleine Fläche)\n // wird eine installierbare Leistung von 0 angenommen\n return MWp < 0 ? 0 : MWp;\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":217,"value":"function calc_energie_pro_zeiteinheit_MWh(ha, pvout) {\n return calc_installierbare_leistung_MWp(ha) * Number(pvout);\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":706,"value":"monate = [\n { monat: \"Januar\", eng_key: \"PVOUT_01_centroid\" },\n { monat: \"Februar\", eng_key: \"PVOUT_02_centroid\" },\n { monat: \"März\", eng_key: \"PVOUT_03_centroid\" },\n { monat: \"April\", eng_key: \"PVOUT_04_centroid\" },\n { monat: \"Mai\", eng_key: \"PVOUT_05_centroid\" },\n { monat: \"Juni\", eng_key: \"PVOUT_06_centroid\" },\n { monat: \"Juli\", eng_key: \"PVOUT_07_centroid\" },\n { monat: \"August\", eng_key: \"PVOUT_08_centroid\" },\n { monat: \"September\", eng_key: \"PVOUT_09_centroid\" },\n { monat: \"Oktober\", eng_key: \"PVOUT_10_centroid\" },\n { monat: \"Dezember\", eng_key: \"PVOUT_11_centroid\" }\n]","pinned":true,"hidden":false,"mode":"ojs"},{"id":806,"value":"flaeche_autobahn_ha = Math.round(10349011869/ 10000)","pinned":true,"hidden":false,"mode":"ojs"},{"id":809,"value":"flaeche_schienen_ha = Math.round(27679550581 / 10000)","pinned":true,"hidden":false,"mode":"ojs"},{"id":812,"value":"perc_autobahn = _.round(\n (flaeche_autobahn_ha / (flaeche_autobahn_ha + flaeche_schienen_ha)) * 100,\n 1\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":817,"value":"perc_schiene = 100 - perc_autobahn","pinned":false,"hidden":false,"mode":"ojs"},{"id":822,"value":"perc_of_total = _.round(\n (total_area / (flaeche_autobahn_ha + flaeche_schienen_ha)) * 100,\n 1\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":68,"value":"total_area = Math.round(\n tb\n .rollup({\n total: (d) => op.sum(d.streifen_flaeche_ha)\n })\n .object().total\n)","pinned":true,"hidden":false,"mode":"ojs"},{"id":468,"value":"total_area_filtered = Math.round(\n tb_gt_treshold\n .rollup({\n total: (d) => op.sum(d.streifen_flaeche_ha)\n })\n .object().total\n)","pinned":true,"hidden":false,"mode":"ojs"},{"id":126,"value":"natura2000_area = total_area - total_area_no_natura2000","pinned":true,"hidden":false,"mode":"ojs"},{"id":119,"value":"total_area_no_natura2000 = Math.round(\n tb\n .filter((d) => !d.natura2000)\n .rollup({\n total: (d) => op.sum(d.streifen_flaeche_ha)\n })\n .object().total\n)","pinned":true,"hidden":false,"mode":"ojs"},{"id":114,"value":"areas_per_land_use = {\n const permanentGrassland = [\"Permanent grassland\"];\n const other = [\n \"Other agricultural areas\",\n \"Small woody features\",\n \"Grapevine\",\n \"Hops\",\n \"Orchard\",\n \"Small woody features\",\n \"Unbekannt\"\n ];\n return {\n Dauergrünland: Math.round(areas_per_crop[\"Permanent grassland\"]),\n Andere: Math.round(\n _.sum(\n areas_per_crop_arr\n .filter((row) => other.includes(row.Kultur))\n .map((row) => row[\"Fläche (ha)\"])\n )\n ),\n Ackerland: Math.round(\n _.sum(\n areas_per_crop_arr\n .filter(\n (row) =>\n !other.includes(row.Kultur) && row.Kultur != \"Permanent grassland\"\n )\n .map((row) => row[\"Fläche (ha)\"])\n )\n )\n };\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":138,"value":"areas_per_land_use_perc = ({\n Dauergrünland: (\n Math.round((areas_per_land_use[\"Dauergrünland\"] / total_area) * 1000) / 10\n ).toLocaleString(),\n Andere: (\n Math.round((areas_per_land_use[\"Andere\"] / total_area) * 1000) / 10\n ).toLocaleString(),\n Ackerland: (\n Math.round((areas_per_land_use[\"Ackerland\"] / total_area) * 1000) / 10\n ).toLocaleString(1)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":56,"value":"areas_per_crop_arr = {\n return Object.keys(areas_per_crop).map((crop) => ({\n Kultur: crop ? crop : \"Unbekannt\",\n \"Fläche (ha)\": areas_per_crop[crop]\n }));\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":107,"value":"areas_per_crop = tb\n .pivot(\"crop_type\", { total: (d) => op.sum(d.streifen_flaeche_ha) })\n .object()","pinned":true,"hidden":false,"mode":"ojs"},{"id":24,"value":"tb = aq.fromArrow(table)","pinned":true,"hidden":false,"mode":"ojs"},{"id":19,"value":"table = pq.readParquet(new Uint8Array(await data_source.arrayBuffer()))","pinned":true,"hidden":false,"mode":"ojs"},{"id":952,"value":"render = (t) => {\n const HEIGHT = 512;\n const WIDTH = 512;\n\n const svg = d3.create(\"svg\").attr(\"width\", WIDTH).attr(\"height\", HEIGHT);\n\n svg.call(t);\n\n svg\n .append(\"rect\")\n .attr(\"width\", WIDTH)\n .attr(\"height\", HEIGHT)\n .style(\"fill\", t.url());\n\n return svg.node();\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":963,"value":"naturschutz_pattern = render(\n textures.lines().size(4).strokeWidth(1).stroke(\"white\").background(\"#7DBB39\")\n)","pinned":true,"hidden":false,"mode":"ojs"},{"id":980,"value":"naturschutz_pattern_png = rasterize(naturschutz_pattern)","pinned":true,"hidden":false,"mode":"ojs"},{"id":879,"value":"// This function converts svg strings into base64 png urls\n// which you can use for creating patterns in mapbox\nfunction rasterize(svgString, size = 512, interpolation = 1) {\n const svg = html`${svgString}`;\n svg.setAttribute(\"width\", size * interpolation);\n svg.setAttribute(\"height\", size * interpolation);\n\n let resolve, reject;\n const promise = new Promise((y, n) => ((resolve = y), (reject = n)));\n const image = new Image();\n image.onerror = reject;\n image.onload = () => {\n const context = DOM.context2d(size * interpolation, size * interpolation);\n for (let i = 0; i < interpolation; i++) {\n for (let j = 0; j < interpolation; j++) {\n context.drawImage(image, i * size, j * size, size, size);\n }\n }\n resolve(context.canvas.toDataURL(\"image/png\"));\n };\n const serializedSVG = new XMLSerializer().serializeToString(svg);\n const base64Data = window.btoa(serializedSVG);\n const imgsrc = `data:image/svg+xml;base64,${base64Data}`;\n image.src = imgsrc;\n\n return promise;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":11,"value":"data_source = FileAttachment(\"lf_an_autobahnen_bahnstrecken_inkl_daten_500.parquet\")","pinned":true,"hidden":false,"mode":"ojs"},{"id":533,"value":"md`### Technischer Anhang - Programmbibliotheken`","pinned":false,"hidden":false,"mode":"ojs"},{"id":584,"value":"maplibregl = require(\"maplibre-gl@2.4.0\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":950,"value":"textures = require(\"textures\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":587,"value":"html`\n(Maplibre Stylesheet)`","pinned":false,"hidden":false,"mode":"ojs"},{"id":603,"value":"pmtiles = import(\"https://unpkg.com/pmtiles@2.7.0/dist/index.mjs?module\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":32,"value":"import { aq, op } from \"@uwdata/arquero\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":655,"value":"import { Swatches, Legend } from \"@d3/color-legend\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":682,"value":"LegendControl = require(\"mapboxgl-legend\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":686,"value":"html`\n(Mapbox-Legend Stylesheet)`","pinned":false,"hidden":false,"mode":"ojs"},{"id":934,"value":"MaplibreGeocoder = {\n const MaplibreGeocoder = await require(\"@maplibre/maplibre-gl-geocoder\");\n\n document.head.appendChild(\n html``\n );\n return MaplibreGeocoder;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1220,"value":"topojson = require(\"topojson\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":13,"value":"pq = {\n const pq = await import(\"https://unpkg.com/parquet-wasm@0.1.1/web.js\");\n // default seems to need to resolve first?\n await pq.default();\n return pq;\n}","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W34","stdlib":"1","titling":"auto"}