Density estimation and adaptive bandwidths: a primer for public health practitioners.  
Jump to Full Text  
MedLine Citation:

PMID: 20653969 Owner: NLM Status: MEDLINE 
Abstract/OtherAbstract:

BACKGROUND: Geographic information systems have advanced the ability to both visualize and analyze point data. While pointbased maps can be aggregated to differing areal units and examined at varying resolutions, two problems arise 1) the modifiable areal unit problem and 2) any corresponding data must be available both at the scale of analysis and in the same geographic units. Kernel density estimation (KDE) produces a smooth, continuous surface where each location in the study area is assigned a density value irrespective of arbitrary administrative boundaries. We review KDE, and introduce the technique of utilizing an adaptive bandwidth to address the underlying heterogeneous population distributions common in public health research. RESULTS: The density of occurrences should not be interpreted without knowledge of the underlying population distribution. When the effect of the background population is successfully accounted for, differences in point patterns in similar population areas are more discernible; it is generally these variations that are of most interest. A static bandwidth KDE does not distinguish the spatial extents of interesting areas, nor does it expose patterns above and beyond those due to geographic variations in the density of the underlying population. An adaptive bandwidth method uses background population data to calculate a kernel of varying size for each individual case. This limits the influence of a single case to a small spatial extent where the population density is high as the bandwidth is small. If the primary concern is distance, a static bandwidth is preferable because it may be better to define the "neighborhood" or exposure risk based on distance. If the primary concern is differences in exposure across the population, a bandwidth adapting to the population is preferred. CONCLUSIONS: Kernel density estimation is a useful way to consider exposure at any point within a spatial frame, irrespective of administrative boundaries. Utilization of an adaptive bandwidth may be particularly useful in comparing two similarly populated areas when studying health disparities or other issues comparing populations in public health. 
Authors:

Heather A Carlos; Xun Shi; James Sargent; Susanne Tanski; Ethan M Berke 
Related Documents
:

21520719  Biosensors for secondary metabolites, two case studies: ochratoxin a and microcystin. 20331949  Venous thromboembolism: a public health concern. 20653969  Density estimation and adaptive bandwidths: a primer for public health practitioners. 16160379  Ontological definition of population for public health databases. 19472059  Finding targets for obesity intervention in urban communities: schoolbased health cent... 21399429  Parallels of radiation and financial risk management: impacts on public acceptance. 12959889  Responding to a bioterrorist attack: environmental investigation of anthrax in new jersey. 24847119  Why public moralities matterthe relevance of socioempirical premises for the ethical d... 9460829  Pesticide poisoning: a major health problem in sri lanka. 
Publication Detail:

Type: Journal Article; Research Support, N.I.H., Extramural Date: 20100723 
Journal Detail:

Title: International journal of health geographics Volume: 9 ISSN: 1476072X ISO Abbreviation: Int J Health Geogr Publication Date: 2010 
Date Detail:

Created Date: 20100813 Completed Date: 20101130 Revised Date: 20140912 
Medline Journal Info:

Nlm Unique ID: 101152198 Medline TA: Int J Health Geogr Country: England 
Other Details:

Languages: eng Pagination: 39 Citation Subset: IM 
Export Citation:

APA/MLA Format Download EndNote Download BibTex 
MeSH Terms  
Descriptor/Qualifier:

Algorithms Humans Population Density* Population Surveillance / methods* Public Health Administration* Texas 
Grant Support  
ID/Acronym/Agency:

1K23AG036934/AG/NIA NIH HHS; AA015591/AA/NIAAA NIH HHS; CA 077026/CA/NCI NIH HHS; K23 AG036934/AG/NIA NIH HHS; P30 CA023108/CA/NCI NIH HHS; R01 AA015591/AA/NIAAA NIH HHS; R01 AA01559107/AA/NIAAA NIH HHS; R01 CA077026/CA/NCI NIH HHS; R01 CA07702614/CA/NCI NIH HHS 
Comments/Corrections 
Full Text  
Journal Information Journal ID (nlmta): Int J Health Geogr ISSN: 1476072X Publisher: BioMed Central 
Article Information Download PDF Copyright ©2010 Carlos et al; licensee BioMed Central Ltd. openaccess: Received Day: 20 Month: 4 Year: 2010 Accepted Day: 23 Month: 7 Year: 2010 collection publication date: Year: 2010 Electronic publication date: Day: 23 Month: 7 Year: 2010 Volume: 9First Page: 39 Last Page: 39 Publisher Id: 1476072X939 PubMed Id: 20653969 DOI: 10.1186/1476072X939 
Density estimation and adaptive bandwidths: A primer for public health practitioners  
Heather A Carlos1  Email: heather.a.carlos@dartmouth.edu 
Xun Shi2  Email: Xun.Shi@dartmouth.edu 
James Sargent1  Email: James.Sargent@dartmouth.edu 
Susanne Tanski1  Email: Susanne.Tanski@dartmouth.edu 
Ethan M Berke1234  Email: Ethan.Berke@TDI.dartmouth.edu 
1Norris Cotton Cancer Center, DartmouthHitchcock Medical Center, Lebanon, NH, USA 

