Google Street View와 딥러닝을 활용한 서울시 녹지 형평성 분석: NDVI와 가로이미지 기반 녹지 산출방법과의 비교를 중심으로

1. 연구의 배경 및 목적

도시 내 녹지는 거주민들의 보행 및 신체활동, 신체 및 정신건강을 증진시키고 열섬 현상을 완화하는 등의 긍정적 효과를 유발한다(Markevych et al., 2017; Thompson et al., 2012; Zhou et al., 2014). 구체적으로 녹지는 보행 및 신체활동을 위한 매력적이고 안전한 장소의 역할을 수행한다(Markevych et al., 2017). 따라서 녹지는 도시민 모두가 공평하게 누려야하는 공공 서비스로 인식되어지고 있다. 하지만 최근 도시공간에서 녹지에 대한 수요가 증가하면서 녹지의 형평성 문제가 지적되고 있다. 즉, 도시 내 녹지가 공간적으로 불균등하게 분포되어 있을 뿐만 아니라, 저소득, 노인 등의 취약계층 주거지와의 관계에 있어서도 불평등하게 분포되어 있다는 것이다(Nesbitt et al., 2019; Oh et al., 2020; Rigolon, 2017; Schwarz et al., 2018).

서울시도 이와 같은 녹지 불평등 문제가 제기되고 있다. 보행권 내에 녹지가 존재하지 않는 지역이 약 27.27km²이며, 이러한 소외지역은 자치구별로 큰 차이를 보이고 있다(서울시, 2015). 2030 서울플랜에서는 이러한 녹지 불평등 문제를 해결하기 위해 녹지 이용의 형평성을 강조하며, 녹지 서비스의 소외지역 해소를 위한 신규 녹지 확충 목표를 수립하였다. 위와 같이 녹지 형평성에 대한 관심이 증대됨에 따라, 서울시의 녹지 소외지역 및 녹지의 불평등 정도를 정확하게 파악하는 것이 중요해질 것으로 판단된다.

이러한 관심으로 인해 서울시를 대상으로 녹지 형평성 분석 및 취약지를 도출해내려는 연구가 이루어지고 있다(고영주 외, 2019; 김미현 외, 2015; 김용국, 2015; Oh et al., 2020). 이들 연구는 주로 위성사진, 토지피복(land cover), 토지이용(land use) 데이터를 활용한 공원 및 녹지면적 변수를 대상으로 근린 수준에서의 형평성 정도를 분석하고 있다. 그러나 위와 같은 기존 녹지 산출방식은 다음과 같은 한계점을 가진다. 첫째, 기존 녹지 산출 방식은 산지 등의 대규모 녹지를 주로 포함함으로써, 도시민들이 일상생활에서 체감할 수 있는 녹지의 양을 과대 산출할 수 있다(Ye et al., 2018). 더욱이 서울시의 경우, 분지지형으로 도시 외곽지역에 다량의 녹지가 분포해 있다. 즉, 사회적 취약계층들이 주로 도시 외곽지역에 거주하는 점을 참고하였을 때, 기존 녹지 산출 방식으로 녹지 형평성을 분석할 시 형평성 정도가 양호하게 산출될 수 있다.

둘째, 전통적 녹지 산출방식은 하늘에서 지상을 바라보며(top-down view) 녹지를 산출하는 방식으로, 평면적이고 공급자 중심적 지표라는 한계점이 존재한다(조용현 외, 2010; Ye et al., 2018). 이러한 한계점으로 기존 녹지 변수는 실제 도시민들이 체감하는 녹지의 양을 제대로 반영하지 못할 것으로 판단된다. 따라서 이러한 한계점을 지닌 기존 지표로 녹지에 대한 불평등 정도 및 취약지를 분석할 시, 실제 거주민들이 체감하는 정도와 괴리가 존재할 수 있다.

한편, 평면적인 기존 녹지 지표를 보완할 수 있는 지표로 일정 지점에 서 있는 사람의 시야 내에서 식생이 보이는 정도를 의미하는 녹시율 개념이 존재한다(조용현 외, 2010; Yang et al., 2009). 이는 시각적 녹지 양을 의미하는 것으로 녹지의 활력 및 입체적인 녹지를 고려함으로써 도시민들이 체감할 수 있는 녹지 양을 대변할 수 있는 변수라고 알려져 있다(Lu et al., 2019; Yang et al., 2009; Ye et al., 2018). 녹시율을 측정하기 위한 기존 방법으로는 현장조사가 있지만, 이는 시간, 금전적 제약이 크다는 한계점이 존재한다. 이러한 한계점을 극복하기 위해 최근 국외 연구를 중심으로 가로경관 빅데이터와 컴퓨터 비전을 이용하여 보행자 관점(eye level)에서의 녹지를 비롯한 근린환경을 측정하려는 시도가 다수 이루어지고 있다(Gong et al., 2018; Ki and Lee, 2021; Nguyen et al., 2019; Rzotkiewicz et al., 2018).

따라서 본 연구의 목적은 다음과 같다. 첫째, 본 연구에서는 서울시를 대상으로 행정동별 기존 녹지 산출방법으로 산출한 녹지 양과 새로운 산출방법인 Google Street View(이하 GSV) 및 컴퓨터 비전을 이용한 가로녹시율을 산출한다. 기존 녹지 산출 방법은 선행연구에서 주로 사용한 Normalized Difference Vegetation Index(NDVI), 토지피복 데이터를 통해 산출한 공원 및 녹지면적 두 가지를 활용한다. 이후 서울시 행정동을 대상으로 산출한 세 가지 녹지변수의 분포를 확인 및 비교한다. 둘째, 본 연구에서는 이들 녹지 변수와 취약계층을 대변할 수 있는 행정동별 인구사회경제 속성과의 공간 회귀분석을 실시한다. 이는 서울시 내 분포되어 있는 녹지와 취약계층 거주지와의 관계를 분석함으로써 각 녹지변수별 불평등 정도를 확인 및 비교하고자 함이다.

1. 녹지 형평성에 관한 연구

도시 내 녹지는 시민들의 건강과 삶의 질에 필수적인 요소이지만, 녹지 분포의 불균형, 불평등성에 대한 문제가 지속적으로 제기되고 있다(김미현 외, 2015). 이는 환경 불평등(environmental inequality) 측면에서 환경 편익, 환경 비용 등이 사회적 소외계층들에게 불리하게 작용한다는 사실에 기반한 것으로(명수정, 2017), 형평적 배분이란 단순히 공간적인 배분의 의미를 넘어, 이를 이용하는 대상의 성별, 연령, 인종, 소득 수준 등 사회경제적 특징에 차별을 받지 않는 것을 의미한다(신지영, 2009). 이러한 환경 형평성 관점에서 도시 내 녹지는 입지의 효율성뿐만 아니라 거주자의 사회경제적 측면에서 형평성이 고려되어야 한다. 또한, 녹지에 대한 불평등은 단순히 도시 환경문제를 넘어서 주거환경의 불평등, 나아가 경제·사회적 불평등의 문제와도 밀접하게 긴밀하게 연관된 중요한 요소라 할 수 있다(고영주 외, 2019).

