Skip to content

Commit

Permalink
[pre-commit.ci] auto fixes from pre-commit.com hooks
Browse files Browse the repository at this point in the history
for more information, see https://pre-commit.ci
  • Loading branch information
pre-commit-ci[bot] committed Sep 25, 2023
1 parent 16263f5 commit fb27bd4
Show file tree
Hide file tree
Showing 5 changed files with 118 additions and 126 deletions.
2 changes: 1 addition & 1 deletion LICENSE.md
Original file line number Diff line number Diff line change
Expand Up @@ -385,7 +385,7 @@ notice like this when it starts in an interactive mode:
This program comes with ABSOLUTELY NO WARRANTY; for details
type `show w'.
This is free software, and you are welcome to redistribute it
under certain conditions; type `show c' for details.
under certain conditions; type`show c' for details.
</pre>

The hypothetical commands `show w` and `show c` should show the appropriate parts of the General Public License. Of course, your program's commands might be different; for a GUI interface, you would use an "about box".
Expand Down
1 change: 0 additions & 1 deletion README.md
Original file line number Diff line number Diff line change
Expand Up @@ -21,7 +21,6 @@ To install fmtm-splitter, you can use pip. Here are two options:
- Latest on PyPi:
`pip install fmtm-splitter`


## Using the Container Image

- fmtm-splitter scripts can be used via the pre-built container images.
Expand Down
93 changes: 42 additions & 51 deletions fmtm_splitter/FMTMSplitterBuildings.py
Original file line number Diff line number Diff line change
Expand Up @@ -16,44 +16,48 @@
# along with fmtm-splitter. If not, see <https:#www.gnu.org/licenses/>.
#

import logging
import argparse
import sys
import logging
import os
import sys
from sys import argv

import psycopg2
from osgeo import gdal

# Instantiate logger
log = logging.getLogger(__name__)

def splitByBuildings(aoi: str, # GeoJSON polygon input file
queries: list, # list of SQL queries
dbd: list, # database host, dbname, user, password
):
"""Split the polygon by buildings in the database using an SQL query"""
dbstring = (f"PG:host={dbd[0]} dbname={dbd[1]} "
f"user={dbd[2]} password={dbd[3]}")
dbshell = psycopg2.connect(host=dbd[0], database = dbd[1],
user = dbd[2], password = dbd[3])

def splitByBuildings(
aoi: str, # GeoJSON polygon input file
queries: list, # list of SQL queries
dbd: list, # database host, dbname, user, password
):
"""Split the polygon by buildings in the database using an SQL query."""
dbstring = f"PG:host={dbd[0]} dbname={dbd[1]} " f"user={dbd[2]} password={dbd[3]}"
dbshell = psycopg2.connect(host=dbd[0], database=dbd[1], user=dbd[2], password=dbd[3])
dbshell.autocommit = True
dbcursor = dbshell.cursor()
dbcursor.execute('DROP TABLE IF EXISTS aoi;')
dbcursor.execute("DROP TABLE IF EXISTS aoi;")
# Add the AOI to the database
log.info(f'Writing {aoi} to database as aoi layer.')
gdal.VectorTranslate(dbstring, aoi, layerName='aoi')
dbcursor.execute('DROP TABLE IF EXISTS project_aoi;'
'CREATE TABLE project_aoi AS (SELECT '
'ogc_fid as fid,wkb_geometry AS geom FROM aoi);'
'ALTER TABLE project_aoi ADD PRIMARY KEY(fid);'
'CREATE INDEX project_aoi_idx '
'ON project_aoi USING GIST (geom);'
'DROP TABLE aoi;')
dbcursor.execute('VACUUM ANALYZE')
log.info(f"Writing {aoi} to database as aoi layer.")
gdal.VectorTranslate(dbstring, aoi, layerName="aoi")
dbcursor.execute(
"DROP TABLE IF EXISTS project_aoi;"
"CREATE TABLE project_aoi AS (SELECT "
"ogc_fid as fid,wkb_geometry AS geom FROM aoi);"
"ALTER TABLE project_aoi ADD PRIMARY KEY(fid);"
"CREATE INDEX project_aoi_idx "
"ON project_aoi USING GIST (geom);"
"DROP TABLE aoi;"
)
dbcursor.execute("VACUUM ANALYZE")
for query in queries:
dbcursor.execute(query)
dbcursor.execute('VACUUM ANALYZE')
log.info('Might very well have completed successfully')
dbcursor.execute("VACUUM ANALYZE")
log.info("Might very well have completed successfully")


if __name__ == "__main__":
# Command Line options
Expand All @@ -64,19 +68,13 @@ def splitByBuildings(aoi: str, # GeoJSON polygon input file
epilog="""
This program splits a Polygon (the Area Of Interest)
examples:
"""
""",
)
p.add_argument("-b", "--boundary", required=True,
help="Polygon AOI GeoJSON file")
p.add_argument("-n", "--numfeatures", default=20,
help="Number of features on average desired per task")
p.add_argument("-v", "--verbose",
action="store_true", help="verbose output")
p.add_argument("-o", "--outfile",
default='fmtm.geojson',
help="Output file from splitting")
p.add_argument("-ho", "--host", help="Database host",
default='localhost')
p.add_argument("-b", "--boundary", required=True, help="Polygon AOI GeoJSON file")
p.add_argument("-n", "--numfeatures", default=20, help="Number of features on average desired per task")
p.add_argument("-v", "--verbose", action="store_true", help="verbose output")
p.add_argument("-o", "--outfile", default="fmtm.geojson", help="Output file from splitting")
p.add_argument("-ho", "--host", help="Database host", default="localhost")
p.add_argument("-db", "--database", help="Database to use")
p.add_argument("-u", "--user", help="Database username")
p.add_argument("-p", "--password", help="Database password")
Expand All @@ -87,41 +85,34 @@ def splitByBuildings(aoi: str, # GeoJSON polygon input file
quit()

# if verbose, dump to the terminal.
formatter = logging.Formatter(
"%(threadName)10s - %(name)s - %(levelname)s - %(message)s"
)
formatter = logging.Formatter("%(threadName)10s - %(name)s - %(levelname)s - %(message)s")
level = logging.DEBUG
if args.verbose:
log.setLevel(level)
else:
else:
log.setLevel(logging.INFO)
ch = logging.StreamHandler(sys.stdout)
ch.setLevel(level)
ch.setFormatter(formatter)
log.addHandler(ch)

# log.debug("debug")
# log.info("info")
# log.info("warning")

# Read in the project AOI, a GeoJSON file containing a polygon
aoi = args.boundary
modulardir = os.path.join(os.path.dirname(__file__),
'fmtm-splitter_osm_buildings')
modulardir = os.path.join(os.path.dirname(__file__), "fmtm-splitter_osm_buildings")
modularsqlfiles = [
"fmtm-split_01_split_AOI_by_existing_line_features.sql",
"fmtm-split_02_count_buildings_for_subsplitting.sql",
"fmtm-split_03_cluster_buildings.sql",
"fmtm-split_04_create_polygons_around_clustered_buildings.sql",
"fmtm-split_05_clean_temp_files.sql"
"fmtm-split_05_clean_temp_files.sql",
]
modularqueries = []
for sqlfile in modularsqlfiles:
with open(os.path.join(modulardir, sqlfile), 'r') as sql:
modularqueries.append(sql.read()
.replace('{%numfeatures%}',
str(args.numfeatures)))
with open(os.path.join(modulardir, sqlfile), "r") as sql:
modularqueries.append(sql.read().replace("{%numfeatures%}", str(args.numfeatures)))
dbdetails = [args.host, args.database, args.user, args.password]
features = splitByBuildings(aoi,
modularqueries,
dbdetails)
features = splitByBuildings(aoi, modularqueries, dbdetails)
93 changes: 47 additions & 46 deletions fmtm_splitter/overpass.py
Original file line number Diff line number Diff line change
@@ -1,81 +1,82 @@
#!/bin/python3

import sys, os
import argparse
import requests
import os
import subprocess
from datetime import datetime

import requests
from osgeo import ogr
import subprocess


def aoiextent(aoifile):
"""
Accept a GeoJSON file, return its extent as a bbox string
"""
indriver = ogr.GetDriverByName('GeoJSON')
"""Accept a GeoJSON file, return its extent as a bbox string."""
indriver = ogr.GetDriverByName("GeoJSON")
indata = indriver.Open(aoifile)
inlayer = indata.GetLayer()
e = list(inlayer.GetExtent())
bboxstring = f'{e[2]},{e[0]},{e[3]},{e[1]}'
return(bboxstring)

bboxstring = f"{e[2]},{e[0]},{e[3]},{e[1]}"
return bboxstring


def query(query_string, overpass_url):
"""
Accept a query in Overpass API query language,
"""Accept a query in Overpass API query language,
return an osm dataset.
"""
try:
response = requests.get(overpass_url,
params={'data': query_string})
response = requests.get(overpass_url, params={"data": query_string})
except:
print("overpass did not want to answer that one\n")
if response.status_code == 200:
print(f'The overpass API at {overpass_url} accepted the query and '\
f'returned something.')
print(f"The overpass API at {overpass_url} accepted the query and " f"returned something.")
return response.text
else:
print(response)
print("Yeah, that didn't work. We reached the Overpass API but "\
"something went wrong on the server side.")
print("Yeah, that didn't work. We reached the Overpass API but " "something went wrong on the server side.")


def dbpush(infile, dbd):
"""
Accept an osm file, push it to PostGIS layers using the Underpass schema
"""
"""Accept an osm file, push it to PostGIS layers using the Underpass schema."""
try:
print(f'Trying to turn {infile} into a PostGIS layer')
style = os.path.join('fmtm_splitter','raw.lua')
pg = ["osm2pgsql", "--create",
"-d",
f"postgresql://{dbd[0]}:{dbd[1]}@{dbd[2]}:{dbd[4]}/{dbd[3]}",
"--extra-attributes", "--output=flex",
"--style", style, infile]
print(pg) # just to visually check that this command makes sense
p = subprocess.run(pg, capture_output=True, encoding='utf-8')
print(f"Trying to turn {infile} into a PostGIS layer")
style = os.path.join("fmtm_splitter", "raw.lua")
pg = [
"osm2pgsql",
"--create",
"-d",
f"postgresql://{dbd[0]}:{dbd[1]}@{dbd[2]}:{dbd[4]}/{dbd[3]}",
"--extra-attributes",
"--output=flex",
"--style",
style,
infile,
]
print(pg) # just to visually check that this command makes sense
p = subprocess.run(pg, capture_output=True, encoding="utf-8")
response = p.stdout
error = p.stderr
print(f'osm2pgsql seems to have accepted {infile} and '\
f'returned {response} \nand\n{error}')
print(f"osm2pgsql seems to have accepted {infile} and " f"returned {response} \nand\n{error}")
return response
except Exception as e:
print(e)


if __name__ == "__main__":
"""return a file of raw OSM data from Overpass API from an input file
of text containing working Overpass Query Language, and push that file
to a PostGIS database as a layer.
"""
p = argparse.ArgumentParser(usage="usage: attachments [options]")
p.add_argument('-q', '--query', help="Text file in overpass query language")
p.add_argument('-b', '--boundary', help="AOI as GeoJSON file")
p.add_argument('-url', '--overpass_url', help='Overpass API server URL',
default="https://overpass.kumi.systems/api/interpreter")
p.add_argument("-ho", "--host", help="Database host",
default='localhost')
p.add_argument("-q", "--query", help="Text file in overpass query language")
p.add_argument("-b", "--boundary", help="AOI as GeoJSON file")
p.add_argument(
"-url", "--overpass_url", help="Overpass API server URL", default="https://overpass.kumi.systems/api/interpreter"
)
p.add_argument("-ho", "--host", help="Database host", default="localhost")
p.add_argument("-db", "--database", help="Database to use")
p.add_argument("-u", "--user", help="Database username")
p.add_argument("-p", "--password", help="Database password")
p.add_argument("-po", "--port", help="Database port",
default='5432')
p.add_argument("-po", "--port", help="Database port", default="5432")

args = p.parse_args()

Expand All @@ -84,16 +85,16 @@ def dbpush(infile, dbd):
(basefilename, extension) = os.path.splitext(args.boundary)
date = datetime.now().strftime("%Y_%m_%d")
dirdate = os.path.join(directory, date)
osmfilepath = f'{dirdate}_{basename_no_ext}.osm'
osmfilepath = f"{dirdate}_{basename_no_ext}.osm"

q = open(args.query)
# TODO get bbox from GeoJSON aoi
bbox = aoiextent(args.boundary)
qstring = q.read().replace('{{bbox}}', bbox)
qstring = q.read().replace("{{bbox}}", bbox)
data = query(qstring, args.overpass_url)
with open(osmfilepath, 'w') as of:
with open(osmfilepath, "w") as of:
of.write(data)
print(f'Wrote {osmfilepath}')
print(f"Wrote {osmfilepath}")

dbdetails = [args.user, args.password, args.host, args.database, args.port]
dblayers = dbpush(osmfilepath, dbdetails)
Loading

0 comments on commit fb27bd4

Please sign in to comment.