코드조각 저장소

PostgreSQL + PostGIS 환경에서 래스터 자료에 대한 함수들을 테스트하면서 결과를 정리해 보았습니다. 래스터 자료 구성이 크게는 두가지로 나눌 수 있는것 같습니다. 1.픽셀 정보 : 크기, 위치( 파일, 데이터 셋 메타데이터 기준 Row, Column[X,Y] 2. Band : Value (Multi, Band1, Band2, ...) https://postgis.net/docs/RT_reference.html Chapter 6. Raster Reference 6.18.1.1. Directory containing many files When GDAL opens a file, GDAL eagerly scans the directory of that file to build a catalog ..

이전 글에서 사용하였던 국토정보 플랫폼에서 수집한 정사영상을 활용할 만한 아이템의 첫 단계로 raster 데이터를 데이터베이스에 적하는 작업 부터 진행해 보겠습니다. raster2pgsql은 postgreSQL + PostGIS과 함께 제공되는 유틸리티입니다. 저는 윈도우 환경에서 진행하였지만 Linux환경도 비슷한것 같습니다. rfriend.tistory.com/439 [Greenplum DB] PostGIS에 raster2pgsql 을 사용하여 raster data import 하기 이번 포스팅에서는 PostgreSQL, Greenplum DB의 PostGIS에 raster2pgsql 유틸리티를 사용하여 raster data를 import하는 방법을 소개하겠습니다. 예제로 사용할 데이터는 'Maste..

작성한 다른 글에서 사용하려고 테스트 데이터를 만드는 과정을 해결하여 정리해 봅니다. 사용하는 데이터셋은 비슷합니다만 GDAL만으로 처리하기에는 실력 부족으로 막혔던 부분을 데이터베이스에 raster데이터를 적재하고 각 포인트의 값을 추출하는 방식으로 해결하였습니다. PointCloud의 자료 형식은 많이 있지만 제일 단순하다고 판단된 [ Text : X Y Z R G B ] 형식으로 원천 데이터를 생성하였고 CloudCompare등을 통해 LAS, PLY등으로 변환하여 사용 할 수 있었습니다. 처음 사용한 자료는 국토정보플랫폼에서 제공하는 공개(DEM) 자료를 사용하였습니다 하지만 이자료는 구역별 높이값을 제공하고 있어서 생성된 결과물에 계단 처럼 층이 저 있었습니다. 그래서 표고데이터 및 등고선 데이..

LineString의 양방으로 지정된 길이만큼 연장한 LineString을 구하는 예제 기울기(벡터)에 따른 양방향 연장 지점에 대한 LineString의 생성 SELECT ST_MakeLine( ST_TRANSLATE(st_pt, sin(az1) * len, cos(az1) * len) ,ST_TRANSLATE(en_pt ,sin(az2) * len, cos(az2) * len) ) FROM ( SELECT st_pt , en_pt , ST_Azimuth(st_pt,en_pt) AS az1, ST_Azimuth(en_pt, st_pt) AS az2 , ST_Distance(a,b) + 1 AS len FROM ( SELECT ST_StartPoint(the_geom) AS st_pt, ST_EndPoi..