chapter09backprojecthoughandradon.pro

来自「IDL可视化工具入门与提高(PDF)及其随书代码」· PRO 代码 · 共 40 行

; Chapter09BackprojectHoughAndRadon.pro
PRO Chapter09BackprojectHoughAndRadon
  imageSize = [64, 64]
  file = FILEPATH('abnorm.dat', SUBDIRECTORY = ['examples', 'data'])
  image = READ_BINARY(file, DATA_DIMS = imageSize)
  displaySize = 4*imageSize
  offset = displaySize/3
  DEVICE, DECOMPOSED = 0
  LOADCT, 0
  houghTransform=HOUGH(image,RHO=houghRadii,THETA=houghAngles,/GRAY)
  WINDOW, 0, XSIZE = displaySize[0] + 1.5*offset[0], $
    YSIZE = displaySize[1] + 1.5*offset[1], TITLE = 'Hough Transform'
  TVSCL, CONGRID(houghTransform, displaySize[0], $
    displaySize[1]), offset[0], offset[1]
  PLOT, houghAngles, houghRadii, /XSTYLE, /YSTYLE, $
    TITLE = 'Hough Transform', XTITLE = 'Theta', $
    YTITLE = 'Rho', /NODATA, /NOERASE, /DEVICE, $
    POSITION = [offset[0], offset[1], displaySize[0] + offset[0], $
    displaySize[1] + offset[1]], CHARSIZE = 1.5
  radonTransform = RADON(image, RHO = radonRadii, $
    THETA = radonAngles, /GRAY)
  WINDOW, 1, XSIZE = displaySize[0] + 1.5*offset[0], $
    YSIZE = displaySize[1] + 1.5*offset[1], TITLE = 'Radon Transform'
  TVSCL, CONGRID(radonTransform, displaySize[0], $
    displaySize[1]), offset[0], offset[1]
  PLOT, radonAngles, radonRadii, /XSTYLE, /YSTYLE, $
    TITLE = 'Radon Transform', XTITLE = 'Theta', $
    YTITLE = 'Rho', /NODATA, /NOERASE, /DEVICE, $
  POSITION = [offset[0], offset[1], displaySize[0] + offset[0], $
    displaySize[1] + offset[1]], CHARSIZE = 1.5
  backprojectHough = HOUGH(houghTransform, /BACKPROJECT, $
    RHO = houghRadii,THETA=houghAngles,NX=imageSize[0],NY=imageSize[1])
    backprojectRadon = RADON(radonTransform, /BACKPROJECT, $
  RHO = radonRadii,THETA=radonAngles,NX=imageSize[0],NY=imageSize[1])
  WINDOW, 2, XSIZE = (3*displaySize[0]), YSIZE = displaySize[1], $
  TITLE = 'Hough and Radon Backprojections'
  TVSCL, CONGRID(image, displaySize[0], displaySize[1]), 0
  TVSCL, CONGRID(backprojectHough, displaySize[0], displaySize[1]), 1
  TVSCL, CONGRID(backprojectRadon, displaySize[0], displaySize[1]), 2
END