Chuyển luồng (kết nối + giá trị) giữa các đa giác

Trong QGIS có hai shapefile biểu thị dữ liệu di chuyển giữa các ô và một lớp bổ sung, xem hình ảnh bên dưới

Di chuyển dữ liệu được xác định bởi:

• Đa giác "LayerA"(hình vuông trong suốt có viền đỏ). Bên cạnh đó, nó cũng liên quan đến các vòng tròn đại diện cho các chuyển động trong các tế bào, được hình dung về vị trí của "LayerA"geocroids.

  

• Lớp polyline "Flows"(mũi tên màu vàng / xám), truyền tải các giá trị thông qua các kết nối giữa các "LayerA"tính năng địa tâm

  

Lớp mục tiêu:

• Đa giác "LayerB"(tính năng màu hoa cà nhạt với viền màu xám đậm).

  

Ngoài ra, tôi đã chuyển "FLUX"và chuyển các giá trị trong các ô từ "LayerA"thành "LayerB"đa giác, xem câu hỏi trước của tôi: Các giá trị được kế thừa giữa các đa giác trong QGIS? . Nó được thực hiện bằng cách sử dụng %các $areatính toán.

Có thể có một giải pháp / cách tiếp cận có ý nghĩa về chuyển / truyền / chuyển đổi các kết nối dòng được đại diện bởi "Flows"và các giá trị của nó từ quan hệ "LayerA"thành quan hệ "LayerB".

Làm thế nào tôi có thể đạt được những kết nối như polylines?

Ngoài ra, các luồng mới sẽ kế thừa một kiểu tương tự "Flows".

Theo yêu cầu, tôi có thể cung cấp một mẫu dữ liệu.

Dòng chảy sẽ tồn tại không phải giữa các tính năng của "LayerA", mà là giữa các tính năng của "LayerB" . Mục đích chính là để đạt được thuộc tính "FLUX"(tức là từ / đến) cho các kết nối giữa "LayerB"có thể như ma trận bảng / Xuất xứ-Điểm đến.

Có một số yêu cầu / tiêu chí cần được tôn trọng:

1. Không có kết nối luồng giữa các phần của tính năng (được chọn màu vàng) trong cùng một ô

2. Không có kết nối giữa cùng một tính năng, ngay cả các bộ phận của nó nằm trong các ô khác nhau

3. Kết nối tồn tại giữa các phần của tính năng "LayerB"(dựa trên "Union"đầu ra) nếu chúng hoàn toàn nằm trong hai "LayerA"tính năng ô riêng biệt

4. Giá"FLUX" trị mới đang truyền tải, sẽ được tính như thể hiện trên hình ảnh bên dưới.

Ví dụ, có một mối liên hệ giữa hai tế bào Ivà II, nơi "FLUX"là 100. Giả sử các giá trị khác, "NEW_FLUX"giữa A'và B''sẽ xung quanh 1.5625. 100chỉ là một ví dụ duy nhất.

Người giới thiệu:

• Các giá trị được kế thừa giữa các đa giác trong QGIS?
• Điều chỉnh các kiểu bản đồ luồng trong QGIS 2
• Tính diện tích tỷ lệ của đa giác trong đa giác của lớp khác
• Mã nguồn mở miễn phí và mã nguồn mở | Dữ liệu di chuyển trong GIS: các vấn đề & ý tưởng
• Bảo tồn hình học với lớp ảo trong QGIS 2

Về mặt lý thuyết, với các Lớp ảo, điều đó là có thể (với các shapefile, quá trình này sẽ kéo dài thêm, nhưng nếu các lớp nằm trong Cơ sở dữ liệu Không gian, tôi nghĩ rằng nó nhanh hơn rất nhiều).

Đây là mã:

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

Đầu ra đồ họa sẽ trông như

Kết quả đã được kiểm tra bằng tay. Sự khác biệt về "FLUX"giá trị là không thể bỏ qua.

Đầu ra cuối cùng sẽ kế thừa các kiểu từ "Flow"và trông giống như

Tôi khuyên bạn nên kiểm tra nó với một vài dữ liệu và nếu mất quá nhiều thời gian cho các tập dữ liệu lớn, hãy thực hiện từng bước các truy vấn ( "inter_ab", "new_flux") và lưu kết quả và thực hiện truy vấn tiếp theo.

Bạn có thể thực hiện nối giữa ba lớp, sau đó tổng hợp theo layerB. Các lớp ảo có thể được sử dụng. Tôi không chắc liệu dữ liệu quan trọng nằm trong layerAhay trong flowlớp .. Đây là một khả năng (chưa được kiểm tra):

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