MATLAB 게시판을 보니, INTERP2 함수를 다음과 같이 사용하는 것으로 나옵니다.

Vq = interp2(X,Y,V,Xq,Yq)

MATLAB을 잘 모르지만, 설명으로 보면, IDL에서 이런 일을 하는 함수는 GRIDDATA 가 아닌가 생각합니다. 위 Matlab 예제와 같은 변수명을 사용한다면,

Vq = GRIDDATA(X, Y, V, XOUT=Xq, YOUT=Yq)

의 형태로 사용할 수 있습니다. GRIDDATA에서 KRIGING 메쏘드를 쓸 수도 있는데, 언급하신 KRIG2D 함수를 사용한다면.

Vq = KRIG2D(V, X, Y, XVALUES=Xq, YVALUES=Yq)

이렇게 사용할 수 있습니다.

1024개의 밴드가 1 column 씩 왼쪽으로 이동되어 있는 상태라는 의미이신가요?
이 경우 맨 왼쪽 컬럼으로 있어야 할 데이터가 맨 오른쪽으로 가 있는 상태인가요? 그렇지 않다면, 맨 뒷쪽 밴드와 맨 앞쪽 밴드는 공유하는 면적이 너무 좁을 것이라고 예상됩니다.

근본적으로는 데이터가 왜 이렇게 생성되는지 파악하면 좋겠지만, 어쨌든 데이터가 어떻게 잘못 저장되어 있는지에대한 파악이라도 확실하다면, 말씀하신 대로 IDL에서 데이터 수정은 가능할 것 같습니다.
이 기능을 ENVI에 넣는 것도 가능합니다. 이를 위해서는 일단 IDL 문법에 익숙할 필요는 있습니다.

개발된 기능을 ENVI에 extension으로 넣는 것도 가능한데,
https://www.harrisgeospatial.com/docs/ToolboxExtensions.html
위 링크의 도움말을 보시면 일단 시작을 하실 수 있을 거라고 생각합니다.

제 생각에는 사용자의 Interaction이 별로 필요하지 않은 데이터 수정 작업 (-1)*(1 column * Band번호) 만큼 이동시키는 거라면 – 밴드번호만큼 왼쪽으로 이동 – 그냥 IDL 프로그램을 만들어 데이터마다 수정하는 게 빠를 것 같기는 합니다. 굳이 ENVI Extension을 만들 필요는 없지 않나 생각합니다.

보여주신 위 코드의 틀에서 어색한 점은 없어 보입니다.
20~30MB 정도는 요즘은 대용량이라고 보기는 어렵습니다. 그 문제는 아닌 듯 합니다.

1) 프로그램을 실행 할 때마다 문제가 발생하는 파일이 달라진다는 말씀이신가요? 예를 들면, 어떤 실행에서는 3월 20일 파일이 비정상적으로 기록되고, 다음 실행에서는 3월 20일은 멀쩡한데, 6월 20일이 비정상적으로 저장되고, 그런다는 의미이신가요? 그렇다면 분명 이상한 상황이기는 합니다.

2) 전체 상황을 잘 모르겠지만, close, 1 을 실행하지 못하면 쓰다 만 파일이 될 가능성이 매우 높습니다.
그런데, 이 경우를 가정해도 다음 반복실행에서 openw, 1 을 실행하면 “이미 열려 있는 파일입니다” 라는 에러가 발생합니다. 혹시 close 대신, /GET_LUN, FREE_LUN을 쓰셨다면 이부분을 살펴보셔야 할 것 같습니다만 올려주신 코드에서는 close를 쓰셨기 때문에 이 문제라고 보기는 좀 어렵죠.
그럼에도, “중간에 작성이 멈춘 것 같다” 라는 말씀은 close, 1을 실행하지 못할 때의 전형적인 상황이어서 일단 이 부분 아닐까 찍어 봅니다.

조금 더 정보가 필요한 상황입니다.

