Le portail francophone de la géomatique

“Introducción a los Sistemas de Información Geográfica con QGIS”

# List all open 3D map canvases
iface.mapCanvases3D()
# [<qgis._3d.Qgs3DMapCanvas object at 0x7f23491b5e10>]
canvas = iface.mapCanvases3D()[0]

# Let's change some settings!
ms = canvas.mapSettings()
ms.setEyeDomeLightingEnabled(True)
ms.setBackgroundColor(QColor('beige'))
ms.setTerrainMapTheme('dtm')
ms.setFieldOfView(110)

# Move the camera to look at specific map coordinates in layer's CRS
cam = canvas.cameraController()
mapPoint = QgsVector3D(-498175.92, -1205400.58, 210)
worldPoint = ms.mapToWorldCoordinates(mapPoint)
cam.setLookingAtPoint(worldPoint, 60, 45, 100)

# Create four new 3D map views
c1 = iface.createNewMapCanvas3D('South View')
c2 = iface.createNewMapCanvas3D('West View')
c3 = iface.createNewMapCanvas3D('North View')
c4 = iface.createNewMapCanvas3D('East View')

# Apply settings to all open 3D map views
for canvas in iface.mapCanvases3D():
	canvas.mapSettings().setEyeDomeLightingEnabled(True)

# Define a camera pose to update the views' cameras
pose = QgsCameraPose()
pose.setCenterPoint(QgsVector3D(0, 210, 0))  # This is in 3D world coordinates
pose.setDistanceFromCenterPoint(100)
pose.setPitchAngle(75)  # Tilt the camera by 75 degrees
pose.setHeadingAngle(0)  # Looking towards North
c1.cameraController().setCameraPose(pose)
pose.setHeadingAngle(90)  # Towards East
c2.cameraController().setCameraPose(pose)
pose.setHeadingAngle(180)  # Towards South
c3.cameraController().setCameraPose(pose)
pose.setHeadingAngle(270)  # Towards West
c4.cameraController().setCameraPose(pose)

# We can set the 3D map views 2D extent to always match the main 2D canvas one
# Our 3D views get updated while zooming/panning the main 2D canvas
canvas = iface.mapCanvas()
canvas.extentsChanged.connect(lambda :c1.mapSettings().setExtent(canvas.extent()))
canvas.extentsChanged.connect(lambda :c2.mapSettings().setExtent(canvas.extent()))
canvas.extentsChanged.connect(lambda :c3.mapSettings().setExtent(canvas.extent()))
canvas.extentsChanged.connect(lambda :c4.mapSettings().setExtent(canvas.extent()))

import remotior_sensus
rs = remotior_sensus.Session(n_processes=2, available_ram=1000)
rs.jupyter().download_interface()

MATCH 
    (t:Trip  {id: "17"})
    <-[:BELONGS_TO]-
    (st:StopTime) 
RETURN st

MATCH 
    (t:Trip {id: "17"})
    <-[:BELONGS_TO]-
    (st:StopTime)
    -[:STOPS_AT]->
    (s:Stop)
RETURN st ,s

from shapely.geometry import Point

QUERY = """MATCH (
    t:Trip {id: "17"})
    <-[:BELONGS_TO]-
    (st:StopTime)
    -[:STOPS_AT]->
    (s:Stop)
RETURN st ,s
ORDER BY st.stopSequence
"""

with driver.session(database="neo4j") as session:
    tx = session.begin_transaction()
    results = tx.run(QUERY)
    df = results.to_df(expand=True)
    gdf = gpd.GeoDataFrame(
        df[['s().prop.name']], crs=4326,
        geometry=df["s().prop.location"].apply(Point)
    )

tx.close() 
m = gdf.explore()
m

QUERY = """
MATCH (t:Trip {id: "17"})
   <-[:BELONGS_TO]-
   (st1:StopTime)
   -[:NEXT_STOP]->
   (st2:StopTime)
MATCH (st1)-[:STOPS_AT]->(s1:Stop)
MATCH (st2)-[:STOPS_AT]->(s2:Stop)
RETURN st1, st2, s1, s2
"""

from shapely.geometry import Point, LineString

def make_line(row):
    s1 = Point(row["s1().prop.location"])
    s2 = Point(row["s2().prop.location"])
    return LineString([s1,s2])

