이상우의 IDL 블로그

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) 데이터의 경우는 보류했었는데요. 오늘은 이 전구 맵 데이터를 대상으로 표출 방법 및 예제를 소개해보고자 합니다. 이번에 예제로 사용할 데이터 파일의 이름은 다음과 같습니다. ..

IDL에서 지도를 표출하는데 있어서 여러가지 투영법(Projection)들이 존재하고 주요 투영법들에 관해서는 제가 이 블로그를 통하여 여러 차례 소개해드린 바 있습니다. 오늘은 지도 투영법들 중에서 GOES-R이라는 투영법을 사용하여 지도를 표출하는 방법 및 예제를 소개해보고자 합니다. 표출에 있어서는 NG 체계에서 MAP 함수를 사용하는 방법을 사용할 것입니다. GOES-R 투영법은 마치 지구로부터 일정 거리만큼 떨어진 상태에서 지구 구체를 바라보는 듯한 느낌으로 표출하는 방식이라는 점에 있어서는 제가 얼마전에 소개했던 NSP(Near Side Perspective) 투영법과 비슷한 측면도 있습니다. 하지만 그 명칭에서 알 수 있듯이 이 투영법은 GOES-R 위성의 위치 정보와 유사한 특성을 갖도록 ..

IDL에서 지도를 표출하는데 있어서 여러가지 투영법(Projection)들이 존재하고 주요 투영법들에 관해서는 제가 이 블로그를 통하여 여러 차례 소개해드린 바 있습니다. 오늘은 지도 투영법들 중에서 Near Side Perspective라는 투영법을 사용하여 지도를 표출하는 방법 및 예제를 소개해보고자 합니다. 표출에 있어서는 NG 체계에서 MAP 함수를 사용하는 방법을 사용할 것입니다. Near Side Perspective 투영법은 마치 지구로부터 일정 거리만큼 떨어진 상태에서 지구 구체를 바라보는 듯한 느낌으로 표출하는 방식인데요. 예전에 제가 DG 체계의 MAP_SET 명령에서 지원되는 Satellite라는 투영법에 관하여 소개하는 게시물을 올린 적이 있는데 그 투영법과 비슷합니다. 그리고 용어..

특정 지역에 대한 2차원적인 데이터가 있을 때, 먼저 그 지역의 지도를 표출하고 그 위에 2차원 데이터를 중첩하는 작업을 해야하는 경우는 꽤 많습니다. 통상적으로 2차원 데이터는 주로 이미지(Image) 또는 등위선(Contour)의 형태로 표출할 수 있습니다. 따라서 지도와 이미지를 중첩하는 경우 그리고 지도와 등위선을 중첩하는 경우로 나눌 수 있습니다. 사실 IDL에서 지도/이미지의 조합으로 중첩 표출을 하는 방법이나 예제는 제가 이 블로그에서 여러번 관련 게시물들을 통하여 소개한 바 있습니다. 반면 지도/등위선의 조합으로 중첩 표출하는 방법에 관해서는 상대적으로 언급이 좀 덜 되었던 측면이 있습니다. 그래서 오늘은 전반적으로 정리하는 의미에서, 지도/이미지의 조합과 지도/등위선의 조합을 함께 소개하..