A linear programming model for preserving privacy when disclosing patient spatial information for secondary purposes

  • Ho-Won Jung1Email author and

    Affiliated with

    • Khaled El Emam2, 3

      Affiliated with

      Contributed equally
      International Journal of Health Geographics201413:16

      DOI: 10.1186/1476-072X-13-16

      Received: 22 January 2014

      Accepted: 7 May 2014

      Published: 29 May 2014



      A linear programming (LP) model was proposed to create de-identified data sets that maximally include spatial detail (e.g., geocodes such as ZIP or postal codes, census blocks, and locations on maps) while complying with the HIPAA Privacy Rule’s Expert Determination method, i.e., ensuring that the risk of re-identification is very small. The LP model determines the transition probability from an original location of a patient to a new randomized location. However, it has a limitation for the cases of areas with a small population (e.g., median of 10 people in a ZIP code).


      We extend the previous LP model to accommodate the cases of a smaller population in some locations, while creating de-identified patient spatial data sets which ensure the risk of re-identification is very small.


      Our LP model was applied to a data set of 11,740 postal codes in the City of Ottawa, Canada. On this data set we demonstrated the limitations of the previous LP model, in that it produces improbable results, and showed how our extensions to deal with small areas allows the de-identification of the whole data set.


      The LP model described in this study can be used to de-identify geospatial information for areas with small populations with minimal distortion to postal codes. Our LP model can be extended to include other information, such as age and gender.


      Health services research Linear programming (LP) De-identified data sets Geographical identifiers HIPAA Privacy Rule


      Patients’ geographical identifiers (e.g., geocodes such as postal/ZIP codes, street addresses and locations on maps) are useful for health research and public health purposes [14]. Geographical identifiers are also fundamental to the practice of spatial epidemiology [5] and are key components of the public health professional’s toolbox [6].

      However, revealing patient data sets, including geographical identifiers, threatens patient privacy if the geographical identifiers can be linked to individuals. In fact, some studies have revealed a threat of re-identification. Sweeney [7] indicated that 87% of subjects could be uniquely identified by their gender, ZIP code and date of birth when linked with other publicly available data, such as voting records. Moskop et al. [8] presented that low-resolution dot maps of diseases published in several medical journals could be used to trace most patients to single addresses. Furthermore, Brownstein et al. [9] also showed that a method of georeferencing and unsupervised classification of the original image could be used to precisely re-identify 26% of 550 patients by using addresses from a presentation quality map and 79% using those from a publication quality map.

      The US Health Insurance Portability and Accountability Act of 1996 (HIPAA) allows the disclosure of personal health information for secondary purposes only if the patients provide authorization (with some exceptions) [10]. If it is not practical to obtain authorization then the data must be de-identified before disclosure [11]. Similar laws exist in Canada where de-identification is required for the disclosure of health information without consent [12, 13].

      According to the HIPAA Privacy Rule, de-identified protected health information (PHI) can be created by one of two ways [10], p.3. The first is the “safe-harbor” method, in which all 18 identifiers, including the five-digit ZIP codes, are removed. Yet, the first three digits of a ZIP code may be included, provided that at least 20,000 people share the same first three digits. The second way is “to have a qualified statistician determine, using generally accepted statistical and scientific principles and methods, that the risk is very small” concerning that such information could be used to identify an individual. The “very small” risk that is used as a threshold for disclosure control depends on the application fields and data users, but has a range of 0.05 to 0.3 of its value [1215]. This study sets a threshold value of 0.2.

      Studies, such as disease mapping or cluster detection in epidemiology, require de-identified data that maximally include the spatial distribution of a disease while complying with a threshold of the re-identification risk. A prevailing method to create de-identified data sets is to aggregate pre-defined areas, such as ZIP codes or counties, into a new area [16]. However, this approach loses useful spatial information while preserving privacy [17]. Furthermore, the level of privacy protection depends on the number of patient records [18]. Another approach uses the deterministic or stochastic function of geographical identifiers [19]. However, this heuristic method cannot quantify the risk to individual privacy and therefore cannot demonstrate that the risk is indeed “very small”.

      Wieland et al. [18] proposed a linear programming (LP) model to create de-identified data sets. The LP model determines the transition probability from an original location of a patient to a new randomized location as a de-identification method. However, it cannot be applied to data sets, including locations with small populations (e.g., the population is smaller than the number of patients). For example, the City of Ottawa has 11,740 postal codes that have a population of more than one. Of these, 98.61% (11,577) of postal codes have a population smaller than the number of patients in our data set (224 patients originated from 161 postal codes). The median population in Ottawa postal codes is 10 people.

      To apply this LP method on real data sets, where small areas will exist, this study revised the previous LP mode to accommodate the case in which some postal codes can have a smaller population than the total number of patients. The results depicted that our revised model can increase the applicability of the LP model in the creation of de-identified data sets.


      Results from WCMB-LP model

      This study solved two Ottawa LP problems using the WCMB-LP model.

      Table 1 shows the range of the maximum re-identification probability and the objective function values. In the table, optimal solutions had the s/ϵ range of 10 to 33, which corresponds to the ϵ value range of 22.4 to 6.79 (where the number of patients, s, is 224). That is, WCMB-LP provided impractical optimal solutions for which the re-identification probability is greater than 1.
      Table 1

      WCMB-LP results for the two Ottawa LP problems

      In equation (3), ν = s N · ϵ = 1 264 , 327 · s ϵ http://static-content.springer.com/image/art%3A10.1186%2F1476-072X-13-16/MediaObjects/12942_2014_Article_594_IEq1_HTML.gif

      Objective function value (unit: meter)

      Nearest 10

      Nearest 30

      s ϵ = 10 ϵ = 22.4 http://static-content.springer.com/image/art%3A10.1186%2F1476-072X-13-16/MediaObjects/12942_2014_Article_594_IEq2_HTML.gif

      ν = 3.783 × 10− 5



      s ϵ = 20 ϵ = 11.3 http://static-content.springer.com/image/art%3A10.1186%2F1476-072X-13-16/MediaObjects/12942_2014_Article_594_IEq3_HTML.gif

      ν = 7.567 × 10− 5



      s ϵ = 30 ϵ = 7.47 http://static-content.springer.com/image/art%3A10.1186%2F1476-072X-13-16/MediaObjects/12942_2014_Article_594_IEq4_HTML.gif

      ν = 11.350 × 10− 5



      s ϵ = 33 ϵ = 6.79 http://static-content.springer.com/image/art%3A10.1186%2F1476-072X-13-16/MediaObjects/12942_2014_Article_594_IEq5_HTML.gif

      ν = 12.485 × 10− 5



      Results from the revised model

      We then solved two Ottawa LP problems formulated on the basis of our revised LP model (Revised-LP). Table 2 shows their optimal solutions for the 224 patients. In the nearest 10 dataset, the revised model provided an optimal solution of less than 0.4 of the maximum re-identification probability across all postal code areas. These results were not considered acceptable for preserving privacy across all postal code areas due to the limited transition postal code areas to the nearest 10 neighbors. As described in the Introduction, this study established a threshold value of 0.2 for the re-identification probability.
      Table 2

      Revised-LP results for the two LP problems

      Re-identification probability

      Objective function value (unit: meter)

      Nearest 10

      Nearest 30

      ϵ = 0.6



      ϵ = 0.5



      ϵ = 0.4



      ϵ =0.3



      ϵ =0.2



      *Note that “Nearest 10” is infeasible when ϵ ≤ 0.4.

      When the transition postal code area was extended to the nearest 30, our LP model provided an optimal solution with patient movement of 1,686.3 meter for 0.2 re-identification probability. In this context, we considered that the transition over the nearest 30 would provide a smaller acceptable re-identification probability than the nearest 10. As expected, patient movement was increased for smaller re-identification probabilities, resulting in a greater loss of patient information. As reference, Groubi solution time in a desktop PC (Windows 7 and Intel Core i5 CPUs with 8G RAM) showed less than 4 seconds for Nearest 10 and 8.23 sec to 90 sec (ϵ =0.2) for Nearest 30.


      Because the area population and latitude and longitude are known for any given postal code, the LP model can generate the optimal transition probability if the number of patients is given. Patient movement in a ZIP code was assumed to follow a multinomial distribution with the transition probability. Thus, two different runs of the same LP problem may provide different patient movements with the same objective function values.

      As we have observed in our empirical studies, a limited number of transition neighbors, such as 10, can render LP models infeasible or impractical in terms of achieving an acceptable re-identification probability. However, increasing the transition neighbors greatly increases the computational burden of the LP problem to obtain the optimal solution, considering the postal codes in a country or region. Thus, it is essential to balance a reasonable number of neighbors with consideration for the LP problem size.


      This study expanded the applicability of the previous LP model regardless of the population across all locations (i.e., postal code areas). Thus, our model can be extended to include other information, such as age and gender. Future research may also include a comparison of the performance of our LP model with that of other methods, such as the previously described aggregation methods.


      An LP model for de-identified data sets

      Wieland et al. [18] introduced an LP model to transform a patient’s spatial identifiers to randomized identifiers in order to create de-identified data sets. In their study, a census block is the only spatial data to be de-identified. Because ZIP codes (postal codesa in Canada) are a common patient residence location indicator [20, 21], this study used ZIP codes as a spatial datum to be de-identified. In order to formulate an LP problem, the following notations are defined:

      A: Set of possible original ZIP codes as identifiers

      B: Set of possible randomized ZIP codes. This could be different from set A

      n i : Population in ZIP code i

      N: Sum of populations across all ZIP codes, i.e., ∑ i ∈ An i  = N

      d ij : Distance between ZIP codes i and j

      s: (Total) number of patients

      ϵ: Probability that any ZIP code from the randomized dataset originating from any specific individual in the underlying population is at most ϵ

      P ij : (Decision variable) Transition probability from an original ZIP code iA to a new ZIP code jB.

      Using an LP solution should ensure that the risk is “very small” while minimizing patient movement in order to reduce substantial information loss. With the notations defined, an LP model named WCMB-LP, where WCMB denotes the first character of each of the four authors’ last names in Wieland et al. [18], can be represented as follows:
      WCMB LP Min i A j B n i N · d ij · P ij subject to j B P ij = 1 , for all i A ; k A n k N . P kj s . P ij 0 , for all i A and j B ; 3 P ij 0 , for all i A and j B . http://static-content.springer.com/image/art%3A10.1186%2F1476-072X-13-16/MediaObjects/12942_2014_Article_594_Equa_HTML.gif

      The objective function in equation (1) minimizes the expected total movement distance of patients, where n i /N denotes a probability that a patient originated from ZIP code i. The constraint in equation (2) specifies that the patients in A should be moved to somewhere in B. This may include self-transition, i.e., patients in a ZIP code may remain there. Constraints in equation (3) implies that “Given the set of s locations comprising the de-identified dataset, the probability that any one of these derived from one specific individual to be at most ϵ. This is guaranteed if the probability that a location from the randomized dataset originated from an arbitrary specific individual is required to be at most ϵ” [18], p. 17612. Further, the transition probability P ij in equation (4) should be greater than or equal to zero. When the decision variable P ij is obtained, patients in ZIP code i are moved to ZIP code j using a multinomial distributionb.

      Three cases can be investigated to improve the understandability of equation (3) as follows:

      [Case 1]: If all ZIP codes include just one person, i.e., n i  = 1 for all i, equation (3) becomes P ij  ≤ ϵ/s.

      [Case 2]: If there is just one ZIP code, there is no transition probability, i.e., equation (3) is reduced to s/N ≤ ϵ

      [Case 3]: If all patients having a randomized ZIP code j are from i, i.e., P kj  = 0 for all k ≠ i, equation (3) becomes s/n i  ≤ ϵ, where n i  = N. This is the same as [Case 2].

      Revised LP model

      For a simple interpretation, equation (3) can be rewritten as follows:
      1 n i · n i N · P ij k A n k N · P kj ϵ s , for all i A and j B . http://static-content.springer.com/image/art%3A10.1186%2F1476-072X-13-16/MediaObjects/12942_2014_Article_594_Equ1_HTML.gif

      Equation (5) is the same equation (5) in WCMB-LP [18], p. 17612, where the first part 1/n i implies the probability that “all individuals in ZIP code i with population n i have an equal chance of having the disease… and the second term is a population-weighted transition probability.”

      Equation (5) can also be rewritten as follows:
      s n i · n i N · P ij k A n k N · P kj ϵ , for all i A and j B . http://static-content.springer.com/image/art%3A10.1186%2F1476-072X-13-16/MediaObjects/12942_2014_Article_594_Equ2_HTML.gif

      In equation (6), the right-hand side ϵ means that all patients in B have the same randomized ZIP code, i.e., a randomized patient list of s patients includes one ZIP code. Contrast to equation (5), s/n i denotes a maximum re-identification probability of s patients with the same randomized ZIP code, assuming its origination from i. In this context, the number of patients cannot exceed the number of people in ZIP code i. That is, s/n i  ≤ 1.

      With these elaborations, equation (3) in WCMB-LP is revised as follows:
      min s n i , 1 . n i N · P ij k A n k N · P kj ϵ , for all i A and j B . http://static-content.springer.com/image/art%3A10.1186%2F1476-072X-13-16/MediaObjects/12942_2014_Article_594_Equ3_HTML.gif
      Rearranging equation (7) and incorporating it into WCMB-LP, we have the following Revised - LP model (our LP model):
      Revised LP Min i A j B n i N · d ij · P ij subject to j B P ij = 1 , for all i A ; k A n k N · P kj min s , n i · P ij 0 , for all i A and j B ; P ij 0 , for all i A and j B . http://static-content.springer.com/image/art%3A10.1186%2F1476-072X-13-16/MediaObjects/12942_2014_Article_594_Equ4_HTML.gif

      Note that equation (8) is the difference between WCMB-LP and our Revised LP.

      Properties of re-identification constraint

      In order to investigate the re-identification constraint further, let v = min(s/, n i /) in equation (8). Then, the constraints in equation (8) can be rewritten as:
      k A n k N · P kj ν · P ij 0 ; k A ; k i n k N · P kj + n i N ν · P ij 0 , for all i A and j B , http://static-content.springer.com/image/art%3A10.1186%2F1476-072X-13-16/MediaObjects/12942_2014_Article_594_Equb_HTML.gif
      where the first part is non-negative because P kj  ≥ 0 for all k and j, and the second part can be represented by the following function g of variable n i :
      g n i = n i N ν = n i N n i = n i N 1 1 ϵ , if s n i 1 i . e . , n i s ; n i N s = 1 N n i s ϵ , if s n i < 1 i . e . , n i > s . http://static-content.springer.com/image/art%3A10.1186%2F1476-072X-13-16/MediaObjects/12942_2014_Article_594_Equ5_HTML.gif
      Function g in equation (9) can be represented as in Figure 1, where g(n i ) has the smallest value of (s/N)(1 − 1/ϵ) when population n i  = s holds and then increases across zero at n i  = s/ϵ.
      Figure 1

      The functional form of equation(9)

      The leftmost region from 1 to s in Figure 1, i.e., n i s, corresponds to the first equation in equation (9). In that region, function g(n i ) always has a negative value and is decreasing because ϵ is assumed to be less than 1. The middle region from s to s/ϵ  (i.e., n i  >s) corresponds to the second equation in equation (9), where function g(n i ) is negative and is increasing. The rightmost region (the second equation in equation (9)) denotes that the function g(n i ) for n i s/ϵ always has a positive value. Thus, the corresponding constraints in equation (8) of our Revised-LP model are always satisfied, i.e., they are redundant because all of its corresponding constraints’ coefficients are nonnegative. A redundant constraint is one that can be left out without changing the model.


      In this study, a data set called Ottawa, which includes only areas with a population of more than one, as in Wieland et al. [18], was applied to both WCMB-LP and our Revised-LP models. Our data set was based on patients’ information in a population of 264,327 children under the age of 18 residing in Ottawa, Canada. Our purpose was to randomize the postal codes of a patient list in CHEO (Children’s Hospital of Eastern Ontario) presenting in the emergency department. The patient list included 224 patients from 126 ZIP codes, in which the number of patients corresponded to 5% of an estimated 4,500 people who visited CHEO in a month during the height of the influenza season. The patients were chosen from a pool of CHEO patient postal codes.

      The area of each postal code was represented by the centroid latitude and longitude. The distance between two postal-code areas was computed by using the Haversine formula [22], which provides the shortest (also termed ‘as-the-crow-flies’ ignoring any hill or great-circle) distance between any two points on a spherical earth from their longitudes and latitudes. Ellipsoidal effects are ignored, but the result is sufficiently accurate for the purpose of the present study.

      Because our data set included 11,740 postal codes, the LP formulation had 137,827,600 variables (i.e., 11,7402) and 137,839,340 constraints [i.e., 11,740 (1 + 11,740)]. In order to reduce the size of this LP problem, transitions from any postal code area were limited to the following two cases: the nearest 10 and 30 postal code areas, i.e., two LP problems with 117,400 and 352,200 variables, and 129,140 and 360,940 constraints, respectively. The two LP problems were solved by using Gurobi 6.5.2 solver [23] with MPL 4.2n modeling language [24].

      Ethics approval for this study was obtained from the CHEO research Institute research ethics board.


      aThis study interchangeably uses term “postal” and “ZIP” codes. However, when we mention data from Canada, the term postal codes are intentionally used.

      bAn R library [25] has a command of generating a multinomially distributed random number in r. WCMB-LP has |A||B| variables and |A| + |A||B| constraints.



      The authors would like to thank the staff members at the Electronic Health Information Laboratory at the CHEO Research Institute. The research was supported by Korea University Business School (2013). This support is gratefully acknowledged.

      Authors’ Affiliations

      Korea University Business School
      Children’s Hospital of Eastern Ontario Research Institute
      Pediatrics, Faculty of Medicine, University of Ottawa


      1. Boulos M: Towards evidence-based, GIS-driven national spatial health information infrastructure and surveillance services in the United Kingdom. Int J Health Geogr. 2004, 3: 1-10.1186/1476-072X-3-1. [http://​www.​ij-healthgeographic​s.​com/​content/​3/​1/​1]PubMed CentralView ArticlePubMed
      2. Cromley EK: GIS and disease. Annu Rev Public Health. 2003, 24: 7-24. 10.1146/annurev.publhealth.24.012902.141019.View ArticlePubMed
      3. Croner CM: Public health, GIS, and the Internet. Annu Rev Public Health. 2003, 24: 57-82.View ArticlePubMed
      4. McLafferty SL: GIS and health case. Annu Rev Public Health. 2003, 24: 25-42. 10.1146/annurev.publhealth.24.012902.141012.View ArticlePubMed
      5. Cassa C, Wieland S, Mandl K: Re-identification of home addresses from spatial locations anonymized by Gaussian skew. Int J Health Geogr. 2008, 7: 45-10.1186/1476-072X-7-45. [http://​www.​ij-healthgeographic​s.​com/​content/​7/​1/​45]PubMed CentralView ArticlePubMed
      6. AbdelMalik P, Boulos M, Jones R: The perceived impact of location privacy: A web-based survey of public health perspectives and requirements in the UK and Canada. BMC Public Health. 2008, 8: 156-10.1186/1471-2458-8-156. [http://​www.​biomedcentral.​com/​1471-2458/​8/​156]PubMed CentralView ArticlePubMed
      7. Sweeney L: k-anonymity: A model for protecting privacy. Int J Uncertainty Fuzziness Knowl Based Syst. 2002, 10: 557-570. 10.1142/S0218488502001648.View Article
      8. Moskop JC, Marco CA, Larkin GL, Geiderman JM, Derse AR: From Hippocrates to HIPAA: Privacy and confidentiality in emergency medicine - Part I: Conceptual, moral, and legal foundations. Ann Emerg Med. 2005, 45: 53-59. 10.1016/j.annemergmed.2004.08.008.View ArticlePubMed
      9. Brownstein J, Cassa C, Kohane I, Mandl K: An unsupervised classification method for inferring original case locations from low-resolution disease maps. Int J Health Geogr. 2006, 5: 56-10.1186/1476-072X-5-56. [http://​www.​ij-healthgeographic​s.​com/​content/​5/​1/​56]PubMed CentralView ArticlePubMed
      10. National Institutes of Health: Research Repositories, Databases, and the HIPAA Privacy Rule. U.S. Department of Human and Health Services, [http://​goo.​gl/​rR28ob]
      11. National Institutes of Health: Dictionary. US Dep Health Hum Serv. [http://​privacyruleandre​search.​nih.​gov/​dictionary.​asp]
      12. Statistics Canada: Therapeutic Abortion Survey.http://​goo.​gl/​v01DsY,
      13. Ministry of the Attorney General: ORDER PO-2037, Appeal PA-010381-1. [http://​goo.​gl/​vuAFFI]
      14. El Emam K: Heuristics for de-identifying health data. IEEE Secur Priv. 2008, 6: 58-61.View Article
      15. Howe H, Lake A, Lehnherr M, Roney : Unique record identification on public use files as tested on the 1994–998 CINA analytic file. North Am Assoc Centr Cancer Registr. 2002, 2002: 2002-[http://​goo.​gl/​nbq6e7]
      16. Fefferman NH, O’Neil EA, Naumova EN: Confidentiality and confidence: Is data aggregation a means to achieve both?. J Public Health Policy. 2005, 26: 430-449. 10.1057/palgrave.jphp.3200029.View ArticlePubMed
      17. Olson KL, Grannis SJ, Mandl KD: Privacy protection versus cluster detection in spatial epidemiology. Am J Public Health. 2006, 96: 2002-2008. 10.2105/AJPH.2005.069526.PubMed CentralView ArticlePubMed
      18. Wieland SC, Cassa CA, Mandl KD, Berger B: Revealing the spatial distribution of a disease while preserving privacy. Proc Natl Acad Sci U S A. 2008, 105: 17608-17613. 10.1073/pnas.0801021105.PubMed CentralView ArticlePubMed
      19. Armstrong MP, Rushon G, Zimmerman DL: Geographically masking health data to preserve confidentiality. Stat Med. 1999, 18: 497-525. 10.1002/(SICI)1097-0258(19990315)18:5<497::AID-SIM45>3.0.CO;2-#.View ArticlePubMed
      20. Ng E, Wilkins R, Perras A: How far is it to the nearest hospital? Calculating distances using the Statistics Canada Postal Code Conversion file. Health reports/Statistics Canada, Canadian Centre for Health Information. 1993, 5: 179-183.
      21. Demissie K, Hanley J, Menzies D, Joseph L, Ernst P: Agreement in measuring socioeconomic status: Area-based versus individual measures. Chronic Dis Can. 2000, 21: 1-7.PubMed
      22. Sinnott RW: Virtues of the haversine. Sky Telescope. 1984, 68: 159-
      23. Gurobi: Gurobi Optimizer Reference Manual. 2009, [http://​www.​gurobi.​com]
      24. MPL: MPL Modeling System. [http://​www.​maximalsoftware.​com/​mplman/​]
      25. R Stats Package: The Multinomial Distribution. 2013, [http://​goo.​gl/​reeZor]


      © Jung and El Emam; licensee BioMed Central Ltd. 2014

      This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.