Eco Life

2023년 첫 포스팅이다. 올해도 많은 분들이 내 블로그를 보고 도움이 되었으면 좋겠다. https://fishingbass.tistory.com/29 [Qgis]웹(web)지도와 벡터파일 좌표계 맞추는 방법 몇일전 어떤 분이 문의를 했던내용이다. 카카오맵에 좌표계가 가지고 있는 벡터파일과 좌표계가 맞지 않아 불어오는 방법에 대한 내용을 문의주셨다. 무료로 공개된 웹 지도를 불러올수 있는 fishingbass.tistory.com 예전 포스팅에서 우리나라 웹 지도를 쉽게 불러올 수 있는 플러그인을 소개했었다. 이번 포스팅에서는 보다 다양한 지도 자료를 이용할 수 있는 플러그인이 있어 소개하고자 한다. 예전 포스팅에서 소개했던 플러그인이다. kakao와 naver, vWorld에서 제공하는 위성 배경지도와..

연구자가 알고 있는 지점, 또는 이미 알고 있는 정보는 내가 작업하고자 하는 공간 내에 몇개의 포인트가 있는지 파악하기 쉽다. 하지만 GBIF 등과 같이 전국 또는 전세계 data와 같이 빅데이터를 이용하는 경우에는 국가별, 지역별로 몇개의 지점이 폴리곤내에 포함되어 있는지 파악할 필요가 있다. 예시로 우리나라 전국에 분포하는 생물자료를 불러왔다고 가정해보자 이 상황에서 도별 출현 지점의 개수를 산출하고자 할 때 사용하면 용이하다. 랜덤포인트 생성은 다음 아래를 참고하면 좋다. https://fishingbass.tistory.com/25 [Qgis]qgis에서 랜덤포인트(Random point) 생성하는 방법, 오류해결 방법 한동안 논문 데이터 정리하느라 많이 바빠서 포스팅을 하지 못했다... 아무튼...

토지이용도 분석 시 buffer를 통해 반경 200m 혹은 연구자들이 판단한 면적내에 도로, 논, 밭 등 어떤 토지를 얼마 만큼의 면적으로 사용하고 있는지 분석할 때에 필요한 분석방법이다. https://fishingbass.tistory.com/22 [Qgis]포인트 벡터에 버퍼(Buffer) 생성하고 면적계산 방법_오류 해결 방법 포인트 벡터를 가지고 buffer를 만드는 방법에 대해서는 많은 분들이 작성해 주셨다. 하지만 내가 필요했던 것은 버퍼 내 하천 또는 기타 등등의 면적을 구하고 싶었다. 여러 방법을 구글링을 통 fishingbass.tistory.com buffer를 생성하고, 면적을 계산하는 방법에 대한 내용은 위의 포스팅 내용과 같다. 토지이용도에서 조사한 정점(point) 반경 200..

8월 중순에 질문 주신 내용에 대한 답을 드리기 위해 작성합니다. 질문 주신 내용은 구글 내 지도에 국내 개발제한구역 레이어 파일을 업로드 하고 싶은데, SHP파일을 KML파일 형식으로 변경 후 업로드를 하고자 하셨습니다. 그러나 구글 내 지도에서는 KML 파일 형식은 5mb 초과 되면 파일을 업로드 할 수 없는데, 방법이 없나요? 라는 질문을 주셨습니다. 파일을 나누어 보는 작업은 질문자 님께서 이미 진행해보셨습니다. 먼저 어떤 방식으로 작업을 하셨는지 모르겠지만 저는 먼저 구글어스pro 프로그램을 이용하여 shp파일을 열어보았습니다. 이후 파일을 나누기 위해 qgis에서 일부 지역의 개발 제한 구역을 나누어 보았습니다. 나누어진 개발제한구역을 KML로 변환하고 용량을 확인해보았습니다. 정말 적은 공간..