2Department of Geography, Dartmouth College, Hanover NH, USA 

3Department of Community and Family Medicine, Dartmouth Medical School, Hanover NH, USA 

4Prevention Research Center at Dartmouth, The Dartmouth Institute for Health Policy and Clinical Practice, Lebanon NH, USA 
From John Snow's Victorian era map of cholera deaths [^{1}] to interactive maps tracking the spread of H1N1 Influenza [^{2}], spatial point patterns have a long and rich history in the public health arena. Disease registries now include geolocation data, which allow detection of clustering (a global tendency) and clusters (a local phenomenon). Public health practitioners focusing on disease prevention use spatial point pattern analysis to quantify social determinants of health (for example, distance to sites of physical activity [^{3}] or to retail outlets [^{4}], discrepancies in access to services by race or ethnicity, or variation in educational attainment).
Geographic information systems (GIS) have advanced the ability to both visualize and analyze these point data. Using GIS, point based maps can easily be aggregated to differing areal units and examined at varying resolutions. This however, creates problems in spatial analysis. In addition to introducing the modifiable areal unit problem (MAUP) [^{5}], where altering the area or shape of an aggregate unit may alter the value within the polygon, any corresponding demographic data must also be available both at the scale of analysis and in the same geographic units. One way to address these issues is to employ kernel density estimation (KDE) techniques rather than geographic aggregation [^{6}^{}^{8}]. KDE is a nonparametric method of extrapolating point data over an area of interest without invoking MAUP or relying on fixed boundaries for aggregation. The density of points is calculated using a specified bandwidth (a circle of a given radius centered at the focal location). This produces a smooth, continuous surface where each location in the study area is assigned a density value, which can then be used as the independent or dependent variable in statistical models. KDE's strength is its ability to provide an estimate of density at any location in the spatial frame (e.g. a geocoded subject or another point of interest), irrespective of arbitrary administrative boundaries.
While various methods exist for calculating KDE surfaces, including some embedded in common GIS software, many public health practitioners and researchers use a static distance for bandwidth patterned after the casecontrol method [^{7},^{9}]. A more indepth discussion on the use of KDE in public health can be found in our prior work [^{10},^{11}] and that of others [^{12}]. When geographic distance (and count of cases) is the primary interest, a static bandwidth KDE appropriately represents the density of a particular attribute, for example understanding how relapse of alcoholism may be predicted in part by proximity to bars or pubs. However, for some health outcomes, a fixed geographic distance is not the appropriate bandwidth. Consider the hypothesis that alcohol outlets (retail alcohol sale locations) are more concentrated in lowincome neighborhoods within a metro area. The problem with using a static bandwidth for each outlet is that we expect the greater density of alcohol outlets in urban areas where the population density is higher than in the suburbs. To the extent that outlet density in poor urban neighborhoods is just a reflection of a higher concentration of people living there, the correlation does not necessarily point to a health disparity.
There are a wide range of analytical methods available to examine spatial point patterns [^{12}] and other researchers have considered the effect of inhomogeneous background populations. Here, inhomogeneous background populations refers not to the population's demographics, rather to the distribution of the source of the event points. For our alcohol outlet example, the background population is the population of the study area, whereas in an analysis of disease cases, only the population at risk is used. Notably, the spatial filtering technique [^{13},^{14}] uses both fixed and adaptive filters/bandwidths to test or map the relationship between a count of cases and a background population while the cluster evaluation permutation procedure [^{15},^{16}] uses an adaptive bandwidth but focuses only on the case count. Alternatively, adding population density as a covariate in a statistical model could address this issue, but a more elegant solution incorporates population density into the outcome variable using a bandwidth that represents the underlying population, rather than a fixed geographical area. An adaptive bandwidth, discussed below, may be preferable when studying issues of population and variations in exposure.
In this paper, we present a number of approaches to density estimation and propose using a KDE method to address uneven population distribution by using an adaptive bandwidth specified by the underlying population. This technique is useful when it is it is important to understand if a density value is just a reflection of the local population or if it may stem from other causes. This methodology was motivated by an analysis of the distribution of alcohol outlets. An illustrative application in that arena is used to compare density methods, but it is also applicable to the analysis of the density of disease, crime, healthcare clinics and other fields where the background population is heterogeneous.
We focus on a suite of density estimation tools: point density, static bandwidth kernel density estimation and adaptive bandwidth kernel density estimation. Density calculations operate on either cases or sites. Cases are event points (e.g. addresses of alcohol outlets or disease cases) whereas sites represent all locations (pixels or each point on a grid) in a study area. Density calculations performed on sites (siteside method) evaluate the density for every location in the study area, whereas the caseside method only looks at case locations and their defined surrounding locations. In order to highlight the differences in the density estimation tools presented here, we limit our discussion to caseside methods of density estimation. Caseside methods are more computationally efficient and in many situations better represent the nature of the application problem. More information about the differences between caseside and siteside methods can be found in Shi [^{11}].
Most generically, a point density function (also called intensity function) defines the number of cases (alcohol outlets, disease cases) per unit area at each location throughout an area of interest. To calculate this density surface, for each case, a "neighborhood" is delineated, usually by defining a search radius (or bandwidth); the number of cases that fall within the neighborhood are divided by the area of the neighborhood; this value is assigned to the neighborhood (Figure 1). The intensity function is expressed as [^{17}]
(1)
where λ(x,y) is the intensity (or point density) at location(x,y), n is the number of events and A is the area of the neighborhood. When neighborhoods overlap, the results are summed to indicate a higher density of cases. The units of λ(x,y) are cases per unit area.
When points are evenly distributed in space, increasing the bandwidth does not have a large impact on the intensity since, as larger neighborhoods are defined, n will likely increase, but so will A. However, increasing the bandwidth does provide a greater smoothing effect (or a more generalized surface), which risks removing meaningful spikes (peaks or valleys) or edges (extent of the influence of a case) from the original data distribution.
Although the point density function is relatively simple and straightforward, it does not convey any information about the spatial configuration of features of interest within the bandwidth. Consider two locations (sites) and one case. Computationally, a site coincident with a case returns the same λ as a site one bandwidth away from the case. While this approach is appropriate for studies which are interested in the number of events per unit area at a specified location (e.g. crime events, residential or population density), in other disciplines there is an expected attenuation with distance (e.g. environmental pollutants which dissipate as they travel from the source). In order to compensate for distance, a density function can incorporate a decay function to assign smaller values to locations which are still in the neighborhood, but more distant from a case. This is the approach employed by kernel density estimation.
Kernel density estimation fits a curved surface over each case such that the surface is highest above the case and zero at a specified distance (the bandwidth) from the case (Figures 2 and 3). In mathematical terms, it is expressed as [^{6}]
where f(x,y) is the density value at location (x,y), n is the number of cases, h is the bandwidth, d_{i }is the geographical distance between case i and location (x, y) and K is a density function (generally a radially symmetric unimodal probability density function) which integrates to one. The units of f(x,y) are cases per unit area.
Static bandwidth kernel density estimation is a technique that is appropriate when geographic distance (and case count) is the primary concern. Since it applies the same geographic extent to each case, static bandwidth KDE does not distinguish the spatial extents of interesting areas [^{11}], nor does it expose patterns above and beyond those due to geographic variations in the density of the underlying population [^{17}].
Because health outcomes involve people, their spatial distribution will often reflect the spatial distribution of the underlying human population. Counts of disease are almost always higher in urban areas than rural areas simply due to the size of the potential exposed population. Likewise, counts of things people use are greater in higher population areas: there are more parks, physical activity sites, and retail outlets in places where more people live. As a result, the density of occurrences should not be interpreted without knowledge of the underlying population distribution [^{7}]. When the effect of the background population is successfully accounted for, differences in point patterns in similar population areas are more discernible; it is generally these variations that are of most interest [^{18}].
A Texas case study illustrates the problem posed by heterogeneous population backgrounds. Figure 4 displays the point data for alcohol outlets while Figure 5 shows the static bandwidth KDE surface for these data. As we would expect, both maps are similar to an image of Texas at night (Figure 6) since they replicate the population distribution. In contrast, Figure 7 shows a density map with the underlying population addressed. An adaptive bandwidth KDE method was used to create this map and it is described below.
Before delving into the adaptive bandwidth methodology, we discuss a data source for the background population. Most GIS based population data are in a polygon format with a population count (or estimate) assigned to each polygon. Depending on the study area and data source, each polygon may be as large as a country, or as small as a city block. These polygons are often irregular shapes and sizes and lack data about how people are geographically dispersed within the polygon. In addition, administrative boundaries may not be consistent with the travel patterns or service utilization of those that live in them. An alternative to polygon based population data is the LandScan™ Global Population Database, which was developed using multiple techniques to disaggregate census counts within an administrative boundary. This worldwide population data product is available on a 30" × 30" latitude/longitude grid (a pixel located in the central United States is approximately 0.65 km^{2}). The advantage of the grid format is that it regularizes the areal unit for population values, unlike administrative boundaries, which vary in size. This makes counts at different locations more spatially comparable and facilitates spatial analysis operations. However, it disconnects the population counts from the related demographic data, which is included in many censuses. Additionally, in urban areas, census blocks may be smaller than the LandScan™ grid and as a result, the larger grid units aggregate the original population counts. A description of the LandScan™ data and the methodology used to create it are described in detail elsewhere [^{19},^{20}].
Geocoded data for the alcohol outlets were obtained from the NAICS (North American Industry Classification System) Association http://www.naics.com. Details on this data and related processing can be found elsewhere [^{4}].
Whereas the static bandwidth kernel density estimation model employs a bandwidth based on a geographic distance, the adaptive bandwidth method uses background population drawn from LandScan™ data to calculate a kernel of varying size for each individual case (which, using the examples above could be an alcohol outlet). This limits the influence of a single case to a small spatial extent where the population density is high as the bandwidth is small [^{10}]. Conversely, in rural areas where population is lower, the kernel is geographically larger and the influence of a single case is greater (Figure 8).
The adaptive method is calculated as follows [^{10}]:
Note that there are several differences between equations 2 and 3. Most notably, in the adaptive method, the bandwidth is represented by P(u,v) which is a function centered on the case located at (u,v) and based on the local population. Additionally, the denominator nh^{2 }is dropped since the output value is not normalized by the geographic area of the kernel (h^{2}). The adaptive bandwidth method bases the influence of a case on the underlying population support, not the area support. There are various choices for the function K; most will not significantly affect the outcome. This study uses a simple form:
(4)
Often the influence of one case will overlap that of another. In this situation, the bandwidth and kernel density estimation calculations are performed separately for each case and then the results are summed (Figure 8, Equation 3).
The expected population parameter determines the extent of the adaptive bandwidth. The expected population normalizes the influence of each case to a certain number of people and thus the bandwidth stops expanding when the expected population is reached. In less populous areas however, the bandwidth could expand beyond the distance where a case may influence health; one can therefore set a limit to the maximum distance of the bandwidth. The maximum distance parameter restricts the bandwidth from expanding further, even if the expected population has not been reached. This may be critical when considering health behaviors influenced by exposures that are beyond a reasonable distances from an individual.
Bandwidth is determined for each individual case by summing the underlying population, starting with the pixel directly under the case and then expanding outward until the expected population is reached. Given reasonable values for population and maximum distance, expected population exerts the most control. In urban areas with high population densities, the expected population limit will often be reached before the maximum distance, thus, by adjusting the expected population, the radius of the influence of the case will diminish. The same is true in rural areas, except that in areas of very low population density the maximum distance may be called into play to limit bandwidth.
The difference between static and adaptive bandwidth KDE methods is best illustrated through visualization. Figure 9 portrays the results from each method applied to alcohol outlets (Figure 9D) in the area surrounding San Antonio, Texas. As mentioned above, the static bandwidth KDE surfaces (Figures 9A and 9E) excel at identifying areas where there are many point sources, but they do not provide a basis for discerning where the point sources are higher or lower than would be expected given the underlying population (Figure 9C). The adaptive bandwidth KDE (Figure 9B and 9F) addresses these issues through utilization of a populationbased bandwidth, allowing for improved detection of neighborhoodlevel differences in exposure, even in areas that have similarly high raw counts of alcohol outlets. This is illustrated in the closeup of San Antonio (Figures 9E and 9F) where the adaptive bandwidth KDE (Figure 9F) shows finegrained variability in the urban center whereas the static bandwidth KDE (Figure 9E) shows little differentiation in alcohol outlet density. This level of analysis is important in associating density of exposure with markers of health disparities, such as poverty, as seen in Figure 9G.
In deciding which method to choose, one needs to consider the research hypothesis. Using alcohol outlets as an example, if the primary concern is distance, a static bandwidth is preferable because it may be better to define the "neighborhood" or exposure risk of each store based on distance. If the primary concern is the "share" of the "service" per person, or differences in exposure across the population, a bandwidth adapting to the population is preferred.
More theoretically, the static bandwidth has a fixed spatial certainty, but varying statistical stability across the area, and thus is more suitable for an application primarily concerning distance. The adaptive bandwidth has a fixed statistical stability but varying spatial certainty. As a result, in a high population density area, it has both high statistical stability (specified by the user) and high spatial certainty (e.g., can better reveal the size of a hot spot), but in a low population density area, its high statistical stability comes with a cost of low spatial certainty.
There are limitations to all methods of spatial analysis, including density estimation, which induces interpolation autocorrelation which may result in over smoothing [^{21}] (this can be controlled by using a global Monte Carlo simulation). Perhaps the greatest limitation is the relatively arbitrary selection of bandwidth limits with both static and adaptive methods. Too large or small a bandwidth poses the risk of over or undersmoothing the original data, respectively. As subsequent analyses are based on this estimated density information, as opposed to original points, a change in bandwidth may have a significant impact on statistical relationships between dependent and independent variables. Methods to estimate appropriate bandwidths are described in more detail elsewhere [^{6},^{17},^{22}]. We recommend, even when applying mathematical models to estimate bandwidth, that one test multiple parameters for bandwidth in a sensitivity analysis. When employing a populationbased adaptive bandwidth, applying a distance limit similar to the example above may be useful when considering the influence of exposure on behavior. In the example of alcohol outlets used above, we placed a 25 km limit on the density calculation if the expected population threshold of 1000 people was not reached. We determined this maximum distance by testing a number of distances as well as choosing a limit based on behavior theory regarding alcohol acquisition and consumption. For a process other than alcohol exposure, a different distance limit may make more sense.
Even though techniques such as static or adaptive bandwidth KDE do not rely on aggregated data or administrative boundaries, issues of data visualization remain. One should use caution when viewing density data at a small scale. This is particularly true when viewing U.S. national maps where variations in density in high population but geographically smaller, northeastern areas are difficult to visualize, leading to visual bias towards the lower population areas of the west. Indeed, display of a density map may not be appropriate at all unless the proper scale is chosen. Quantitative data from the KDE may be better reported in a tabular presentation. Finally, the analyses in this study were conducted using a program created by one of the authors, as opposed to commercially available software. With recognition of multiple methods of KDE and its use on the rise, we expect that these approaches will appear in common GIS and spatial statistics software in the near future.
Researchers in the health sciences should be aware that multiple approaches to density estimation exist. Kernel density estimation is a useful way to consider exposure at any location within a spatial frame, irrespective of administrative boundaries. The ability for the researcher to analyze data easily at multiple scales reduces the risk of misinterpretation of results due to the MAUP. Utilization of an adaptive bandwidth may be particularly useful in comparing two similarly populated areas when studying health disparities.
The authors declare that they have no competing interests.
HC conceived of the study, performed the analysis, and drafted the manuscript. XS designed and implemented the population based adaptive bandwidth method and assisted with manuscript preparation. JS assisted with data interpretation and manuscript preparation, ST assisted with data interpretation and manuscript preparation. EB conceived of the study, supervised analyses, and drafted the manuscript. All authors read and approved the final manuscript.
This work was supported by the National Institutes of Health (AA015591 and CA 077026). Ethan Berke is supported by 1K23AG036934.
References
Snow J,On the mode of communication of choleraYear: 18552London,: J. Churchill  
Tracking Swine Flu Cases WorldwideThe New York TimesYear: 2009http://www.nytimes.com/interactive/2009/04/27/us/20090427fluupdategraphic.html  
Maroko AR,Maantay JA,Sohler NL,Grady KL,Arno PS,The complexities of measuring access to parks and physical activity sites in New York City: a quantitative and qualitative approachInt J Health GeogrYear: 200983410.1186/1476072X83419545430  
Berke EM,Tanski S,AlfordTeaster J,Shi X,Sargent J,Retail Alcohol Density and Poverty in Urban U.S. Census Tracts: A Geographic AnalysisAmerican Journal of Public HealthYear: 2010 in press .  
Openshaw S,Ecological fallacies and the analysis of areal census dataEnviron Plan AYear: 198416173110.1068/a16001712265900  
Silverman BW,Density estimation for statistics and data analysisYear: 1986London; New York: Chapman and Hall  
Bithell JF,An application of density estimation to geographical epidemiologyStatistics in MedicineYear: 1990969170110.1002/sim.47800906162218172  
Guagliardo MF,Spatial accessibility of primary care: concepts, methods and challengesInt J Health GeogrYear: 20043310.1186/1476072X3314987337  
Bithell JF,A classification of disease mapping methodsStatistics in MedicineYear: 2000192203221510.1002/10970258(20000915/30)19:17/18<2203::AIDSIM564>3.0.CO;2U10960848  
Shi X,A Geocomputational Process for Characterizing the Spatial Pattern of Lung Cancer Incidence in New HampshireAnnals of the Association of American GeographersYear: 20099952153310.1080/00045600902931801  
Shi X,Selection of bandwidth type and adjustment side in kernel density estimation over inhomogeneous backgroundsInternational Journal of Geographical Information ScienceYear: 20102464366010.1080/13658810902950625  
Gatrell AC,Bailey TC,Diggle PJ,Rowlingson BS,Spatial Point Pattern Analysis and Its Application in Geographical EpidemiologyTransactions of the Institute of British GeographersYear: 19962125627410.2307/622936  
Rushton G,Lolonis P,Exploratory spatial analysis of birth defect rates in an urban populationStat MedYear: 19961571772610.1002/(SICI)10970258(19960415)15:7/9<717::AIDSIM243>3.0.CO;209132899  
Tiwari C,Fisher PFUsing Spatially Adaptive Filters to map Late Stage Colorectal cancer Incidence in IowaDevelopments in spatial data handling: 11th International Symposium on Spatial Data HandlingYear: 2005Berlin; New York: Springer676 xix.  
Talbot TO,Kulldorff M,Forand SP,Haley VB,Evaluation of spatial filters to create smoothed maps of health dataStat MedYear: 2000192399240810.1002/10970258(20000915/30)19:17/18<2399::AIDSIM577>3.0.CO;2R10960861  
Turnbull BW,Iwano EJ,Burnett WS,Howe HL,Clark LC,Monitoring for clusters of disease: application to leukemia incidence in upstate New YorkAm J EpidemiolYear: 1990132S1361432356825  
Waller LA,Gotway CA,Applied spatial statistics for public health dataYear: 2004Hoboken, N.J.: John Wiley & Sons full_text.  
Brunsdon C,Estimating probability surfaces for geographical point data: An adaptive kernel algorithmComputers & GeosciencesYear: 199521877894  
Dobson JE,Bright EA,Coleman PR,Durfee RC,Worley BA,LandScan: A global population database for estimating populations at riskPhotogrammetric Engineering and Remote SensingYear: 200066849857  
LandScan™Global Population DatabaseOak Ridge, TN; Oak Ridge National Laboratoryhttp://www.ornl.gov/landscan/  
Shi X,Evaluating the uncertainty caused by Post Office Box addresses in environmental health studies: A restricted Monte Carlo approachInternational Journal of Geographical Information ScienceYear: 20072132534010.1080/13658810600924211  
Scott DW,Multivariate density estimation: theory, practice, and visualizationYear: 1992New York: Wiley  
World Stable LitesData Processing by Noaa's National Geophusical Data Center, Dmsp Data Collected by the US air Force Weather Agencyhttp://www.ngdc.noaa.gov/dmsp/download_Night_time_lights_9495.html 
Figures
Article Categories:

Previous Document: Development of stable reporter system cloning luxCDABE genes into chromosome of Salmonella enterica ...
Next Document: Socioeconomic inequalities in the use of outpatient services in Brazil according to health care need...