R에서 포인트 패턴 분석 - 분석 영역 마스킹

이전 기사( R에서 포인트 패턴 분석 )에서는 spatstat 패키지를 이용한 포인트 패턴의 클러스터링 식별과 마크 첨부 점 과정을 소개했습니다.

그 결과, 문문의 유적은 조몬의 유적 근처에 위치하는 것으로 나타났습니다.
이 결과에 약간 불만을 느끼는 것은, 유적 입지 그 자체가 지형적 제약을 받기 때문에, 매크로인 시점에서는 유적의 입지는 집중해 버린다고 생각되는 것입니다. 분석의 전제가 이상하다는 것입니다. 유적이 존재하는 영역에서 분석 영역을 마스킹 그래서 유적이 입지하는 영역을 추출하고 분석 영역을 한정하기로 합니다. 다음 단계로 진행합니다. 유적에서 반경 1000m의 버퍼 다각형 만들기 버퍼 다각형 결합 버퍼 다각형의 범위 내 영역으로 분석 수행 반경 1000m의 버퍼 다각형 우선은 전회 분석한 곳까지 데이터 읽기와 처리를 진행합니다. ### 패키지 로드 library(tidyverse) library(sf) 라이브러리 (GISTools) library(spatstat) library(ggthemes) library(maptools) ### 데이터 로드 map<-st_read(dsn="data/",layer="kokudo_gyousei_utm54") point<-st_read(dsn="data/",layer="site") ### 찰문과 조몬으로 분류한 새로운 필루드(class) 추가 point<-point%>% mutate (class = case_when ( grepl("찰문",period) == TRUE ~ "sat", grepl("조몬",period) == TRUE ~ "jou" ))%>% filter(class=="jou" | class=="sat") %>% dplyr::select(class) point_df<-point%>%as.data.frame() #point를 데이터 프레임으로 변환 coord<-st_coordinates(point)%>%as.data.frame() #point 객체에서 xy 좌표 검색 point_join<-bind_cols(point_df,coord) #point_df와 coord 결합 point_join2<-point_join%>%filter(!grepl("NaN",geometry)) #NaN 제거(왜 2행 있었다) 버퍼 원을 만듭니다. 객체 "point"는 sf 클래스입니다. 조금 길지만 spatstat 패키지로 분석 영역을 정의하는 owin이라는 유형의 객체로 변환하기 전까지의 단계입니다. st_buffer (dist = 중심점으로부터의 거리, endCapStyle = 버퍼 모양) 함수로 반경 1000m의 버퍼 원 생성 st_union() 함수로 모든 버퍼 원을 결합 as("Spatial")로 sp 클래스로 변환 as("owin")로 spatstat 패키지의 owin 객체로 변환 buf<-point%>%st_buffer(dist=1000,endCapStyle = "ROUND")%>% st_union()%>%as("Spatial")%>%as("owin") st_buffer()로 작성한 반경 1000m의 버퍼 원. 스케일이 너무 작기 때문에 점만 볼 수 있습니다.

st_union()으로 버퍼 원을 결합. 이것을 sp 객체 형식을 통해 owin 객체로 변환합니다.

ppp 객체 만들기

ppp() 함수의 인수 window에 owin 객체로 변환한 결합 버퍼 원(buf)을 지정합니다.
이제 buf를 대상 영역으로 하는 ppp 객체가 됩니다.

point_p<-ppp(point_join2$X,point_join2$Y,
    window=buf,
    marks=factor(point_join2$class))

대상 영역이 설정되어 있으므로 밀도도는 대상 영역만 계산됩니다.

point_p%>%density()%>%plot()

엔벨로프 플롯 출력

### cross-L関数によるエンベローププロット
point_p%>%envelope(Lcross,nism=99)%>%plot()

직사각형의 대상 영역을 설정했을 때에 비해 처리 시간이 대폭 증가합니다.
3시간 정도 걸렸습니다.


대상 영역이 좁은 가운데는 조몬과 찰문의 유적 분포는 서로 무관하게 분포하는 것 같습니다만, 대상 영역이 25km정도가 되면 조몬으로부터 멀어진 곳에 찰문 유적이 분포하는 경향이 있는 것 같습니다.

조몬 유적이 밀집하는 이른바 시모카이안(남지부 지구나 혜산 지구)에 찰문 유적이 거의 발견되지 않은 것 등이 영향을 미치고 있을지도 모릅니다.