※ 좌표계 : EPSG:5186
using gdal:ClipRasterByMaskLayer
- 전체 데이터
- 픽셀 단위 확대
using native:PixelsToPoint
- 전체 데이터
- 포인트 단위 확대
import geopandas as gpd dem_point = gpd.read_file('D:/2022 공간빅데이터 경진대회/한반도/11_dem_point(5186).shp')
dem_point
dem_point.crs
- 초기 좌표계 : EPSG:5186
dem_point = dem_point.to_crs(5179)
dem_point = dem_point.reset_index().rename(columns={'index':'id'})
create table dem_point ( id int, value numeric, geometry public.geometry(geometry,5179), -- 'sbd' schema에서 'public.geometry' 타입 사용 uploaded_at timestamp default current_timestamp ); create index index_dem_point_on_geometry on dem_point (geometry);
import psycopg2 from sqlalchemy import create_engine # geometry type 송수신용 라이브러리 conn = psycopg2.connect( 'host=localhost port=5432 dbname=postgres user=postgres password=postgres', options='-c search_path=sbd' ) conn.set_session(autocommit=True) cur = conn.cursor() gconn = create_engine( 'postgresql://postgres:postgres@localhost:5432/postgres' )
cur.execute( 'delete from dem_point' )
dem_point.to_postgis( 'dem_point', gconn, schema='sbd', if_exists='append' )
DEM 자료를 활용해 얻을 수 있는 경사도 데이터 생성 후 업로드
using native:Slope
※ 범례 : 파랑(경사 완만함) ~ 빨강(경사 가파름)
DEM 자료를 활용해 얻을 수 있는 경사 방향 데이터 생성 후 업로드
using native:Aspect
※ 범례 : 파랑(≥0˚) ~ 빨강(<360˚)