Have issue when using RDLS to perform limited-angle cone beam reconstruction
Nobody has claimed this yet.
Assessment
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Newbie friendliness
- 25/100
- Issue type
- Bug
- Clarity
- Needs clarification
- Activity status
- Stale
- Tech stack
- python
- Domain
- computer-vision
Research direction
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.
Description
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/360984)
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).
- Dominant language
- Cuda
- Stars
- 249
- Forks
- 37
- PR merge metrics
- No merged PRs in 30d
Contributor guide
No contributing guide indexed for this repository
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
More from llnl/LEAP
-
Difficulty 1/5 Under an hour Newbie friendliness 88/100
-
Difficulty 3/5 1-2 days Newbie friendliness 52/100
-
Difficulty 5/5 Over a week Newbie friendliness 35/100
-
Difficulty 3/5 1-2 days Newbie friendliness 38/100
-
Difficulty 1/5 Under an hour Newbie friendliness 15/100
Similar issues
-
Difficulty 1/5 Under an hour Newbie friendliness 92/100
open-compass/VLMEvalKit#1695 ·
-
Difficulty 1/5 1-3 hours Newbie friendliness 92/100
autowarefoundation/autoware_universe#13413 ·
-
Difficulty 2/5 1-3 hours Newbie friendliness 90/100
-
Difficulty 2/5 1-3 hours Newbie friendliness 78/100
alicevision/AliceVision#2180 · 1 comment ·
-
lk_single_level's det.abs() < f64::EPSILON degeneracy guard is on the wrong scale after #130/#138 Openbug help wanted video-module
Difficulty 2/5 1-3 hours Newbie friendliness 84/100