다각형 간 흐름 전송 (연결 + 값)

QGIS에는 셀과 하나의 추가 레이어 간 이동 데이터를 나타내는 두 개의 shapefile이 있습니다 (아래 이미지 참조).

다음에 의해 정의 된 데이터 이동 :

• 다각형 "LayerA"(빨간색 윤곽선이있는 투명 사각형). 또한 그것은 "LayerA"지오 센터의 위치에 시각화 된 세포 내 움직임을 나타내는 원과 관련이 있습니다.

  

• 폴리 라인 레이어 "Flows"(노란색 / 회색 화살표), "LayerA"지형의 중심 중심 간의 연결을 통해 값을 전달

  

대상 레이어 :

• 다각형 "LayerB"(진한 회색 윤곽선이있는 연한 라일락 피처).

  

또한 "FLUX"셀 내에서 폴리곤 "LayerA"으로 값 을 이미 전송 하고 이동했습니다 "LayerB". 이전 질문 : QGIS에서 폴리곤간에 상속 된 값? . 그것은 사용하여 수행되었다 %의 $area계산.

전사 의미 용액 / 방법이있을 수 / 송신 / 변형 유동 연결 로 표현 "Flows"의 관계에서 그 값 "LayerA"의 관계로 "LayerB".

폴리 라인으로 이러한 연결을 어떻게 달성 할 수 있습니까?

또한 새 흐름은와 유사한 스타일을 상속합니다 "Flows".

요청에 따라 데이터 샘플을 제공 할 수 있습니다.

흐름의 기능 사이에하지 존재하지 않을 것이다 "LayerA", 그러나 의 기능 사이 "LayerB" . 주요 목표는 가능한 테이블 / 원래 목적지 매트릭스 "FLUX"간의 연결을위한 속성 (즉, 출발 / "LayerB"도착 ) 을 달성하는 것입니다 .

준수해야 할 몇 가지 요구 사항 / 기준 이 있습니다.

1. 동일한 셀에서 피처 부품 (노란색으로 선택) 사이에 흐름 연결이 없습니다.

2. 동일한 피처가 서로 다른 셀에 있어도 동일한 피처간에 연결이 없습니다

3. 기능 이 완전히 두 개의 개별 셀 기능 내에있는 경우 기능의 일부 "LayerB"( "Union"출력 기반) 간에 연결이 존재 "LayerA"합니다.

4."FLUX" 전달되는 새로운 값은 아래 이미지와 같이 계산됩니다.

예를 들어, 두 셀 간의 연결이 I하고 II, "FLUX"된다 100. 다른 값을 가정하면 "NEW_FLUX"사이 A'와 사이에 B''있습니다 1.5625. 100하나의 예일뿐입니다.

참고 문헌 :

• QGIS에서 다각형 사이에 상속 된 값?
• QGIS 2에서 플로우 맵 스타일 조정
• 다른 레이어의 다각형 내에서 다각형의 비례 영역 계산
• 무료 및 오픈 소스 GIS 산란 | GIS의 운동 데이터 : 문제 및 아이디어
• QGIS 2에서 가상 레이어로 지오메트리 유지

이론적으로 가상 레이어를 사용하면 가능합니다 (셰이프 파일을 사용하면 프로세스가 훨씬 오래 걸리지 만 레이어가 공간 데이터베이스에 있으면 훨씬 빠릅니다).

여기 코드가 있습니다 :

WITH inter_ab AS ( 
--create intersection between LayerA and LayerB 
SELECT LayerA.id || '_' || LayerB.FLAECHEID AS id, 
LayerA.id AS id_a, 
ST_AREA(LayerA.geometry) AS area_a, 
LayerB.FLAECHEID AS id_b, 
ST_INTERSECTION(LayerB.geometry, LayerA.geometry) AS geom 
FROM LayerA, LayerB 
WHERE ST_INTERSECTION(layerB.geometry, layerA.geometry) IS NOT NULL 
),

--calculation of the new flux value 
new_flux AS (SELECT t1.id_b AS origine, 
t2.id_b AS dest, 
SUM(Flows.flux * ST_AREA(t1.geom) / t1.area_a * ST_AREA(t2.geom) / t2.area_a) AS value  
FROM inter_ab t1, inter_ab t2, flows 
-- no connection between the same feature 
WHERE t1.id <> t2.id 
-- rule 1 
AND t1.id_a <> t2.id_a 
-- rule 2 
AND t1.id_b <> t2.id_b 
-- get flow data 
AND flows.origine = t1.id_a 
AND flows.dest = t2.id_a 
GROUP BY t1.id_b, t2.id_b
)

--create flows between original layerB features
SELECT new_flux.origine,
new_flux.dest,
new_flux.value AS flux,
make_line(ST_CENTROID(t3.geometry), ST_CENTROID(t4.geometry)) AS geom --ST_MakeLine under postGIS
FROM LayerB t3,
LayerB t4,
new_flux
WHERE t3.FLAECHEID = new_flux.origine
AND t4.FLAECHEID = new_flux.dest

그래픽 출력은 다음과 같습니다

결과는 수동으로 테스트되었습니다. "FLUX"값 의 차이 는 무시할 수 있습니다.

최종 출력은 "Flow"다음과 같은 스타일을 상속받습니다.

몇 가지 데이터로 테스트하는 것이 좋으며 큰 데이터 세트에 너무 오래 걸리면 쿼리를 단계별로 실행하고 ( "inter_ab", "new_flux") 결과를 저장하고 다음 쿼리를 실행하십시오.

세 계층간에 조인을 수행 한 다음로 집계 할 수 layerB있습니다. 가상 레이어를 사용할 수 있습니다. 중요한 데이터가 레이어 에 있는지 layerA또는 flow레이어 에 있는지 확실하지 않습니다 . 여기에 (예상치 않은) 가능성이 있습니다.

SELECT b.id, b.geometry, sum(a.myVar)
FROM layerB b
LEFT JOIN flow f
   ON ST_Intersects(ST_EndPoint(f.geometry),b.geometry)
 JOIN layerA a
   ON ST_Intersects(ST_StartPoint(f.geometry),a.geometry)
GROUP BY b.id