ds format in GREIT algorithm

[JeanRintoul]

I've noticed that running JAC and Back Projection algorithm result in ds format which can be overlayed on the mesh as shown below:

Is there anyway to obtain ds in GREIT algorithm in the same format to make plotting it consistent across algorithms? Currently it outputs an n by n array.

[JeanRintoul]

Also, please let me know if this is your preferred mode of communication, or whether there is a forum or other for comments like this.

[liubenyuan]

For greit you can use pcolorfast from matplotlib,

# suppose xg, yg are the grids along x- and y- axis, eit_im is the greit image
xv = xg[0]
yv = yg[:, 0]
eit_im = eit_im.reshape(xg.shape)
ax.pcolorfast(xv, yv, eit_im, cmap=your_cmap)

Sure, it's ok to communicate via the issues, currently, we do not have a forum yet.

[JeanRintoul]

Great - that works and should do the trick, although the different format still means different plotting methods need to be used to create a similar plot across algorithms. Perhaps this is unavoidable.

[liubenyuan]

yes. In our group, GREIT has been used for pixel-based image analysis of EIT applications.