with driver.session(database="neo4j") as session:
    tx = session.begin_transaction()
    results = tx.run(QUERY)
    df = results.to_df(expand=True)
    gdf = gpd.GeoDataFrame(
        df[['s1().prop.name']], crs=4326,
        geometry=df.apply(make_line, axis=1)
    )

tx.close() 
gdf.explore(m=m)

MATCH (t:Trip {id: "17"})
   <-[:BELONGS_TO]-
   (st1:StopTime)
   -[:NEXT_STOP]->
   (st2:StopTime)
MATCH (st1)-[:STOPS_AT]->(s1:Stop)
MATCH (st2)-[:STOPS_AT]->(s2:Stop)
RETURN st1.departureTime AS time1, 
   st2.arrivalTime AS time2, 
   s1.location AS geom1, 
   s2.location AS geom2, 
   duration.inSeconds(
      time(st1.departureTime), 
      time(st2.arrivalTime)
   ).seconds AS traveltime

pip install neo4j

from neo4j import GraphDatabase

URI = "neo4j://localhost"
AUTH = ("neo4j", "password")

with GraphDatabase.driver(URI, auth=AUTH) as driver:
    driver.verify_connectivity()

QUERY = "MATCH (p:Stop) RETURN p.name AS name, p.location AS geom"

records, summary, keys = driver.execute_query(
    QUERY, database_="neo4j",
)

for rec in records:
    print(rec)

import geopandas as gpd
import numpy as np

with driver.session(database="neo4j") as session:
    tx = session.begin_transaction()
    results = tx.run(QUERY)
    df = results.to_df(expand=True)
    df = df[df["geom[].0"]>0]
    gdf = gpd.GeoDataFrame(
        df['name'], crs=4326,
        geometry=gpd.points_from_xy(df['geom[].0'], df['geom[].1']))
    print(gdf)

tx.close() 

LOAD CSV with headers 
FROM "file:///stops.txt" 
AS row 
CREATE (:Stop {
   stop_id: row["stop_id"],
   name: row["stop_name"], 
   location: point({
    longitude: toFloat(row["stop_lon"]),
    latitude: toFloat(row["stop_lat"])
    })
})

load csv with headers from
'file:///agency.txt' as row
create (a:Agency {
   id: row.agency_id, 
   name: row.agency_name, 
   url: row.agency_url, 
   timezone: row.agency_timezone, 
   lang: row.agency_lang
});

load csv with headers from
'file:///routes.txt' as row
match (a:Agency {id: "rigassatiksme"})
create (a)-[:OPERATES]->(r:Route {
   id: row.route_id, 
   shortName: row.route_short_name,
   longName: row.route_long_name, 
   type: toInteger(row.route_type)
});

load csv with headers from
'file:///stops.txt' as row
create (s:Stop {
   id: row.stop_id, 
   name: row.stop_name, 
   location: point({
      latitude: toFloat(row.stop_lat), 
      longitude: toFloat(row.stop_lon)
   }),
   code: row.stop_code
});

load csv with headers from
'file:///trips.txt' as row
match (r:Route {id: row.route_id})
create (r)<-[:USES]-(t:Trip {
   id: row.trip_id, 
   serviceId: row.service_id,
   headSign: row.trip_headsign, 
   direction_id: toInteger(row.direction_id),
   blockId: row.block_id,
   shapeId: row.shape_id
});

:auto
load csv with headers from
'file:///stop_times.txt' as row
CALL { with row
match (t:Trip {id: row.trip_id}), (s:Stop {id: row.stop_id})
create (t)<-[:BELONGS_TO]-(st:StopTime {
   arrivalTime: row.arrival_time, 
   departureTime: row.departure_time,
   stopSequence: toInteger(row.stop_sequence)})-[:STOPS_AT]->(s)
} IN TRANSACTIONS OF 10 ROWS;