선행연구들은 다수의 사례지역을 대상으로 녹지 형평성을 분석하였으며, 분석결과로 사회적 소외계층들이 주로 거주하는 지역에 녹지의 양이 적거나, 접근성이 좋지 않아 녹지에 대한 불평등이 존재한다는 것을 실증하고 있다(고영주 외, 2019; 김미현 외, 2015; Nesbitt et al., 2019; Oh et al., 2020; Rigolon, 2017; Schwarz et al., 2018).

우선, 김미현 외(2015)는 서울시를 대상으로 공원을 위계별로 구분하여, 위계별 공원 유치권 기준을 통해 공원의 접근권역과 사회적 소외계층인 노인, 기초생활 수급자, 장애인 등의 분포를 비교하였다. 분석결과, 서울시 내에서 지역별로 공원 접근권역에서의 차이가 크게 나타난다고 보고하고 있으며, 특히 강서구, 영등포구, 용산구 등에서 공원 접근도가 낮은 지역이 밀집해있다고 보고하고 있다.

이와 유사한 연구로 서울시를 대상으로 진행한 Oh et al. (2020) 연구 또한 공원을 위계별로 세분화하여 공원 접근 가능 면적을 행정동별로 분석하였고, 다양한 사회적 소외계층의 분포와 비교를 하였다. 분석결과, 장애인, 노인, 어린이, 기초생활 수급자 계층은 공원의 위계별 종류에 따라 형평성 정도가 다르게 나타났으나, 외국인 계층의 경우 모든 위계의 공원에서 형평성 정도가 낮게 나타났다.

고영주 외(2019)의 연구는 서울시 중구, 성동구, 동대문구를 대상으로 지표 온도, 녹지 면적 등의 다양한 환경 형평성 지표를 활용하였고, 녹지 면적에서는 이전 연구와는 다르게 위성사진을 활용한 NDVI, 가로수 숲 지붕 비율(tree canopy ratio)을 변수로 활용하였다. 연구의 결과로 소외계층들이 주로 거주하는 지역에 NDVI는 낮게, 지표온도는 높게 나타나, 이들 자치구의 환경 형평성이 낮은 것으로 분석되었다. 추가적으로 주택유형 별로 녹지 면적의 차이가 큰 것으로 분석되었는데, 구체적으로 아파트 단지의 녹지 양이 단독 및 다가구 밀집지역에 비해 약 2배 이상 높게 나타났다. 이는 상대적으로 저소득층들이 주로 거주한다고 볼 수 있는 다가구 주택에 녹지가 적은 것으로 나타나, 녹지 측면에서 환경 불평등이 존재한다고 판단될 수 있다.

해외에서도 녹지 형평성에 대한 연구가 활발하게 이루어졌다. 대표적으로 Nesbitt et al.(2019) 연구는 미국의 10개 도시를 대상으로 녹지 형평성을 분석하였다. 본 연구는 녹지 분포 특성상 특정 지역에 밀집해있다는 점을 들어, 공간적 자기상관의 가능성을 언급하였고, 이를 제어할 수 있는 공간통계 모형을 활용하였다. 분석결과, 도시별로 차이는 존재하지만 학력이 낮고, 소득이 낮은 집단이 주로 거주하는 지역에 환경 불평등 정도가 심각한 것으로 나타났다.

추가적으로 Rigolon(2017) 연구는 공원 및 녹지의 속성을 세분화하여 양, 접근성, 질 세 가지 측면에서 녹지 형평성을 분석하였다. 공원의 질 측면에서는 공원의 안전성(safety)이 대리 변수로 설정되었다. 분석결과, 사회적 소외계층들이 주로 거주하는 지역은 녹지의 양, 질 측면에서는 열악하지만, 접근성 측면에서는 양호한 것으로 나타나 기존 연구들과 다소 상반된 결과를 보고하고 있다. 추가적으로 미국을 대상으로 한 연구들은 인종간의 이슈가 존재하는 특성으로 사회적 소외계층을 저소득, 노인, 교육 수준이 낮은 집단뿐만 아니라, 유색인종까지도 포함하고 있다(Nesbitt et al., 2019; Rigolon, 2017; Schwarz et al., 2018).

이와 같이 녹지 불평등이 나타나는 이유는 환경 불평등적인 관점에서 해석되고 있다. 일례로 저소득계층의 경우, 주택 구입에 있어 시장경쟁으로 인해 질 높은 녹지 인프라를 포함한 쾌적한 환경자원에 대해 낮은 접근성을 지니게 된다(명수정, 2017). 이를 뒷받침하는 결과로 다양한 사례지를 대상으로 주택가격은 녹지에 대한 접근성 및 녹지의 양과 긍정적인 관계를 지니는 것으로 보고되고 있다(Fu et al., 2019; Park et al., 2017; Zhang and Dong, 2018). 따라서 녹지 접근성을 포함한 녹지 인프라가 좋은 근린에 위치하는 주택의 경우, 주택 가격이 높기 때문에 사회·경제적 소외계층이 질 높은 녹지 인프라가 형성되어 있는 근린에 거주하는 것이 어렵게 된다. 추가적인 녹지 불평등 이유로 소외계층들이 주로 거주하는 지역의 지방 정부는 부족한 세원으로 인해 공원 및 녹지를 충분하게 공급 및 관리하지 못한다는 점도 존재한다(Chen et al., 2020; Gordon-Larsen et al., 2006; Oh et al., 2020).

이러한 맥락에서 녹지 형평성 관련 연구는 환경 불평등 관점에서 사회·경제적 소외계층의 분포와 녹지의 분포를 통해 도출된 녹지 취약지를 대상으로 정부 차원에서 녹지에 대한 공급의 중요성을 강조하고 있다(고영주 외, 2019; 김미현 외, 2015; 명수정, 2017; Gordon-Larsen et al. 2006; Oh et al., 2020; Rigolon, 2017)

2. Google Street View를 활용한 연구

최근 빅데이터와 딥러닝에 대한 관심이 증대됨에 따라 도시계획 분야에서도 이를 이용하여 근린 및 가로환경을 측정하려는 시도가 다수 이루어지고 있다(Rzotkiewicz et al., 2018). 근린 및 가로환경을 측정하기 위해 Google Street View, Tencent Street View 등 다양한 민간업체에서 제공하는 가로경관 이미지를 활용하고 있으며, 이러한 방법론은 기존 측정방법인 현장조사에 비해 정확도 및 효율성이 높다고 언급되고 있다(Gong et al., 2018; Rzotkiewicz et al., 2018). GSV를 활용하여 주로 측정하고 있는 요소로는 천공개폐율(Gong et al., 2018), 가로수 음영 정도(Li and Ratti, 2018), 수공간이 보이는 정도(Helbich et al., 2019), 가로녹시율(Chen et al., 2020; Fu et al., 2019; Li et al., 2015; Lu, 2018; Ye et al., 2018; Zhang and Dong, 2018) 등이 존재한다.

