 Research
 Open Access
 Published:
Computeraided analysis of micromorphological structure of porous membranes
BioMedical Engineering OnLinevolume 17, Article number: 68 (2018)
Abstract
Background
The paper presents an approach and computeraided method of numerical evaluation of the quality of porous membranes used as scaffolds for cultivation of chondrocytes in regeneration of biological tissues.
Materials
The scanning electron microscope (SEM) images of 300× and 1000× magnification presenting the sections of artificial polyvinylpyrrolidone membranes obtained in two alternative production processes are examined.
Theory and methods
There is presented a combined morphological and statistical method of the assessment of artificial membranes’ porosity, based on computeraided segmentation and analysis of the size and shape of pores. Theoretical backgrounds of description pores as irregular objects in discrete 3dimensional space are presented. The parameters characterizing the quality of pores: pores irregularity coefficient and pores density are defined. The quality of the examined specimens of materials is characterized by the size (mean 2dimensional section areas) of pores. The main concept presented in the paper is the extraction of lacking information concerning the third dimension of pores from the 2dimensional SEM images of their sections. Two approaches to evaluation of the parameters characterizing pores on the basis of computeraided analysis of their crosssections are proposed: (1) based on statistical extension of geometrical data and (2) based on analysis of brightness profiles. The corresponding methods are based on the assumption of isotropy of the examined porous materials. The results of automatic measurements of the areas of pores, lengths of their chords and recording the brightness profiles along fixed lines crossing the analyzed images are illustrated by examples.
Conclusions
Practical usefulness of the proposed methods to evaluation of the quality of porous membranes consists in their ability to be used in case if alternative methods for some reasons cannot be used.
Background
Numerous physiological processes taking place in living organisms: absorption of nutritive products and gases, excretion of waste products etc. are based on the phenomenon of diffusion of various types of material particles through porous membranes. Similar processes play a substantial role in numerous areas of medical therapy based, in particular, on artificial organs: artificial lungsheart, artificial kidney, encapsulated drugs or living cells, cultivation of artificial tissues etc. [1,2,3,4,5]. Porous membranes application to water or air purification indirectly influences also human health [6,7,8]. In all the abovementioned cases porous membranes play the role of filters separating some biochemical media or of scaffolds for seeding and cultivation of tissues. The mechanisms of separation or of settling particles (cells) depend on biophysical, mechanical and microstructural properties of the porous materials. Examples of biophysical properties are electric and heatconductivity, dielectric constant, pH level, etc. Mechanical properties denote specific gravity, Young module, break or fractureresistance, etc. The notion of microstructural properties concerns spatial distribution, density, size and shape of the components the investigated material consists of, in particular—of pores and/or of the walls that separate them.
It was shown in several former papers [9,10,11,12,13] that computer technology combined with adequately chosen microscopic imaging methods may be an effective tool for material’s porosity analysis and evaluation. The 2D crosssections of the samples of porous materials provide interesting but highly intricate images of magnified pores as cavities and channels penetrating the material, as shown in Fig. 1.
In most cases the pores are of different size and highly irregular shape. However, the notion of irregularity needs some explanation. We call regular curves (contours) or 3D shapes which can be precisely described by simple analytical functions or geometrical objects. Weakly irregular are contours or 3D shapes which can be described as parameterized compositions of finite number of regular curves (contours) or 3D shapes. Curves (contours) or 3D shapes are strongly irregular if they are neither regular nor weakly irregular. The class of strongly irregular curves (contours) and 3D shapes can be divided into the subclasses of random and chaotic objects. Random curves (contours) or 3D shapes are assumed to be able to be described by probabilistic models [14, 15]. Chaotic objects can be described by fractals [16] or any other, more sophisticated, theoretical tools. Below, it is assumed that pores with sufficient accuracy can be described as randomshape objects. Therefore, no geometrical model is suitable for their accurate description. The statistical geometry or morphology based methods [17, 18] seem to be more suitable for numerical analysis of this type of structures. However, considering images of pores as instances of specific 2D random fields seems to be more rational. On the other hand, porosity is a 3D rather than a 2D property of a rigid body. Therefore, it arises a nontrivial problem of obtaining information about the 3D structure of a porous material by analysis of images presenting its 2D crosssections. The problem is relatively simple if the samples of the investigated material can be cut into regularly distanced parallel slices, as proposed in [12, 19] Unfortunately, this cannot be done, e.g. in the case of high fragility of the samples. In such case, we are faced with a more sophisticated statistical decision making problem to be solved by using advanced computerbased methods. In this paper some concepts concerning overcoming the difficulty of 3D morphological structures description not based on analysis of stacks, but of single 2D images, are presented. The missing information about the 3rd dimension of the examined porous materials’ samples can (at least partially) be extracted from the brightness level distribution in the crosssection images or should be deduced under the assumption of 3dimensional morphological isotropy of the samples. Below, both methods will be shortly described. However, statistical porosity characteristics describe only some aspects of the investigated materials’ utility in their biomedical applications. The relationships between measurable morphological parameters and largescale properties of porous materials are still an open problem requiring additional observations and experiments and it is not considered here.
The analysis of irregular structures (in particular, of pores) is not a new problem. Its origins date the Buffon’s needle problem which initiated the development of stochastic geometry [20]. Roughly speaking, this area of investigations concerns the properties of random compositions of regular geometrical objects. This is not quite the problem of highly irregular shapes investigation arising in biomedical engineering, material engineering, geophysics, geomorphology, etc. [11, 14, 17, 21]. However, each of the disciplines mentioned here uses its specific experimental data acquisition methods and this causes that the corresponding data processing methods are not directly applicable in other application areas. For example, echosounding methods of geophysical data acquisition and SEM imaging of pores provide formally different types of experimental data, despite the fact that in both cases they concern less or more irregular morphological structures.
The paper is organized as follows. The “Materials” section shortly presents the materials used as a source of image data subjected to computer analysis. The “Theory and methods” section consists of four subsections. The “Preliminary image processing” subsection presents initial operations aimed at enhancement the quality of SEM images before their analysis. The “Basic assumptions” subsection describes some discrete geometry objects used in computer analysis of pores. The “Characteristics of irregular object” subsection contains definitions of basic parameters used to describe irregular shapes of pores. The “Rough evaluation of 3D objects’ porosity” subsection is devoted to presentation of proposed two approaches to evaluation of the volumes of pores on the basis of analysis of their 2D sections: 1st based on statistical extension of geometrical data and 2nd based on analysis of brightness profiles. “Conclusions” contain remarks summarizing the work.
Materials
The materials used in this study have been provided by the Laboratory of Semipermeable Membranes and Bioreactors of the Nalecz Institute of Biocybernetics and Biomedical Engineering, Polish Academy of Sciences. They had the form of series of images presenting the crosssections of samples of porous membranes produced in the Laboratory. The membranes were produced by an inversephase method using the Polillactide (PLLA), and polyvinylpyrrolidone (PVP) as basic input components [22]. Several variants of technological processes were used in order to select the most effective ones considering the highest membrane’s quality and the highest production process stability. The samples of investigated materials were cut and covered by a 7‒10 nm coat of gold and in this form they have been used as specimens for analysis in a scanning electron microscope (SEM) type Hitachi TM1000 using a 15 kV acceleration voltage. The crosssection images in 300× and/or 1000× magnifications in electronic form of (approximately) 1100 × 1280 pix TIFF standard visual objects were then provided for computer analysis.
Theory and methods
Preliminary image processing
All calculated morphological characteristics of pores are based on the distances expressed in pixels (pix), areas expressed in pix^{2} etc., recalculated to metric units. In order to make the parameters evaluated in various samples of material mutually comparable it is necessary to present the images in standardized size, magnification and scale (in pix/μm). Moreover, in order to outline the pores more exactly the contrast of the original SEM images should be enhanced. Typical image enhancement methods are described in [23, 24].
Figure 2 shows a typical crosssection of a pore: original SEM image (a), its brightness histogram (b) and the image enhanced by its histogram extension to the left side of brightness scale (c). This operation makes the borders and inner details of pores better visible, however, it changes the original relationship between the brightness level and the depth of the visible pores insight.
Basic assumptions
Morphological analysis of pores is aimed at investigation of their size, shapes and spatial distribution. As such, it is based on some topological and metrical notions. Original SEM images can be considered as objects imbedded in Euclidean (2 or 3dimensional) space, while their representatives subjected to computer analysis are objects of discrete (2D or 3D) spaces. The difference between the two types of formal models is particularly substantial if smallsize objects are taken into consideration. This is caused by the fact that distances between any pair of points in Euclidean space can be expressed by any nonnegative real number while in discrete space defined on a rectangular grid of points they are limited to enumerable sets of values, depending on an assumed type of connectivity of adjacent points. Pairs of different points can be called connective when their absolute (Manhattan) distances d: (a) satisfy the condition d = 1 or (b) satisfy the condition d ≤ n where n is the dimensionality of the space. As a consequence, in 2D space two types of connectivity: 4connectivity and 8connectivity can be established, as it follows from the below shown tables of points (⋇) adjacent to a central point (0) (Fig. 3).
Similarly, in a 3D space a 6connectivity and 26connectivity can be established, as illustrated by the tables (Fig. 4).
An open discrete linesegment can be defined as a minimal linearly ordered set of points Λ = [p_{ 0 }, p_{ 1 }…, p_{ κ },…, p_{ K }] such that p_{ 0 } and p_{ K } are connective to exactly one, while all the others to exactly two other points of the set. A closed discrete line segment (a contour) C is defined similarly excepting that it does not contain any point connective to less than two other points of the segment (p_{ K } is connective both to p_{ K−1 } and p_{ 0 }).
A finite discrete 2D or 3D object can be defined as a finite set Ω = {p_{ κ }}, κ = 1,2,…,K, such that each point p_{ k } is connective to at least one other point of Ω. The surface of Ω is defined as a maximal subset S_{ Ω }⊆ Ω such that each its point p_{ κ } is connective to at least one point not belonging to Ω. An object Ω is called thin if S_{ Ω }≡ Ω. Evidently, contours and surfaces of objects are thin objects themselves. The abovedefined discrete objects have also a property of compactness, i.e. any pair of their points can be linked by a line totally belonging to the object.
A discrete length of a line segment Λ or of a contour C, measured in pixels [pix], is defined as the number of points constituting them. Similarly, a discrete area of a surface S_{ Ω }, measured in [pix^{2}] and a discrete volume of an object, measured in [pix^{3}], are defined as numbers of points that constitute them.
The discrete 3D objects will be described in a lefthanded system of coordinates (x, y, z), adequate to be presented on a computer screen, as shown in Fig. 5.
The higherconnective models suit, in general, better to approximation of continuous shapes then the lowconnective models. However, evaluation of distance, area or volume measures in discrete models of geometrical objects by counting the numbers of pixels belonging to the objects may lead to substantial differences concerning the respective Euclidean measures of the objects. This should be taken into account in morphological analysis of pores.
Characteristics of irregular objects
The abovedefined geometrical notions will be used to describe some basic properties of pores. In particular, the distance and/or area measures can be used as basic data to evaluate morphological characteristics of irregular shapes of natural objects. A 2D contour of a pore’s crosssection can be roughly characterized by its length L, area S and regularity coefficient:
taking value η=1 in circles and decreasing to 0 if the shape of contour becomes more irregular.
Similarly, a regularity coefficient of a 3D object can be defined as:
where V denotes the volume of the object,
and γ = 1 in spheres and it is decreasing to 0 if the shape of a 3D object (pore) becomes more irregular.
The abovegiven formulae suit well if the input data L, S and/or V can be exactly measured. This takes place, if a 2D object has been outlined, as shown in Fig. 6 where a crosssection of a microcapsule has been outlined by an ellipse.
If the larger D and smaller d diameters of the ellipse are parallel to the screen axes (as it, approximately, takes place in Fig. 6) then the Euclidean lengths of the diameters are given by the formulae:
where D and d are (both in case of 4 and 8connectivity) expressed by the number of pixels and δ denotes a scale factor given in metric units (e.g. in μm) per pixel. Then the wellknown formulae for ellipses can be used:
and, from (1),
However, the situation differs in a more general case when the distances and areas in a discrete space 2D (both, in the 4 and 8connectivity case) are expressed by the numbers of pixels constituting the objects (lines or contours’ inside). In such case the numerical values of contour lengths or of areas not only differ from those in Euclidean space but also depend on the type of connectivity and on the size and spatial orientation of the objects. This is show in Table 1 where four circles of different radii have been approximated by their discrete representatives and their regularity coefficients η have been calculated.
It follows from the Table that assumption of 8connectivity of the discrete space suits to describe irregular contours better than 4connectivity. Moreover, in both cases the regularity coefficient η is better suited for characterization of larger than of small objects. In any case, the type of connectivity should be a priori strongly established in order to make the results of porosity evaluation mutually comparable. However, it still remains the problem, how to evaluate the values of S and V of 3D objects observable only in 2D space without a possibility to measure them directly in all three dimensions.
Rough evaluation of 3D objects’ porosity
The differences between the discrete and continuous measures of distance, area, or volume may be less substantial if morphological parameters based on them are used only to comparative analysis of porous materials. In such case the easier measurable discrete (8connective) values L_{ i } and S_{ i } of outlined pores (i enumerating the pores in the image) within the frame of an examined SEM image can be evaluated. For them the regularity coefficients η can be calculated. The following parameters are widely used as basic porosity characteristics [13, 25]:

(a)
Mean pores’ volume m_{ V },

(b)
Standard deviation of pores’ volume ∆_{ V },

(c)
Total pores’ density:
$$\rho = \frac{{V_{\sum } }}{V}$$(7)
where V_{ Σ } denotes a total volume of pores in the volume V of the analyzed specimen. The density ρ may vary between 0 and 1.
In order to calculate the abovelisted parameters the volumes V_{ i } of pores in the observation area should be evaluated. In [26] a method of statistical evaluation of V_{ i } based on the assumption of isotropy of the morphological structure of pores and of their quasispherical shapes was described. Below, two other approaches to solution of the problem are presented.
Approach based on statistical extension of geometrical data
This approach is particularly suitable to analysis of pores whose shapes can roughly be predicted by taking into consideration the method of porous material production. In particular, porosity can be caused by leaving free areas in a squeezed bunch of fibers. In such case the structure of pores is anisotropic: irregular (with holes) in a plane perpendicular to the dominating direction of fibers and lengthened along the axes of fibers. Example of such porous structure is shown in Fig. 7.
In Fig. 7a brightness values varying mostly near two extreme levels: the highest level corresponding to the crosssections of fibers and the lowest, dark level corresponding to the pores are remarkable. This suggests that the pores density ρ in this case can be given approximately by the formula
where S denotes the area of the image frame and S_{ Σ } stands for total area of all white spots in the frame corresponding to the crosssections of the fibers.
Another situation arises in the case of isotropic type of porosity. Evaluation of porosity can then be based on the following assumptions.
Assumption 1
The morphological structure of the examined porous material is isotropic in 3D; therefore, the probability distributions of the lengths of parallel chords crossing the pores do not depend on the direction of the chords.
Assumption 2
The crosssections of the pores in (y, z) plane are of elliptic form with the major and minor axes ordered in parallel to the y and z coordinate axes.
Assumption 3
The pore’s crosssections in the (x, y) plane are of compact irregular form.
Assumption 3 seems to suit better to characterization of irregular shapes of pores than this based on the concept of effective pore diameter proposed in [27,28,29].
Figure 8 presents a typical contour of a pore, as visible in (x, y) plane, divided into horizontal strips of equal breadths δ.
The lengths of soobtained horizontal chords are denoted by l_{ i }. According to Assumption 2, the partition of the contour into strips corresponds to cutting the 3D pores into horizontal disks of a constant thickness. The disks then can be approximated by flat cylinders of elliptic base and constant height. The chord length l_{ i } can be taken as an axis of the elliptic base of the cylinder, the other (perpendicular) axis should be directed along the (invisible) coordinate z. In order to calculate the volume of the cylinder:
the length g_{ i } of the second axis should be evaluated. For this purpose Assumption 1 justifying the proposed statistical extension approach will be used. For any given, ith horizontal chord a set of crossing it vertical equally distanced parallel chords of the pore’s contour should be found and their lengths h_{ 1 }, h_{ 2 },…, h_{ J } should be measured, as shown in Fig. 9. For the sake of simplicity, the distances between the vertical chords can also be fixed equal ε, ε > 0, in particular, ε = δ.
Thus, it is assumed that for any given (ith) horizontal chord a conditional probability distribution (cpd) of the length of the assigned to it vertical chords w(hl_{ i }) can be taken into consideration and the same cpd is in effect also for the length z of horizontal chords crossing the ith chord and directed perpendicularly to the visualization plane.
Let us consider the set Φ_{ i } = {h_{ 1 }, h_{ 2 },…,h_{ J }} of measured lengths of the vertical chords. Two approaches to evaluate the missing value g_{ i } on the basis of Φ_{ I } can be taken into account:

(a)
Based on the mean value:
$$g_{i} = M_{i} = \frac{1}{J}\sum\limits_{j = 1}^{J} {h_{j} } ;$$(10) 
(b)
Based on the highest conditional probability.
This approach is reasonable only in case of large values of J (e.g. J ≥ 10). For this purpose let us find the minimal and maximal h_{ max } values h_{ min }, h_{ max } ∈ Φ_{ i }. Then, for a chosen value K, 0 < K < J the length of an auxiliary subinterval:
should be calculated. For k, 1 ≤ k ≤ K, the numbers n_{ κ } of h_{ j } ∈Φ satisfying the inequalities:
are then defined and should be calculated.
Let κ* be such that n_{ κ* } = max[n_{ 1 }, n_{ 2 },…, n_{ K }]. If n_{ κ* } is unique then we put:
Otherwise, if more than one subinterval satisfies the condition of maximization n_{ κ } the number K should be reduced by 1 and the procedure should be repeated. In a marginal case, when K = 1 is reached, all h_{ j } have approximately the same value h* and g_{ i } = h* as a solution can be accepted.
Finally, the values g_{ i } given by (10) or (13) can be used in formula (9) for calculation of the volumes v_{ i } and
can be used in formula (7) for calculation of the pores density ρ.
Approach based on analysis of brightness profiles
In this section porous materials whose pores have inner shape similar to this shown in Fig. 10a are considered.
The visualized pores are not only of various sizes and shapes but also of various brightness caused by their different depths and inner profiles. The profiles recorded along several horizontal lines crossing the image, shown in Fig. 10b, 1–4, suggest a possibility of extraction of information concerning the depths of pores in the third dimension. Therefore, the below presented method will be based on the following:
Assumption 4
The brightness levels in 2D visualization of pores are linearly proportional to the metric depths of pores with negative coefficient of proportionality.
Moreover, the former Assumptions 1 and 3 hold as well. Taking into account that the recorded brightness levels β are given by the integers satisfying the inequality 0 ≤ β ≤ 255, the following relation between β and the corresponding depth ζ of a pore can be established:
the function being plotted in Fig. 11. Here D_{ max } denotes the maximal expected depth of pores, 255 − ε is the lowest brightness level below of which the pixels are qualified as belonging to pores;
is a scaling coefficient.
The constants D_{ max } and ε should be established on the basis of image analysis. Figure 12 shows the histograms of brightness levels β recorded along the lines 1–4 in the image in Fig. 10a.
The brightness scale of 0–255 units has been divided into subintervals by 10 units enumerated from 0 to 26. Each bar in the histogram represents the number of pixels whose brightness levels belong to the given subinterval.
In all histograms domination of the 26th bars over the several preceding bars is remarkable. This means that in the given case ε = 10 in the formulae (15, 16) can be used.
At the next step, the profiles shown in Fig. 10b should be taken into consideration. On the level 255 − ε the lengths D_{ α } of subintervals crossing the areas where 0 ≤ β <255 − ε should be recorded.
The value:
then can be used in the formulae (15, 16). Finally, if ζ_{ ν } ν = 1,2,…, N, denote the calculated depths of pores and N is a total number of pixels in the analyzed image then the value V_{Σ} in formula (7) will be given as a sum:
because each pixel in the analyzed image represents a cuboid of 1 × 1 × ζ_{ ν } [pix] size. For similar reason, it is
which can be inserted into the formula (7).
Conclusions
Artificial porous membranes have a wide spectrum of applications in various areas of medicine. Porosity of membranes directly influences their utility as filters separating some types of particles from a fluid matter or as scaffolds for cultivation of biological tissues. However, the notion of porosity contains several properties directly connected with morphological microstructure of pores. The size, density and shape of pores can be mentioned as the basic porosity parameters. Evaluation of the parameters is thus an auxiliary but still important problem to be solved in order to improve the quality of membranes and to increase their utility. However, the problem is not trivial because of irregularity of the shapes of pores and lack of access to some characterizing them data. In this paper some concepts concerning the methods of artificial membranes’ porosity evaluation based on computeraided analysis of microscope images of their crosssections are presented. The proposed methods are the result of a longterm cooperation between two working groups in the Nalecz Institute of Biocybernetics and Biomedical Engineering PAS: first (headed by Ph.D. Przytulska), specialized in computeraided image processing methods and second (headed by prof., D.Sc. Chwojnowski) specialized in artificial membranes design and applications. The methods are based on combined morphological and statistical approach and on some assumptions concerning the statistical isotropy of the materials under investigation. Unlike other methods of 3D pores analysis (e.g. the Focused Ion Beam or Serial Block Face methods) assuming easy access to deep sections of the specimens of examined materials, the proposed methods can be recommended in case of pores analysis based on single SEM images of 2D sections of the materials. It is shown that all the parameters needed to practically use the proposed methods can be obtained by relatively simple numerical analysis of images: measuring distances, areas, brightness profiles, calculation of histograms and their parameters, etc. The necessary procedures are available in various statistical or image processing libraries of programs [23, 24]. However, in order to make our research more effective an original specialized software MeMoExplorer™ has been designed and it is still under development [27].
References
 1.