call apoc.periodic.iterate('match (t:Trip) return t',
'match (t)<-[:BELONGS_TO]-(st) with st order by st.stopSequence asc
with collect(st) as stops
unwind range(0, size(stops)-2) as i
with stops[i] as curr, stops[i+1] as next
merge (curr)-[:NEXT_STOP]->(next)', {batchmode: "BATCH", parallel:true, parallel:true, batchSize:1});

Communauté	Ressources	L'emploi	Le marché	L'association	Fils RSS
Forum Blogs Wiki Planet Sites amis	Agenda Biblio Fiches AMO-CNIG Paris PCGIS Téléchargements	Offres d'emploi Banque de CV Formations Démarche Métiers	Geo-Entreprises Geo-Communiqué Appels d'offres	L'association Qui sommes-nous ? Nous aider Bienfaiteurs Partenaires Historique Statuts	Tout le site Job Géomatique Services GeoBlogs ...	Faire un don Mentions légales

Le portail francophone de la géomatique

Canaux

Planet OSGeo

Publié: 22 February 2024, 9:00pm CET

Publié: 21 February 2024, 12:30pm CET

Publié: 20 February 2024, 8:17pm CET

Publié: 19 February 2024, 1:00am CET

Pièce jointe: [télécharger]

Publié: 18 February 2024, 1:00am CET

Publié: 15 February 2024, 8:09pm CET

Publié: 15 February 2024, 1:00am CET

Publié: 13 February 2024, 4:19am CET

Publié: 13 February 2024, 3:29am CET

Pièce jointe: [télécharger]

Publié: 8 February 2024, 1:00am CET

Publié: 7 February 2024, 1:55pm CET

Publié: 6 February 2024, 5:00pm CET

Publié: 6 February 2024, 12:03pm CET

Publié: 6 February 2024, 1:00am CET

Publié: 5 February 2024, 9:00am CET

Publié: 5 February 2024, 12:09am CET

Publié: 4 February 2024, 3:42pm CET

Publié: 4 February 2024, 12:16am CET

Publié: 3 February 2024, 10:54am CET

Publié: 1 February 2024, 11:07pm CET

Publié: 1 February 2024, 6:45pm CET

Publié: 1 February 2024, 6:41pm CET

Publié: 1 February 2024, 1:00am CET

Publié: 30 January 2024, 11:23am CET

Publié: 29 January 2024, 12:00pm CET

Publié: 29 January 2024, 12:00pm CET

Publié: 27 January 2024, 2:38pm CET

Publié: 26 January 2024, 3:38pm CET

Publié: 25 January 2024, 2:08pm CET

Publié: 25 January 2024, 7:56am CET

Publié: 24 January 2024, 1:00am CET

Publié: 23 January 2024, 12:25pm CET

Publié: 20 January 2024, 12:46pm CET

Publié: 18 January 2024, 3:04pm CET

Publié: 18 January 2024, 12:24pm CET

Publié: 18 January 2024, 9:15am CET

Publié: 17 January 2024, 6:19pm CET

Publié: 1 January 2024, 12:00am CET

Publié: 31 December 2023, 8:53pm CET

Publié: 31 December 2023, 2:34am CET

Publié: 30 December 2023, 1:38pm CET

Publié: 29 December 2023, 4:58pm CET

Publié: 28 December 2023, 2:36pm CET

Publié: 28 December 2023, 2:35pm CET

Publié: 28 December 2023, 2:25pm CET

Publié: 27 December 2023, 3:29pm CET

Publié: 27 December 2023, 3:00pm CET

Publié: 26 December 2023, 8:34pm CET

Publié: 26 December 2023, 8:33pm CET

Publié: 26 December 2023, 8:32pm CET

Publié: 26 December 2023, 3:20pm CET

Publié: 23 December 2023, 11:00pm CET

Publié: 22 December 2023, 4:29pm CET

Publié: 22 December 2023, 1:30pm CET

Publié: 21 December 2023, 10:52am CET

Publié: 20 December 2023, 8:23pm CET

Publié: 19 December 2023, 5:00pm CET

Publié: 19 December 2023, 9:12am CET

Publié: 18 December 2023, 12:14pm CET

Publié: 18 December 2023, 1:00am CET

Publié: 15 December 2023, 9:59pm CET

Publié: 15 December 2023, 5:00pm CET

Publié: 15 December 2023, 5:00pm CET

Publié: 14 December 2023, 5:00pm CET

Publié: 13 December 2023, 10:45pm CET

Publié: 12 December 2023, 5:00pm CET

Publié: 12 December 2023, 9:35am CET

Publié: 11 December 2023, 5:53pm CET

Publié: 10 December 2023, 3:12pm CET

Publié: 9 December 2023, 10:33pm CET

Publié: 8 December 2023, 5:42pm CET

Publié: 7 December 2023, 7:18pm CET