Nguyen et al.(2019)의 연구는 GSV를 활용하여 미국 전역을 대상으로 고속도로, 녹지 등의 근린환경을 파악하였으며, 이를 이용하여 County 단위에서의 건강 수준과의 관련성을 분석하였다. 해당 연구는 근린환경 데이터가 세밀하게 구축되어 있지 않은 비도시 및 교외 지역의 근린환경 측정 시 GSV가 유용하게 활용될 수 있음을 언급하였고, GSV를 통해 산출한 근린환경과 건강 수준과의 유의한 관계성을 보고하였다. 그러나 넓은 공간적 범위로 인한 GSV 구득의 한계로 다소 세밀하지 못한 근린환경 측정에 그쳤다는 점에서 한계점이 존재한다.

보다 세밀한 근린환경을 측정한 연구로 Helbich et al.(2019)의 연구는 가로 이미지를 활용하여 보행자 관점에서 시각적으로 보이는 수공간 비율, 녹지 비율을 측정하였고, 이들과 거주민의 정신건강과의 관계를 분석하였다. 분석결과, 정신건강과 GSV로 산출한 수공간 비율, 녹지 비율 사이의 관계는 유의한 것으로 도출되었지만, 기존 녹지 및 수공간 산출 방법인 NDVI, Normalized Difference Water Index(NDWI)는 유의하지 않은 관계로 도출되었다. 추가적으로 GSV로 산출한 수공간 비율과 NDWI, GSV로 산출한 녹지 비율과 NDVI의 상관관계를 각각 분석한 결과, 두 방법 모두 녹지 및 수공간을 측정하는 방법임에도 각 변수 간의 유의한 관계가 도출되지 않았다. 이와 같은 결과는 동일한 근린환경 요소를 측정함에도 세부적인 요소의 포함 여부, 측정 방식(top-down view, eye level)에 따라 측정값이 크게 차이날 수 있음을 시사한다.

이와 유사한 연구 결과로 전통적 방식으로 산출한 녹지변수는 보행 및 신체활동과 일관되지 않은 결과를 보이는 반면 GSV를 통해 산출한 가로녹시율은 보행 및 신체활동과 밀접한 관련을 지닌다고 보고되고 있다(Ki and Lee, 2021; Lu et al., 2019). 이에 대한 이유로 GSV를 통해 근린환경을 측정하는 방식은 보행자가 일상생활 중에 접할 수 있는 녹지를 보행자 시점에서 산출하며, 기존 녹지 산출 방식으로 포함하지 못하는 벽면녹화, 정원, 관목 등의 소규모 녹지를 포함하고, 녹지의 활력까지 고려할 수 있기 때문이라 언급되고 있다. 이러한 특징으로 가로녹시율은 거주민의 실제 녹지 체감 정도를 잘 대변할 수 있다는 장점을 지닌 변수라고 볼 수 있다(Helbich et al., 2019; Ki and Lee, 2021; Lu et al., 2019; Ye et al., 2019).

GSV를 활용하여 녹지 형평성을 분석한 연구들은 최근 소수 진행되었다(Chen et al., 2020; Li et al., 2015; Li and Ratti, 2018). 일례로 Li et al.(2015)의 연구는 미국 Hartford를 대상으로 GSV와 Color band 방법론을 활용하여 가로녹시율을 산출하였다. 해당 연구는 대상지 내 block group별로 가로녹시율의 차이가 크며, 이러한 가로녹시율은 인구사회경제 속성과 밀접한 관련이 있다고 보고하고 있다. 구체적으로 저소득, 유색인종, 고등교육을 받지 못한 거주민이 많은 지역일수록 가로녹시율이 낮은 것으로 나타나, 녹지 불평등이 존재한다고 언급하고 있다.

이와 유사한 방법으로 Li and Ratti(2018) 연구는 GSV를 활용하여 가로수 음영 정도(tree shade provision)를 산출하였다. 이를 활용하여 토지이용별로, 인구사회경제 특성별로 음영 정도의 차이를 확인하였고, 이전 연구와 유사하게 사회적 취약계층들이 주로 거주하는 지역의 가로수 음영 정도가 낮게 나타났다. 그러나 이 두 연구에서 활용한 Color band 방법론은 색(RGB)을 기반으로 이미지 내 요소를 분류하는 방법으로, 녹색의 인공물(차량, 벽, 건물 등)이 녹지로 분류될 수 있는 한계점이 존재한다(Lu, 2018).

Chen et al.(2020)의 연구는 광저우, 선전 등 중국 9개 도시를 대상으로 Tencent Street View를 활용해 가로녹시율을 산출하였고, 이와 지역별(n=32) 총 GDP, 1인당 GDP, 재정수입 등과의 상관분석을 실시하였다. 분석결과, 지역의 경제적 지표와 가로녹시율 간의 관계는 강한 양(+)의 관계로 나타났다. 이는 도시의 녹지 공급의 주체가 지방정부이기 때문에 재정상황이 양호한 지역일수록 녹지에 대한 공급이 충분히 이루어졌기 때문이다. 하지만 해당 연구는 넓은 공간적 범위로 GSV 구득 지점이 다소 세밀하지 못한 점(300×300m)과 지역별 재정상황에 초점을 맞췄다는 점에서 본 연구와 차이가 있다고 볼 수 있다.

정리해보면, 최근 GSV 및 딥러닝을 활용하여 근린환경을 측정하려는 시도가 다수 이루어지고 있다. 이러한 접근방식은 녹지 측정 시 보행자 레벨에서의 녹지 측정이 가능하며, 기존 GIS 데이터로 구득이 어려운 상세한 녹지 형태 및 활력까지 고려할 수 있다. 이와 같은 장점에 기반하여 GSV로 측정한 녹지 변수는 기존 측정 방식의 변수보다 개인의 행태와 밀접한 관련을 지니며, 실제 거주민이 느끼는 녹지 체감 정도를 잘 대변할 수 있다고 보고되고 있다.

3. 선행연구 한계점 및 연구의 차별성

