<div dir="ltr">I am trying to extract a subtile from a vtkimage using the following code:<div><br></div><div><div>def GetTile(self,image,x1,y1,x2,y2):</div><div>    i = vtk.vtkImageData().NewInstance()</div><div>    i.SetDimensions(x2-x1+1,y2-y1+1,1)</div><div>    i.AllocateScalars(vtk.VTK_UNSIGNED_CHAR,3)</div><div>    for x in range(x1,x2+1):</div><div>      for y in range(y1,y2+1):</div><div>        i.SetScalarComponentFromDouble(x-x1,y-y1,0,0,image.GetScalarComponentAsDouble(x,y,0,0))</div><div>    return i</div><div><br></div><div><br></div><div>x1,y1 being the upper left corner of the subimage</div><div>x2,y2 being the lower right corner of the subimage</div><div><br></div><div>and it appears to work when I view the subtile using:</div><div><div>    data = self.GetTile(myvtkimage,0,0,256,256).GetPointData()<br></div><div>    data = numpy_support.vtk_to_numpy(data.GetArray(0))</div><div>    d = data[:,0]</div><div>    dd = d.reshape(256,256)</div><div>    iii = sitk.GetImageFromArray(dd)</div><div>    sitk.WriteImage(iii,"/home/erich/subtile.jpg")</div></div><div><br></div><div>But, they don't appear the same because a segmentation program that works fine on the original image doesn't work on the subtiled image (even using the same original dimensions)</div><div><br></div><div>Any thoughts?</div><div><br></div><div><br></div><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>==========================================================</div><div>Erich Bremer, M.Sc.</div><div>Director for Cyberinfrastructure</div><div>Department of Biomedical Informatics</div><div>Stony Brook Medicine</div><div>Tel. : 1-631-444-3560</div><div>Fax  : 1-631-444-8873</div><div>Cell : 1-631-681-6228</div><div><a href="mailto:erich.bremer@stonybrook.edu" target="_blank">erich.bremer@stonybrook.edu</a></div><div>Office Location/Mailing Address</div><div>HSC, L3: Room 119</div><div>Stony Brook, NY 11794-8330</div><div><br></div></div></div></div></div>
</div></div>