- This topic has 5개 답변, 3명 참여, and was last updated 8 years, 11 months 전에 by
.
-
글
-
-
안녕하세요, 항상 게시판을 보며 많은 것을 배워가고 있습니다.
감사합니다.이번엔 제가 질문이 있어서 게시판을 이용하게 되었네요.
현재 제가 사용하는 data는 정지궤도위성 적외영상 아시아 영역과 모델장입니다.
위성영상 data는 1544*1934이고, 아시아 영역만 data가 있고 나머지는 nan값입니다.
위성영상의 격자는 scanline을 따라 grid가 만들어져서 위경도가 균일하지 않습니다.
동시에 사용하고자 하는 모델장은 해상도가 0.5도*0.5도로, data는 720*361 입니다.
저는 모델장을 위성영상 도메인으로 확장시켜 사용하려고 하는데 여기서 자꾸 부딪히네요.
몇주정도 고민하다가 질문 올려봅니다.제가 시도해본 방법은 다음과 같습니다.
1. rebin이나 congrid를 사용
– rebin은 정수배가 아니라서 안되고, congrid는 확장은 되지만 불규칙적인 위성의 격자에 맞춰 확장되지 않을 것이라서, 지푸라기 잡는 심정으로 해보긴했습니다만 안되더군요.
2. sht_cat은 linear한 데이터라고 오류가 뜨네요. (griddata도 sph_scat과 비슷한 regular data에서 irregular data로 맞추는거라 그런지 시도해도 잘안됩니다 ㅠㅠ)
3. 그래서 현재 그냥 720*361번 for문을 돌려서 모델장의 각 위경도 격자에 들어오는 모든 위성영상 격자에 모델값을 할당하는 코드를 짜서 임시로 쓰고 있습니다만… 문제점은 해상도가 낮은 모델장의 네모난 위경도안에 들어올 시에 위성영상 위경도에 모델값을 주다보니 grid는 고해상도인데 내삽이 안되서 계단처럼 값이 할당되고, 무엇보다 for문이라 시간이 너무 오래걸린다는 점이네요. 해결할 수 있는 방법이 있을까요? -
완벽한 해결책이 될지는 모르겠으나, 비슷한 이슈에 관하여 제 IDL 블로그에 관련 게시물을 올린 적이 있습니다. 이 내용을 일단 참조해보시면 어떨까 합니다. 링크는 아래에 있습니다.
-
-
-
griddata의 이름이 (irregular를 2차원배열로 만드는 것이어서) 안될 것처럼 보이지만, 2차원 배열을 irregular처럼 사용하고, XOUT, YOUT 키워드를 이용해서 알고자 하는 지점(들)의 X,Y좌표를 명시하면 해당 좌표값만 계산해 냅니다. 예를 들어, 우리나라 전역의 기상예보격자 기온 같은 것이 있다면(동네예보), 이것을 입력자료로 써서, 원하는 위경도(들)의 기온을 유추(내삽)할 수 있다는 것입니다.
도움말을 보시면 확인되는데, Minimum Curvature 알고리즘과 Quintic 알고리즘은 무조건 격자형으로만 리턴을 하고, 나머지 알고리즘(IDW, KRIGING, Nearest Neighbour 등등)은 모두 XOUT, YOUT 키워드를 사용할 수 있습니다.
-
-