선행연구 고찰을 통한 선행연구의 한계점 및 본 연구의 차별성은 다음과 같다. 첫째, 녹지 형평성에 관한 실증 연구는 다수 이루어졌지만, 대부분의 연구가 NDVI, 토지이용 및 토지피복 데이터를 활용한 공원 및 녹지면적 변수를 형평성 분석을 위한 변수로 활용하고 있다. 이들 지표는 하늘에서 지상을 내려다보며(top-down view) 측정하는 방식으로 평면적이고 공급자 중심적이라는 한계점이 존재한다. 이로 인해 기존 녹지 지표는 거주민이 일상생활 중에 실제로 인식하는 입체적인 녹지의 양을 대변하지 못한다. 더욱이, 본 연구 대상지인 서울시의 경우 분지 지형으로 도시 외곽에 다량의 산지를 비롯한 대규모 녹지가 위치해 있다. 이러한 지형적인 특성과 외곽지역에 주로 거주하는 사회적 소외계층들의 분포 특성을 고려하였을 때, 기존 녹지 지표로 형평성을 분석할 시 녹지의 형평성 정도가 상대적으로 양호하게 산출될 수 있다. 그러나 외곽지역에 위치한 대규모 산지는 일상생활과 밀접하게 관련 있는 녹지로 보기 힘들기 때문에 거주민이 일상생활 중에 직접적으로 체감할 수 있는 녹지의 양을 과대산출 할 수 있다는 한계점이 존재한다. 따라서 본 연구에서는 GSV 이미지와 컴퓨터 비전을 활용하여 거주민의 실질적 녹지 체감 정도를 대변할 수 있는 가로녹시율 변수를 산출하여 녹지 형평성을 분석하였다는 점에서 차별성이 있다.

둘째, 녹지 형평성을 분석한 선행연구들은 단일 녹지 변수를 활용하여 형평성 정도를 분석하였다. 하지만 도시 내 녹지는 다양한 형태로 존재하며, 각 녹지 변수가 포함할 수 있는 세부 녹지 특성에 있어 차이가 나타날 것으로 예상된다. 즉, 단일한 녹지 변수로 형평성을 분석할 시, 도시 내 다양한 형태로 존재하는 녹지를 포괄적으로 다루지 못하며, 이는 실질적인 녹지의 불평등 정도를 확인하지 못할 수 있다. 일례로, 기존 녹지 측정 방식인 토지이용 및 토지피복 데이터만을 활용하여 녹지 분포를 분석할 시, 가로변에 위치하는 소규모 녹지, 사설 녹지 등을 포함하지 못하게 되며, NDVI 또한 해상도의 문제로 소규모의 녹지를 세밀하게 식별하지 못한다는 한계점이 존재한다. 이와 반대로 가로녹시율의 경우, 가로변에 존재하는 소규모 녹지, 사설 녹지 등을 포함할 수 있으며, 이들 녹지의 활력까지 고려하여 산출될 수 있다. 그러나 이는 가로변의 이미지만을 활용하기 때문에 전통적인 녹지 산출방법에 포함되는 대규모 녹지, 공원 등의 면적인 녹지를 포함하지 못한다는 특성이 있다. 즉, 본 연구가 활용하고자 하는 NDVI, 녹지 및 공원면적, 가로녹시율 각 변수가 포함하는 세부 녹지 특성은 상이하며, 이로 인해 각 녹지 변수의 분포 및 형평성 정도에 있어 차이가 나타날 것으로 예상된다. 따라서 본 연구에서는 위와 같은 세 가지 녹지 변수를 활용하여 형평성 정도를 분석하며, 변수 간 그 차이를 규명한다는 점에서 차별성이 있다.

셋째, GSV를 활용하여 산출된 가로녹시율을 대상으로 녹지 형평성을 분석한 소수 연구는 가로 이미지에서 녹지를 추출할 시, Color band 방법을 활용하고 있다. 본 방법론은 RGB 값만을 활용하여 이미지 내 속성을 분류하는 것이다. 즉, Color band를 활용할 시 녹색의 인공물(차량, 벽, 건물 등)이 녹지로 분류되는 경우가 존재하며(Lu, 2018), 이는 해당 지점의 가로녹시율을 과대 산출할 수 있다는 한계점이 존재한다. 이에 반해 컴퓨터 비전의 한 방법론인 의미론적 분할(semantic segmentation) 기법은 라벨링(labeling)된 이미지 데이터 셋을 활용한 학습과정을 통해 이미지 내에 존재하는 물체의 색 뿐만 아니라, 모형, 분포 등까지 고려하여 이미지 내 요소를 분류해주는 방법이다. 즉, 의미론적 분할기법은 Color band 방법론보다 높은 정확성으로 기존 연구의 한계점을 극복할 수 있을 것으로 보인다. 따라서 본 연구에서는 이미지 내 요소를 분류해주는 발전된 방법론이라 할 수 있는 의미론적 분할 기법을 활용하였다는 점에서 차별성이 존재한다.

1. 연구의 범위

본 연구는 가로 이미지 및 사회적 소외계층 데이터 구득의 용이성 측면에서 연구의 공간적 범위를 서울특별시로 설정하였다. 본 연구의 중점 변수인 가로녹시율을 산출하기 위한 GSV 이미지 데이터의 경우, 다른 도시들에 비해 서울시를 대상으로 조밀하게 구축되어 있다. 추가적으로 서울시는 본 연구의 분석단위인 행정동별 다양한 측면의 사회적 소외계층과 관련된 데이터 구득이 용이하다. <그림 1>에서 확인할 수 있는 것처럼 서울시는 분지지형으로, 외곽지역을 중심으로 상당한 녹지가 분포되어 있는 것을 확인할 수 있다.

Figure 1.

Study area

본 연구의 시간적 범위는 2018년 봄을 대상으로 하였다. 이에 대한 이유로 GSV는 한 시점에 서울시 전역의 가로경관이 구축되어 있는 것이 아니라 장기간에 걸쳐 시기별로 이미지 데이터가 구축되고 있는데, 2018년 봄에 구축된 사진이 가장 많기 때문에 이를 시간적 범위로 설정하였다(<표 2> 참고). 따라서 NDVI를 산출하기 위한 위성사진도 해당 시점에 맞춰 구득하였다.

2. 연구 질문 및 분석 방법

2) 분석과정

본 연구는 위와 같은 연구 질문을 검증하기 위해 분석과정은 크게 데이터 구축 단계와 통계분석 단계로 구분하여 진행된다(<그림 2> 참고). 데이터 구축 단계에서는 본 연구의 종속변수인 녹지 변수 구축을 위해 위성영상(Landsat 8), GSV 이미지, 서울시 생태현황도 자료를 활용하여 행정동별 NDVI, 가로녹시율, 1인당 공원 및 녹지면적 변수를 구축한다. 추가적으로 녹지 분포와 사회적 소외계층 거주지 간의 분포 비교를 위해 행정동별 인구·사회·경제 속성 변수를 구축한다.

Figure 2.

Research framework

