Kali ini kita akan mencoba mengunduh data dari OpenStreetMap berupa Polygon Bangunan. Kita akan coba gunakan Metode XMax dan Xmin, serta Y Max dan Y Min.
Pertama kali saudara Buka Jupyter Notebook nya dan Gunakan Script Berikut.
import os
import requests
import geopandas as gpd
from shapely.geometry import Polygon
data_folder = "data"
os.makedirs(data_folder, exist_ok=True)
def fetch_osm_buildings_bbox(xmin, ymin, xmax, ymax):
print(f"📥 Mengunduh data bangunan dalam BBOX ({xmin}, {ymin}, {xmax}, {ymax})...")
query = f"""
[out:json];
(
way["building"]({ymin},{xmin},{ymax},{xmax});
);
out body; >; out skel qt;
"""
url = "http://overpass-api.de/api/interpreter"
response = requests.get(url, params={"data": query})
if response.status_code == 200:
data = response.json()
if "elements" in data and len(data["elements"]) > 0:
print(f"✅ Data berhasil diunduh dalam BBOX ({xmin}, {ymin}, {xmax}, {ymax}).")
return data
else:
print(f"⚠️ Tidak ada data bangunan yang ditemukan dalam BBOX tersebut.")
else:
print(f"❌ Gagal mengambil data dalam BBOX ({xmin}, {ymin}, {xmax}, {ymax}).")
return None
def convert_buildings_to_gdf(osm_data):
if not osm_data or "elements" not in osm_data:
print("⚠️ Data kosong atau tidak valid.")
return None
buildings = []
nodes = {node["id"]: (node["lon"], node["lat"]) for node in osm_data["elements"] if node["type"] == "node"}
for element in osm_data["elements"]:
if element["type"] == "way" and "nodes" in element:
coordinates = [nodes[node_id] for node_id in element["nodes"] if node_id in nodes]
if len(coordinates) > 2 and coordinates[0] == coordinates[-1]:
buildings.append({
"geometry": Polygon(coordinates),
"building": element.get("tags", {}).get("building", "unknown"),
"name": element.get("tags", {}).get("name", "unknown"),
"levels": element.get("tags", {}).get("building:levels", "unknown"),
"material": element.get("tags", {}).get("building:material", "unknown")
})
if not buildings:
print("⚠️ Tidak ada bangunan yang ditemukan dalam data OSM.")
return None
return gpd.GeoDataFrame(buildings, geometry="geometry", crs="EPSG:4326")
xmin, ymin, xmax, ymax = 107.60235, -6.91000, 107.63879, -6.88949
osm_data = fetch_osm_buildings_bbox(xmin, ymin, xmax, ymax)
if osm_data:
gdf = convert_buildings_to_gdf(osm_data)
if gdf is not None and not gdf.empty:
file_path = os.path.join(data_folder, "buildings_bbox.geojson")
gdf.to_file(file_path, driver="GeoJSON")
print(f"📁 Data bangunan berhasil disimpan di: {file_path}")
else:
print(f"⚠️ Tidak ada data yang dapat disimpan dalam BBOX tersebut.")
print("🎉 Semua proses selesai!")
Pada xmin, ymin, xmax, ymax = , bisa teman-teman gunakan export yang disediakan oleh OpenStreetMap. Contohnya sebagai berikut:
lalu lakukan download, dan setelah berhasil data berada di folder data/buildings_bbox.geojson
untuk memastikan data yang tersedia itu benar-benar ada, gunakan code dibawah:
import geopandas as gpd
import folium
from folium.plugins import MarkerCluster
from IPython.display import display
file_path = "data/buildings_bbox.geojson"
gdf = gpd.read_file(file_path)
gdf = gdf.to_crs(epsg=4326)
m = folium.Map(location=[-6.9175, 107.6191], zoom_start=12, tiles="CartoDB dark_matter")
folium.GeoJson(
gdf,
name="Bangunan",
style_function=lambda feature: {
"fillColor": "yellow",
"color": "orange",
"weight": 1,
"fillOpacity": 0.5
},
tooltip=folium.GeoJsonTooltip(fields=["name", "building"], aliases=["Nama", "Tipe Bangunan"])
).add_to(m)
folium.LayerControl().add_to(m)
display(m)
Hasilnya sebagai berikut:
Demikianlah hasil unduhan yang telah berhasil dilakukan. Apabila data yang didownload lebih besar, disarankan untuk menggunakan Aplikasi GIS Desktop.
