안녕하세요. 논문 준비 중인 대학원생입니다. 

이번에는 다름이 아니고, 저번에 올렸던 질문에 대한 궁금증이 생겨서 비슷한 내용으로 다시 올리게 되었습니다.


제가 사용하는 프로그램은 Arcgis pro 2.8.1 버전입니다.


들어가기에 앞서 간단히 제가 사용하는 데이터 자료를 설명드리자면, 중국 동북지역의 DEM(30m) 파일을 merge하여 사용하고 있습니다.

여쭙고자 하는 본 내용은 polygon 파일의 면적과 둘레를 어떻게 알 수 있는지에 관한 것입니다.


1. 보란색 polygon이 둘레와 면적을 구하고자 하는 대상입니다


사진1.JPG

보시면 Shape_Length가 둘레일텐데, 직접 폴리곤을 따라서 구한 길이 값과 차이가 큽니다. Shape_Length:2229m로 제가 구한 길이는 3,018m로 차이가 납니다. 이 오류를 잡고 싶습니다. 


2. 면적도 마찬가지 방법으로 진행해봤습니다.


Inked사진2_LI.jpg


면적도 제가 직접 따라서 그린 면적값과 다르게 나옵니다. Shape_Area에서는 239,920 평방미터로 나왔고, 제가 직접 측정한 수치는 408,807 평방미터로 나와서 역시 차이가 큽니다.


3. 이러한 문제점을 잡기 위해 나름대로 Add Geometry Attributes 기능을 사용했습니다.


사진3.JPG


이런 식으로 해당된 폴리곤을 넣고 Geodesic Area(면적), Geodesic length of perimeter(길이)를 실행했습니다.


4. 그 결과는 아래와 같습니다.


사진4.JPG

원래 나오던 데이터 수치랑 별반 차이가 없게 결과가 나왔는데, 해결 방법이 필요할 것 같습니다.


제가 측정해야 될 폴리곤 수가 50개가 넘어서 일일이 폴리곤 라인을 따라서 그리자니 시간이 부족하고, 무엇보다도 제가 라인을 따라서 그린 값은 정확한

데이터 값이 아니기 때문에 정확하게 폴리곤의 면적과 둘레를 구하는 방법을 알고 싶습니다.




추가적으로 저번 질문에 감사하게도 다른 분께서 약간의 힌트를 주신적이 있습니다만 이해하지 못했습니다.


그분 말씀으로는 해당 레이어의 속성테이블에 들어가서 필드계산기를 열어서 "SHAPE_LENGTH" = '$length' 함수를 이용해 길이와 면적을 설정해주면

된다고 하셨는데, 구체적으로 어떻게 함수를 써야 길이와 면적이 계산되는지 모르겠습니다.

필드계산기까지 열어봤는데 이 이상 어디에 어떤 식으로 함수식을 풀어야 할지 알고 싶습니다. 


그리고 이외에 다른 방법도 있는지 알고 싶습니다.

Inked필드계산_LI.jpg


profile

mosaictailor

2021.08.09 10:05
*.70.80.32

지금 사용하고 계신 좌표계가 뭔지는 모르겠지만, 

Add Geometry에서 사용하신 좌표계(EPSG 3857)은 타일맵 서비스에 특화된 좌표계이므로 면적, 거리 측정에 적합하지 않습니다.


각 좌표계마다 거리, 방향, 면적에 있어 왜곡이 발생하는데 위 좌표계는 이를 최소화하는 방향으로 설계된 좌표계가 아닙니다.

UTM 등을 사용하면 조금 더 근사치에 해당하는 결과를 얻을 수 있을 것으로 보입니다.



profile

sincerekjh

2021.08.09 14:32
*.152.224.28

Define Projection에서 UTM 들어가서 WGS 1984 UTM Zone 51N 혹은 52N으로 설정하면 아래 첨부파일처럼

동남아시아로 dem 파일이 설정됩니다. 혹시 UTM에서 어떤 걸 참고해야될까요?


캡처.JPG






좌표 정보는 이걸 참고했습니다.


좌표계 정리.jpg

첨부
profile

sincerekjh

2021.08.09 14:35
*.152.224.28

현재 Dem 파일에서 사용중인 좌표는 Korea_2000_Transverse_Mercator 입니다

profile

mosaictailor

2021.08.09 16:50
*.70.80.32