통계분석에서는 다음과 같은 세부적인 두 단계로 구성된다. 첫째, 본 연구는 행정동별 구축된 세 가지 녹지 변수의 분포를 비교하기 위해 행정동별 녹지 분포 분석 및 Hot Spot 분석을 실시하여 각 녹지가 집중되어 있는 지역(High-High, HH), 녹지 양이 적거나 없는 밀집 지역(Low-Low, LL)을 확인한다. 추가적으로 세 가지 녹지 변수 분포의 차이를 비교한다.

두 번째 단계에서 공간통계 모형 활용을 검토하기 위해 공간적 종속성(spatial dependency)을 확인하고, 이를 계량화할 수 있는 공간적 자기상관(spatial autocorrelation) 지수인 전역적 모란지수(Global Moran’s I)를 산출한다. 이후 전역적 모란지수 값을 통해 공간적 자기상관이 존재한다고 판단되면, 이를 제어할 수 있는 공간통계모형을 활용한다. 최종적으로 사회적 소외계층들 간의 분포와 녹지변수별 분포 간의 공간회귀분석을 실시하여 서울시 녹지 형평성을 분석하고, 이를 기반으로 정책적 시사점을 도출한다.

3. 전통적 녹지 측정 변수

본 연구에서는 전통적 녹지 측정 변수로 토지피복 데이터를 활용한 녹지 및 공원면적, 위성사진을 활용한 NDVI를 측정하였다. 산출 시 본 연구의 분석단위인 행정동별 1인당 공원 및 녹지 면적, 평균 NDVI로 산출하였다.

2) NDVI

NDVI는 식생의 밀도 및 건강상태를 간접적으로 대변하는 대표적인 식생 지표로(정명희 외, 2012), 파장대의 반사특성을 이용하여 산출된다. 이를 산출하기 위해 USGS에서 제공하는 Landsat 8 위성영상을 활용하였다. 본 연구의 시간적 범위인 2018년 봄(3~5월)에 구득 가능한 위성사진이 3장 존재하나, 영상들이 찍힌 날짜 중 운량이 가장 적은 날짜인 3월 6일 위성사진을 활용하여 식 (1)과 같이 NDVI를 산출하였다(<표 1>, <그림 3> 참고).

$$

(1)

Table 1.

Description of Landsat 8 data

Figure 3.

Distribution of NDVI

NDVI 값은 0~1 값을 지니게 되며, 1에 가까울수록 식생이 상대적으로 많으며, 활동성이 높은 것으로 해석할 수 있다. <그림 3>의 NDVI 분포를 참고하면, <그림 1>의 공원 및 녹지면적과 유사하게, 큰 규모의 산지가 위치한 서울시 외곽지역과 산지 및 공원 등이 위치한 내부지역에 높게 나타나는 것을 확인할 수 있다.

본 연구에서는 30×30m 격자별로 산출되는 NDVI 값을 활용하여 격자 중심점을 통해 행정동별 평균 NDVI를 산출하였다. 특기할 만한 점으로 행정동별 NDVI 산출 시, 한강을 포함하는 격자는 제외되었다.

4. 가로녹시율 변수

1) GSV 구득 과정

본 연구에서는 가로녹시율을 산출하기 위해 가로경관 빅데이터 셋인 GSV 파노라마 이미지(panorama image)를 구득하였다. 파노라마 이미지는 좌우 360°, 상하 180°의 가로 경관을 담고 있는 이미지로(<그림 4> 위 참고), 해당 지점에서의 가로녹시율을 산출하는 데 용이하다고 판단된다.

Figure 4.

Examples of GSV panorama image (a) and preprocessed image (b)

그러나 일부 연구에서는 파노라마 이미지 특성상 왜곡현상이 존재하고, 이를 활용하여 근린환경을 추출할 시 근린환경이 과대 및 과소 산출될 수 있다는 점을 지적하고 있다(Tsai and Chang, 2013; Li and Ratti, 2018; Yin and Wang, 2016). 일례로 Tsai and Chang(2013)은 파노라마 이미지 왜곡현상이 이미지 상단 부분(하늘), 이미지 하단 부분(지면)에 심각하다고 언급하고 있다. 이와 반대로 파노라마 이미지의 중앙부분은 왜곡현상이 적기 때문에 해당 부분이 적절한 Point of Interest(POI)가 될 수 있다고 언급하고 있다. 또한, Xia et al.(2021)는 동일한 지점을 대상으로 파노라마 이미지와 일반 이미지(static image) 각각을 대상으로 가로녹시율을 산출하여 그 값을 비교하였다. 분석결과, 파노라마 이미지는 구득이 용이하여(open-access) 활용성이 높고, 직관적이며, 보행자의 녹지 인지를 잘 대변할 수 있기 때문에 가로녹시율 산출에 적합한 이미지임을 보고하고 있다.

일부 연구에서는 파노라마 이미지의 한계를 극복하기 위해 파노라마 이미지를 재투영시키거나(Li and Ratti, 2018), 파노라마 이미지 중 왜곡현상이 적은 일부분을 잘라 근린환경을 산출하고 있다(Yin and Wang, 2016). 구체적으로 Yin and Wang(2016)의 연구는 파노라마 이미지 중 왜곡현상이 적은 중앙부분을 잘라 시각적 개방성(하늘 비율)을 산출하였다. 해당 연구는 파노라마 이미지의 일부분을 자른 영역이 보행자의 시점과 더 일치한다고 언급하고 있다.

따라서 본 연구에서는 파노라마 이미지의 왜곡 문제를 해결하기 위해 왜곡현상이 상대적으로 적으며, 보행자의 시각과 좀 더 밀접하다고 판단되는 부분을 잘라 분석에 활용하였다. 구체적으로 <그림 4a>는 파노라마 원본 이미지로 선행연구에서 언급한 것처럼, 상단부분과 하단 부분의 왜곡현상이 심한 것을 확인할 수 있다. 이후 Yin and Wang(2016) 연구에서 적용한 방법으로 이미지의 상단 및 하단 부분을 자른 <그림 4b>의 경우, 왜곡현상이 상대적으로 감소한 것을 확인할 수 있다. 따라서 본 연구에서는 파노라마 이미지의 왜곡현상을 해결하기 위해 이와 같은 방법을 활용하였으며, 가로녹시율 산출 시 전처리된 이미지(그림 4b)를 활용하였다.

GSV 이미지 구득을 위해 서울시 내 도로망을 기준으로 50m 간격별로 지점을 생성하였고(n= 210,746), 각 지점별로 파노라마 이미지를 구득하였다. 그러나 가로녹시율의 경우, 년도 및 계절에 따라 변할 수 있기 때문에 시점이 다른 사진을 활용할 시 특정 지역의 가로녹시율이 과대 및 과소 추정될 수 있는 가능성이 존재한다. 따라서 본 연구에서는 Google Street View meta-data API를 활용하여 각 지점에서 촬영된 사진이 찍힌 년도, 계절을 분석하였고, 그 결과는 <표 2>와 같다.

