PyQGIS와 함께 동일한 크기의 다각형을 생성합니까?

다음 단계에서 AtlasCreator에 사용할 다각형을 선을 따라 만들고 싶습니다.

ArcMap에는 Strip Map Index Features 라는 도구가 있습니다 .

이 도구를 사용하여 다각형의 높이와 너비 (예 : 8km x 4km)를 선택하고 선을 따라 자동으로 생성 / 회전 할 수 있습니다.

각 다각형의 생성 된 속성 중 하나는 나중에 Atlas Generator에서 북쪽 화살표를 회전해야하는 회전 각도입니다.

누구나 QGIS에서 / pyQGIS 로이 작업을 해결하는 방법을 알고 있습니까? 사용자 정의 플러그인 내에서 사용할 수있는 잔디 또는 SAGA 알고리즘 또는 prossinging-toolbox 모델도 좋습니다.) Edit1 : 인쇄 범위뿐만 아니라 다각형을 사용해야하는 일종의 개요 맵으로 모든 다각형 / 익스텐트.

Edit2 : QGIS와 별도로 소프트웨어를 설치할 필요없이 QGIS 플러그인에서 사용할 수 있는 PyQGIS 솔루션을 찾고 있는데 현상금을 제공하고 있습니다 (PostGIS / Oracle과 같은 RDBMS 없음)

재미있는 질문! 그것은 내가 직접 시도하고 싶었던 일이므로 갈 수 있습니다.

PostGRES / POSTGIS에서 다각형 세트를 생성하는 함수를 사용하여이 작업을 수행 할 수 있습니다.

필자의 경우 철도를 나타내는 하나의 기능 (MULTILINESTRING)이있는 테이블이 있습니다. 미터 단위로 CRS를 사용해야하며 osgb (27700)를 사용하고 있습니다. 4km x 2km '페이지'를 수행했습니다.

여기에서 결과를 볼 수 있습니다 ... 녹색 물건은 도로 네트워크이며, 철도 주변의 1km 버퍼에 클리핑되어 다각형의 높이에 잘 맞습니다.

기능은 다음과 같습니다.

CREATE OR REPLACE FUNCTION getAllPages(wid float, hite float, srid integer, overlap float) RETURNS SETOF geometry AS
$BODY$
DECLARE
    page geometry; -- holds each page as it is generated
    myline geometry; -- holds the line geometry
    startpoint geometry;
    endpoint geometry;
    azimuth float; -- angle of rotation
    curs float := 0.0 ; -- how far along line left edge is
    step float;
    stepnudge float;
    currpoly geometry; -- used to make pages
    currline geometry;
    currangle float;
    numpages float;
BEGIN
    -- drop ST_LineMerge call if using LineString 
    -- replace this with your table.
    SELECT ST_LineMerge(geom) INTO myline from traced_osgb; 
    numpages := ST_Length(myline)/wid;

    step := 1.0/numpages;
    stepnudge := (1.0-overlap) * step; 
    FOR r in 1..cast (numpages as integer)
    LOOP
        -- work out current line segment

        startpoint :=  ST_SetSRID(ST_Line_Interpolate_Point(myline,curs),srid);
        endpoint :=  ST_SetSRID(ST_Line_Interpolate_Point(myline,curs+step),srid);
        currline := ST_SetSRID(ST_MakeLine(startpoint,endpoint),srid);

        -- make a polygon of appropriate size at origin of CRS
        currpoly := ST_SetSRID(ST_Extent(ST_MakeLine(ST_MakePoint(0.0,0.0),ST_MakePoint(wid,hite))),srid);

        -- then nudge downwards so the midline matches the current line segment
        currpoly := ST_Translate(currpoly,0.0,-hite/2.0);

        -- Rotate to match angle
        -- I have absolutely no idea how this bit works. 
        currangle := -ST_Azimuth(startpoint,endpoint) - (PI()/2.0) + PI();
        currpoly := ST_Rotate(currpoly, currangle);

        -- then move to start of current segment
        currpoly := ST_Translate(currpoly,ST_X(startpoint),ST_Y(startpoint));

        page := currpoly;

        RETURN NEXT page as geom; -- yield next result
        curs := curs + stepnudge;
    END LOOP;
    RETURN;
END
$BODY$
LANGUAGE 'plpgsql' ;

이 기능 사용

다음은 예입니다. 4km x 2km 페이지, epsg : 27700 및 10 % 겹침

