- This topic has 4개 답변, 2명 참여, and was last updated 7 years, 8 months 전에 by
.
-
글
-
-
안녕하세요. 성균관대 조경학과 박사과정 임중빈입니다.
다름이 아니라 IDL로 ENVI band math 기능을 반복적으로 활용하는 방법에 대해 문의드리고 싶습니다.
제가 인터넷에서 검색을 해 봤으나 해결하지 못해서 이렇게 문의드립니다. 죄송합니다.
제가 가지고 있는 영상이 0~255 DN 값을 갖는 영상입니다.
그래서 ((x*0.0078125)-1) 수식을 적용하여 -1부터 1까지의 DN 값을 갖는 영상으로 변환하고 싶습니다.
파일명들은 JAPAN831221-30_NDVI.tif, JAPNA831211-20_NDVI.tif, JAPAN831201-10_NDVI.tif, etc. 로
JAPAN**(년도)****(날짜)-**(10,20,30 일 합성)_NDVI.tif 의 규칙으로 반복이 됩니다.
이 때 IDL 명령어를 어떻게 작성해야 할까요?
제가 인터넷 찾아보고 작성해본 명령어는 아래와 같습니다.
e=envi()
cd, ‘C:\downloads\AVHRR\AVHRR_NDVI\AVHRR_NDVI\1983’files1983=file_search(‘JAPAN83*NDVI.tif’, count=ct1983)
for i = 0, ct1983-1 do $
Raster=e.openraster(Files1983[i])
OUTFILE = Files1983[i]+’R.tif’
fid = ENVIRasterToFID(raster)
ENVI_File_Query, fid, DIMS=dims
ENVI_Doit, ‘Math_Doit’, $
FID = [fid, fid], $
DIMS = dims, $
POS = [0], $
EXP = ‘(b1*float(0.0078125))-float(1)’, $
OUT,NAME = OutFileenvi_batch_exit
end
오류가 발생하고 진행이 되지를 않았습니다. 어떻게 수정해야 할까요?
그리고 년도별 영상의 평균 영상을 제작하고 싶은데요
file_search와 같이 자동으로 년도를 검색하고, 그 파일들의 평균 파일을 구하는 명령어가 있을까요?
제 질문을 요약하면 다음과 같습니다.
Q1. JAPAN******-**_NDVI.tif 의 이름을 갖는 모든 영상의 DN값 변환
Q2. 년도별 NDVI 평균값 계산 반복
저도 계속 찾아보겠습니다. 감사합니다.
-
안녕하세요, 이제 온라인에서도 뵙게되었네요 ㅎㅎ 반갑습니다.
위 내용과 관련하여 에러 메시지도 알려주셨으면 어떤 문제인지 더 정확히 파악할 수 있을 것 같은데요, 우선은 알려주신 코드로 실행해본 후 몇가지 문제에 대해 말씀드립니다.– 먼저 for 문 이후 begin 이 실행되지 않았습니다. i 가 순차적으로 증가할 때마다 outfile까지 처리가 되어야 할 것 같은데요, 위 경우라면 e.openraster가 한번에 처리되고 이후, outfile 지정하는 부분으로 넘어가는데 문제가 생깁니다.
– out,name은 단순 오타인지, 정말로 저렇게 사용하신 것인지 모르겠습니다. out_name으로 쓰셔야 합니다.
– 각 영상을 호출하고 단일 영상에서 각 픽셀에 대한 값을 변경하는 것이므로 fid는 1개이면 됩니다.
심각한 문제는 아니지만 조금 더 말씀드리면,
– outfile 이름도 저 상태로라면 다소 이상하게 저장됩니다. files1983[i]에서 이미 .tif가 들어가 있기 때문에 “원본.tifR.tif”가 될 것 같습니다.
– 0.0078125는 이미 실수형이기 때문에 float으로 감싸지 않으셔도 됩니다.위 내용으로 수정한 예시 코드를 추가합니다.
compile_opt idl2
e=envi()files1983=file_search(‘JAPAN83*NDVI.tif’, count=ct1983)
for i = 0, ct1983-1 do begin
Raster=e.openraster(Files1983[i])
temp_name = strmid(files1983[i], 0, 19)
OUTFILE = temp_name+’_R’fid = ENVIRasterToFID(raster)
ENVI_File_Query, fid, DIMS=dims
ENVI_Doit, ‘Math_Doit’, FID = [fid], DIMS = dims, POS = [0], EXP = ‘b1*0.0078125-1.’, OUT_NAME = OutFile
endfor -
-
-