https://upload.wikimedia.org/wikipedia/commons/e/ed/Utm-zones.jpg


UTM zone 코드는 위 링크 참조하시면 됩니다.

지금 Polygon으로 작업해서 면적 연산중이신걸로 보이는데, Polygon 생성할 때 위의 UTM 좌표에 해당하는 자료로 생성하시면 될 것 같습니다.

profile

sincerekjh

2021.08.09 19:16
*.152.224.28

감사합니다. 덕분에 UTM Zone 51T와 52T가 각각 필요하다는 걸 알게 되었습니다.

다만 Define projection에서 UTM에서 UTM Zone 51T 및 52T가 보이지 않습니다. 

혹시 다운받는 사이트가 있나요? 구글링을 못하는 건지 방법을 모르겠습니다...ㅜㅜ


제가 한 방법은 1~3번으로 이루어져 있습니다.

1.투영좌표계에 들어가서

2.UTM을 클릭, WGS 1984에 들어가면

3.Northern Hemisphere이 나옵니다. 역시 클릭하면, WGS 1984 UTM Zone 1N~60N 혹은 WGS 1984 Complex UTM Zone 20N~30N까지만 나옵니다. 

혹시 UTM Zone 51T 혹은 52T가 나오게 할 해결 방법이 있을까요?

profile

mosaictailor

2021.08.09 19:53
*.147.161.49
중국 동북부 지역이면 wgs84 51N 52N 사용하시면 될 것 같아보이는데 중국 데이터는 다뤄본 적이 없어서 확실친 않습니다. 위 프로젝션으로 테스트 해보시면 좋을 것 같습니다.
profile

sincerekjh

2021.08.09 20:04
*.152.224.28

죄송하지만, 위에서 댓글 달았던 것처럼 wgs84 51N을 사용하면 이렇게 나오는데, 해결 방법이 있을까요?844dbda4be81c370c03f70ef9bae8851.jpg

첨부
profile

mosaictailor

2021.08.10 14:16
*.102.128.89
현재 올려주신 사진이 갖고 계신 dem 인거죠?
dem을 utm으로 define 하라는 의미는 아닙니다.

define은 본 자료가 구축될 혹은 유통될 당시의 좌표계를 기준으로 설정하는 것이기 때문에 취득하셨을 당시의 좌표계를 유지하시면 됩니다.

위에서 언급한 utm은 자료의 면적 연산이 되는 기준 자료. 즉, 4번 그림에서 언급하신 폴리곤의 좌표가 utm이어야 한다는 의미입니다.

그래서 새로 만드는 폴리곤의 껍데기를 utm 51n 또는 52n으로 하시거나, 지금 사용하시는 폴리곤을 project 도구를 이용해서 utm으로 변경한 뒤 폴리곤을 그리고, geometry를 확인해보시면 될 것 같습니다.

define projection 도구는 좌표를 변경하는 도구가 아니고 좌표계를 정의해주는 도구입니다.
profile

sincerekjh

2021.08.10 19:17
*.152.224.28

dem 파일을 utm으로 define 하는 것이 아닌 polygon 파일을 utm으로 좌표계 설정하라는 말씀 이해했습니다. 

하다보니까 끝도 없이 궁금한게 꼬리에 꼬리를 물고 생기네요...제가 이쪽으로 배움이 부족해서 이해하는데 많이 느렸습니다.

마지막 질문이 될 것 같은데요. 

제가 가진 최초의 dem 파일 좌표계가 WGS 1984인데요. contour를 만들어서 지형데이터를 기반으로 polygon을 만들어서 면적과 둘레를 연산하고 있습니다. 여기서 궁금한게 생겼습니다. 위에 답변 주신대로 polygon의 좌표계는 utm 51n 또는 52n으로 한다면 dem 좌표계와 다른데 그래도 상관없다는 말씀이시죠? polygon 파일을 만들때 contour를 고려하고 있습니다. 그렇다면 polygon이 위치한 contour 또한 좌표계를 제대로 설정해야 적합한 면적이나 둘레가 계산될텐데, 그러면 contour의 좌표를 utm 51n이나 52n으로 설정해야 하나요? 아니면 WGS 1984로 설정해야하나요?

profile

mosaictailor

2021.08.10 19:42
*.70.80.32


가장 최선의 방법은 아닐 것으로 인지하고 있지만, 