select st_asEwkt(getallpages) from getAllPages(4000.0, 2000.0, 27700, 0.1);

이것을 실행 한 후 PgAdminIII에서 csv 파일로 내보낼 수 있습니다. 이것을 QGIS로 가져올 수 있지만 레이어에 대해 CRS를 수동으로 설정해야 할 수도 있습니다. QGIS는 EWKT의 SRID를 사용하여 레이어 CRS를 설정하지 않습니다.

베어링 속성 추가

이것은 postgis에서 더 쉬울 것입니다. QGIS 표현식으로 수행 할 수 있지만 코드를 작성해야합니다. 이 같은...

create table pages as (
    select getallpages from getAllPages(4000.0, 2000.0, 27700, 0.1)
);

alter table pages add column bearing float;

update pages set bearing=ST_Azimuth(ST_PointN(getallpages,1),ST_PointN(getallpages,2));

경고

그것은 약간 해킹되었으며 하나의 데이터 세트에서만 테스트 할 수있는 기회를 가졌습니다.

베어링 속성 업데이트에서 어떤 두 정점을 선택해야하는지 100 % 확실하지 않습니다 query. 실험해야 할 수도 있습니다.

왜 현재 선분과 일치하도록 다각형을 회전시키기 위해 그러한 복잡한 공식을 수행해야하는지 모르겠다. ST_Rotate ()에서 ST_Azimuth ()의 출력을 사용할 수 있다고 생각했지만 그렇지 않은 것 같습니다.

다른 솔루션이 있습니다. 그리고 이것은 간단한 폴리 라인과 선택된 여러 엔티티와 함께 ​​작동 할 수 있습니다

블록 다이어그램 :

1. 매개 변수

   1. 생성 방향을 선택하고 인덱스를 읽습니다 (왼쪽에서 오른쪽으로, 북쪽에서 남쪽으로 ...)
   2. 객체 크기 설정

   shape = (4000,8000) # (<width>,<length>)

   3. 중첩 계수 정의 (기본적으로 10 %?)
2. 초기화
   1. 폴리선 정렬 (시작점과 끝점 비교) 순서는 방향 선택에 따라 달라집니다.> 정점 순서 지정 피처 클래스

3. OrderNode에서 루프

   1. 앵커로 첫 번째 포인트를 만들

   2. 각 정점에 대해 dict x, y, id에 추가하고 벡터를 계산하십시오.

   3. 중첩을 줄임으로써 (길이 및 벡터 방향에 따라) 다각형 생성 (10 % / 2)> 5 % 왼쪽 다각형
   4. 선행 정점이 다각형을 벗어나거나 벡터 길이가 모양보다 길 경우 중지
   5. 이전의 좋은 솔루션으로 다각형을 생성하고 마지막으로 좋은 위치로 앵커 포인트 설정
   6. 새로운 루프를 수행하고 dict x, y, id를 재설정하여 다음 다각형 객체를 생성하십시오.

이 제안이 실제로 명확하지 않거나 언급되지 않은 경우 변경할 수 있습니다.

Steven Kays는 pyqgis로 대답합니다. 스크립트를 실행하기 전에 레이어에서 선을 선택하십시오. 스크립트는 라인 병합을 지원하지 않으므로 multilinestring을 가진 레이어에서는 작동하지 않습니다.

#!python
# coding: utf-8

# https://gis.stackexchange.com/questions/173127/generating-equal-sized-polygons-along-line-with-pyqgis
from qgis.core import QgsMapLayerRegistry, QgsGeometry, QgsField, QgsFeature, QgsPoint
from PyQt4.QtCore import QVariant


