이상우의 IDL 블로그

IDL에서 DEM(Digital Elevation Model) 자료 파일을 읽고 처리 및 표출까지 해보는 방법을 소개하고 있습니다. 먼저 지난 회차 게시물에서는 한반도 지역에 대한 ASTER DEM V003 자료 파일을 예제로 사용하면서 이 파일을 읽고 이미지의 형태로 기본적인 표출만 하는 과정까지 진행해보았습니다. 오늘은 DEM 데이터를 이미지의 형태로 한반도 지도상에 중첩하여 표출하는 작업을 진행해보기로 합니다. 사용할 DEM 자료 파일은 지난 회차에서 사용했던 것과 동일합니다. 일단 READ_TIFF 함수를 사용하여 이 파일을 읽어서 배열로 가져오는 과정까지는 지난번과 동일하게 진행합니다. file = 'ASTGTMV003_N36E126_dem.tif'data = READ_TIFF(file, OR..

DEM(Digital Elevation Model)은 수치표고모델이라고도 하는데 식생 및 인공구조물 등을 제외한 지형 자체만의 표고(Elevation) 값을 뜻합니다. 그리고 DEM 자료는 경도 및 위도 방향으로 일정한 간격으로 표고 값을 측정한 자료로서 원격탐사, 기상, 지질 등 다양한 분야에서 널리 사용되는 자료이기도 합니다. 주로 인공위성들에 의하여 측량되는 이러한 데이터는 그 종류가 다양하며 관측 위성의 종류 및 관측의 해상도에 따라 세분화됩니다. 이러한 자료들은 일반에게 공개적으로 배포되는 경우도 많지만 그렇지 않은 경우도 있습니다. 물론 IDL에서도 DEM 자료를 읽고 처리 및 표출하는 것이 가능합니다. 그래서 이번 기회에 그러한 방법을 예제와 함께 소개해보고자 합니다. 먼저 어떤 DEM 자료..

IDL의 NG 체계에서 지원되는 그래픽 함수들은 대부분 MARGIN 또는 POSITION 속성을 지원합니다. 이러한 속성들은 그래픽 창 내에서 개별 그래픽 개체의 여백 또는 위치를 설정하는 역할을 합니다. MAP 함수의 경우도 마찬가지로 이러한 속성들을 지원합니다. 다만 MAP 함수의 경우는 MARGIN 또는 POSITION 속성을 사용하여 여백 또는 위치를 설정하는데 있어서 다른 그래픽 함수들과는 약간 다른 특성이 있기 때문에 유의가 필요합니다. 오늘은 이와 관련된 내용을 소개해보기로 합니다. 먼저 다음과 같이 WINDOW 함수로 그래픽 창을 띄우고 MAP 함수로 지도를 표출하는 예제를 살펴봅시다. win = WINDOW(DIMENSION=[600, 600], /NO_TOOLBAR) m = MAP('G..

IDL에서 지도를 표출하면서 국경선을 함께 표시해야 하는 경우가 종종 있습니다. 그러한 표출을 위해서는 NG 체계 기준으로는 MAPCONTINENTS 함수에서 /COUNTRIES 키워드를 사용하면 됩니다. 이러한 방법을 사용하여 동남아시아 지역을 국경선과 함께 표시하는 예제를 보면 다음과 같습니다. win = WINDOW(DIMENSIONS=[600, 600], /NO_TOOLBAR) limit = [0, 90, 50, 140] m = MAP('Geographic', LIMIT=limit, $ FILL_COLOR='light blue', ASPECT_RATIO=0, $ LABEL_POSITION=0, FONT_SIZE=11, CLIP=0, $ MARGIN=0.1, /CURRENT) mc = MAPCONT..

LCC(Lambert Conformal Conic) 투영법 기반의 지도를 표출하는 방법에 관해서는 제가 예전에도 관련 게시물들을 통하여 몇번 다뤄본 적이 있습니다. 또한 LCC 투영법 기반으로 격자 분포를 하는 2차원 데이터를 지도와 함께 중첩 표출하는 예제들도 관련 게시물들을 통하여 소개했던 적이 있습니다. 따라서 LCC 투영법 기반의 지도를 표출하는 방법에 관한 기본적인 내용을 굳이 반복할 생각은 없지만, 추가적으로 다뤄보고 싶은 관련 주제가 있어서 오늘 소개해보고자 합니다. IDL의 NG 체계에서 LCC 투영법으로 지도를 표출할 때 약간 골치아프게 다가오는 이슈가 있습니다. 바로 지도의 외곽 부분에 대한 처리인데요. LCC 투영법 자체의 특성상 경도 및 위도 격자선이 곡선 또는 기울어진 선이 됩니다..

지난 회차에 이어서 오늘은 2)에 관한 내용, 즉 지도상에 2차원 데이터를 중첩 표출하면서 대륙 부분에서만 데이터가 보이도록 처리하는 방법을 소개해보겠습니다. 지난 회 게시물에서 잠시 언급했듯이 그 반대의 경우인 1)은 구현 방법이 꽤 간단하지만 2)는 좀 복잡합니다. 1)의 경우는 지도와 2차원 데이터의 중첩 처리를 한 이후에 대륙경계선을 표시하기 위하여 MAPCONTINENTS 함수를 사용할 때 FILL_COLOR 속성을 사용함으로써 대륙 부분을 색상으로 채우는 방식의 처리만 하면 됩니다. 꽤 간단하죠. 그런데 2)의 경우는 그냥 단순하게 생각해보면 MAPCONTINENTS 함수에서 대륙 대신 해양을 색상으로 채우는 기능만 있으면 간단하지않을까하는 생각을 가져볼 수 있습니다. 그렇게만 된다면야 매우 ..

