GIS

개요이 포스팅의 원문이 되는 글은 사실 3년 전에 쓰여졌습니다. 당시에는 원래 쓰던 언어도 바꾸고(C++ -> Python) 일하는 환경도 바꾸고 다루는 데이터도 바뀌는 등 큰 변화를 겪고 있었어요. 근 1년 간은 Raster를 집중적으로 다뤘습니다만 원래는 Vector를 중점적으로 다뤘습니다. 그 때의 추억도 소환할 겸, 복기도 하고, 다시 공부하는 느낌으로 포스팅을 쓰게 되었습니다.오늘 포스팅은 좀 어려운 내용들을 다룹니다. 여러분들은 포스팅을 보면서 아래의 표정을 짓게 될 겁니다.네, 이해합니다. 석사 과정 당시에 공부한 내용들이 이 포스팅을 통해 빛을 발하는군요.(+) 2024.08.05 기준 원문을 기반으로 작성한 글입니다. 추후 python 코드로 좀 더 쉽게 이해할 수 있게 예시를 덧붙이겠습..

GIS 엔지니어로 수많은 좌표계를 다루지만, 보통은 Mercator projection이 기반인 좌표계를 주로 다룹니다.그런데 말입니다,어느날 혜성과도 같이 이 GEOS projection이라는 녀석이 절 찾아왔습니다.악몽의 시작이었죠.1년은 시달린 것 같습니다.얘는 이전에 썼던 글인 GeoTransform과 GeoReference에서georeference를 처리하는 방법에 추가로 필요한 요소가 있습니다.그리고 얼마 전에 제가 드디어 방법을 찾았습니다.오늘 글의 숨은 주제는 정답이 없는 문제를 푸는 과정이기도 합니다.작은 실마리 찾기다른 사례를 보면서 역엔지니어링하기메뉴얼이 부족한 경우 배경 지식을 활용해서 값들을 계산해서 구하기역엔지니어링한 경우와 내 해결방법을 교차 검증하기개요포스팅의 목적GEOS p..

GIS 데이터를 다루다 보면 이 두 단어는 상당히 많이 보실 겁니다.특히 제일 머리 아플 때는 이미지 좌표계와 공간 좌표계 간 원점 위치의 차이일 겁니다. 오늘은 이 헷갈리는 점을 짚고 넘어가려 합니다.개요포스팅의 목적gdal의 geotransform의 개념과 구성 요소를 알아보자gdal에서는 어던 좌표계 체계를 채택하고 있는지 알아보자이미지를 georeference하려면 어떻게 해야할까? 이미지 좌표계와 공간 좌표계 간의 차이에 주의하자.개념이미지 좌표계? 공간 좌표계?이미지 좌표계위 그림은 우리가 알고리즘 문제를 풀 때 많이 접하는 배열의 형태입니다. 가로(width)를 column, 세로(height)를 row로 대응합니다.원점은 좌상단부터 시작합니다.이미지 좌표계에서도 위와 같은 형태로 배열을 취..

python으로 공간 데이터 시각화를 하는 방법은 두 가지가 있습니다.matplotlibplotly과학자들은 matplotlib으로 그래프를 그리는 것이 익숙합니다.공간 데이터 시각화에도 matplotlib 기반으로 많이 접근합니다.plotly는 지도를 그려놓고 그 위에서 공간 데이터를 시각화합니다.오늘 할 이야기는 matplotlib을 기반으로 공간 데이터를 시각화하는 얘기, 그 중에서도 Basemap을 사용하는 방법에 대해 알아보고자 합니다.왜 Basemap이냐?matplotlib기반 공간 정보 시각화 방법에는 cartopy와 basemap이 있습니다.이때까지는 cartopy를 많이 썼습니다.제가 matplotlib 기반으로 그래프는 많이 그려봤지만 cartopy는 그 명성은 익히 들었지만 거의 사용..

목적 배경 공간 정보 데이터는 맞는데, 그렇다고 이를 바로 지도에 표출하려니 데이터가 다른 분야에서 다루는 데이터라서 공간 정보 데이터의 형식에 맞춰줘야 할때 기상 데이터들은 공간에서 일어나는 현상을 다룬 데이터라 위치 정보들이 있지만, 우리가 생각하는 EPSG:4326 혹은 EPSG:3857 같은 GIS 계에서 흔히 쓰는 좌표계를 쓰지 않는다 lambert conformal coordinate 각 픽셀마다 좌표 정보가 있다지만, 다음의 두 가지 형식으로 저장된다 픽셀 마다 좌표 정보가 있다. 각 픽셀들은 lcc 좌표계 내에서 등간격으로 존재하고 표기는 위경도로 표기된다 상대적인 범위 내에 격자 형식으로 좌표를 기술할 수도 있다. 아래 그림은 기상청에서 쓰는 격자 좌표계를 기술한 문서의 일부 기상 데이터..