def getAllPages(layer, width, height, srid, overlap):
    for feature in layer.selectedFeatures():
        geom = feature.geometry()
        if geom.type() <> QGis.Line:
            print "Geometry type should be a LineString"
            return 2
        pages = QgsVectorLayer("Polygon?crs=epsg:"+str(srid), 
                      layer.name()+'_id_'+str(feature.id())+'_pages', 
                      "memory")
        fid = QgsField("fid", QVariant.Int, "int")
        angle = QgsField("angle", QVariant.Double, "double")
        attributes = [fid, angle]
        pages.startEditing()
        pagesProvider = pages.dataProvider()
        pagesProvider.addAttributes(attributes)
        curs = 0
        numpages = geom.length()/(width)
        step = 1.0/numpages
        stepnudge = (1.0-overlap) * step
        pageFeatures = []
        r = 1
        currangle = 0
        while curs <= 1:
            # print 'r =' + str(r)
            # print 'curs = ' + str(curs)
            startpoint =  geom.interpolate(curs*geom.length())
            endpoint = geom.interpolate((curs+step)*geom.length())
            x_start = startpoint.asPoint().x()
            y_start = startpoint.asPoint().y()
            x_end = endpoint.asPoint().x()
            y_end = endpoint.asPoint().y()
            # print 'x_start :' + str(x_start)
            # print 'y_start :' + str(y_start)
            currline = QgsGeometry().fromWkt('LINESTRING({} {}, {} {})'.format(x_start, y_start, x_end, y_end))
            currpoly = QgsGeometry().fromWkt(
                'POLYGON((0 0, 0 {height},{width} {height}, {width} 0, 0 0))'.format(height=height, width=width))
            currpoly.translate(0,-height/2)
            azimuth = startpoint.asPoint().azimuth(endpoint.asPoint())
            currangle = (startpoint.asPoint().azimuth(endpoint.asPoint())+270)%360
            # print 'azimuth :' + str(azimuth)
            # print 'currangle : ' +  str(currangle)

            currpoly.rotate(currangle, QgsPoint(0,0))
            currpoly.translate(x_start, y_start)
            currpoly.asPolygon()
            page = currpoly
            curs = curs + stepnudge
            feat = QgsFeature()
            feat.setAttributes([r, currangle])
            feat.setGeometry(page)
            pageFeatures.append(feat)
            r = r + 1

        pagesProvider.addFeatures(pageFeatures)
        pages.commitChanges()
        QgsMapLayerRegistry.instance().addMapLayer(pages)
    return 0

layer = iface.activeLayer()
getAllPages(layer, 500, 200, 2154, 0.4)

두 개의 답변 (게시 시점)은 독창적이고 잘 설명되어 있습니다. 그러나 이것에 대해 매우 간단하지만 효과적인 해결책이 있습니다 (강을 기반으로 한 임의의 북쪽 방향이 아닌 전통적인 방식으로 모든지도를 북쪽으로 정렬한다고 가정). 회전을 원한다면 가능하지만 약간 더 복잡합니다 (아래 참조).

먼저 여기 내 게시물을 살펴보십시오 . 이를 통해 Atlas에 대한지도 적용 범위를 만드는 방법을 알 수 있습니다. 원하는 방법은 사용 방법의 '워크 플로 2'입니다. 정점 또는 길이로 선형 형상을 분할하고 원하는만큼 형상을 버퍼링합니다. 버퍼링하는 양은 부분적으로 오버랩을 지시하지만 (아래 참조) 더 중요한 것은 영역이있는 피처를 만듭니다. 여러 플러그인을 사용하여 선을 분할 할 수 있지만 GRASS v.split.length 및 v.split.vert가 좋은 옵션입니다 (Processing Toolbox에서 사용 가능).

Map Composer에서 Atlas Generation을 활성화하고 버퍼링 된 레이어를 선택한 후 항목 탭으로 다시 전환하고 맵 객체를 선택하십시오. 'Controlled by Atlas'를 확인하고 사용 사례에서 Margin around 기능을 선택합니다. 이렇게하면지도 간의 겹침을 제어 할 수 있습니다 (또는 고정 배율을 선호 할 수도 있습니다).

작곡가의 상단 툴바에있는 미리보기 아틀라스 버튼을 사용하여 아틀라스를 미리 볼 수 있으며 얼마나 많은 페이지가 생성되는지 확인할 수 있습니다. 모든 페이지를 단일 PDF 또는 별도의 파일로 내보내도록 선택할 수 있습니다.

선을 따라 맵을 회전 시키려면 맵 작성기 항목 속성에 회전 필드가 있습니다. 표현식을 설정해야합니다 (회전 상자 오른쪽에있는 작은 버튼 사용). 변수를 옵션으로 선택한 다음 편집을 선택하십시오. 표현식 빌더가 팝업되고 아틀라스 피쳐의 지오메트리 또는 필드에 액세스 할 수 있습니다. 그런 다음 피쳐의 회전에 따라 맵을 회전하는 급행을 구축 할 수 있습니다 (각 선분의 시작점과 끝점과 약간의 삼각점을 사용하여 베어링을 계산할 수 있음). 동일한 프로세스 또는 반복 계산 된 변수를 사용하여 동일한 프로세스를 반복하여 북쪽 화살표를 회전시킵니다.