You could regionprops and ask for MajorAxesLength to get "equivalent" diameter (an ellipse would be fitted to the data.)
I suspect that might be accurate enough for your purpose. If not then bwtraceboundary(), or regionprops pixellist . Then pdist(). And do notaverage the distances: take themaximum distance to get the diameter.
Best Answer