개요 해당 문서의 대상 독자(이하 서술 조건의 대부분에 속하신다면 대상이십니다) 급하게 ArcGIS로 데이터를 처리해야 함 javascript보다는 python이 편함 ArcGIS에 데이터는 올렸는데 편집이 필요하다. 근데 나 맥 써서...ArcGIS Pro로 대응을 할 수 없다. conda에서 offline 설치 가능하게 환경 설정 > 명령어 공식 문서에는 The conda utility will pull all the arcgis package dependencies from the installed set of Anaconda libraries instead of searching the internet. 라고 되어있지만 cli로 테스트해보니까 제대로 import 되지 않아 포기 docker ins..

개요 qgis 3.22버전에서 새롭게 annotation layer라는 layer가 새로 생겼다. qgis에서 annotation layer가 포함된 qgs를 열면 다음과 같이 보인다. raster와 qgs가 한 쌍으로 움직인다. 이 layer는 벡터 데이터가 아니며, qgis project file(qgs)에 xml로 도형을 기록한다 어떻게 annotation의 geometry를 추출할 수 있는가? 참고 자료 qgis annotation layer란? 3.22에는 annotation layer을 대응할 수 있는 tool bar가 있는 것으로 보이나, 내가 쓰는 QGIS 3.26 버전에서는 찾을 수 없었다. (수정 : 22.11.18) macos QGIS 3.28 버전 기준, View > Toolbar ..

전제 및 가정 지오레퍼런싱(georeference)을 하기 위해서는 다음 세 가지 정보가 꼭 필요하다. 좌표계 정보(EPSG 코드가 명시되어 있으면 가장 좋다.) 한국에서 쓰는 좌표계는 몇 가지가 있는데, 중부 원점(5186)을 가장 많이 쓰고 조선소의 경우 동부 원점으로 기술된 경우가 있다. 한국 내에서 영상의 지리 정보를 기술할 때 표준과 상관없이 기술하는 경우가 대다수이므로, 이럴 때는 다음을 확인해 보자. 한국에서 정의된 좌표계는 좌표계의 원점을 (0,0)으로 뒀을 때, 이를 좌표계 상에서 동쪽으로 얼마, 북쪽으로 얼마 이동한 가상 원점을 사용한다. 예를 들어 동쪽으로 200,000m 북쪽으로 600,000m 이동했다고 적혀있으면 epsg:5186에 해당하는 중부 원점이다. 좌표계에서의 영상의 좌..

GIS에서 흔하게 볼 수 있는 데이터 포멧은 shp, kml, geojson이 있는데 구글어스에 뭔가를 올려서 보고싶다면 kml을 피할 수 없다. xml 문서답게 kml은 그 트리 구조부터 attribute까지 사람을 헷갈리게 하는 경우가 많다. 그래서 이번에 특정 조건 하에 kml 파일을 대량 생산하면서 배운 것들을 여기에 정리하려 한다. 1. 쓰기 이전에 코드를 받아 써서 simplekml으로 파일을 썼는데 주의할 점이 있다. 루프당 kml을 생성한다면 반드시 kml 문서 instance를 초기화했는지 확인하자. 그렇지 않다면 나처럼 눈덩이처럼 불어버린 kml을 마주할 수 있게 된다. ''' - key(str) : 문서 식별자 - input_df(dataframe) : kml로 ..

PostGIS에 저장되어있는 Polygon에 대해 가장 긴 변의 길이를 구하고 싶다는 요청이 들어왔다. 이에 구글링하면서 레퍼런스를 찾았고 그 결과에 대해 좀 더 덧붙여서 기록하려고 한다. 원본 레퍼런스 : https://stackoverflow.com/questions/7595635/how-to-convert-polygon-data-into-line-segments-using-postgis SELECT ST_AsText( ST_MakeLine(sp,ep) ) FROM -- extract the endpoints for every 2-point line segment for each linestring -- call two points as start point and end point from each..