汇总重叠线段,测量有效长度
我有一张road_events桌子:create table road_eventsevent_id number road_id number year number from_meas number(10,2) to_meas number(10,2) total_road_length number );insert into road_events (event_id,road_id,year,from_meas,to_meas,total_road_length) values (1,1,2020,25,50,100);insert into road_events (event_id,road_id,year,from_meas,to_meas,total_road_length) values 100);insert into road_events (event_id,road_id,year,from_meas,to_meas,total_road_length) values (3,1,1980,0,25,100);insert into road_events (event_id,road_id,year,from_meas,to_meas,total_road_length) values (4,1,1960,75,100,100);insert into road_events (event_id,road_id,year,from_meas,to_meas,total_road_length) values (insert into road_events (event_id,road_id,year,from_meas,to_meas,total_road_length) values (6、2、2000、10、30、100);insert into road_events (event_id,road_id,year,from_meas,to_meas,total_road_length) values (7、2、1975、30、60、100);insert into road_events (event_id,road_id,year,from_meas,to_meas,total_road_length) values (8、2、1950、50、90、100);insert into road_events (event_id,road_id,year,from_meas,to_meas,total_road_length) values (9、3、2050、40、90、100);insert into road_events (event_id,road_id,year,from_meas,to_meas,total_road_length) values (insert into road_events (event_id,road_id,year,from_meas,to_meas,total_road_length) values ;insert into road_events (event_id,road_id,year,from_meas,to_meas,total_road_length) values ;insert into road_events (event_id,road_id,year,from_meas,to_meas,total_road_length) values (13,5、1985、50、70、300);insert into road_events (event_id,road_id,year,from_meas,to_meas,total_road_length) values (14、5、1985、10、50、300);insert into road_events (event_id,road_id,year,from_meas,to_meas,total_road_length) values (15、5、1965、1、301、300);commit;select * from road_events; EVENT_ID ROAD_ID YEARFROM_MEAS TO_MEAS TOTAL_ROAD_LENGTH---------- ---------- ---------- ---------- ---------- ----------------- 5 5 5 我想选择代表每条路上最新工作的事件。
因为事件 通常只涉及道路的一部分 。这意味着我不能简单地选择每条路的最新事件。我只需要选择最近不重复的 活动里程
。
可能的逻辑(按顺序):
我不想猜测如何解决这个问题,因为它最终可能会带来更大的伤害(类似于
XY问题)。另一方面,它可以提供对问题本质的洞察,所以它来了:
选择每条道路的最新事件。我们将调用最近的事件:event A。
如果event A 是>= total_road_length,这就是我所需要的。算法到此结束。
否则,获得下一个时间顺序事件(event B),其范围与范围不同event A。
如果event B重叠范围event A,则仅得到event B不重叠的部分。
重复步骤3和4,直到事件总长度为= total_road_length。或者,当路上没有更多的活动时,停下来。问题:
我知道这是一项艰巨的任务,但 该怎么办?
这是经典 线性参考 问题。如果我能把线性引用操作作为查询的一部分,会很有帮助。
结果将是:
EVENT_ID ROAD_ID YEARTOTAL_ROAD_LENGTH EVENT_LENGTH---------- ---------- --------------------------- ------------ 5 300
解决方案:
我的主要DBMS是Teradata,但这在Oracle也可以使用。
WITH all_meas AS ( -- get a distinct list of all from/to points SELECT road_id,from_meas AS meas FROM road_events UNION SELECT road_id,to_meas FROM road_events )-- select * from all_meas order by 1,2 ,all_ranges AS ( -- create from/to ranges SELECT road_id,meas AS from_measLead(meas) Over (PARTITION BY road_id ORDER BY meas) AS to_meas FROM all_meas) -- SELECT * from all_ranges order by 1,2,all_event_ranges AS ( -- now match the ranges to the event ranges SELECT ar.*re.event_idre.yearre.total_road_lengthar.to_meas - ar.from_meas AS event_length -- used to filter the latest event as multiple events might cover the same range Row_Number() Over (PARTITION BY ar.road_id,ar.from_meas ORDER BY year DESC) AS rn FROM all_ranges ar JOIN road_events re ON ar.road_id = re.road_id AND ar.from_measre.from_meas WHERE ar.to_meas IS NOT NULL )SELECT event_id,road_id,year,total_road_length,Sum(event_length)FROM all_event_rangesWHERE rn = 1 -- latest year onlyGROUP BY event_id,road_id,year,total_road_lengthORDER BY road_id,year DESC;如果您需要返回实际承保范围from/to_meas(如编辑之前的问题所示),可能会更复杂。第一部分是一样的,但是查询可以在不聚合的情况下返回相同的event_id相邻行(例如,事件3:0-1和1-25):
SELECT * FROM all_event_rangesWHERE rn = 1ORDER BY road_id,from_meas;要合并相邻行,还需要两个步骤(使用标准方法,标记组第一行并计算组号):
WITH all_meas AS ( SELECT road_id,from_meas AS meas FROM road_events UNION SELECT road_id,to_meas FROM road_events )-- select * from all_meas order by 1,2 ,all_ranges ASSELECT road_id,meas AS from_measLead(meas) Over (PARTITION BY road_id ORDER BY meas) AS to_meas FROM all_meas)-- SELECT * from all_ranges order by 1,2,all_event_ranges AS ( SELECT ar.*re.event_idre.yearre.total_road_lengthar.to_meas - ar.from_meas AS event_lengthRow_Number() Over (PARTITION BY ar.road_id,ar.from_meas ORDER BY year DESC) AS rn FROM all_ranges ar JOIN road_eventsre ON ar.road_id = re.road_id AND ar.from_measre.from_meas WHERE ar.to_meas IS NOT NULL )-- SELECT * FROM all_event_ranges WHERE rn = 1 ORDER BY road_id,from_meas,adjacent_events AS( -- assign 1 to the 1st row of an event SELECT t.*CASE WHEN Lag(event_id) Over(PARTITION BY road_id ORDER BY from_meas) = event_id THEN ELSE 1 END AS flag FROM all_event_ranges t WHERE rn = 1 )-- SELECT * FROM adjacent_events ORDER BY road_id,from_meas ,grouped_events AS ( -- assign a groupnumber to adjacent rows using a Cumulative Sum over 0/1 SELECT t.*Sum(flag) Over (PARTITION BY road_id ORDER BY from_meas ROWS Unbounded Preceding) AS grp FROM adjacent_events t)-- SELECT * FROM grouped_events ORDER BYroad_id,from_measSELECT event_id,road_id,year,Min(from_meas),Max(to_meas),total_road_length,Sum(event_length)FROM grouped_eventsGROUP BY event_id,road_id,grp,year,total_road_lengthORDER BY 2,Min(from_meas);
页:
[1]