Modelling of oedemous limbs and venous ulcers using partial differential equations
Identifieur interne : 001C18 ( Pmc/Corpus ); précédent : 001C17; suivant : 001C19Modelling of oedemous limbs and venous ulcers using partial differential equations
Auteurs : Hassan Ugail ; Michael J. WilsonSource :
- Theoretical Biology & Medical Modelling [ 1742-4682 ] ; 2005.
Abstract
Oedema, commonly known as tissue swelling, occurs mainly on the leg and the arm. The condition may be associated with a range of causes such as venous diseases, trauma, infection, joint disease and orthopaedic surgery. Oedema is caused by both lymphatic and chronic venous insufficiency, which leads to pooling of blood and fluid in the extremities. This results in swelling, mild redness and scaling of the skin, all of which can culminate in ulceration.
We present a method to model a wide variety of geometries of limbs affected by oedema and venous ulcers. The shape modelling is based on the PDE method where a set of boundary curves are extracted from 3D scan data and are utilised as boundary conditions to solve a PDE, which provides the geometry of an affected limb. For this work we utilise a mixture of fourth order and sixth order PDEs, the solutions of which enable us to obtain a good representative shape of the limb and associated ulcers in question.
A series of examples are discussed demonstrating the capability of the method to produce good representative shapes of limbs by utilising a series of curves extracted from the scan data. In particular we show how the method could be used to model the shape of an arm and a leg with an associated ulcer.
We show how PDE based shape modelling techniques can be utilised to generate a variety of limb shapes and associated ulcers by means of a series of curves extracted from scan data. We also discuss how the method could be used to manipulate a generic shape of a limb and an associated wound so that the model could be fine-tuned for a particular patient.
Url:
DOI: 10.1186/1742-4682-2-28
PubMed: 16078992
PubMed Central: 1198260
Links to Exploration step
PMC:1198260Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Modelling of oedemous limbs and venous ulcers using partial differential equations</title><author><name sortKey="Ugail, Hassan" sort="Ugail, Hassan" uniqKey="Ugail H" first="Hassan" last="Ugail">Hassan Ugail</name><affiliation><nlm:aff id="I1">School of Informatics, University of Bradford, Bradford BD7 1DP, UK</nlm:aff></affiliation></author><author><name sortKey="Wilson, Michael J" sort="Wilson, Michael J" uniqKey="Wilson M" first="Michael J" last="Wilson">Michael J. Wilson</name><affiliation><nlm:aff id="I2">Department of Applied Mathematics, University of Leeds, Leeds LS2 9JT, UK</nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">16078992</idno><idno type="pmc">1198260</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC1198260</idno><idno type="RBID">PMC:1198260</idno><idno type="doi">10.1186/1742-4682-2-28</idno><date when="2005">2005</date><idno type="wicri:Area/Pmc/Corpus">001C18</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">001C18</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Modelling of oedemous limbs and venous ulcers using partial differential equations</title><author><name sortKey="Ugail, Hassan" sort="Ugail, Hassan" uniqKey="Ugail H" first="Hassan" last="Ugail">Hassan Ugail</name><affiliation><nlm:aff id="I1">School of Informatics, University of Bradford, Bradford BD7 1DP, UK</nlm:aff></affiliation></author><author><name sortKey="Wilson, Michael J" sort="Wilson, Michael J" uniqKey="Wilson M" first="Michael J" last="Wilson">Michael J. Wilson</name><affiliation><nlm:aff id="I2">Department of Applied Mathematics, University of Leeds, Leeds LS2 9JT, UK</nlm:aff></affiliation></author></analytic><series><title level="j">Theoretical Biology & Medical Modelling</title><idno type="eISSN">1742-4682</idno><imprint><date when="2005">2005</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><sec><title>Background</title><p>Oedema, commonly known as tissue swelling, occurs mainly on the leg and the arm. The condition may be associated with a range of causes such as venous diseases, trauma, infection, joint disease and orthopaedic surgery. Oedema is caused by both lymphatic and chronic venous insufficiency, which leads to pooling of blood and fluid in the extremities. This results in swelling, mild redness and scaling of the skin, all of which can culminate in ulceration.</p></sec><sec sec-type="methods"><title>Methods</title><p>We present a method to model a wide variety of geometries of limbs affected by oedema and venous ulcers. The shape modelling is based on the PDE method where a set of boundary curves are extracted from 3D scan data and are utilised as boundary conditions to solve a PDE, which provides the geometry of an affected limb. For this work we utilise a mixture of fourth order and sixth order PDEs, the solutions of which enable us to obtain a good representative shape of the limb and associated ulcers in question.</p></sec><sec><title>Results</title><p>A series of examples are discussed demonstrating the capability of the method to produce good representative shapes of limbs by utilising a series of curves extracted from the scan data. In particular we show how the method could be used to model the shape of an arm and a leg with an associated ulcer.</p></sec><sec><title>Conclusion</title><p>We show how PDE based shape modelling techniques can be utilised to generate a variety of limb shapes and associated ulcers by means of a series of curves extracted from scan data. We also discuss how the method could be used to manipulate a generic shape of a limb and an associated wound so that the model could be fine-tuned for a particular patient.</p></sec></div></front><back><div1 type="bibliography"><listBibl><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Theor Biol Med Model</journal-id><journal-title>Theoretical Biology & Medical Modelling</journal-title><issn pub-type="epub">1742-4682</issn><publisher><publisher-name>BioMed Central</publisher-name><publisher-loc>London</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">16078992</article-id><article-id pub-id-type="pmc">1198260</article-id><article-id pub-id-type="publisher-id">1742-4682-2-28</article-id><article-id pub-id-type="doi">10.1186/1742-4682-2-28</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research</subject></subj-group></article-categories><title-group><article-title>Modelling of oedemous limbs and venous ulcers using partial differential equations</article-title></title-group><contrib-group><contrib id="A1" corresp="yes" contrib-type="author"><name><surname>Ugail</surname><given-names>Hassan</given-names></name><xref ref-type="aff" rid="I1">1</xref><email>h.ugail@bradford.ac.uk</email></contrib><contrib id="A2" contrib-type="author"><name><surname>Wilson</surname><given-names>Michael J</given-names></name><xref ref-type="aff" rid="I2">2</xref><email>mike@maths.leeds.ac.uk</email></contrib></contrib-group><aff id="I1"><label>1</label>School of Informatics, University of Bradford, Bradford BD7 1DP, UK</aff><aff id="I2"><label>2</label>Department of Applied Mathematics, University of Leeds, Leeds LS2 9JT, UK</aff><pub-date pub-type="collection"><year>2005</year></pub-date><pub-date pub-type="epub"><day>3</day><month>8</month><year>2005</year></pub-date><volume>2</volume><fpage>28</fpage><lpage>28</lpage><ext-link ext-link-type="uri" xlink:href="http://www.tbiomed.com/content/2/1/28"></ext-link><history><date date-type="received"><day>11</day><month>5</month><year>2005</year></date><date date-type="accepted"><day>3</day><month>8</month><year>2005</year></date></history><permissions><copyright-statement>Copyright © 2005 Ugail and Wilson; licensee BioMed Central Ltd.</copyright-statement><copyright-year>2005</copyright-year><copyright-holder>Ugail and Wilson; licensee BioMed Central Ltd.</copyright-holder><license license-type="open-access" xlink:href="http://creativecommons.org/licenses/by/2.0"><p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"></ext-link>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p><pmc-comment>
Ugail
Hassan
h.ugail@bradford.ac.uk
Modelling of oedemous limbs and venous ulcers using partial differential equations
2005 Theoretical Biology and Medical Modelling 2(1): 28-. (2005) 1742-4682(2005)2:1<28> urn:ISSN:1742-4682
</license></permissions><abstract><sec><title>Background</title><p>Oedema, commonly known as tissue swelling, occurs mainly on the leg and the arm. The condition may be associated with a range of causes such as venous diseases, trauma, infection, joint disease and orthopaedic surgery. Oedema is caused by both lymphatic and chronic venous insufficiency, which leads to pooling of blood and fluid in the extremities. This results in swelling, mild redness and scaling of the skin, all of which can culminate in ulceration.</p></sec><sec sec-type="methods"><title>Methods</title><p>We present a method to model a wide variety of geometries of limbs affected by oedema and venous ulcers. The shape modelling is based on the PDE method where a set of boundary curves are extracted from 3D scan data and are utilised as boundary conditions to solve a PDE, which provides the geometry of an affected limb. For this work we utilise a mixture of fourth order and sixth order PDEs, the solutions of which enable us to obtain a good representative shape of the limb and associated ulcers in question.</p></sec><sec><title>Results</title><p>A series of examples are discussed demonstrating the capability of the method to produce good representative shapes of limbs by utilising a series of curves extracted from the scan data. In particular we show how the method could be used to model the shape of an arm and a leg with an associated ulcer.</p></sec><sec><title>Conclusion</title><p>We show how PDE based shape modelling techniques can be utilised to generate a variety of limb shapes and associated ulcers by means of a series of curves extracted from scan data. We also discuss how the method could be used to manipulate a generic shape of a limb and an associated wound so that the model could be fine-tuned for a particular patient.</p></sec></abstract></article-meta></front><body><sec><title>1 Introduction</title><p>Oedema, commonly known as tissue swelling, is associated with a range of causes such as venous disease, trauma, infection, joint disease, orthopaedic surgery and removal of the lymph nodes. Oedema and associated venous ulcers occur on mainly on the leg and the arm. It can be a painful, embarrassing and costly disorder [<xref ref-type="bibr" rid="B1">1</xref>,<xref ref-type="bibr" rid="B2">2</xref>]. It occurs widely in the general population, especially from late middle age, in diabetics and in immobile patients [<xref ref-type="bibr" rid="B3">3</xref>-<xref ref-type="bibr" rid="B5">5</xref>]. Apart from the tissue swellings the ulcers themselves can typically range in size from around 0.5 cm to 10 cm across, and are of variable depth [<xref ref-type="bibr" rid="B6">6</xref>,<xref ref-type="bibr" rid="B7">7</xref>]. Fig. <xref ref-type="fig" rid="F1">1</xref> shows an example of an oedemous leg infected with venous ulcers.</p><fig position="float" id="F1"><label>Figure 1</label><caption><p>An example of an oedemous leg infected with oedema venous ulcers.</p></caption><graphic xlink:href="1742-4682-2-28-1"></graphic></fig><p>An important task, during the treatment of oedema and venous ulcers, is the measurement of the amount of oedema as well as the area and volume of the ulcer wounds. This is because without an accurate and objective means of measuring changes in the size or shape of ulcers, it is difficult or impossible to evaluate the efficiency of the available therapies properly. Therefore, a prerequisite for this development is a reliable method of measuring ulcers. There exist a variety of measurement methods none of which is ideal. At present direct contacting measurements are widely used but they are not accurate, carry a risk of infection and are, to say the least, uncomfortable for the patient. For example, conventional techniques for measuring the area and volume of wounds depend on making physical contact with the wound, for example by drawing around the periphery on an acetate sheet or by making an alginate cast of the wound [<xref ref-type="bibr" rid="B6">6</xref>]. There is currently significant interest in developing non-invasive measurement systems using optical methods such as 'structured light' (a technique that projects stripes on to a surface and infers the shape from changes in the linearity of the reflected stripe) [<xref ref-type="bibr" rid="B8">8</xref>] or stereo-photogrammetry. The availability of high-resolution 3D digital cameras, increasing computing power and the development of software techniques for manipulating three-dimensional information have benefited this area. However, equipment associated with these sorts of measurement methods is not often portable and is often costly, thus making it prohibitive for routine medical use.</p><p>The aim of this paper is to show how it is possible to develop a system for measuring the shape and size of limbs and venous ulcers by means of utilising an economical mathematical model. In particular, one of the outcomes we hope to achieve from this work is a technique with potential for clinical use. For this, a small number of key measurements of limbs (with minimal possible contact with the limb and the associated ulcer), made using readily available instruments such as callipers and tape measures, can be input to a computer program. The program will then be able to reconstruct a good estimate of the limb shape and dimensions. It is believed that such a technique will provide a cheap, efficient, non-invasive instrument for measuring the degree of oedema and consequently enabling various treatment plans to be evaluated.</p><p>At present there exists a wide variety of methods that can be utilised to generate the geometry of limbs affected by oedema and venous ulcers. These include boundary based methods such as polygon based design [<xref ref-type="bibr" rid="B9">9</xref>], extrusions and surface of revolution [<xref ref-type="bibr" rid="B10">10</xref>] and polynomial patches [<xref ref-type="bibr" rid="B11">11</xref>]; procedural modelling such as implicit surfaces [<xref ref-type="bibr" rid="B12">12</xref>] and fractals [<xref ref-type="bibr" rid="B13">13</xref>]; and volumetric models such as constructive solid geometry [<xref ref-type="bibr" rid="B14">14</xref>] and subdivision [<xref ref-type="bibr" rid="B15">15</xref>]. Many of these techniques, especially polygon based design and polynomial patches, would be very appropriate for limb shape reconstruction, although they may not be ideally suited for the problem we address here. For example, conventional spline patches would require a large array of control points and weights in order to represent a realistic shape of a limb and associated wounds.</p><p>In this initial stage of the work we are concerned with developing efficient techniques in order to perform two important tasks. They are: the generation of smooth surfaces resembling the surface data obtained from a 3D scanner; and, once a smooth surface is obtained, manipulation of the geometry so as to obtain a good representation of the limb shape for any given patient. To do this we utilise real data from a series of surface scans provided by a medical partner namely, the Department of Medical Physics and Vascular Surgery of Bradford Teaching Hospitals National Health Services Trust (BTHNHST), UK, with whom we work closely on these problems. The department of Medical Physics at BTHNHST acquired the surface data using multiple-camera photogrammetry with a DSP400 system from 3dMD Ltd. This commercial technology has been widely used for acquiring medical images, especially in the USA, and captures data in a few milliseconds. The surface resolution (i.e. the separation of data points) is approximately 2 mm with a positional accuracy of approximately 0.2 mm. When developing our PDE based techniques for modelling human limbs, which are affected by oedema and venous ulcers, our medical partner has two aims. Firstly they require a compact and smooth surface representation of their captured data. Secondly, and rather more importantly, they require a modelling tool that would enable them to manipulate the shape of a given limb so as to provide a good representative limb shape of any given patient.</p><p>In this paper we utilise the so called PDE method [<xref ref-type="bibr" rid="B16">16</xref>-<xref ref-type="bibr" rid="B18">18</xref>] to address the problem. A positive feature of the PDE method is that it can define surfaces in terms of a small set of design variables [<xref ref-type="bibr" rid="B19">19</xref>], instead of many hundreds of control points. In broad terms this is because its boundary-value approach means that PDE surfaces are defined by data distributed around just their boundaries, instead of data distributed over their surface area, e.g. control points. Thus, a PDE model, when changed by altering the values of its design parameters, remains continuous; there is no need for a designer to intervene in order to close up any holes that might appear at patch boundaries. In the present context, this means that PDE surfaces can be made to adapt to changes in the shape of the limb and the associated wounds.</p></sec><sec><title>2 PDE Surfaces</title><p>A PDE surface is a parametric surface patch <inline-graphic xlink:href="1742-4682-2-28-i1.gif"></inline-graphic>, defined as a function of two parameters <italic>u </italic>and <italic>v </italic>on a finite domain Ω (⊂) <italic>R</italic><sup>2 </sup>by regarding the function <inline-graphic xlink:href="1742-4682-2-28-i2.gif"></inline-graphic> as a mapping of a point in Ω to a point <inline-graphic xlink:href="1742-4682-2-28-i1.gif"></inline-graphic> in the physical space. The shape of the surface patch is usually determined by specifying a set of boundary data at the edge of (∂)Ω. Typically the boundary data are specified in the form of <inline-graphic xlink:href="1742-4682-2-28-i1.gif"></inline-graphic> and a number of its derivatives on (∂)Ω. Hence, by casting the surface generation as a boundary value problem, the surface <inline-graphic xlink:href="1742-4682-2-28-i1.gif"></inline-graphic> is regarded as a solution of an elliptic PDE.</p><p>Various elliptic PDEs could be used; the ones we utilise for this work are based on the biharmonic and triharmonic equations, namely,</p><p><inline-graphic xlink:href="1742-4682-2-28-i3.gif"></inline-graphic></p><p>and</p><p><inline-graphic xlink:href="1742-4682-2-28-i4.gif"></inline-graphic></p><p>Also, periodic boundary conditions are very often considered. Assuming we are working with the above two elliptic PDEs, we require them to satisfy a set of 2<italic>N </italic>conditions, where <italic>N </italic>is 2 in the case of Equation (1) and <italic>N </italic>3 in the case of Equation (2). The general form of these conditions can then be written as,</p><p><bold>X</bold>(0, <italic>v</italic>) = <bold>f</bold><sub>1</sub>(<italic>v</italic>), (3)</p><p><bold>X</bold>(<italic>u</italic><sub><italic>i</italic></sub>, <italic>v</italic>) = <bold>g</bold><sub><italic>i</italic></sub>(<italic>v</italic>), <italic>i </italic>= 2 ... 2<italic>N </italic>- 1 (4)</p><p><bold>X</bold>(1, <italic>v</italic>) = <bold>f</bold><sub>2<italic>N</italic></sub>(<italic>v</italic>), (5)</p><p>where <bold>f</bold><sub>1</sub>(<italic>v</italic>) in Equation (3) and <bold>f</bold><sub>2<italic>N</italic></sub>(<italic>v</italic>) in Equation (5) are function conditions specified at <italic>u </italic>= 0 and <italic>u </italic>= 1 respectively. The conditions <bold>X</bold>(<italic>u</italic><sub><italic>i</italic></sub>, <italic>v</italic>) = <bold>g</bold><sub><italic>i</italic></sub>(<italic>v</italic>) in Equation (4) can take the form either</p><p><bold>X</bold>(<italic>u</italic><sub><italic>i</italic></sub>, <italic>v</italic>) = <bold>f</bold><sub><italic>i </italic></sub>for 0 <<italic>u</italic><sub><italic>i </italic></sub>< 1, <italic>i </italic>= 2 ... 2<italic>N </italic>- 1, (6)</p><p>or</p><p><inline-graphic xlink:href="1742-4682-2-28-i5.gif"></inline-graphic></p><p>In simpler terms the above conditions imply that for a PDE surface patch of order 2<italic>N</italic>, we can specify two function conditions, as given in Equations (3) and (5), that should be satisfied at the edges (at <italic>u </italic>= 0 and <italic>u </italic>= 1) of the surface patch, and a number of function or derivative conditions, as given in Equation (4), amounting to 2<italic>N </italic>– 2 conditions that the PDE should also satisfy.</p><sec><title>2.1 Solution of the PDEs</title><p>There exist many methods for solving Equations (1) and (2) ranging from analytic solutions to sophisticated numerical methods. The problems we address in this paper involve modelling of human limbs, which are essentially closed and cylindrical, and therefore the broad range of shapes encountered can be incorporated by solving the chosen PDEs with periodic conditions. Note here periodic conditions imply that for the <italic>v </italic>parameter the condition, <inline-graphic xlink:href="1742-4682-2-28-i6.gif"></inline-graphic>, is satisfied. Thus, for the work described here, we restrict ourselves to periodic conditions and obtain a closed form analytic solution of Equations (1) and (2).</p><p>Choosing the parametric region to be 0 ≤ <italic>u </italic>≤ 1 and 0 ≤ <italic>v </italic>≤ 2π, and assuming that the conditions given in Equations (3), (4) and (5) are periodic functions, we can use the method of separation of variables and spectral approximation [<xref ref-type="bibr" rid="B20">20</xref>] to write down the analytic solution of Equations (1) and (2) as,</p><p><inline-graphic xlink:href="1742-4682-2-28-i7.gif"></inline-graphic></p><p>where <inline-graphic xlink:href="1742-4682-2-28-i8.gif"></inline-graphic> is a polynomial function and <inline-graphic xlink:href="1742-4682-2-28-i9.gif"></inline-graphic> and <inline-graphic xlink:href="1742-4682-2-28-i10.gif"></inline-graphic> are exponential functions. The specific forms of <inline-graphic xlink:href="1742-4682-2-28-i11.gif"></inline-graphic> and <inline-graphic xlink:href="1742-4682-2-28-i10.gif"></inline-graphic> for the case of Equations (1) can be found in [<xref ref-type="bibr" rid="B17">17</xref>] and for the case of Equations (2) can be found in [<xref ref-type="bibr" rid="B21">21</xref>].</p><p>The main point to bear in mind regarding the above solution method is that it enables one to represent a set of general periodic conditions in terms of a finite <italic>M </italic>Fourier series, where <italic>M </italic>is typically taken to be ≤ 10, whilst the term <inline-graphic xlink:href="1742-4682-2-28-i10.gif"></inline-graphic>, which acts as a correction term, enables the conditions to be satisfied exactly. Detailed discussions of this solution method can be found in [<xref ref-type="bibr" rid="B20">20</xref>].</p></sec><sec><title>2.2 Methods of Generating PDE Surfaces</title><p>In this section we discuss a series of examples, showing the various methods by which PDE surfaces can be generated where the PDEs are chosen to be Equations (1) and (2) and the conditions are taken in the format described in Equations (3), (4) and (5).</p><p>As a first example we show how a fourth order PDE surface is generated where all the conditions are taken to be function conditions. Fig. <xref ref-type="fig" rid="F2">2(b)</xref> shows the shape of a surface generated by the fourth order PDE where the conditions are specified in terms of the curves shown in Fig. <xref ref-type="fig" rid="F2">2(a)</xref>. In particular, the conditions are such that: <inline-graphic xlink:href="1742-4682-2-28-i12.gif"></inline-graphic> and <inline-graphic xlink:href="1742-4682-2-28-i13.gif"></inline-graphic>. Since we are taking four function conditions to solve the fourth order PDE, all the curves in this case lie on the resulting surface. Thus, in this particular case the resulting PDE surface is a smooth interpolation between the given set of functional conditions.</p><fig position="float" id="F2"><label>Figure 2</label><caption><p>The shape of a surface generated by the fourth order <bold>PDE </bold>where the conditions are all taken to be function conditions (<bold>a</bold>) The conditions defined in the form of curves in 3-space. (<bold>b</bold>) The resulting surface shape.</p></caption><graphic xlink:href="1742-4682-2-28-2"></graphic></fig><p>The next example shows how a fourth order PDE surface is generated when the conditions are taken to be a mixture of function conditions and derivative conditions. Fig. <xref ref-type="fig" rid="F4">4(b)</xref> shows the shape of a surface generated by the fourth order PDE where two function boundary conditions and two derivative boundary conditions are specified in terms of the curves shown in Fig. <xref ref-type="fig" rid="F4">4(a)</xref>. In particular, the boundary conditions are chosen such that: <inline-graphic xlink:href="1742-4682-2-28-i14.gif"></inline-graphic> and <inline-graphic xlink:href="1742-4682-2-28-i15.gif"></inline-graphic>, where <italic>s </italic>is a scalar. In this case the surface patch generated as a solution to the fourth order PDE contains the boundary curves <italic>c</italic><sub>1 </sub>and <italic>c</italic><sub>4 </sub>whilst it does not necessarily contain the curves <italic>c</italic><sub>2 </sub>and <italic>c</italic><sub>3</sub>. A typical scenario where a surface of this nature is required would be a blend design where the derivative boundary curves can be adjusted to produce a smooth blend surface that bridges between two primary surfaces.</p><fig position="float" id="F4"><label>Figure 4</label><caption><p>The shape of a surface generated by the fourth order <bold>PDE </bold>where the boundary conditions are taken to be both positions and derivatives (<bold>a</bold>) The boundary conditions defined in the form of curves in 3-space. (<bold>b</bold>) The resulting surface shape.</p></caption><graphic xlink:href="1742-4682-2-28-4"></graphic></fig><p>Fig. <xref ref-type="fig" rid="F3">3(b)</xref> shows the shape of a surface generated by the sixth order PDE where the conditions are all taken to be positions specified in terms of the curves shown in Fig. <xref ref-type="fig" rid="F3">3(a)</xref>. In particular, the boundary conditions are such that <inline-graphic xlink:href="1742-4682-2-28-i16.gif"></inline-graphic> and <inline-graphic xlink:href="1742-4682-2-28-i17.gif"></inline-graphic>. As in the first example of the fourth order case, since we are taking six function conditions to solve the sixth order PDE, all the curves in this case lie on the resulting surface. Thus, the resulting PDE surface is a smooth interpolation between the six prescribed curves.</p><fig position="float" id="F3"><label>Figure 3</label><caption><p>The shape of a surface generated by the sixth order <bold>PDE </bold>where the conditions are all taken to be positions (<bold>a</bold>) The conditions defined in the form of curves in 3-space. (<bold>b</bold>) The resulting surface shape.</p></caption><graphic xlink:href="1742-4682-2-28-3"></graphic></fig><p>As a final example we show how a sixth order PDE surface is generated where the boundary conditions are taken to be a mixture of function boundary conditions and derivative conditions (both first and second order). Fig. <xref ref-type="fig" rid="F5">5(b)</xref> shows the shape of a surface generated by the sixth order PDE where two function boundary conditions, two first order derivative boundary conditions and two second order derivative boundary conditions are specified in terms of the curves shown in Fig. <xref ref-type="fig" rid="F5">5(a)</xref>. In particular, the boundary conditions are chosen such that <inline-graphic xlink:href="1742-4682-2-28-i18.gif"></inline-graphic> and <inline-graphic xlink:href="1742-4682-2-28-i19.gif"></inline-graphic>. where <italic>s </italic>and <italic>t </italic>are scalars. As in the example of fourth order case shown in Fig. <xref ref-type="fig" rid="F4">4</xref> the surface generated in this case contains the curves <italic>c</italic><sub>1 </sub>and <italic>c</italic><sub>6 </sub>whilst it does not necessarily contain the rest of the curves. Again this type of surface shape can be utilised in blend design where higher order continuity is desired in producing a smooth blend surface that bridges between two primary surfaces. As one can see from the format of these derivative condition definitions, the derivative conditions are all defined using simple finite difference schemes. The curves defining the derivative conditions provide an intuitive shape manipulation tool in that the shape of the surface closely follows the shape of the boundary conditions.</p><fig position="float" id="F5"><label>Figure 5</label><caption><p>The shape of a surface generated by the sixth order <bold>PDE </bold>where the boundary conditions are taken to be both positions and derivatives (<bold>a</bold>) The boundary conditions defined in the form of curves in 3-space. (<bold>b</bold>) The resulting surface shape.</p></caption><graphic xlink:href="1742-4682-2-28-5"></graphic></fig><p>The above examples demonstrate how PDE surfaces of order four and six can be utilised to generate surface shapes, which are applicable to a wide variety of design scenarios. Thus, the basic idea here is to generate a series of curves (both function and derivative) that can be utilised to define the boundary conditions for the chosen PDE. As seen in the examples, the resulting surface shape can always be intuitively predicted from the shapes of the chosen curves.</p></sec></sec><sec><title>3 Modelling of Limbs and Ulcers</title><p>In this section we discuss the shape modelling of human limbs affected by oedema and venous ulcers. In what follows, we discuss two examples of shape modelling of human limbs namely modelling of an arm shape and modelling of a leg shape with an ulcer. We utilise a mixture of PDEs of order four and six in order to model the surface shapes in question. In order to generate a representative smooth PDE surface shape, we extract a series of curves along the profile of the geometric model.</p><p>Fig. <xref ref-type="fig" rid="F6">6</xref> shows a typical surface scan data set provided by the medical partner where in this particular case the data set corresponds to an arm shape. Note that scan data are only available for half the surface. In order to generate a representative PDE surface shape, we extract a series of curves along the profile of the geometric model. To do this, first we import the geometric model into an interactive graphical environment through which we can examine and interact with the model. The geometric definitions of the scan data are provided in .obj file format where the 3D polygonal data with connectivity information are readily available. This enables us to display the model as well as compute the normal curvature distribution across the surface. A series of regions on the scan data model are then manually identified based on changes in the surface curvature. These regions are then utilised to determine the number of surface patches required to produce a good representative model of the limb in question. In determining the number of PDE surface patches required the aim is to reduce the number of patches that need to be utilised to produce a good representative geometric model with given accuracy. Once the number of surface patches required is decided the appropriate number of curves for each surface patch is extracted from the scan geometry data. To do this we create a series of free-form cubic spline curves within the interactive environment. The spline curves are then projected on to the scan geometry at the positions where the PDE curves are to be extracted. Note that the surface data obtained in this case do not naturally give us curves that are periodic. Thus, in this case, for each curve extracted, a series of fictitious points is added to each curve in order to make the curve periodic. The surfaces are then generated using the analytic solution described previously, where the surface is generated for the region 0 ≤ <italic>v </italic>≤ π which forms the portion of the curves extracted from the scanned data.</p><fig position="float" id="F6"><label>Figure 6</label><caption><p>Scanned surface data of an arm.</p></caption><graphic xlink:href="1742-4682-2-28-6"></graphic></fig><sec><title>3.1 Example 1: Modelling of an Arm Shape</title><p>As a first example, we discuss the modelling of the shape of a human arm. Fig. <xref ref-type="fig" rid="F6">6</xref> shows the scan surface data corresponding to an arm shape provided by the medical partner. Fig. <xref ref-type="fig" rid="F7">7(a)</xref> shows a series of curves extracted from the original scan surface data. Fig. <xref ref-type="fig" rid="F7">7(b)</xref> shows the arm shape generated using PDE surfaces. In particular the shape is generated as a combination of two fourth order patches and a single sixth order surface patch. i.e. the curves <italic>c</italic><sub>5</sub>, <italic>c</italic><sub>6</sub>, <italic>c</italic><sub>7 </sub>and <italic>c</italic><sub>8 </sub>and <italic>c</italic><sub>8</sub>, <italic>c</italic><sub>9</sub>, <italic>c</italic><sub>10 </sub>and <italic>c</italic><sub>11 </sub>form boundary conditions for two fourth order surface patches with the common boundary at <italic>c</italic><sub>8 </sub>whereby all the conditions are taken to be position conditions. The curves <italic>c</italic><sub>1</sub>, <italic>c</italic><sub>2</sub>, <italic>c</italic><sub>3 </sub>and <italic>c</italic><sub>5 </sub>form a sixth order surface patch where <italic>c</italic><sub>1</sub>, <italic>c</italic><sub>2</sub>, <italic>c</italic><sub>4</sub>, <italic>c</italic><sub>5 </sub>are taken to be four position conditions and the differences between <italic>c</italic><sub>2</sub>, <italic>c</italic><sub>3 </sub>and <italic>c</italic><sub>5</sub>, <italic>c</italic><sub>4 </sub>are taken to be two first order derivative boundary conditions. The value of the parameter <italic>s </italic>is taken to be 0.34.</p><fig position="float" id="F7"><label>Figure 7</label><caption><p>The arm shape generated using <bold>PDE </bold>surfaces by means of utilising curves extracted from scanned data (<bold>a</bold>) The extracted curves. (<bold>b</bold>) The resulting surface shape generated using two fourth order patches and a sixth order patch.</p></caption><graphic xlink:href="1742-4682-2-28-7"></graphic></fig></sec><sec><title>3.2 Example 2: Modelling of a Leg Shape with a Venous Ulcer</title><p>As a second example we discuss the modelling of the shape of a leg infected with a venous ulcer. Fig. <xref ref-type="fig" rid="F8">8</xref> shows the scan surface data, corresponding to the infected leg with a venous ulcer. As in the previous example, in order to generate a representative smooth surface shape, we first extract a series of curves along the leg and the associated wound. Fig. 9(a) shows a series of curves extracted from the original scan data.</p><fig position="float" id="F8"><label>Figure 8</label><caption><p>Scanned surface data of a leg infected with an ulcer.</p></caption><graphic xlink:href="1742-4682-2-28-8"></graphic></fig><p>In order to create a smooth shape that closely resembles the geometry of the leg, we utilise two sixth order patches to generate the main portion of the leg. Thus, the curves <italic>c</italic><sub>1</sub>, <italic>c</italic><sub>2</sub>, <italic>c</italic><sub>3</sub>, <italic>c</italic><sub>4</sub>, <italic>c</italic><sub>5 </sub>and <italic>c</italic><sub>6 </sub>form the position condition for a sixth order surface patch where the surface patch passes through these curves. The other surface patch is generated using the curves <italic>c</italic><sub>6</sub>, <italic>c</italic><sub>7</sub>, <italic>c</italic><sub>8</sub>, <italic>c</italic><sub>9 </sub>and <italic>c</italic><sub>10 </sub>where the curve <italic>c</italic><sub>6 </sub>is common to both surface patches. Moreover, for the later surface patch the curves <italic>c</italic><sub>6</sub>, <italic>c</italic><sub>7 </sub>and <italic>c</italic><sub>9</sub>, <italic>c</italic><sub>10 </sub>form four position boundary conditions and the differences between the curves <italic>c</italic><sub>7</sub>, <italic>c</italic><sub>8 </sub>and <italic>c</italic><sub>9</sub>, <italic>c</italic><sub>8 </sub>form two first order derivative boundary condition thus ensuring a smooth geometry transition between the foot and the leg. The parameter <italic>s </italic>is taken to be 0.12.</p><p>To generate the wound shape on the leg, we define a curve on the PDE leg surface that closely resembles the edge of the wound. This curve, marked as <italic>c</italic><sub>11 </sub>as shown in Fig. <xref ref-type="fig" rid="F9">9(a)</xref>, is generated using the (<italic>u</italic>, <italic>v</italic>) parameter space of the corresponding the PDE surface. Next the surface portion corresponding to the interior of the curve <italic>c</italic><sub>11 </sub>is trimmed out. This trimming process is again carried out using the (<italic>u</italic>, <italic>v</italic>) parameter space as described in [<xref ref-type="bibr" rid="B17">17</xref>]. Fig. <xref ref-type="fig" rid="F10">10(a)</xref> shows the main leg surface with the trim.</p><fig position="float" id="F9"><label>Figure 9</label><caption><p>The main leg shape generated using <bold>PDE </bold>surfaces by means of utilising curves extracted from scan data (<bold>a</bold>) The extracted curves (including the wound). (<bold>b</bold>) The resulting surface shape corresponding to the main shape of the leg, generated using two sixth order patches and a sixth order patches.</p></caption><graphic xlink:href="1742-4682-2-28-9"></graphic></fig><fig position="float" id="F10"><label>Figure 10</label><caption><p>Leg and ulcer geometry. (<bold>a</bold>) A portion of the surface is trimmed out using a curve resembling the edge of the ulcer wound (<bold>b</bold>) The complete leg geometry with the ulcer wound.</p></caption><graphic xlink:href="1742-4682-2-28-10"></graphic></fig><p>Once the appropriate trimming is carried out, a separate fourth order patch resembling the shape of the wound is generated where the curve <italic>c</italic><sub>11 </sub>which lies on the main leg surface is utilised as one of the four position boundary conditions. Fig. <xref ref-type="fig" rid="F10">10(b)</xref> shows the complete leg shape along with the ulcer wound.</p><p>Both the examples discussed above show how PDE surfaces of low order (i.e. order 4 and 6 in this case) can be utilised to generate good representative shapes using little information from the scan data. One could argue that a single PDE surface of higher order can be equally well suited to generating a single surface patch through a given number of curves. However, from the min-max principle for elliptic PDEs it is well known that PDEs of higher order (i.e. orders above 6) are difficult to control. Choosing lower order PDEs to generate the surface therefore makes sense. It is also noteworthy that the parameters <italic>s </italic>and <italic>t </italic>and the difference between the corresponding position and derivative curves enable both the size and the direction of the derivative boundary conditions at the edge of a given surface patch to be controlled. The derivative boundary conditions are used to control the smoothness of the blend between two surface patches. Such a tool cannot be deployed to reduce the number of curves used and hence the number of surface patches utilised to model the complete limb.</p></sec></sec><sec><title>4 Conclusion</title><p>This paper describes how the PDE method can be utilised to model a wide variety of geometries of limbs affected by oedema and venous ulcers. The shape modelling is based on solving a PDE subject to a set of curves extracted from 3D scan data providing the shape of the affected limbs. For this work we utilise a mixture of fourth order and sixth order PDEs, depending on the accuracy and continuity requirements for obtaining a good representative shape of the limb and associated ulcers in question.</p><p>In this work we are concerned with developing efficient techniques in order to undertake two important tasks. They are: the generation of smooth surfaces closely resembling the surface data obtained from a 3D scanner; and once a smooth surface is obtained, manipulation of the geometry so as to provide a good representative limb geometry shape for any given patient. Thus, the prime aim of the technique we discuss here is to generate a good representative shape of the limb quickly from the scanned data and to be able to manipulate that shape efficiently. It is noteworthy that the process of PDE geometry generation from the scan data is currently carried out manually. We are currently working on developing a methodology for automating this process. We have shown examples that clearly demonstrate the ability of PDE shape modelling techniques to generate a variety of limb shapes and associated venous ulcers. The geometry models themselves are flexible in terms of their manipulation capabilities, i.e. the manipulation of geometry can be carried out via the manipulation of curves defining the surface.</p><p>Our future direction in this work is to define a shape parameterisation tool for limbs where a set of shape parameters can be associated with the curves. Such shape parameterisation can then be utilised to fine tune a given generic limb model to suit to a handful of data measured from a given patient's limb. This will enable one to develop efficient non-invasive techniques for measuring various properties (such as surface area and volume) of oedema and venous ulcers.</p></sec></body><back><ack><sec><title>Acknowledgements</title><p>The authors wish to thank Dr. R.G. Cameron and Dr. W. Gardner of the Department of Medical Physics and Vascular Surgery at Bradford Teaching Hospitals National Health Services Trust of UK for fruitful discussions and supplying 3D scan data of limbs.</p></sec></ack><ref-list><ref id="B1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bosanquet</surname><given-names>N</given-names></name></person-group><article-title>Costs of Venous Ulcers: From Maintenance Therapy to Investment Programmes</article-title><source>Phlebology</source><year>1992</year><volume>7</volume><fpage>44</fpage><lpage>46</lpage></citation></ref><ref id="B2"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Doherty</surname><given-names>D</given-names></name><name><surname>Ross</surname><given-names>F</given-names></name><name><surname>Yeo</surname><given-names>L</given-names></name><name><surname>Uttley</surname><given-names>J</given-names></name></person-group><article-title>Leg Ulcer Management in an Integrated Service</article-title><source>The South Thames Evidence Based Practice (STEP) Project Report (6)</source><year>2000</year><publisher-name>Kingston University</publisher-name></citation></ref><ref id="B3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Prasad</surname><given-names>A</given-names></name><name><surname>Ali-Khan</surname><given-names>A</given-names></name><name><surname>Mortimer</surname><given-names>PS</given-names></name></person-group><article-title>Leg Ulcers and Oedema: A Study Exploring the Prevalence, Aetiology and Possible Significance of Oedema in Leg Ulcers</article-title><source>Phlebology</source><year>1990</year><volume>5</volume><fpage>181</fpage><lpage>187</lpage></citation></ref><ref id="B4"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Fletcher</surname><given-names>A</given-names></name></person-group><person-group person-group-type="editor"><name><surname>Harding KG, Leaper DL, Turner TD</surname></name></person-group><article-title>Common Problems of Wound Management in the Elderly</article-title><source>First European Conference on the Advances in wound management</source><year>1992</year><publisher-name>London: Macmillan</publisher-name><fpage>25</fpage><lpage>29</lpage></citation></ref><ref id="B5"><citation citation-type="other"><article-title>Lymphoedema Framework, Towards a National Framework for Lymphoedema Management</article-title><ext-link ext-link-type="uri" xlink:href="http://www.lymphoedemaframework.org"></ext-link></citation></ref><ref id="B6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Plassmann</surname><given-names>P</given-names></name><name><surname>Melhuish</surname><given-names>JM</given-names></name><name><surname>Harding</surname><given-names>KG</given-names></name></person-group><article-title>Methods of Measuring Wound Size: A Comparative Study</article-title><source>WOUNDS: A Compendium of Clinical Research and Practice</source><year>1994</year><volume>6</volume><fpage>54</fpage><lpage>61</lpage></citation></ref><ref id="B7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Elder</surname><given-names>D</given-names></name><name><surname>K</surname><given-names>G</given-names></name></person-group><article-title>Venous Disease: How to Heal and Prevent Chronic Leg Ulcers</article-title><source>Geriatrics</source><year>1995</year><volume>50</volume><fpage>30</fpage><lpage>36</lpage><pub-id pub-id-type="pmid">7635325</pub-id></citation></ref><ref id="B8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Krouskop</surname><given-names>TA</given-names></name><name><surname>Baker</surname><given-names>R</given-names></name><name><surname>Wilson</surname><given-names>MS</given-names></name></person-group><article-title>A Noncontact Wound Measurement System</article-title><source>Journal of Rehabilitation Research and Development</source><year>2002</year><volume>39</volume><fpage>337</fpage><lpage>346</lpage><pub-id pub-id-type="pmid">12173754</pub-id></citation></ref><ref id="B9"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Hoschek</surname><given-names>J</given-names></name><name><surname>Lasser</surname><given-names>D</given-names></name></person-group><source>Computer Aided Geometric Design</source><year>1993</year><publisher-name>A K Peters</publisher-name></citation></ref><ref id="B10"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Mortenson</surname><given-names>ME</given-names></name></person-group><source>Geometric Modelling</source><year>1985</year><publisher-name>New York: Wiley-Interscience</publisher-name></citation></ref><ref id="B11"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Farin</surname><given-names>G</given-names></name></person-group><source>Curves and Surfaces for Computer Aided Geometric Design, A Practical Guide</source><year>2001</year><publisher-name>Morgan-Kaufmann</publisher-name></citation></ref><ref id="B12"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Bajaj</surname><given-names>C</given-names></name><name><surname>Blinn</surname><given-names>J</given-names></name><name><surname>Bloomenthal</surname><given-names>J</given-names></name><name><surname>Cani-Gascuel</surname><given-names>M</given-names></name><name><surname>Rock-wood</surname><given-names>A</given-names></name><name><surname>Wyvill</surname><given-names>B</given-names></name><name><surname>Wyvill</surname><given-names>G</given-names></name></person-group><source>Introduction to Implicit Surfaces</source><year>1997</year><publisher-name>Morgan-Kaufmann</publisher-name></citation></ref><ref id="B13"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Szeliski</surname><given-names>R</given-names></name><name><surname>Terzopoulos</surname><given-names>D</given-names></name></person-group><article-title>From Splines to Fractals</article-title><source>Proceedings of the 16th annual conference on Computer Graphics and Interactive Techniques</source><year>1989</year><publisher-name>ACM New York</publisher-name><fpage>51</fpage><lpage>60</lpage></citation></ref><ref id="B14"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Rappoport</surname><given-names>A</given-names></name><name><surname>Spitz</surname><given-names>S</given-names></name></person-group><article-title>Interactive Boolean Operations for Conceptual Design of 3-d Solids</article-title><source>Proceedings of the 24th Annual Conference on Computer Graphics and Interactive Techniques</source><year>1997</year><publisher-name>ACM New York</publisher-name><fpage>269</fpage><lpage>278</lpage></citation></ref><ref id="B15"><citation citation-type="book"><person-group person-group-type="author"><name><surname>DeRose</surname><given-names>T</given-names></name><name><surname>Kass</surname><given-names>M</given-names></name><name><surname>Truong</surname><given-names>T</given-names></name></person-group><article-title>Subdivision Surfaces in Character Animation</article-title><source>Proceedings of SIGGRAPH 98</source><year>1998</year><publisher-name>Addison-Wesley</publisher-name><fpage>85</fpage><lpage>94</lpage></citation></ref><ref id="B16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bloor</surname><given-names>MIG</given-names></name><name><surname>Wilson</surname><given-names>MJ</given-names></name></person-group><article-title>Generating Blend Surfaces Using Partial Differential Equations</article-title><source>Computer-Aided Design</source><year>1989</year><volume>21</volume><fpage>165</fpage><lpage>171</lpage><pub-id pub-id-type="doi">10.1016/0010-4485(89)90071-7</pub-id></citation></ref><ref id="B17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ugail</surname><given-names>H</given-names></name><name><surname>Bloor</surname><given-names>MIG</given-names></name><name><surname>Wilson</surname><given-names>MJ</given-names></name></person-group><article-title>Techniques for Interactive Design Using the PDE Method</article-title><source>ACM Transactions on Graphics</source><year>1999</year><volume>18</volume><fpage>195</fpage><lpage>212</lpage><pub-id pub-id-type="doi">10.1145/318009.318078</pub-id></citation></ref><ref id="B18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Du</surname><given-names>H</given-names></name><name><surname>Qin</surname><given-names>H</given-names></name></person-group><article-title>Direct Manipulation and Interactive Sculpting of PDE Surfaces</article-title><source>Computer Graphics Forum (Proceedings of Eurographics 2000)</source><year>2000</year><volume>19</volume><fpage>261</fpage><lpage>270</lpage><pub-id pub-id-type="doi">10.1111/1467-8659.00418</pub-id></citation></ref><ref id="B19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ugail</surname><given-names>H</given-names></name><name><surname>Wilson</surname><given-names>MJ</given-names></name></person-group><article-title>Efficient Shape Parameterisa-tion for Automatic Design Optimisation using a Partial Differential Equation Formulation</article-title><source>Computers and Structures</source><year>2003</year><volume>81</volume><fpage>2601</fpage><lpage>2609</lpage><pub-id pub-id-type="doi">10.1016/S0045-7949(03)00321-3</pub-id></citation></ref><ref id="B20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bloor</surname><given-names>MIG</given-names></name><name><surname>Wilson</surname><given-names>MJ</given-names></name></person-group><article-title>Spectral Approximations to PDE Surfaces</article-title><source>Computer-Aided Design</source><year>1996</year><volume>28</volume><fpage>145</fpage><lpage>152</lpage><pub-id pub-id-type="doi">10.1016/0010-4485(95)00060-7</pub-id></citation></ref><ref id="B21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kubeisa</surname><given-names>S</given-names></name><name><surname>Ugail</surname><given-names>H</given-names></name><name><surname>Wilson</surname><given-names>M</given-names></name></person-group><article-title>Interactive Design Using Higher Order PDE's</article-title><source>The Visual Computer</source><year>2004</year><volume>20</volume><fpage>682</fpage><lpage>693</lpage><pub-id pub-id-type="doi">10.1007/s00371-004-0261-3</pub-id></citation></ref></ref-list></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Sante/explor/LymphedemaV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 001C18 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 001C18 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Wicri/Sante
|area= LymphedemaV1
|flux= Pmc
|étape= Corpus
|type= RBID
|clé= PMC:1198260
|texte= Modelling of oedemous limbs and venous ulcers using partial differential equations
}}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:16078992" \
| HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \
| NlmPubMed2Wicri -a LymphedemaV1
| This area was generated with Dilib version V0.6.31. |  |