28

2007-Nov

Spatial Analyst를 이용한 보간(IDW와 같은 지역보간) 절차

작성자: 김한국 IP ADRESS: *.181.253.30 조회 수: 29306

김수영님께서 질문하신 내용인데 다른 분들에게도 필요할 것 같아서 정리해 봤습니다.


 


1. 점데이터를 올립니다. 당연히 점데이터에는 보간할 데이터가 숫자형으로 들어있어야 합니다.



 


2. Spatial Analyst에서 Interpolate to Raster->Inverse Distance Weighted...를 선택합니다.


(보간에는 여러가지 방법이 있습니다. 각 방법마다 특성이 다르구요... 분석가의 경험과 노하우를 이용하여 가장 알맞는 방법을 선택하여햐 합니다.)


 


Inverse Distance Weighted는 IDW라고합니다. 뜻을 헤아려보면... 거리에 가중 값을 주는데 가중값을 역으로 준다는 것입니다. 무슨 말이냐면 점으로부터 거리가 가까울 수록 점의 값에 가까워진다는 지극히 평범한 내용입니다. ^^



 


3. 아래와 같은 대화상자가 나타납니다.


여기서 Z value field라는 것이 보간할 값을 입력하는 곳입니다.


다시말해 서울시 대기오염 측정 자료가 있는데 그 중에서 일산화탄소 량을 보간하고 싶다고 하면 일상화탄소 값이 들어 있는 필드를 선택하여야 합니다.



 


4. 아래보면 결과가 나왔습니다.


그림을 보면 점데이터가 분포한 EXTENT에 해당하는 범위에만 보간이 되어 있습니다.


당연한 결과입니다.


 


다음 점의 value를 파악하고 보간을 해나가야 하는데 가장 변두리에 있는 녀석의 경우 보간을 할 점이 없으니 더 이상 보간을 못하는 것입니다.


이 부분을 좀더 확대해서 뻗어나가게 하고 싶다면 꽁수는 있습니다.


네 귀퉁이에 임의의 점 4개를 더 만들어 주고 다시 보간을 하면 해당 점까지 확장이 되어 보간이 되겠죠? 하지만 이 방법은 측정값이 들어가는 것이 아니고 임의로 입력한 값이기 때문에 상당히 위헙합니다.



 


5. 그림이 너무 평범해 보입니다.


좀더 변화의 추이를 더 잘 파악하도록 해보겠습니다.



 


6. 등급을 나누는 방법을 선택할 수 있습니다.


일반적으로 등간격으로 관찰합니다.



 


7. 한번 Quantile로 바꿔보겠습니다.(통계책을 뒤져보면 Quantile이 무엇인지 알 수 있습니다.)



 


8. 달라졌습니다.(보기방법이 달라진 것이지 내부에 데이터가 변경된 것은 아닙니다.)



 


9. 확대해 보겠습니다.



 


 


이러한 보간 기법은 GIS에서 가장 많이 쓰는 기법 중 하나입니다.


좀 직설적으로 얘기하면 보간방법을 한번도 활용해 보지 않았다고 한다면 GIS를 단순 뷰어정도로만 활용했다고 볼 수 있습니다. ^^


 


강우관측소 데이터가 있다면 이러한 보간 방법을 이용하여 전국 단위의 강우맵을 만들수가 있습니다.


그렇다면 홍수예측이나 토사유출량을 구할 때 유용하게 사용할 수 있습니다.


 


하지만 대부분의 엔지니어링 회사에서는 티센폴리곤을 이용하여 강우량을 예측하고 있습니다.(관공서에서 조례가 티센을 이용하도록 나와 있으니 어쩔수 없지만... ^^)


티센 폴리곤보다 이러한 보간방법이 훨씬 정확한데도 말입니다.


 


 

profile

Map for your Life ~~~~~~ !

biz-gis.com 운영진입니다.

profile

애인있어요

2007.11.28 14:01
*.228.188.222

와우!! 감사합니다.


정말 제가 필요로 하던 내용이었고,


특히 평범한 도면 속에 숨어있는 차이점을 이처럼 약간의 테크닉만으로도 눈에 확띄게 바꿔놓을 수 있다는 점에서


더욱 고수의 손길이 느껴지는군요....


제 데이터에 바로 적용할 수 있을 것 같습니다.


감사합니다.

List of Articles
번호 제목 글쓴이 날짜 조회 수
64 좌표계 정의가 되지 않은 캐드 파일(dxf, dwg) 좌표부여하기... 김한국 2011-07-22 56436
63 GRS80 좌표계 변환결과 300m 정도 오차나는 경우 file + 3 김한국 2011-05-24 49289
62 엑셀에서 shp 파일 만들기 - 이상희님 자료입니다. file + 4 김한국 2011-04-12 48550
61 래스터 디스플레이 할때 Too many unique values (>65536) 에러 해결법 file + 2 김한국 2011-04-01 30933
60 ArcMap 매크로(VBA) 기능이 비활성화 되는 문제 해결법 (ArcMap Visual Basic Editor macro is disabled) + 1 김한국 2010-08-09 37967
59 ArcMap 테이블의 복잡한 필드 계산(field calculator)시 방법 file + 3 김한국 2010-06-18 39573
58 DEM의 표면면적(3D면적) 구하기 file + 2 김한국 2010-04-14 38856
57 래스터 데이터(토지피복 등) 폴리곤별 포함면적 계산 file + 5 김한국 2010-01-12 36707
56 [re] 수리수문, 풍향풍속과 같은 현황을 화살표를 활용하여 어떻게 표현하나요?? + 1 박세진 2009-09-22 38532
55 수리수문, 풍향풍속과 같은 현황을 화살표를 활용하여 어떻게 표현하나요?? file + 3 eric~~~~ 2009-09-21 38745
54 경위도(Degree->DMS) 변환설명 이기우 2009-06-10 40265
53 CAD MAP에서 속성별로 레이어만들기 + 6 강세경 2009-06-08 43492
52 Raster Calculator에서 삼각함수(sin, tan, cos) 사용법 + 1 김한국 2009-06-04 41027
51 Cad→ Shape변환 관련자료입니다. file + 10 조덕연 2009-05-07 39923
50 지적도(폴리곤)에서 중첩된 영역 추출하는 방법 + 3 김한국 2009-03-17 42559
49 아크맵상에서 shp파일 만드는 방법입니다. + 4 김영우 2009-01-29 45631
48 아크GIS 9.2 에서 dxf 파일을 불러오면 한글이 깨짐니다 + 5 정익선 2009-01-19 36346
47 데이타(data) 폴더 이동시 레이어가 안보일때 해결방법은 없나요? + 5 김도환 2008-12-10 33065
46 지번 19자리 설명 + 2 김한국 2008-11-06 44559
45 PolyLine을 Polygon으로 변환할때 속성정보까지 같이 변환 file + 4 송영기 2008-10-10 38638