Have issue when using RDLS to perform limited-angle cone beam reconstruction

Nobody has claimed this yet.

Hi Kyle:
Recently I want to use iterative algorithm to perform limited-angle cone beam reconstruction. The specfication is as follows:

proj size is (984 ,256, 821) # view, row, column
proj = proj[656:984, :,:] #choose 328 view from 984 view
self.numRows = 256
self.numCols = 821
self.sod = 625.61
self.sdd = 1097
self.pixelSizex = 1.09
self.pixelSizey = 0.8224 * 1097 / 821
self.start_angle = 180
self.end_angle = self.start_angle + ang_num
self.numAngles = round(ang_num/360*984)
self.voxelsize = 0.9
self.angles = np.linspace(self.start_angle, self.end_angle, self.numAngles, dtype=np.float32)
leapct.set_conebeam(self.numAngles, self.numRows, self.numCols, self.pixelSizex, self.pixelSizey, 0.5*(self.numRows-1), 0.5*(self.numCols-1), self.angles, self.sod, self.sdd)
leapct.set_curvedDetector()
leapct.set_volume(img_size, img_size, self.numRows, voxelWidth = self.voxelsize, voxelHeight = self.voxelsize) #numX, numY, numZ, voxelWidth=None, voxelHeight=None
rec = copy.deepcopy(init)
rec = rec.squeeze(1)
filters = filterSequence()
filters.append(TV(leapct, delta=0.01/100.0, p=1.0))
leapct.ASDPOCS(proj,rec,iter_num,50,10,filters)
leapct.RDLS(proj,rec,iter_num,filters)

And the ASDPOCS is works fine(left image), but the RDLS result is not good contain lots of noise(right).


Start by reproducing the reported Python API calls with the supplied projection dimensions, limited-angle slice, cone-beam geometry, and TV filter. Compare the RDLS result with ASDPOCS using the same reconstruction and iteration settings. The issue names no source file or test; done would require identifying why RDLS produces noise for this setup and adding a verified fix or explanation.

Written by the indexing model from the issue text.