arbitraryrotation.pro

IDL开发源码

PRO arbitraryRotation
; 选择文件并合作使用己知维数读入数据
file=FILEPATH('m51.dat',SUBDIRECTORY=['examples','data'])
image =READ_BINARY (file,DATA_DIMS =[340,440])
; 准备显示器并装载一个黑白颜色表
DEVICE ,DECOMPOSED =0,RETAIN =2
LOADCT,0
;创建窗口显示原始图像
WINDOW ,0 ,XSIZE=340 ,YSIZE=440
TVSCL,image
; 将新图像顺时针旋转 33 度并将图像缩小 50%及用内插法防止数组边界溢出
arbitraryImg=ROT(image,33,.5,/INTERP,MISSING=127)
; 显示旋转的图像
WINDOW ,1 ,XSIZE=340,YSIZE=440
TVSCL,arbitraryImg
END