헷갈리실 것 같아서 조금 더 쉽고 정확한 방법으로 답변드립니다.


1. UTM51n 또는 UTM52n을 기준으로 Polygon을 생성합니다. (이하 Polygon1)

2. 면적 연산을 위한 영역을 Polygon1에 생성합니다.

3. Add Geometry Attributes 도구를 이용해서 면적을 연산/비교합니다.

-> 이상이 없을 경우 4번으로 이동. 이상이 있다면 1번으로 돌아가 좌표계를 바꾸어보는 것을 추천합니다.


4. Polygon1을 DEM과 동일한 좌표계로 변경하기 위해서 Project 도구를 실행합니다.

4.1. Input에는 Polygon1을 Output에는 원하는 이름을 지정합니다. (이하 Polygon2)

4.2. Output Coordinate System에는 DEM과 동일한 좌표계를 입력합니다. (위의 내용에 따르면 WGS1984, 즉 EPSG 4326으로 판단됨)


5. Contour 등의 작업 시 Polygon2를 사용합니다.



// 요지는 면적 계산용과 그렇지 않은 Polygon을 구분하여 사용하시는 것으로 설명하였습니다.

// 반드시 분리하여 사용할 필요는 없을 수 있으나, 혹시 모를 오류에 대비하여 구분하는 것으로 설명하였습니다.


잘 해결 되시면 좋겠네요.


profile

sincerekjh

2021.08.10 20:49
*.152.224.28

말씀대로 polygon1과 polygon2로 나눠서 작업 중입니다.

알려주신대로 하다가 궁금해서 샘플로 몇 개만 뽑아서 원보고서와 비교해보았는데 수치가 거의 일치하게 나와서 기쁘네요!!

거의 기초지식이 없는 상황에서 선생님 덕분에 잘 해결됐습니다!

정말 정말 감사합니다!!

