[Rtk-users] ITK ERROR: Geometry is empty error when generating DRRs

[Suranga W]

Hi Simon,

Many thanks. I have another few questions though.

1.) How can I pass different gantry angles if I use the aforementioned
AddProjection method, which lacks the gantryAngle parameter? What I want is
to generate DRR images for angles 36, 72, 108, and so on.

2.) If I have the field of view origin from a corresponding kV X-ray dicom
image (let's say for gantry angle zero Field of View Origin at detector is
[116.853836060, 4.566089153] in pixels see
<https://dicom.innolitics.com/ciods/enhanced-xa-image/enhanced-xa-image-multi-frame-functional-groups/52009229/00189432/00187030>),
how can I use them to produce DRRs with aforementioned machine
configurations in AddProjection method ?



On Mon, Dec 12, 2022 at 8:43 AM Simon Rit <simon.rit at creatis.insa-lyon.fr>
wrote:

> Hi,
> You might have missed an earlier message
> WARNING: In
> /home/srit/src/itk/itk/Modules/Remote/RTK/src/rtkThreeDCircularProjectionGeometry.cxx,
> line 210
> DataObject (0x55d436d7ee10): Failed to AddProjection
> The problem is that the two vectors describing the projection orientation
> are not unit vector. You can correct this with:
> geometry.AddProjection([-60.0, 59.0 + sid, 34.0],[-60.0 + 115.,
> 59.0+sid-sdd, 34.0],[1./math.sqrt(2.),1./math.sqrt(2.),0.],[0.,0.,1.])
> Simon
>
> On Fri, Dec 9, 2022 at 1:43 PM Suranga W <
> isurusuranga.wijesinghe at gmail.com> wrote:
>
>> Hi,
>>
>> I attempted to generate DRRs by adding the geometry configuration to the
>> ThreeDCircularProjectionGeometry object, as shown in the attached code
>> (python based code). However, I encountered a runtime error as below
>> (please look at the end of the email).
>>
>> On the other hand, how can I pass the relevant gantry angle to it for
>> generating DRRs at different angles if I add geometry configuration using
>> the addProjection method as mentioned in my code?
>>
>> What I mean is that if I use the code below to set machine parameters,
>> I'm not able to add the corresponding gantry angle for generating DRRs.
>>
>> *geometry.AddProjection([-60.0, 59.0 + sid, 34.0],[-60.0 + 115.,
>> 59.0+sid-sdd, 34.0],[1.,1.,0.],[0.,0.,1.])*
>>
>> where [-60.0, 59.0, 34.0] represents the patient's isocenter.
>>
>> However, if I use the code below, I can use the gantry angle but I
>> cannot set the machine geometry configuration.
>>
>> *geometry.AddProjection(sid, sdd, gantryAngle, projOffsetX, projOffsetY,
>> outOfPlaneAngle, inPlaneAngle, sourceOffsetX, sourceOffsetY)*
>>
>> I attached a CT volume for your further reference.
>>
>> Please assist me in resolving this issue.
>>
>> ---------------------------------------------------------------------------RuntimeError                              Traceback (most recent call last)
>> Cell In [29], line 45     43 rei.SetInput(0, constantImageSource.GetOutput())     44 rei.SetInput(1, CT)---> 45 rei.Update()     47 Dimension = 3     48 OutputPixelType = itk.UC
>> RuntimeError: D:\a\im\include\rtkProjectionsRegionConstIteratorRayBased.hxx:88:
>> ITK ERROR: Geometry is empty, cannot determine iterator type.
>>
>>  drr_projection_with_diff_angles.txt
>> <https://drive.google.com/file/d/1Ags89Lpczvrk_rdHpGSqDcBrashzNSZr/view?usp=drive_web>
>>  sample_CT.nii.gz
>> <https://drive.google.com/file/d/1SGh0cjspFwTj2v8qTLXtcCfbUvyUENS5/view?usp=drive_web>
>>
>>
>>
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.creatis.insa-lyon.fr/pipermail/rtk-users/attachments/20221212/8a4e53c4/attachment.htm>