보여주신 3개의 문장 중, 첫번째 문장에 있는 hdf_sd_select() 함수는 IDL 함수로 IDL만 있으면 실행할 수 있습니다. sd_file 은 파일을 상징하는 번호이고, 0는 그 파일 안에서 SD 데이터셋 번호입니다. HDF 파일 안에는 여러개의 SD 데이터셋이 들어 있을 수 있습니다.
아마도 sd_file 변수는 HDF_SD_START() 함수의 리턴값일 것입니다. 소스코드에서 약간 앞쪽 부분을 찾아 보세요.

이후에 나오는 두 문장은 ENVI라는 영상처리 소프트웨어에서 지원하는 프로시저입니다. ENVI는 IDL로 만들어져서 IDL에서 기능 호출이 됩니다. 이를 ENVI API라고 하는데, ENVI가 설치되어 있는 IDL 환경에서 사용할 수 있습니다.

envi_write_envi_file, ji, outname=aa, r_fid=bb
과 같이 프로시저 형태로 사용되는 것이 맞을텐데요, ji 배열(일반적으로 영상입니다)을 aa라는 파일 이름으로 저장하라는 의미입니다. 프로시저 이름처럼, ENVI 표준 포맷으로 저장하고요, 저장된 이후 이 파일을 가리키기 위한 번호를 bb 변수로 받겠다는 의미입니다. (r_fid : Return File ID) bb 안에 어떤 숫자가 들어갈지는 프로그램이 실행될 때마다 다르지만 어쨌든 이 때 생성되는 파일을 가리킵니다.

envi_file_query, aa, ns=n1, nl=n2, nb=n3
과 같은 형태로 씌였을 것 같은데요, 프로시저 이름대로, ENVI 파일(보통 영상입니다, ENVI는 영상처리 플랫폼이니까요)의 속성을 알아내는 프로시저입니다. aa는 File ID입니다. 앞 문장의 r_fid 같은 것으로 받아내는 값이죠. 이후 ns= 이후에 나오는 변수에 # of Samples(가로크기), nl= 이후에 나오는 변수에 # of Lines(세로크기), nb= 이후 나오는 변수에 # of Bands (밴드 개수)를 받아냅니다.

N=7
data=bytarr(N, N)
data[0:*:2]=1
print, data
data=data*255

sz=N*60
img=rebin(data, sz, sz, /SAMPLE)
window, xsize=sz, ysize=sz
tv, img
end

일단 위 예제는 N이 홀수일 경우에만 정상적으로 보입니다.  프로그램 구조를 일단 단순화하는 것이 좋으므로 짝수일 경우는 배제하였습니다.

위 예제에서 print로 출력되는 부분은 다음과 같습니다.

   1   0   1   0   1   0   1
   0   1   0   1   0   1   0
   1   0   1   0   1   0   1
   0   1   0   1   0   1   0
   1   0   1   0   1   0   1
   0   1   0   1   0   1   0
   1   0   1   0   1   0   1

보이는 대로 일단 배열 자체를 체크 플래그로 만들고, 곱하기 255를 해서 흑색0과 백색255로 나누고 마지막으로 rebin을 통해 확대한 이미지로 변환하는 것이 이 예제의 흐름입니다.

처음 올려주신 예제가 repeat를 사용하기 위한 연습이라면 제가 잘못 답변을 드린 것인데요, 처음 올려주신 예제처럼 하시면 X와 Y가 함께 증가하게 되므로 대각선 방향만 원하시는 대로 생성될 수밖에 없을 것 같습니다.

아직 준비중인 사이트인데 첫 손님으로 오셨네요. 올리신 글이 어느날 없어지거나 할 수도 있는데요, 준비중이어서 그런 것이니 널리 양해 부탁드립니다. 일부러 지우거나 하는 그런 일은 없을 것입니다.

예상치 못했지만, 환영합니다.