Table 2.

GSV image collection date in Seoul

분석결과, 2018년에 촬영된 이미지는 전체 이미지 중 84.9%로 가장 많은 것으로 나타났으며, 2018년 144,902장의 이미지 중에서는 봄에 찍힌 이미지가 138,930장으로 가장 많은 것으로 나타났다. 따라서 본 연구에서는 행정동별 가로녹시율의 대표성을 확보하기 위해 가능한 많은 사진을 구득할 수 있는 2018년 봄에 찍힌 사진만을 활용하였고, 138,930장에 해당한다. 이는 행정동별 평균적으로 약 327.6장의 이미지가 존재하는 수치이며, 분석단위별 대표성 있는 가로녹시율을 산출하기에 적합한 수치라고 판단된다.

2) 의미론적 분할 기법

본 연구는 이미지 내 녹지 요소를 추출하기 위해 컴퓨터 비전의 한 기법인 의미론적 분할 기법을 활용하였으며, 이는 차량, 건물, 벽, 하늘, 녹지, 사람 등의 이미지 내에 존재하는 요소를 픽셀 단위에서 구분해주는 기법이다(<그림 5, 6> 참고). 의미론적 분할 모형으로는 FCN8s, DeepLabv3+, U-Net, ReSeg, PSPNet 등이 있다. 이와 같은 모형 중 본 연구는 구글이 2015년 개발한 이후 현재까지 지속적으로 발전하여 우수한 성능을 보이는 Deep-Labv3+ 모형을 활용하였다(Chen et al., 2018). 해당 모형은 Wang and Vermeulen(2020) 연구에서 근린환경을 측정하는 데 활용된 바 있다. 전술한 바와 같이, 의미론적 분할 기법은 대용량의 라벨링된 이미지 데이터 셋과 딥러닝을 통해 분할 모형이 학습(training)되고 이를 활용하여 GSV 이미지를 분할한다.

Figure 5.

Example of semantic segmentation (1) Green View Index: 18.85

Figure 6.

Example of semantic segmentation (2) Green View Index: 10.86

따라서 본 연구에서는 GSV와 유사한 가로경관 이미지를 포함하고 있는 Cityscape 데이터 셋을 활용하여 학습이 완료된 모형을 이용하였다(Cordts et al., 2016).¹⁾ 학습이 완료된 모형의 경우, Cityscape 데이터 셋을 대상으로 정확도가 0.907로 나타나 GSV 이미지를 분할하는 데 적합한 것으로 볼 수 있다.²⁾ 해당 모형은 도로, 보행로, 빌딩, 벽, 펜스, 녹지, 하늘, 보행자, 차량 등 총 19개의 항목을 분류할 수 있으나, 본 연구에서는 중점 변수인 녹지만을 분석에 활용하였다.

5. 변수 선정

환경정의에 관한 선행연구는 특정 인구계층, 인종에 집중되어 환경 불평등이 발생한다는 것을 전제로 하고 있다. 이들 연구는 인구계층을 크게 사회적 요인, 경제적 요인으로 구분하고 있으며, 인종, 소득, 교육 수준, 종교 등이 이에 해당한다. 구체적으로 녹지 형평성을 분석한 선행연구에서는 취약계층을 연령별(노인, 미성년자), 경제 수준별(저소득), 학력별(저학력자), 인종별(유색인종)로 구분하고 있다(Li et al., 2015; Nesbitt et al., 2019; Rigolon, 2017; Schwarz et al., 2018).

그러나 이들 연구의 대상지 대부분이 미국 도시인 점을 참고하였을 때, 본 연구의 대상지인 서울은 인종별 이슈가 상대적으로 크지 않고, 행정동별 학력, 직접적인 경제 수준 데이터(소득 데이터)를 구득할 수 없는 한계점이 존재한다. 이러한 이유로 국내 사례지를 대상으로 녹지 형평성을 분석한 연구는 연령별(노인, 미성년자), 경제 수준별(기초생활 수급자, 주택가격, 다세대 주택), 기타(독거노인, 장애인, 외국인)의 변수를 활용하고 있다(고영주 외, 2019; 김미현 외, 2015; Oh et al., 2020).

따라서 본 연구에서는 취약계층을 대변할 수 있는 변수들 중 인종 및 학력을 제외하였고, 경제 수준의 경우 소득의 대리변수로 활용되는 행정동별 아파트 가격을 활용하였다. 행정동 아파트 가격 산출 시, 국토교통부 실거래가를 활용하였으며, 행정동별 충분한 표본을 확보하기 위해 2018년 1년 동안 거래된 81,388건을 대상으로 행정동별 평균 아파트 가격을 산출하였다.

추가적으로 본 연구에서는 경제수준을 대리할 수 있는 변수인 기초생활 수급자 비율, 다세대·연립주택 비율을 추가하였다. 다세대·연립주택 비율 변수는 직접적인 경제 수준으로 보기에는 한계가 존재하나, 상대적으로 다른 주택 유형에 비해 주택 가격이 낮다는 점, 이로 인해 취약계층들이 주로 거주할 수 있다는 점을 참고하여(고영주 외, 2019), 간접적으로 취약계층을 대변할 수 있는 변수라 판단된다.

최종적으로 본 연구가 선정한 변수는 <표 3>과 같으며, 취약계층 변수는 연령별(65세 이상, 19세 이하, 독거노인), 경제 수준별(주택 가격, 기초생활 수급자, 다세대·연립주택) 이상 6개 변수이며, 추가적으로 인구밀도 변수를 고려하였다.

Table 3.

Description of variable

6. 공간 통계 모형

본 연구에서는 전역적 모란지수를 통해 공간적 자기상관을 확인하고 자기상관이 존재할 시, 이를 제어하기 위한 공간 통계모형을 활용한다. 대표적인 공간통계 모형인 공간 시차 모형을 활용하며, 공간 시차 모형, 공간 오차 모형은 각각 식 (3), 식 (4)와 같이 표현할 수 있다(Anselin, 2006).

• W: 공간가중행렬, X: 설명변수, λ, β: 계수,
• ε: 공간 의존된 오차항, μ: 오차항

두 모형 중 최종모형 선택은 Anselin(2006)이 제시한 것처럼 Lagrange Multiplier 테스트인 LM-Lag와 LM-Error 테스트를 통해 모형의 적합도인 Akaike Info Criterion(AIC), Schwarz Criterion(SC) 값을 활용하여 판단된다. 이와 같은 분석을 위해 공간적 자기상관을 제어할 수 있는 모듈을 포함하고 있는 GeoDa 프로그램을 활용하였다.

1. 녹지 변수별 분포 비교