Chang TMS, Macintosh FC, Mason SG. Semi permeable aqueous microcapsules: I. Preparation and properties. Can J Pharmacol. 1966;44:115–28.
 2.
Tanaka Y, Yamaoka H, Nishizawa S, Nagata S, Ogasawara T, Asawa Y, et al. The optimization of porous polymeric scaffolds for chondrocyte/atelocollagen based tissueengineered cartilage. Biomaterials. 2010;31:4506–16.
 3.
Wu CH, Lee FK, Kumar SS, Ling QD, Chang Y, Chang Y, et al. The isolation and differentiation of human adiposederived stem cells using membrane filtration. Biomaterials. 2012;33:8228–39.
 4.
Kinasiewicz A, Dudziński K, Chwojnowski A, Weryński A, Kawiak J. Threedimensional culture of hepatocytes on spongy polyethersulfone membrane developed for cell transplantation. Transplant Proc. 2007;39(9):2914–6.
 5.
Hollister SJ. Porous scaffold design for tissue engineering. Nat Mater. 2005;4:518–24.
 6.
Lee A, Elam JW, Darling SB. Membrane materials for water purification: design, evelopment, and application. Environ Sci Water Res Technol. 2016;2(1):17–42.
 7.
Shannon MA, Bohn PW, Elimelech M, Georgiadis JG, Mariñas BJ, Mayes AM. Science and technology for water purification in the coming ecades. Nature. 2008;452(7185):301–10.
 8.
