말씀하신 것처럼 NaN값이 있을 때는 흑백으로 이상하게 저장되는 것 같습니다.

아마 IDL> i.Save, ‘[파일이름].pdf’ 으로 저장하셨을 텐데

/BITMAP 키워드를 사용하면 해결되는 것 같습니다.

/BITMAP 키워드에 대한 설명은 IDL 도움말을 보시면 될 것 같고요…;;

eps 파일은 확인해보지 않았지만 pdf 파일은 이 방법으로 해결되는 것 같습니다.

올려주신 내용을 확인했습니다.
우선 INSIDE 함수는 Coyote 라이브러리를 사용하셨을 것으로 예상하여 해당 함수를 사용하였습니다.

확인 결과 값의 차이는 잘 모르겠으나 일부 영역에서 NaN 값이 나오는 것은 확인하였습니다.
이 부분은 아마도 내삽을 할 때 NaN 값의 영향을 많이 받는 부분은 전부 NaN 값으로 계산되는 것이 아닐까 추측합니다.
그래서 일부 영역은 값이 있었으나 주위의 NaN 값의 영향을 받아서 NaN 값으로 나온 것 같습니다.

해결 방법은 아니지만, 다른 방법을 제안하자면 아래와 같습니다.
작성하신 것과는 반대로 영역 내의 자료만 추출하여 내삽을 진행하는 방법을 사용하였습니다.
다만 이렇게 하면 외삽 부분이 발생하기 때문에, 다시 영역을 잘라 주어야 하며,
값의 일부는 외삽값이 됩니다.

지난 번과 마찬가지로 flag 옵션을 통해 살펴보시면 될 것 같습니다.
장단점을 살펴보자면, flag = 0 일 때는 NaN 값의 영향을 받는 모든 영역이 NaN이 됩니다(실제로 값이 있던 격자도). 하지만 값의 왜곡은 없어 보입니다.
flag = 1 일 때는 추출하고자 하는 영역에 대한 값을 구할 수 있습니다. 하지만 외삽된 값이 사용됩니다.

다만 값이 차이나는 경우는 저도 확인은 못해봤습니다(자료의 영향이 있을 수도 있을 것 같습니다.)

=============================================

PRO TEST_TRIGRID2

flag = 0

lat = [20:50:1]
lon = [120:150:1]

lonx = [128,138,140,146,142,130,128]
laty = [34,38,46,44,35,31,34]

lat_interpol = TRANSPOSE(REBIN(lat,31,31))
lon_interpol = REBIN(lon,31,31)

lat_cit = REFORM(lat_interpol,31*31)
lon_cit = REFORM(lon_interpol,31*31)

dtec_cit = HANNING(31,31)

IF flag EQ 0 THEN BEGIN
index = WHERE(INSIDE(lon_cit, lat_cit, lonx, laty) EQ 0 )
dtec_cit_div = dtec_cit
dtec_cit_div[index] = !VALUES.F_NAN
ENDIF ELSE BEGIN
index = WHERE(INSIDE(lon_cit, lat_cit, lonx, laty) EQ 1 )
dtec_cit_div = dtec_cit[index]
lon_cit_div = lon_cit[index]
lat_cit_div = lat_cit[index]
ENDELSE

IF flag EQ 0 THEN BEGIN
TRIANGULATE, lon_cit, lat_cit, tr
grid = TRIGRID(lon_cit, lat_cit, dtec_cit_div, tr, [0.1,0.1], [120,20,150,50],XGRID=xgrid,YGRID=ygrid)
ENDIF ELSE BEGIN
TRIANGULATE, lon_cit_div, lat_cit_div, tr
grid = TRIGRID(lon_cit_div, lat_cit_div, dtec_cit_div, tr, [0.1,0.1], [120,20,150,50],XGRID=xgrid,YGRID=ygrid,MISSING = !VALUES.F_NAN)
xgrid_interpol= REFORM(REBIN(xgrid,301,301),301L*301L)
ygrid_interpol= REFORM(TRANSPOSE(REBIN(ygrid,301,301)),301L*301L)
index = WHERE(INSIDE(xgrid_interpol, ygrid_interpol, lonx, laty) EQ 0)
grid[index] = !VALUES.F_NAN
ENDELSE

index = WHERE(INSIDE(lon_cit, lat_cit, lonx, laty) EQ 1)
gridSample = grid[0:300:10,*]
gridSample = gridSample[*,0:300:10]
gridSample = gridSample[index]
dataSample = dtec_cit[index]

diff = ABS(gridSample – dataSample)
PRINT, ‘Maximum different: ‘, MAX(diff[WHERE(FINITE(diff) EQ 1)])

i1 = IMAGE(dtec_cit,lon,lat,LAYOUT=[2,1,1],MIN_VALUE = 0,MAX_VALUE = 1,TITLE = ‘Origin’)
p1 = POLYLINE(lonx,laty,’r’,TARGET=i1,/DATA)
i2 = IMAGE(grid,xgrid,ygrid,LAYOUT=[2,1,2],MIN_VALUE = 0,MAX_VALUE = 1,TITLE = ‘Interpolation’,/CURRENT)

END

• 이 답변은 mwkim에 의해 8 years 전에 수정됐습니다.
• 이 답변은 mwkim에 의해 8 years 전에 수정됐습니다.