IDL에서 지도를 표출하고 그 위에 2차원 데이터를 중첩하여 표출하는 방법 및 예제들은 제가 이 블로그를 통하여 여러 차례 소개해드린 바 있습니다(관련 링크). 그런데 오늘 소개하고자하는 내용은 이와 같이 중첩 표출을 하면서 대륙 또는 해양에 해당되는 부분만 보이도록 하는 방법에 관한 것입니다. 예를 들면 지도상에 2차원 데이터를 등위선(Contour)의 형태로 중첩표출을 하되, 왼쪽은 대륙을 가리고 해양 부분만 보이도록 처리하는 경우 또는 반대로 해양을 가리고 대륙 부분만 보이도록 처리한 경우입니다. 물론 등위선 대신 이미지가 될 수도 있을 것입니다. 어쨌든 실제로 이와 같은 표출을 하고싶은 경우가 분명히 있습니다. 그래서 이번 기회에 그 방법을 소개해보고자 합니다. 소개할 순서는 다음과 같습니다. 1..

IDL에서 표출된 지도에서 라벨 문자들이라면 당연히 경위도 격자들에 대응되는 경위도 값들을 뜻합니다. 지도를 표출하게 되면 어떤 투영법에서든간에 경위도 격자선 및 경위도 값 문자들이 함께 표시되는 것이 기본입니다. 오늘은 이렇게 지도상에서 표시되는 라벨 문자들에 대하여 그 형식(Format)을 제어하는 방법을 관련 예제들과 함께 소개해보고자 합니다. 일단 지도부터 먼저 표출해놓고 시작합니다. win = WINDOW(DIMENSIONS=[800, 700], /NO_TOOLBAR) limit = [33, 121, 41, 131] m = MAP('Lambert Conformal Conic', LIMIT=limit, $ STANDARD_PAR1=30, STANDARD_PAR2=60, $ CENTER_LONGIT..

GOES 16 및 17은 미국 NASA/NOAA에서 운용중인 기상관측용 정지궤도 인공위성들입니다. GOES 16은 GOES-R 또는 GOES East라고도 부르며 경도 위치가 75.2W이고, GOES-17은 GOES-S 또는 GOES West라고도 부르며 경도 위치는 137.2W입니다. 따라서 아메리카 대륙 및 대평양 지역에 대한 기상 연구에 있어서 GOES 16 및 17의 데이터가 자주 사용됩니다. 오늘은 이 데이터 파일을 읽고 표출하는 예제를 살펴보고자 합니다. 어차피 데이터의 종류가 매우 많지만 그 중에서 GOES 17의 ABI(Advanced Baseline Image)에서 관측되는 복사량(Radiance) 데이터를 예제 데이터로 사용해보겠습니다. GOES 16/17의 데이터를 받을 수 있는 경로..

천리안 2A 위성의 영상 데이터를 표출하는 방법과 예제를 소개하는 게시물들을 제가 지난 해(2020)에 2회에 걸쳐서 게재했던 바 있습니다. 천리안 위성 2A호(GK2A)의 데이터 표출 [1] 천리안 위성 2A호(GK2A)의 데이터 표출 [2] 여기서는 적외선 영상인 IR 11.2 um 데이터를 담고 있는 netCDF 포맷의 파일을 사용하였습니다. 그리고 영상 데이터의 영역 커버리지가 전구(Global), 동아시아, 한반도 등 세가지 종류로 제공되는데, 동아시아 및 한반도 영역 데이터를 대상으로 하였고 전구 맵(Global Map) 데이터의 경우는 보류했었는데요. 오늘은 이 전구 맵 데이터를 대상으로 표출 방법 및 예제를 소개해보고자 합니다. 이번에 예제로 사용할 데이터 파일의 이름은 다음과 같습니다. ..