Jackson EA, Hillmyer MA. Nanoporous membranes derived from Block copolymers: from drug delivery to water filtration. ACS Nano. 2010;4(7):3548–53.
 9.
Zeman L, Denault L. Characterization of microfiltration membranes b y image analysis of electron micrographs. Part I. Method development. J Membr Sci. 1992;71(3):221–31.
 10.
Chwojnowski A, Przytulska M, Wierzbicka D, Kulikowski JL, Wojciechowski C. Membrane’s porosity evaluation by computeraided analysis of SEM images, a preliminary study. Biocybern Biomed Eng. 2012;32(4):65–75.
 11.
Przytulska M, Kulikowski JL, Lewińska D, Grzeczkowicz M, KupikowskaStobba B. Computeraided image analysis for microcapsules’ quality assessment. Biocybern Biomed Eng. 2015;35(4):342–50.
 12.
Sundaramoorthia G, Hadwigera M, BenRomdhanea M, Behzadb AR, Madhavanc P, Nunesc SP, Pereira S. 3D membrane imaging and porosity visualization. Ind Eng Chem Res. 2016;55:3689–95.
 13.
Al Marzooqi FA, Bilad M, Mansoor B, Arafat HA. A comparative study of image analysis and porometry techniques for characterization of porous membranes. J Mater Sci. 2016;51(4):2017–32.
 14.
