GIS/래스터 데이터 분석
학습 목표
[편집]- 좌표 등록 과정
- 아래 세 가지가 유사한 개념(용어)임.
- 좌표등록 : image registration
- 기하보정 : geometric correction. 예 : 영상을 보면 실제와 다르게 길이 쭈글쭈글하게 찍힐 수 있는데 이걸 펴줌.
- 지리좌표화(참조) 과정 : georeferencing
- 실제 원격 탐사 자료를 활용하여 기하보정하고, 벡터자료와 통합하기 위한 분석과정에 대해 배운다.
래스터 자료의 지리좌표화
[편집]래스터를 실세계화. 행렬 형태의 래스터 좌표계를 지도 좌표계로 변경(참조)하는 것.
유사 용어
- 영상 대 지도의 편위 수정(Image to map rectification)
- 영상 대 영상 등록(image to image registration)
지상기준점
[편집]ground control point(GCP). 영상에서 식별 가능하며 지도상에 정확히 위치되는 점.
원격탐사 자료 기하보정에 필수.
♣과정 알기
- GCP 쌍 추출
- 수학적 모델 결정(예 : 어핀 변환(Affine transformation)
- 모델의 기하변환계수 추출
- 기하보정
기하보정 종류
[편집]영상 대 지도 보정(Image to map registration)
[편집]- 영상이 기하보정 대상(slave)
- 지도가 참조자료(master)
영상의 직접 식별 가능한 점 좌표을 어떤 함수 f를 통해 지도 좌표(x, y)로 바꾸는 게 영상 대 지도 보정.
이 두 좌표를 GCP 쌍이라 함.
기하보정 위해 GCP 수집하는 방법(지도 좌표 (x, y))
- 하드카피 형태 평면지도
- 수치평면지도
- 수치정사사진 : 이미 좌표가 등록이 된 영상.
- GPS이용 현장 수집.
영상 대 영상 등록(image to image registration)
[편집]영상 좌표를 이미 좌표가 등록되어 있는 영상에 등록.
동일 지역 두 영상에서 동일한 물체들이 서로 같은 위치에 나타나도록 회전, 변환.
기하보정 과정
[편집]♣
- 입력화소 좌표와 동일 위치 지도 좌표 사이 기하학적 관계 식별(공간 내삽(spatial interpolation))
- 공간 내삽으로 보정된 출력화소의 밝기값을 결정하는 방법 결정(강도 내삽(intensity interpolation))
공간 내삽
[편집]어핀 변환(Affine transformation)
대표적인 좌표변환식. =일차 선형변환.
영상 간 x, y 방향 이동, 축척(scale) 변화, 휨(nonorthogonality), 회전(rotation) 고려.
영상 내 왜곡(오차) 제거 가능. but 왜곡 원인 분석은 안 함.
위 여섯 계수(변환 계수)들에 대한 비선형 결합식을 선형결합식으로 표현하는 것.
영상 좌표를 x, y라 하고, 지도 좌표를 X, Y라 하면
GCP 1개 당 방정식 2개, 미지수 6개
GCP 3개면 방정식 6개, 미지수 6개. 풀이 가능.
그러나 GCP는 오차를 포함하고 있으므로 4개 이상 GCP를 이용하여 최소제곱법으로 푼다.
A가 정사각행렬이 아니기 때문에 역행렬을 구할 수 없다. 따라서 최소제곱법으로
다항식 차수가 증가하면 GCP 모델링이 정교해지지만, GCP 존재하지 않는 먼 거리의 지역은 오차가 오히려 커지기 때문에 1차식으로 하는 게 좋다.
순방향 매핑(forward mapping)
[편집]GIS 분야에서는 사용하나, 원격탐사 자료(래스터 자료) 공간내삽에는 쓰지 않음.
영상좌표를 변환해줬을 때 지도좌표 격자에 딱 맞지 않고 애매하게 걸리는 문제 발생.
역방향 매핑(backward mapping)
[편집]어핀변환을 거꾸로 만드는 것.(지도좌표가 영상좌표로 변환되도록) 역시 딱 격자에 맞진 않는데, DN값을 가지고 결정하여 빈칸없이 지도좌표계를 생성 가능.
강도 내삽과 결합하여 효과적으로 기하보정 가능.
강도 내삽
[편집]map에서 영상으로 변환했는데 그 위치에 화소가 없는 경우 재배열하는 것.
종류
- 최근린 내삽법(nearest neighbor) : 가장 가까운 화소 밝기값으로 보간.
- 공일차 내삽법(bilinear interpolation) : 2개의 직교하는 방향에 있는 4개 화소 밝기값의 거리가중 평균.
- 직접 표로 계산해보기.
- 입방회선법(bicubic interpolation, cubic convolution)
- 가장 가까운 16개 화소 밝기값 거리 가중 평균.
- 계산 방식은 동일.
오차 계산
[편집]♣♣
각 GCP에 대해 RMS(Root Mean Square) 오차 계산.
- ' 붙은 좌표 : 변환된 영상 좌표
- ' 없는 좌표 : 지도 좌표
영상에서 GCP의 정확도를 나타냄. 어떤 GCP가 큰 오차 내는지 알 수 있다.
응용
[편집]- 영상 모자이킹 : 여러 장 영상을 이어붙여서 파노라마 사진처럼 광대한 지역의 연속적인 전경 만듦.
- 히스토그램 매칭 : 날씨에 따른 색상 및 명암 차이 보정하는 것.
- 형상화(Feathering) : 모자이크된 영상 사이 경계선을 섞거나 흐릿하게 해줌. 사람 눈에 경계선이 잘 인식되지 않게 함.