| from math import pi | |
| line_lyr = QgsProject.instance().mapLayersByName('line')[0] | |
| line_geom = line_lyr.getGeometry(1) | |
| results_lyr = QgsVectorLayer(f'Point?crs={line_lyr.sourceCrs().authid()}', 'temp', 'memory') | |
| results_provider = results_lyr.dataProvider() | |
| results_provider.addAttributes([QgsField('id', QVariant.Int), | |
| QgsField('dist_along_line', QVariant.Double), | |
| QgsField('dist_from_previous_point', QVariant.Double), |
| from collections import defaultdict | |
| # Initialize layer variables | |
| streets_lyr = QgsProject.instance().mapLayersByName('streets')[0] | |
| origin_lyr = QgsProject.instance().mapLayersByName('schools')[0] | |
| destination_lyr = QgsProject.instance().mapLayersByName('police_stations')[0] | |
| od_matrix = QgsProject.instance().mapLayersByName('od_pairs')[0] | |
| # Initialize variables for analysis | |
| origin_ids = QgsExpression('array_agg("fid")').evaluate(origin_lyr.createExpressionContext()) |
| create table route_lines as ( | |
| select r.*, s.geom | |
| from ( | |
| select route_id, route_short_name, route_long_name, route_color, min_headway_minutes | |
| from routes | |
| group by route_id | |
| ) r | |
| join ( | |
| select route_id, shape_id | |
| from trips |
Batch command, all GeoJSON in folder —> PostGIS
FOR %%f IN (*.geojson) DO (
ogr2ogr -f PostgreSQL PG:"dbname={dbname} user={user} password={pass}" %%f
)
Batch command, Numbered files —> PostGIS
Run as: to_postgis.bat [base_filename] [file_extenstion] [num_files]. Assumes first file is "[base_filename]_0.[file_extension]".