Moran PAP. The probabilistic basis of stereology. In: Proceedings of symposium on statistical and probabilistic problem in metallurgy. Seattle, Washington, 4–6 August, 1971, ed. by Nicholson—Special Supplement to Advanced Applied Probability. 1972, p. 69–91.
 15.
Harding EF, Kendall DG, editors. Stochastic geometry: a tribute to the memory of Rollo Davidson. New York: Wiley; 1974.
 16.
Mandelbrot BB. The fractal geometry of nature. New York: W.H. Freeman and Co.; 1983.
 17.
Vulfson NI. Methods of stereology in geophysics. Leningrad: Gidrometeoizdat; 1989 (in Russian).
 18.
Giardina CR, Dougherty ER. Morphological methods in image and signal processing. Englewood Cliffs: Prentice Hall; 1987.
 19.
Nunes SP, Peinemann KV. Membrane technology: in the chemical industry. New York: Wiley; 2006.
 20.
Sylvester JJ. On a funicular solution of Buffon’s needle problem. Acta Math. 1890;14:185–205.
 21.
Armand NA, Krapivin VF, Mkrtshyan FA. Methods of radiophysical data processing in natural environment investigations. Moscow: Nauka; 1987 (in Russian).
 22.
Dudziński K, Chwojnowski A, Gutowska M, Płończak M, Czubak J, Łukowska E, Wojciechowski C. Three dimensional polyethersulphone scaffold for chondrocytes cultivation—the future supportive material for articular cartilage regeneration. Biocybern Biomed Eng. 2010;30(3):65–76.
 23.
