ST_Affine을 사용하여 PostGIS에서이 작업을 수행 할 수 있습니다 . 임의의 점을 중심으로 회전하는 기능이 ST_Rotate 에 추가되었습니다 . PostGIS 2.0 용 에 .

PostGIS 1.5 또는 그 이전 버전과 같은 이전 버전이있는 경우 다음 기능을 추가 할 수 있습니다.

CREATE OR REPLACE FUNCTION st_rotate(geometry, double precision, geometry)
  RETURNS geometry AS
'SELECT ST_Affine($1,  cos($2), -sin($2), 0,  sin($2),  cos($2), 0, 0, 0, 1, ST_X($3) - cos($2) * ST_X($3) + sin($2) * ST_Y($3), ST_Y($3) - sin($2) * ST_X($3) - cos($2) * ST_Y($3), 0)'
  LANGUAGE sql IMMUTABLE STRICT
  COST 100;
COMMENT ON FUNCTION st_rotate(geometry, double precision, geometry) IS 'args: geomA, rotRadians, pointOrigin - Rotate a geometry rotRadians counter-clockwise about an origin.';

CREATE OR REPLACE FUNCTION st_rotate(geometry, double precision, double precision, double precision)
  RETURNS geometry AS
'SELECT ST_Affine($1,  cos($2), -sin($2), 0,  sin($2),  cos($2), 0, 0, 0, 1,    $3 - cos($2) * $3 + sin($2) * $4, $4 - sin($2) * $3 - cos($2) * $4, 0)'
  LANGUAGE sql IMMUTABLE STRICT
  COST 100;
COMMENT ON FUNCTION st_rotate(geometry, double precision, double precision, double precision) IS 'args: geomA, rotRadians, x0, y0 - Rotate a geometry rotRadians counter-clockwise about an origin.';

ST_Rotate 에서 예제 를 사용하여 x , y 주위의 지오메트리를 회전하는 데 사용하는 방법에 대한 아이디어를 얻으십시오.중심점 (공통 중심)을 포함하여 지점 .

우리는 모두 수학을 좋아하기 때문에 위 함수의 변환 행렬은 다음과 같이 표현됩니다.

[ cos(θ)  | -sin(θ)  ||  x0 - cos(θ) * x0 + sin(θ) * y0 ]
[ sin(θ)  |  cos(θ)  ||  y0 - sin(θ) * x0 - cos(θ) * y0 ]

여기서 θ 는 원점에 대한 시계 반대 방향 회전입니다. x0 은 의 동쪽 / 경도, y0 은 북쪽 / 위도입니다. 이 수학은 모든 아핀 변환 도구에 적용될 수 있습니다 .

qgsAffine 도구를 사용하려면 매트릭스 값이 흐르는 위치를 이해해야합니다. 사전 계산을 수행하려면 올바른 스프레드 시트 템플릿이 필요합니다. qgsAffine 대화 상자는 다음과 같습니다.

              X   Y
            +---+---+
      Scale | a | e |
            +---+---+
   Rotation | d | b |
            +---+---+
Translation | c | f |
            +---+---+

어디:

• a : cos (θ)
• b : -sin (θ)
• 씨 : x0-cos (θ) * x0 + sin (θ) * y0
• 디 : 죄 (θ)
• 이자형 : cos (θ)
• 에프 : y0-sin (θ) * x0-cos (θ) * y0

예를 들어, 다각형을 시계 방향으로 42 ° S, 174 ° E 주위로 30 ° 회전하려면 다음은 스프레드 시트에 입력 된 내용입니다.

• x0 = 174
• y0 = -42
• θ = -30도 또는 -0.523598776 라디안

그런 다음 스프레드 시트의 결과를 오른쪽 상자에 복사 / 붙여 넣기하십시오. 대화 상자에서의 탭 순서 사용 :

• a : 0.866025404
• d : -0.5
• c : 44.31157974
• e : 0.866025404
• b : 0.5
• f : 81.37306696

PostGIS의 동일한 예는 다음과 같습니다.

SELECT ST_Rotate(geom, -30*pi()/180, 174.0, -42.0)

답변