List of Articles
번호 제목 글쓴이 날짜 조회 수
공지 질문올리 실때 조금만 질문 내용을 조금만 더 자세히 올려주시면 답변받기 수월합니다. + 11 김한국 2008-01-04 240834
7926 QGIS 일련번호 + 2 벼리맘 2021-09-08 199
7925 XY좌표만 입력되어있는 CSV파일을 쉐이프파일로 변환하고 싶습니다 ㅠㅠ... + 1 뇽뇽뇽 2021-09-08 74
7924 gis로 우수흐름표기 하는방법좀 알려주세요 ㅠㅠ + 3 용이아디 2021-09-06 105
7923 qgis 에서 수치지형도2.0과 지적도를 겹쳤는데 dxf 내보내기시 따로노는 현상.. file + 1 망치맨 2021-09-06 67
7922 ArcGIS에서 미지의 TM 좌표계를 GRS80-UTM-K, EPSG: 5181로 좌표변환 해주고 싶습니다. file + 2 RENUS 2021-09-05 128
7921 Arcgis질문사항 : 격자를 입힌 shp파일 - 점 갯수만큼 격자에 표시하고 싶습니다. file + 1 san 2021-09-03 58
7920 수치표고모델 SRTM과 ASTER 중에서 SRTM의 정확도가 더 높다고하는데 혹시 SRTM(Ver 3.0)을 어디서 얻을 수... sincerekjh 2021-09-03 29
7919 지오코딩 공개버전 파일로 주소 올릴때 지도에서 보기 기능 file kyi1221 2021-09-02 41
7918 속성테이블 도구창 문의 file + 1 gunseok 2021-09-02 39
7917 biz-gis 내에서 연도별 데이터를 알아볼 수 있을까요? file chelestinoo 2021-09-02 17
7916 Arcgis Pro 내에서, polyline을 기준으로 Polygon을 자르는법이 어떻게 될까요..조언 부탁드립니다. file + 1 전진 2021-09-02 39
7915 표고분석 및 경사분석시 필요한 자료가 무었이 있나요? + 1 watergis 2021-09-01 73
7914 [과거 임상도] 자료를 구하고자 합니다.! + 1 jane27 2021-08-31 89
7913 q gis 심볼 추가 히몬 2021-08-31 29
7912 DEM 파일에서 고도자료 추출하는 방법 도와주세요..윽엑 2021-08-31 55
7911 서울시 DEM 데이터 어디서 구할 수 있을까요? (해상도 1mX1m) + 1 GIS공부하는학생 2021-08-31 67
7910 1 file 폴짝 2021-08-30 30
7909 정사영상에 자꾸 선이 갑니다. file 김시락 2021-08-29 56
7908 전국 지하철 역사 Shapefile과 부산시 경계면 Shapefile의 좌표가 자꾸 어긋납니다(Project) file RENUS 2021-08-27 52
7907 관로 shp를 cad로 도면화 할때 절점이 없는 관로는 선으로 표기되는데 모두 폴리선으로 도면화 할수 없나요? 행방불명 2021-08-26 33
7906 service area 폴리곤 내 포인트 수 아는 법 + 1 peace 2021-08-24 73
7905 qgis 내 심볼에는 등간격으로 나눠서 등급별로 색을 부여하는 기능이 있는데,, 필드 계산기 내에서도 값을 ... file 소라빵 2021-08-24 68
7904 군사보호구역 shp 파일은 어디서 구할수있을까요 구알라구알라 2021-08-23 67
7903 EPSG 3857에 대해서 궁금해요 file 신난다햄토리 2021-08-23 78
7902 qgis 거리 계산 라인 만들기 우하하 2021-08-22 72
7901 Qgis에서 플러그인>glass기능 관련 yoonja01 2021-08-20 33
7900 안녕하세요,, 외주관련해서 여쭤보고 싶습니다 + 2 도라에몽박 2021-08-20 143
7899 QGIS 몇가지 여쭤봅니다. + 1 아요티 2021-08-19 109
7898 Building Pyramids 무시해도 되나요?? file 다이라트스 2021-08-19 33
7897 인접한 지역간 경계를 데이터화 할수있을까요? 마도임 2021-08-18 24
7896 동일한 공간 위치의 레스터별로 주성분 분석 및 상관분석?을 어떻게 할까요..? WM 2021-08-17 26
7895 고속도로 차로수 변경지점 추출방법에 대한 질문 쥐린이 2021-08-17 31
7894 캐드 일반파일에 특정 좌표를 넣고싶습니다 악악악으아아악 2021-08-13 43
7893 이동경로 추적 질문입니다! file 쥐린이 2021-08-12 41
7892 qgis 텍스트20자리인데 새로운필드를 10자리부터 15자리까지로 만들고싶어요 + 1 니목 2021-08-12 44
7891 여러개 SHP 파일을 래스터로 변환할 때 색상 통일하는 방법 없나요? file + 3 다이라트스 2021-08-12 90
7890 국가공간정보포털에서 제공하는 지적도와 다음지도가 중첩이 안 됩니다. 도움 부탁드립니다. file + 1 환경인 2021-08-10 100
7889 arcgis 위치 조인 질문입니다! 도와주세요ㅠㅠ file + 4 빠방 2021-08-10 125
7888 - 빠방 2021-08-10 24
7887 ARCGIS에서 MAP을 이미지파일로 EXPORT할떄 타일로 출력가능한가요? file 민지렁이 2021-08-10 37
7886 테이블에서 마우스로 클릭하지 않고 폴리곤 확인하는 방법 있나요? file + 1 다이라트스 2021-08-10 165
7885 좌표계 관련 질문! WGS 1984 UTM Zone 51T 혹은 52T를 어디서 찾을 수 있습니까? sincerekjh 2021-08-09 31
7884 두 자료가 겹치지 않아요ㅜㅜ 좌표계도 통일시켰는데.. file + 4 구뜨 2021-08-09 183
7883 tif파일을 shp파일로 변환하는 방법 좀 알려주세요 마르꼬 2021-08-08 55
» 안녕하세요. Arcgis pro를 이용하여 Polygon의 둘레와 면적을 산출하려고 합니다. 제가 하는 방법이 맞는지... file + 11 sincerekjh 2021-08-06 137
7881 아크gis 작업한 지번라벨들을 qgis로 그모양 그대로 불러올 수 있나요? QGISER 2021-08-06 45
7880 레이어를 이미지 내보내기 하면 이미지가 사라집니다. file 흥구네 2021-08-06 32
7879 레이어 추가가 안됩니다 ㅠㅠ(사진첨부) file + 1 나플라 2021-08-06 64
7878 arcgis 솔루션컨텐츠 다운로드 문의 비트 2021-08-06 23
7877 표고분석도 질문입니다 file + 1 cobe 2021-08-05 98