Russ JC. The image processing handbook. 6th ed. Boca Raton: CRC Press; 2011.
 24.
Pitas I. Digital image processing algorithms and applications. New York: A WileyInterscience Publication; 2013.
 25.
Tung KL, Chang KS, Wu TT, Lin NJ, Lee KR, Lai JY. Recent advances in the characterization of membrane morphology. Curr Opin Chem Eng. 2014;4:121–7.
 26.
Przytulska M, Kulikowski J, Wasyłeczko M, Chwojnowski A, Piętka D, Evaluation of 3D morphological structure of porous membranes based on computeraided analysis of their 2D images. In: XII scientific conference “membranes and membrane processes in environmental protection” MEMPEP 2018 13th–16th of June 2018, Zakopane, Poland (submitted to presentation).
 27.
Shea FH, Tungb KL, Konga LX. Calculation of effective pore diameters in porous filtration membranes with image analysis. Robot Comput Integr Manuf. 2008;24:427–34.
 28.
MeyMarom A, Katz MG. Measurement of active pore size distribution of microporous embrane—a new approach. J Membr Sci. 1986;27(2):119–30.
 29.
Alhichri HS, Kamel M. Virtual circles: a new set of features for fast image registration. Pattern Recognit Lett. 2003;24(9–10):1181–90.
Authors’ contributions
JLK carried out the detailed theoretical derivation, drafted and revised the manuscript; MP collected the data and carried out the experimental analysis, checked the theoretical derivation and revised the manuscript; AC supervised the theoretical development and revised the manuscript. All authors read and approved the final manuscript.
Acknowledgements
We would like to express our thanks to Prof. R. Koprowski for encouraging us to prepare the manuscript, the Reviewers for their valuable critical remarks and recommendations as well as to MSc. P. Mazur and MSc. R. Maksymowicz for their help to improve English style of the paper.
Competing interests
The authors declare that they have no competing interests.
Ethics approval and consent to participate
Not applicable.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Author information
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Porous membranes
 Porosity evaluation
 Image processing
 Membrane’s quality assessment