GIS/래스터 데이터 분석

< GIS

학습 목표

편집
  • 좌표 등록 과정
    • 아래 세 가지가 유사한 개념(용어)임.
    • 좌표등록 : image registration
    • 기하보정 : geometric correction. 예 : 영상을 보면 실제와 다르게 길이 쭈글쭈글하게 찍힐 수 있는데 이걸 펴줌.
    • 지리좌표화(참조) 과정 : georeferencing
  • 실제 원격 탐사 자료를 활용하여 기하보정하고, 벡터자료와 통합하기 위한 분석과정에 대해 배운다.

래스터 자료의 지리좌표화

편집

래스터를 실세계화. 행렬 형태의 래스터 좌표계를 지도 좌표계로 변경(참조)하는 것.

유사 용어

  • 영상 대 지도의 편위 수정(Image to map rectification)
  • 영상 대 영상 등록(image to image registration)

지상기준점

편집

ground control point(GCP). 영상에서 식별 가능하며 지도상에 정확히 위치되는 점.

원격탐사 자료 기하보정에 필수.

♣과정 알기

  1. GCP 쌍 추출
  2. 수학적 모델 결정(예 : 어핀 변환(Affine transformation)
  3. 모델의 기하변환계수 추출
  4. 기하보정

기하보정 종류

편집

영상 대 지도 보정(Image to map registration)

편집
  • 영상이 기하보정 대상(slave)
  • 지도가 참조자료(master)

영상의 직접 식별 가능한 점 좌표 을 어떤 함수 f를 통해 지도 좌표(x, y)로 바꾸는 게 영상 대 지도 보정.

이 두 좌표를 GCP 쌍이라 함.

기하보정 위해 GCP 수집하는 방법(지도 좌표 (x, y))

  • 하드카피 형태 평면지도
  • 수치평면지도
  • 수치정사사진 : 이미 좌표가 등록이 된 영상.
  • GPS이용 현장 수집.

영상 대 영상 등록(image to image registration)

편집

영상 좌표를 이미 좌표가 등록되어 있는 영상에 등록.

동일 지역 두 영상에서 동일한 물체들이 서로 같은 위치에 나타나도록 회전, 변환.

기하보정 과정

편집

  1. 입력화소 좌표와 동일 위치 지도 좌표 사이 기하학적 관계 식별(공간 내삽(spatial interpolation))
  2. 공간 내삽으로 보정된 출력화소의 밝기값을 결정하는 방법 결정(강도 내삽(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) : 모자이크된 영상 사이 경계선을 섞거나 흐릿하게 해줌. 사람 눈에 경계선이 잘 인식되지 않게 함.