IDL图像旋转间接法双线性内插法

 

 

pro test

  ;***1. 输入旋转文件，可以修改
  file = 'C:\Users\ASUS\Desktop\作业\遥感程序设计\ch06-2作业图像旋转idl程序\jz512-rgb.tif'
  ;***1. 输入旋转角度，可以修改 
  angle = 20   

  ;***2. 读取 tiff 格式图像
  tm_data = READ_TIFF(file)     
  dims = size(tm_data, /DIMENSIONS)     ;图像数据大小，行，列 
  print,dims
  ;***3. 显示图像 
  ;定义显示窗口,和数据同样大的窗口
  window, 0, xsize = dims[1], ysize = dims[2], $  
  TITLE = 'Landsat TM data of Jiaozuo' 
  ;显示图像,注意加入order关键字 ，图像变正
  TV, tm_data,TRUE=1,/order 

  pi=3.14159265

  ;***4. 定义旋转矩阵,
  ;还记得旋转矩阵的使用吗？？？ 
  ;[x',y']=|cosA, sinA| * [x,y] 
  ;        |-sinA,cosA|

  Rot=[[cos(angle*pi/180),sin(angle*pi/180)],   $
       [-1*sin(angle*pi/180),cos(angle*pi/180)]]  
  
  ;***5. 定义四个角点坐标，[0,0]为原点 
  BT1=[0,0]             ;左上角点坐标
  BT2=[0,dims[1]-1]     ;右上角点坐标
  BT3=[dims[2]-1,0]     ;左下角点坐标
  BT4=[dims[2]-1,dims[1]-1] ;右下角点坐标
  
  BT1_rot=Rot           ;旋转左上角点坐标
  BT2_rot=Rot#BT2       ;旋转右上角点坐标
  BT3_rot=Rot#BT3       ;旋转左下角点坐标
  BT4_rot=Rot#BT4       ;旋转右下角点坐标
  
  x_min=min([BT1_rot[0],BT2_rot[0],BT3_rot[0],BT4_rot[0]])  ;旋转后图像最小x值
  y_min=min([BT1_rot[1],BT2_rot[1],BT3_rot[1],BT4_rot[1]])  ;旋转后图像最小y值
  x_max=max([BT1_rot[0],BT2_rot[0],BT3_rot[0],BT4_rot[0]])  ;旋转后图像最大x值
  y_max=max([BT1_rot[1],BT2_rot[1],BT3_rot[1],BT4_rot[1]])  ;旋转后图像最大y值
  
  width=x_max-x_min+1   ;旋转后图像宽度
  height=y_max-y_min+1  ;旋转后图像高度
  
  ;***6. 新建MxNx3的空矩阵，用于存储旋转后彩色图像
  img_rot=intarr(3,height,width)   
  
  ;***7. *********采用双线性内插法旋转图像********
  for i = 0, width-1 do begin
      for j = 0, height-1 do begin
      
          x_rot=i+x_min
          y_rot=j+y_min
          
          ; 坐标变换
          x = (x_rot*cos(angle*pi/180))+(y_rot*sin(angle*pi/180))
          y = (y_rot*cos(angle*pi/180)) - (x_rot*sin(angle*pi/180))
          
          ; 确定四个领域像素的坐标和权重
          x1=floor(x)
          x2=x1+1
          y1=floor(y)
          y2=y1+1
          
          wx2=x-x1
          wy2=y-y1
          wx1=1-wx2
          wy1=1-wy2
          
          ; 边界情况的处理
       if x1 gt 0 and x2 le dims[1]-1 and y1 gt 0 and y2 le dims[2]-1 then begin
      
          ; 四个领域像素取值
          p11=tm_data[*,x1,y1]
          p12=tm_data[*,x1,y2]
          p21=tm_data[*,x2,y1]
          p22=tm_data[*,x2,y2]
          
          ; 双线性内插计算 
          for k=0,2 do begin
              img_rot[k,j,i]=wx1*wy1*p11[k]+wx1*wy2*p12[k]+wx2*wy1*p21[k]+wx2*wy2*p22[k]
          endfor
       endif
     endfor
 endfor
  
  ;***8. 显示旋转后图像
  dimsrot= size(img_rot, /DIMENSIONS)
  window, 1, xsize = dimsrot[1], ysize = dimsrot[2], $ 
  TITLE = 'Rotation of Landsat TM data Jiaozuo'   
  TV, img_rot, TRUE=1   
        
  ;***9.保存旋转图像为tiff格式
  tif_rot_file = 'C:\Users\ASUS\Desktop\作业\遥感程序设计\ch06-2作业图像旋转idl程序\jz512-rgb-rotation-bilinear.tif'
  WRITE_TIFF, tif_rot_file, img_rot
  
  print,'end pro'
End

  ;***8. 显示旋转后图像
  dimsrot= size(img_rot, /DIMENSIONS)
  window, 1, xsize = dimsrot[1], ysize = dimsrot[2], $ 
  TITLE = 'Rotation of Landsat TM data Jiaozuo'   
  TV, img_rot, TRUE=1   
        
  ;***9.保存旋转图像为tiff格式
  tif_rot_file = 'C:\Users\ASUS\Desktop\ch06-2作业图像旋转idl程序\jz512-rgb-rotation-bilinear.tif'
  WRITE_TIFF, tif_rot_file, img_rot
  
  print,'end pro'
End