안녕하세요. Eco Life 블로그의 주인장 bioteco 입니다. 요새 많은 분들이 메일과 댓글로 문의를 주셔서 질문에 대한 답을 드리는 Q&A 카테고리를 만들었습니다. 여러분께서 문의주시는 답변은 Q&A 카테고리에 게재될 예정이니 참고해주시기 바랍니다. 저도 여러분들의 질문에 답을 드리기 위해 포스팅하는 것이 공부가 되는 것 같습니다. 부담 없이 문의해 주시면 제가 아는 선에서 답변을 드리겠습니다. ps. 개인적으로 출장이 많아 바쁜 경우가 많습니다. 최대한 빠른 답변 드리도록 노력하겠습니다.

R에서 그래프를 작성하고 출력할 때에 해상도가 낮아 논문이나 보고서에 넣기 어려울 때가 있다. 이럴 때 R에서 해상도를 높여주는 작업을 할 수 있다. 각 코드의 기능을 살펴보자. ggsave( filename, plot = last_plot(), device = NULL, path = NULL, scale = 1, width = NA, height = NA, units = c("in", "cm", "mm", "px"), dpi = 300, limitsize = TRUE, bg = NULL) filename: 생성할 파일이름 width, height, units: 이미지 크기(단위: in, cm, mm, px) plot: last_plot()이 기본값 dpi: 저장할 해상도 device: 파일 형식(eps..

오랜만에 Qgis에 관련된 포스팅을 작성하는 것 같다. 개인적으로 이런저런 바쁜일들이 많았다.. 이번 포스팅은 Qgis에서 폴리곤 벡터를 자르는 방법에 관한 내용을 작성하고자 한다. 먼저 벡터를 자르기 위해 폴리곤 벡터를 준비해 준다. (위성영상 사진은 필요 없다.) 기본적으로 폴리곤 편집모트 켜고 끄기는 아실 것이라고 생각한다. 모르시는 분들이라면.. 왼쪽하단 레이어에 수정하고자 하는 벡터를 선택하고 위쪽의 연필 모양의 편집모두 켜고 끄기를 눌러주면 된다. 편집모드가 활성화 되면 아래쪽의 편집도구들이 활성화 되는 것을 알 수 있다. 여러 편집도구들 중에서 피처 분할이라는 탭을 선택하면, 작업창에 보이는 마우스 포인트가 변화하는 것을 확인 할 수 있다. 자르고자 하는 폴리곤의 위치를 파악한다. 필자는 위..

Qgis를 하다보면 누군가 작성해준 point 벡터를 사용하는 경우가 많다. 하지만 내가 추가로 포인트 벡터를 만들어야 하는 상황이 반드시 생기게 될 것이다. https://fishingbass.tistory.com/59?category=434156 [Qgis]포인트(point) 벡터 쉽게 생성하는 방법_google earth pro를 이용하자 gis 작업을 하다보면, 포인트 벡터를 사용할 때가 많다. 포인트를 생성할 수 있는 방법은 굉장히 많이 있다. https://fishingbass.tistory.com/19?category=434156 [Qgis]사진파일 좌표를 이용한 포인트 벡터 만들. fishingbass.tistory.com 지난번 포스팅에서는 google earth에서 저장한 데이터를 qg..

https://fishingbass.tistory.com/14?category=434156 [Qgis]위성영상을 이용한 NDVI(정규식생지수) 분석 NDVI= Normalized Difference Vegetation Index NDVI 분석을 위해서는 위성 영상자료가 필요함. Landsat8외에도 다양한 영상자료가 존재하나 이번 포스팅에서는 Landsat8을 이용한 NDVI 분석을 실시하고자 함. La.. fishingbass.tistory.com 예전에 포스팅 했던 NDVI(정규식생지수) 분석에서 약간의 오류가 발생할 수 있다는 걸 공부하면서 찾았다. 다양한 논문에서 보다 면밀한 분석이 이루어 질 수 있는(?) 방안이 있어 이번 포스팅을 작성하고자 한다. 예전의 포스팅에서 이용했던 위성영상은 Land..

이번 포스팅에서는 R에서 히스토그램(histogram)을 그리는 방법에 대해 이야기하고자 한다. #ggplot2와 readxl이 없는 사람만 수행 install.packages("ggplot2") install.packages("readxl") #ggplot2와 readxl 불러오기 library(ggplot2) library(readxl) 먼저 install.packges()를 통해 ggplot2와 readxl을 다운받아주고 library()를 통해 각 패키지를 불러와 준다. ggplot2의 경우 R을 하시는 분들이라면 많이 들어보시고, 사용하시는 패키지이며, 정말 많은 분들이 소개하고 작성해 놓은 글들이 많다. 자신에게 맞은 스타일을 정리해 놓는 것이 포인트라고 생각한다. #테마를 지정(많은 스타일이..

gis 작업을 하다보면, 포인트 벡터를 사용할 때가 많다. 포인트를 생성할 수 있는 방법은 굉장히 많이 있다. https://fishingbass.tistory.com/19?category=434156 [Qgis]사진파일 좌표를 이용한 포인트 벡터 만들기_importphotos 나도 qgis를 접하기 전에 조사지점의 좌표를 google earth를 이용해서 point를 만들고 저장해 놓곤했다. 여간 귀찮은 일이 아니다... 단순히 몇 지점만 찍는다면 간단할 수 도 있지만 지점이 수십개, 수 fishingbass.tistory.com https://fishingbass.tistory.com/25?category=434156 [Qgis]qgis에서 랜덤포인트(Random point) 생성하는 방법, 오류해결..

gis를 사용하다보면 지도상에 데이터 시각화가 필요할 때가 있다. 다른 여러 방법도 있지만, qgis에서 파이차트를 이용한 데이터 시각화 방법이 있어 소개하고자 한다. 행정경계를 굳이 불러올 필요는 없지만 대부분 행정경계나 기타 작업을 하고나서 시각화 작업을 하기때문에 예시로 행정경계를 불러왔다. 예시로 만들어본 포인트 벡터. 포유류, 양서 파충류, 육상곤충, 어류, 수서곤충 등... 생각나는대로 분류군을 적어보았다.. (숫자는 나름 종수를 생각해서 아무렇게나 입력했다..) 행정경계 위에 포인트벡터가 잘 찍혀있다. 노란색 박스로 표시한 레이어다이어그램 탭을 눌러 들어가준다. 레이어다이어그램 탭을 잘 선택했다면 다음과 같은 화면을 볼 수 있다. 상단의 다이어그램 없음 탭을 눌러준다. 다이어그램 없음 탭을 ..

국가공간정보포털 환경공간정보서비스 등의 포털에서 제공하는 다양한 자료들이 있다. 그중에서 토지이용현황도는 파일 크기도 크고 워낙 많은 정보가 들어있어서 쉽게 다루기가 힘들다. 이번에 원하는 자료를 쉽게 찾는 방법을 찾게 되어 포스팅하게 되었다. 먼저 한반도 내륙은 도별로 자료가 따로 제공되어 있고.. 토지이용현황도를 병합하기에는 너무 오랜시간이 걸리기 때문에 가장 한눈에 쉽게 볼 수 있는 제주도를 예시로 하고자 한다. 토지이용현황도를 불러오면 제주도를 꽉차게 덮을 만큼 많은 정보가 들어있다. 다양한 정보들 중에서 내가 추출하고자 하는 벡터(폴리곤 등)의 정보를 먼저 파악해야 한다. 객체 식별 도구를 이용해서 원하는 폴리곤을 눌러보자. 식별결과를 확인해보면 PNU, JIBUN, BCHK, SGG... 등의..

공간 분석을 하다보면 래스터의 픽셀크기를 맞추어야 하는 경우가 발생한다. 예를 들면 픽셀크기가 작아서 대용량의 래스터파일을 이용하기가 힘든경우(?)에 이용하면 좋을 듯하다. 물론 픽셀크기가 작으면 작을 수록 선명도가 높아지고, 지점에 대한 공간정보가 정확하게 나타나겠지만, 파일의 용량은 점점 더 올라가게 된다. ^^; 행정경계를 잘 보이게 하고 싶어서 폴리곤을 넣었지만, 픽셀을 수정하기 위해서는 래스터파일만 있으면 된다. 예시로 가지고 온 픽셀의 크기이다. 파일마다 픽셀의 크기가 다르기 때문에 위의 사진과 픽셀 크기가 같지 않을 수 있다. 래스터의 픽셀이 다를 경우 저 픽셀 크기를 맞추어 주어야 하는 경우가 있다. 오늘은 그 방법에 대해서 포스팅하고자 한다. 공간처리툴박스에 '재투영' 이라고 검색하면, ..

몇일전 MaxEnt 분석을 진행하면서 subsampling이 필요한 상황이 발생했다. 꼭 이경우가 아니더라도 조사지점의 선정에 도엽이 필요한 경우가 있는데, 이러한 상황에서 그리드(grid)의 생성이 필요한 경우가 종종 있는 것 같다. 먼저 그리드를 생성할 행정경계 자료와 생물이 서식하고 있는 지점들의 좌표가 입력되어 있는 벡터를 불러온다. 공간처리 툴박스에 그리드를 검색해준다. 그리드를 검색하면 벡터생성 탭에 그리드 생성을 선택해준다. 그리드 생성탭에 들어오면 다음과 같은 화면이 생성된다. 근데 여기서 발생하는 문제가 있다. dgree 옆에 경고? 문구 같은게 있어서 확인해보니, 좌표계가 '도'로 입력되면 안된다고 한다. 이번에 지정한 좌표계가 4326이었는데, 이게 문제가 되었던 것 같았다. 바로 벡..

몇일전 어떤 분이 문의를 했던내용이다. 카카오맵에 좌표계가 가지고 있는 벡터파일과 좌표계가 맞지 않아 불어오는 방법에 대한 내용을 문의주셨다. 무료로 공개된 웹 지도를 불러올수 있는 방법이 많이 있지만, 나는 우연히 플러그인에서 쉽게 지도를 불러올 수 있는 방법을 찾아내었다. 플러그인에 korea를 검색하면 다음과 같이 TMS for korea 라고 뜰 것이다. 해당 탭을 누르고 플러그인 설치를 누른다. 플러그인이 오류없이 잘 설치되었다면 웹 이라는 탭을 누르면 사진과 같이 국내에서 제공하는 지도들이 뜬다. 이 방법이 다양한 웹지도를 불러오는데 있어서는 가장 편했던 것 같다. 아마 나에게 문의를 주셨던 분의 문제가 다음과 같은 것이라고 생각한다. 사진에서 보면 알겠지만 벡터파일의 위치가 이상한 것을 알 ..

생태계쪽에서 일하는 사람들, 실험실에 있는 사람들은 알겠지만.... 요즘 정말 바쁘게 출장을 다니고 있다. 핑계아닌 핑계로 바쁜 출장스케줄로 인해서 블로그에 글을 작성하지 못했다. 갑자기 주말에 시간이 나게 되어서 이것저것 하다보니 qgis에서 배경색상을 바꾸는 방법에 대해 갑자기 궁금해졌다. 해결방법은 정말 간단하다. 나도 처음에 기본색상인 하얀색을 가지고 편집하고, 논문, 보고서 등에 이용했다. qgis책을 보던 도중 검은색 배경을 보고나서 "아... 검은색으로 표현할 수도 있구나." 라는 생각이 들어서 이에 대한 방법을 찾아보았다. 윈도우 환경에서도 동일하다. 프로젝트 탭을 누르고, 속성을 찾아서 들어가 준다. 프로젝트 속성의 일반설정에 들어가면 사진에서 보다싶이 배경색상을 지정할 수 있는 탭이 있..

기존에 사용하던 삼성 노트북이 10년도 넘었다... 고등학교 3학년때 EBS를 듣겠다고 부모님께 부탁해서 구매했던 추억이 담긴 노트북... 발열이 너무 심하고, 무엇보다 CPU성능과 32비트로 다양한 통계프로그램을 사용하지 못하는 경우가 많아져서 고민고민 끝에 맥북을 구매하기로 결정하게 되었다. 맥북을 결정하기전에 엄~~청 고민했던 것들을 뽑아보면... 1. 어렵고, 적응하기 힘들것 같다. 2. 윈도우에서 사용하던 프로그램이 호환되지 않을 것이다. 이 2가지가 가장 큰 고민들이었던 것 같다. 결론부터 말하자면... case by case 인것 같다. 나 같은 경우에는 GIS와 R을 많이 사용하기 때문에 사전에 호환이 가능한지 알아봤고 내가 맥으로 오는데에는 큰 문제가 없었던 것 같다. 서론이 길었다.....

2021/01/21 - [공부/Qgis] - [Qgis]qgis에서 랜덤포인트(Random point) 생성하는 방법, 오류해결 방법 [Qgis]qgis에서 랜덤포인트(Random point) 생성하는 방법, 오류해결 방법 한동안 논문 데이터 정리하느라 많이 바빠서 포스팅을 하지 못했다... 아무튼... 오늘은 qgis에서 랜덤포인트를 생성하는 방법에 대해 작성해보고자 한다. 검색창에 random이라고 만 작성해도 벡터 fishingbass.tistory.com 랜덤포인트 생성하는 방법에 대해서 이전 글에서 설명을 했다. 그렇다면, 이들의 좌표를 확인해야할 필요가 있는 상황이 올텐데... 랜덤 포인트를 생성하거나, 임의로 포인트 벡터를 생성하면, 위의 사진처럼 좌표계가 출력되어 있지 않는다. 지금부터 각..

한동안 논문 데이터 정리하느라 많이 바빠서 포스팅을 하지 못했다... 아무튼... 오늘은 qgis에서 랜덤포인트를 생성하는 방법에 대해 작성해보고자 한다. 검색창에 random이라고 만 작성해도 벡터생성 그룹에 라인을 따라 랜덤 포인트 생성 레이어 경계 안에 랜덤포인트 생성 범위 안에 랜덤포인트 생성 폴리곤 내부에 랜덤포인트 생성 이 있는 것을 확인할 수 있다. 이 중에 본인이 필요한 툴을 이용해서 랜덤포인트를 만들 수 있다. 나는 이번 포스팅에서 폴리곤 내부에 랜덤포인트 생성을 해보고자 한다. 작업창에 폴리곤 파일 하나를 불러온다. 예를 들면... 한반도 경계, 혹은 하천경계, 시도 경계... 등등 범위를 나타내는 폴리곤 파일을 올려두면된다. 샘플링 전략은 2가지가 있는데, 나는 포인트 수를 적용했고,..

엑셀자료 또는 CSV 파일을 불러오다 보면 행이름이나 열이름이 원치않게 변경되는 경우가 있다. 예를들어... 엑셀에서 이렇게 정리했다고 가정해보자. 이걸 R에서 불러오면.... ...1이라는 열이름이 새로 생성되고, 행이름을 St.1, St.2 ... 로 하고 싶었지만, 1,2,3,4,5 로 지정되어 있다. 이걸 바꾸기 위해서는 다음의 코드를 이용하면 된다. data = read_xlsx("c:/example.xlsx") name = data[,1] frame = data.frame(data, row.names = name$data...1.) frame = frame[,-1] frame 실행하면.... > data = read_xlsx("c:/example.xlsx") #엑셀자료를 불러오고... New ..