본 연구의 첫 번째 연구 질문을 검증하기 위해 행정동별 1인당 공원 및 녹지면적, NDVI, 가로녹시율 분포를 확인하였고 그 결과는 <그림 7>과 같다. 우선 전통적인 녹지 변수인 1인당 공원 및 녹지면적, NDVI의 경우, 상대적으로 유사한 분포 양상을 보이고 있다. 두 변수 모두 서울시 외곽지역을 중심으로 높게 나타났으며, 이들 지역이 Hot Spot(HH)으로 도출됐다. 구체적으로 북한산, 도봉산 등이 위치한 서울시 북부지역, 관악산이 위치한 남서부 지역, 대모산, 청계산 등이 위치한 남동부 지역에 대규모 산지로 인해 녹지가 다량으로 분포하는 것으로 확인됐다.

Figure 7.

Distribution and Hot spot analysis by each green variable (Top: park & green area, Middle: NDVI, Bottom: green view index)

추가적으로 두 녹지변수의 경우, 서울시 내부 지역 중 일부분이 Cold Spot(LL)으로 나타났는데, 여의도를 중심으로 한 영등포구, 마포구, 구로구, 동작구, 성동구 일부 지역, 강남구 북부 지역이 이에 해당된다.

이와 반대로 본 연구에서 산출한 가로녹시율 변수의 경우, 전통적인 녹지 변수와 비교했을 때 상이한 분포를 보이고 있다 (<그림 7> 아래). 우선 전통적인 변수에서 외곽지역에 대규모 산지로 인해 녹지가 많이 분포하게 나타난 것에 비해, 가로녹시율은 이러한 현상이 줄어든 것으로 나타났다. 구체적으로 전통적 녹지 변수에서 Hot Spot으로 도출된 북부, 남서부, 남동부 지역의 Hot Spot이 작아지거나 Hot Spot이 아닌 지역으로 도출되었다. 더욱이, 전통적인 녹지 변수로 산출할 시 녹지가 많이 분포한 지역으로 나타난 광진구 동쪽(아차산), 관악구(신림 일대, 관악산 주변) 등은 오히려 가로녹시율 분포에서는 Cold Spot으로 나타났다. 이러한 차이를 보이는 이유로 기존 녹지 변수는 서울시 외곽지역에 위치한 대규모 산지 및 면적인 녹지를 주로 포함하는 것에 비해, 가로녹시율은 이보다 가로변에 존재하는 소규모 및 생활 녹지 등을 중심으로 포함하기 때문이다.

추가적으로 가로녹시율 분포에서 두드러지는 점은 강남구, 서초구를 중심으로 가로녹시율이 높게 분포하는 지역이 밀집해 있다는 것이다. 구체적으로 전통적인 녹지 변수에서는 강남구, 서초구 남쪽에 위치한 청계산, 대모산 일부 지역만 Hot Spot으로 도출된 것에 비해, 가로녹시율은 이들 자치구의 북쪽 지역까지 광범위한 Hot Spot이 도출되었다. 또한, Hot Spot으로 도출되진 않았지만, 잠실 일대, 청담 및 압구정동 일대까지 가로녹시율이 높게 분포한 것으로 확인되었다. 이와 같은 가로녹시율의 분포는 추후 진행될 사회적 소외계층과 가로녹시율의 분포의 관련성 분석 시, 소득이 높은 계층이 거주하는 근린에 가로녹시율이 높게 나타날 수 있는 경향성을 암시한다.

2. 공간적 자기상관 검토

이전 단계에서 각 녹지 변수의 분포를 확인한 결과, 세 녹지 변수 모두 녹지가 특정 지역에 공간적으로 밀집해있거나, 그 반대로 존재하지 않거나 적은 지역들이 밀집해 있는 것으로 나타났다. 즉, 이는 녹지분포에서 공간적 자기상관이 존재할 수 있다는 점을 암시한다.

따라서 본 연구에서는 이를 확인하기 위해 공간적 자기상관의 정도를 나타내는 지표인 전역적 모란지수(Global Moran’s I)를 세 가지 녹지변수별로 산출하였고 그 결과는 <그림 8>과 같다. Moran’s I 값은 -1~1 값을 지니며, 부호에 상관없이 그 값이 클수록 공간적 자기상관 정도가 높다고 해석할 수 있다. 분석결과, 세 가지 녹지변수 모두 공간적 자기상관이 강하게 유의한 것으로 나타났다(p＜0.01). 변수 별 차이로는 NDVI, 가로녹시율, 1인당 공원 및 녹지 면적 순으로 자기상관성이 높게 나타났다. 이러한 결과는 세 가지 녹지 변수에 공간적 종속성이 존재하며, 일반 회귀분석이 아닌 공간 통계모형을 활용하여 공간적 자기상관을 제어해주어야 함을 시사한다.

Figure 8.

Global Moran’s I value by each green variable

3. 기초통계량

본 연구의 기초통계량은 <표 4>와 같다. 본 연구의 종속변수인 세 가지 녹지 변수를 살펴보면 행정동별 그 값이 상이한 것을 알 수 있다. 일례로, 가로녹시율 평균은 약 10.1 이지만, 최소 2.7(광진구 중곡 1동), 최대 30.6(노원구 상계 2동)으로 행정동별 불균등하게 분포되어 있다. NDVI 또한 최소 0.121(마포구 대흥동), 최대 0.432(관악구 신림동)로 나타났다.

Table 4.

Descriptive statistics of variable

서울시 행정동 고령자 및 19세 이하 평균 비율은 각각 14.3%, 15.5%지만, 이 변수 또한 행정동별로 큰 차이를 보이는 것으로 나타났다. 구체적으로 고령자, 19세 이하 비율이 가장 높은 행정동은 약 30.7%, 33.7%이다. 경제적 수준을 대리하는 변수 중 하나인 주택 가격 변수는 1m²당 평균 827만 원으로 나타났으나, 가장 주택 가격이 높은 행정동은 서초구 반포본동으로 1m²당 2,940만 원인 것으로 확인되었다. 특기할 만한 점으로 본 연구에서는 2018년 1년간 거래된 실거래가 자료를 활용하였으나, 4개 행정동(삼청동, 능동, 가회동, 고덕 2동)에서는 아파트가 존재하지 않거나, 실거래가 이루어지지 않아 분석에 활용되지 않았다. 즉, 본 연구의 최종 표본은 420개이다.

4. 공간통계 모형

