2011년 3월 6일 일요일

Creating a slope using data in Geographic coordinates

경위도 좌표계를 가진 DEM으로 경사도 분석시 오류 대처 방안입니다.


▣ 원문 출처

▣ DEM 출처 및 샘플데이터
 - ASTER GDEM : http://www.gdem.aster.ersdac.or.jp/
 - 이 블로그 내 소개: ASTER Global Digital Elevation Map
 - 샘플데이터
  : 다운로드한 DEM Merge 후 제주도 부분 clip
  : ASTER_GDEM_JEJU.zip

▣ 샘플 DEM
 - Columns and Rows: 3148, 1892
 - Coordinate Reference System: GCS_WGS_1984(EPSG:4326)
 - Min X: 126.121410883696
 - Min Y: 33.094139773836
 - Max X: 126.99585532814
 - Max Y: 33.6196953293916
 - CellSize = 0.00027777778

▣ DEM -> Hillshade

▣ DEM ->Slope : Z factor = 1

▣ DEM -> Slope : Z factor = 0.000003

▣ DEM -> Slope : Z factor = 0.000008834378
 1. Determine what the middle latitude of the area of interest is.
 : middle latitude = (Minimum Y + Maximum Y) / 2.0
            = (33.094139773836 + 33.6196953293916) / 2.0
            = 33.3569175516138

2. Convert that degree value to radians. 
 : 1 degree = 0.0174532925 radians = PI / 180.0
 : input latitude in radians = 33.3569175516138 * (PI / 180.0)
                                     = 0.5821880395919463711183838467741
                         
3. Use the value in radians in the following equation:
 : Z factor = 1.0 / (113200 * cos(input latitude in radians)) 
 : cos(input latitude in radians) = 0.99994837648629650508861336534359
 :  ∴ Z factor = 1.0 / (113200 * cos(0.5821880395919463711183838467741)) 
               = 0.0000088343783231345251334838953833445