본 연구에서는 가로녹시율 및 전통적인 녹지 변수의 형평성 정도를 분석하기 위해 대표적인 공간통계 모형인 SEM, SLM 모형을 활용하였고, 이들 모형과 비교를 위해 OLS 모형 결과를 함께 기재하였다.³⁾ 본 연구에서는 1인당 공원 및 녹지면적, NDVI, 가로녹시율 세 녹지변수를 산출하여 분석하였으나, 1인당 공원 및 녹지면적 모형의 경우, 다른 전통적 녹지 변수라 할 수 있는 NDVI에 비해 설명력이 매우 낮게 나타났으며(OLS 모형 기준 0.183), 대부분의 사회·경제적 취약계층 변수와 유의한 관계를 지니지 않는 것으로 나타났다. 따라서 전통적 녹지변수 중 NDVI 결과만을 본문에 기재하였고, 이와 새로운 녹지변수인 가로녹시율 간의 비교를 실시하였다. 우선 <표 5>의 가로녹시율 분석결과를 살펴보면 OLS 모형은 Moran’s I이 유의하게 나타나, 공간적 자기상관이 존재하는 것을 재차 확인하였다. 그러나 SEM, SLM 모형의 경우 Moran’s I 값이 유의하지 않은 것으로 나타났고, 이는 두 모형이 공간적 자기상관을 제어하는 데 적합한 모형임을 의미한다.

Table 5.

Result of Green View Index model

SEM, SLM 두 모형 중 최종 모형 선택 기준은 LR test 값이 높고, 상대적 적합도를 의미하는 AIC, SC 값이 낮은 모형이다. 두 모형 중 SEM 모형이 LR test 값이 높으며, AIC, SC 값이 낮은 것으로 나타나 SEM을 최종 모형으로 선정하였다. 추가적으로 SEM 모형이 SLM에 비해 높은 R² 값을 보이는 것으로 확인되었다.

SEM 모형에서는 65세 이상 인구 비율, 19세 이하 인구 비율, 주택 가격, 기초생활 수급자 비율, 다세대 및 연립주택 비율, 인구 밀도가 유의한 변수로 나타났다. 65세 이상, 19세 이하 인구 비율 두 변수 모두 양(+)의 관계로 나타났고, 이는 이들 인구 집단이 많이 거주하는 행정동에 높은 가로녹시율이 분포되어 있다고 해석할 수 있다. 즉, 노인, 미성년자 경우, 가로녹시율의 형평성 정도가 양호한 수준이라 볼 수 있다.

그러나 경제적 수준과 관련된 소외계층 변수인 주택 가격, 기초 생활 수급자 비율, 다세·연립주택 비율은 가로녹시율과 각각 양(+), 음(-), 음(-)의 관계로 나타났다. 즉, 주택 가격이 저렴하고, 기초생활 수급자가 많으며, 다세대·연립주택 비율이 높은 행정동에 가로녹시율이 낮다고 해석할 수 있다. 이는 특정 사회·경제적 집단의 거주지에 녹지의 분포가 불평등하게 나타날 것이라는 가설을 뒷받침 하는 결과라 볼 수 있다. 선행연구에서 지적하는 녹지에 대한 환경 불평등이 서울시 가로변 녹지에도 적용되는 것으로 가로변 녹지 인프라가 좋은 근린에 높은 주택 가격이 형성되어 있기 때문에(Fu et al., 2019; Zhang and Dong, 2018), 저소득층들의 경우 가로녹지 인프라가 충분한 근린에 거주하지 못하는 현상이 반영된 것으로 보인다.

전통적 녹지변수 중 하나인 NDVI 모형의 결과는 <표 6>과 같다. NDVI 또한 가로녹시율 모형과 마찬가지로 OLS 모형에서는 공간적 자기상관이 제어가 되지 않은 것을 확인할 수 있다. 그러나 이전과는 상이하게 두 공간통계 모형 중 SEM 모형에서만 공간적 자기상관이 완전히 제어되었다. 추가적으로 AIC, SC, LR test 값 모두 SEM 모형이 상대적으로 더 적합함을 보이고 있기 때문에, SEM 모형을 최종 모형으로 설정하였다.

Table 6.

Result of NDVI model

우선 소외계층 중 인구특성을 대변하는 변수인 65세 이상, 19세 이하 모두 양(+)의 관계로 나타나, 가로녹시율 모형과 동일하게 연령에 관한 소외계층들에 있어서 형평성 정도가 양호하다고 볼 수 있다. 경제적 능력과 관련된 취약계층 변수인 주택 가격, 기초생활 수급자 비율, 다세대 및 연립주택 비율 변수의 경우, 각각 음(-), 양(+), 양(+)의 관계로 유의한 것으로 나타나 모든 변수에 있어 가로녹시율 모형과 반대의 결과를 나타내고 있다. 이는 주택가격이 저렴하고, 기초생활 수급자 및 다세대·연립주택이 많은 지역에 NDVI가 높다는 것으로 NDVI로 측정한 녹지 분포는 경제적 취약계층과의 관계에 있어서도 형평성이 양호한 것으로 볼 수 있다.

이처럼 NDVI 모형은 경제적 소외계층 측면에서 가로녹시율 모형과 상반된 결과를 보이고 있다. 앞서 언급한 것처럼 NDVI는 연구 대상지인 서울시 외곽에 주로 위치하는 대규모 산지 등을 주로 포함함으로써 외곽지역을 중심으로 높게 나타난다(<그림 7> 참고). 물론 NDVI는 가로변 식생, 소규모 녹지 등을 일부 포함할 수 있지만, 상대적으로 낮은 해상도(30×30m)로 인해 이러한 녹지를 세밀하게 파악하는 것에는 한계가 있다. 즉, NDVI의 분포 특성은 도시 외곽지역에 주로 거주하는 사회·경제적 소외계층들의 분포와 유사하기 때문에 녹지 형평성 정도가 양호한 것으로 산출되게 된다.

이를 좀 더 구체적으로 확인하기 위해 본 연구에서는 서울시 행정동별 경제적 소외계층의 분포를 확인하였다(<그림 9> 참고). 기초생활 수급자, 다세대 주택 비율 모두 서울시 외곽지역에 높게 분포하는 것을 확인할 수 있다. 이와 반대로 주택 가격 변수는 동남부(강남구, 서초구 남쪽)를 제외한 서울시 외곽지역은 낮게 형성되어 있고, 강남구, 서초구, 송파구를 중심으로 매우 높게 형성되어 있다. 즉, 사회·경제적 소외계층은 도시 외곽지역을 중심으로 주로 거주하고, 이러한 외곽지역을 중심으로 전통적 녹지 측정변수가 높게 나타나는 점으로 인해 NDVI 지표로 녹지 형평성을 분석할 시 양호하게 나타나는 것으로 판단된다.

Figure 9.

Distribution of socio-economic vulnerable group (Top: housing price, Middle: basic livelihood security recipient rate, Bottom: multiplex housing rate)

본 연구 결과를 종합해보면 가로녹시율 모형에 한해서 녹지에 대한 불평등이 존재하는 것으로 나타났다. 즉, 대규모 공원 등의 면적인 녹지 측면에서는 불평등이 존재하지 않지만, 일상생활 중에 쉽게 접할 수 있는 가로변 녹지에 대한 형평성이 열악하다고 해석될 수 있다.