Spatio-temporal analysis of COVID-19 incidence rate using GIS: a case study-Tehran metropolitan, Iran.

COVID-19 has been distinguished as a zoonotic coronavirus, like SARS coronavirus and MERS coronavirus. Tehran metropolis, as the capital of Iran, has a high density of residents that experienced a high incidence and mortality rates which daily increase the number of death and cases. In this study, the IDW (Inverse Distance Weight), Hotspots, and GWR (Geography Weighted Regression) Model are used as methods for analyzing big data COVID-19 in Tehran. The results showed that the majority of patients and deaths were men, but the death rate was higher in women than in men; also was observed a direct relationship between the area of the houses, and the infected rate, to COVID-19. Also, the results showed a disproportionate distribution of patients in Tehran, although in the eastern regions the number of infected people is higher than in other districts; the eastern areas have a high population density as well as residential land use, and there is a high relationship between population density in residential districts and administrative-commercial and the number of COVID-19 cases in all regions. The outputs of local R2 were interesting among patients and underlying disorders; the local R2 between hypertension and neurological diseases was 0.91 and 0.79, respectively, which was higher than other disorders. The highest rates of local R2 for diabetes and heart disease were 0.67 and 0.55, respectively. From this study, it can be concluded the restrictions must be considered especially, in areas densely populated for all people.

Introduction

COVID-19 has been distinguished as a zoonotic coronavirus, like SARS coronavirus and MERS coronavirus (Dong et al., 2020; Gao et al., 2020a; Zheng et al., 2020). Tehran metropolis, the capital of Iran, has a high density of residents, which makes up approximately 20% of the country's population, disproportionately distributed throughout the city. This disproportion causes high incidence and mortality rates where a daily increase in the number of death and cases can be observed.

Disproportionate distribution of different land uses, including residential, commercial, administrative, recreational, green, and sports areas have led to many problems, such as traffic, air pollution, noise pollution, etc. Because of the economic imbalance caused by this disproportionate land use, it can be said for this case the residents of the northern areas have a far better economic rank as well as a better climate than the southern residents of the city (Zarandi et al., 2021). While the nature of this Virus is unique and strange, it is a very high spread infectious illness. Because of the geographic spread of the Virus, city officials, epidemiologists, and experts need to react and plan a strategy to the infection outbreak fast and suitably (Esri, 2011). Today, Pandemic monitoring is based and centralized on the potentials of geographic information systems (GIS) (Esri, 2020). This software gives various beneficial opportunities for improving the quality of the controls over the pandemic. It is done through the GIS abilities on improving monitoring and identifying paths, and it can predict the location of the infected people (Boulos & Geraghty, 2020; Esri, 2011; Koch & Koch, 2005). The World Health Organization (WHO) has consistently applied spatial analysis to manage and control disease prevalence, during past years (Esri, 2020). The monitoring and controlling of SARS, Ebola, and Zika and other vector-borne diseases are thriving experiences by application of GIS systems (Esri, 2020). The knowledge around the location of the patient is very notable in the process of decision-making (Esri, 2011). At the beginning of the COVID-19 outbreak, some countries began researching infected people and the rate of mortality (WHO, 2020). Using GIS systems to do this, South Korea became one of the best countries in the world to detect patients with the COVID-19 virus in the early stages of this outbreak (Mollalo et al., 2020). Like other countries that have used GIS and have succeeded in controlling this disease, we are using this tool to give a better spatial view to planners in the field of control.

This study aims to show spatial modeling on COVID-19 incidence and mortality rate, COVID-19 mobility pattern and aims to identify the relationship between COVID-19 and comorbidity, population density, and land use.

Literature review

Some studies have been conducted on various aspects related to the COVID-19 outbreak in Iran. For example, a study investigated effective climatology parameters on the COVID-19 Outbreak in Iran For this study, the main parameters included the number of infected people with COVID-19, population density, intra-provincial movement, infection days at the end of the study period, average temperature, average precipitation, humidity, wind speed, and average solar radiation was investigated. The Partial Correlation Coefficient (PCC) and Sobol’-Jansen methods have been used for analyzing the effect and correlation of variables with the COVID-19 spreading rate. The results of the sensitivity analysis showed that the population density and intra-provincial movement have a direct relationship with the infection outbreak. Conversely, areas with low values of wind speed, humidity, and solar radiation exposure had a high rate of infection that supported the virus's survival. The provinces such as Tehran, Alborz, Mazandaran, Qom, and Gilan are more susceptible to SARS-COV-2 infection because of high population density, intra-provincial movements, and high humidity rate in comparison with Southern provinces (Mohsen Ahmadi et al., 2020).

Another study investigated epidemiological characteristics of coronavirus disease 2019 (COVID-19) patients at a hospital in Iran. In their retrospective study, data related to the epidemiological characteristics of COVID-19 patients admitted to Baqiyatallah Hospital in Tehran, Iran, from February 19, 2020, to April 15, 2020, have been analyzed. During this period, 12,870 patients were referred to the hospital emergency department, of which 2968 were hospitalized with a COVID-19 diagnosis. The majority of cases were in the age group of 50 to 60 years old. The male-to-female ratio was 1.93:1. A total of 239 deaths happened among all cases for an overall CFR of 1.85%. Out of all patients, 10.89% had comorbidity. Their results showed that diabetes, chronic respiratory diseases, hypertension, cardiovascular diseases, chronic kidney diseases, and cancer were the most common comorbidities with 3.81, 2.02, 1.99, 1.25, 0.60, and 0.57%, respectively. Male gender (OR = 1.45, 95% CI: 1.08–1.96), older age (OR = 1.05, 95% CI: 1.04–1.06), and having underlying diseases (OR = 1.53, 95% CI: 1.04–2.24) were significantly associated with mortality (Nikpouraghdam et al., 2020). A study was conducted in Shiraz city that studied 113 COVID-19 positive cases admitted to hospitals affiliated to Shiraz University, Iran, from February 20 to March 20. Their study results showed that the mean age was 53.75 years, and 71 (62.8%) were males. The most typical symptoms at onset were fever (67: 59.3%), cough (73: 64.6%), and fatigue (75: 66.4%). Laboratory data showed a significant correlation between lymphocyte count (P-value = 0.003), partial thromboplastin time (P-value = 0.000), international normalized ratio (P-value = 0.000) with the severity of the disease. The most common abnormality in chest CT scans was ground-glass opacity (77: 93.9%), followed by consolidation (48: 58.5%) (Shahriarirad et al., 2020).

There are a limited number of GIS-based spatial modeling studies since the initial outbreak of COVID-19. Some studies indicated various applications of GIS for the precise description of the COVID-19 spread, spread probability, recognition health care priority places about vulnerable communities, including palliative, elderly, and comorbidity patients, mobility pattern, draw solutions (Arab-Mazar et al., 2020; Boulos & Geraghty, 2020; Gao et al., 2020b; Gibson & Rush, 2020; Lakhani, 2020; Mollalo et al., 2020; Rezaei et al., 2020; Sarwar et al., 2020).

Material and methods

Study area

It extends from 900 to 1800 m above sea level, with a decrease from north to south. This height decrease caused different climates to form in different areas of the Tehran metropolis (Alibakhshi et al., 2020). In 2019, the Tehran population was 8.964 million with a median age of 30.8 years in 22 districts, with it almost entirely distributed in the central, eastern, and northern regions. Moving from the north to the south, the quality of life and economic situations decrease (Zarandi et al., 2021).

Data required for the research process

Population data and land-use information were obtained from the Statistics Center of Iran, satellite images, and Tehran municipality, respectively. The data of patient's homes area obtained with phone calls to patients and their families. The phone calls were used only to complete the incomplete information in the file registered in the medical centers. In this study, we were contacted with all cases registered in medical centers, that 80% cooperated in completing the information.

Tehran has 118 hospitals, of which 49 are government hospitals, 55 are private hospitals, and 14 are hospitals belonging to the Armed Forces, called Armed Forces Hospitals.

The study is only related to hospitalized cases in the Tehran hospitals, and the information on non-hospital cases is unclear. The main source of information is the deputy healthcare system, while patient's incomplete information, such as test results or outcomes, has been completed utilizing other systems and registration systems (patients have adapted to each other using the national code and dossier registration number). The data is also related to January to mid-September. To obtain dispersion maps related to the COVID-19 virus in the Tehran metropolis, first, the patient's location of residence was registered as a unique location code in a spreadsheet next to the address column with the addresses named as the same location code register on Google Earth. Next, the patient's location of residence on Google Earth was called to ARC GIS software and the linguistic data were assigned the appropriate weight. For example, concerning patients' discharge status, were given to death cases, weight 2, hospitalization cases weight 1, and discharge cases weight 0. also those who had the chronic disease were assigned weight 1 and those who did have not were given weight 0.

Methods

The methods were used in GIS, including IDW, Hotspots as well as normalized distribution. Hotspots and Getis-Ord-Gi statistics were used to obtain hot and cold population points in the neighboring areas of each district from different areas of each region in the Tehran metropolis. Inverse Distance Weight (IDW) was used to distribute patients' condition based on deaths, hospitalizations and discharges, distribution of patients in terms of sex, distribution of patients in families, and distribution of patients in the family, as well as distribution based on area of residence location.

In this study, the GWR model was used to assess the relationship between patients infected to COVID-19 and peoples of infected to COVID-19 with chronic diseases (people with comorbidity).

IDW and hotspots model

The IDW model assumes a weight for each measurement point based on the distance between that point and the position of the point with unknown value. These weights are then controlled by the weighting power. Larger powers reduce the effect of points farther than the estimated point, and smaller powers distribute the weights more evenly between adjacent points. However, this method only takes into account their distance, regardless of the position and arrangement of the points. The points have the same distance from the estimation, the point has the same weight. The weight factor value is calculated using the following equation (Gong et al., 2014):where λi is the weight of station I, Di is the station distance to the unknown point, and α is the weighting power.

The Getis-Ord-Gi statistic was also used to distribute hot and cold spots spatially when a complication is known as a hotspot with both itself and those adjacent to it that have large amounts. When the Getis-Ord-Gi GI statistic is 1, 2, 3, the confidence level is 99%, 95%, 90%, respectively. The Getis-Ord-Gi statistic is calculated as follows (Gemmer et al., 2004):

Geographically weighted regression (GWR) model

The geographically weighted regression model (GWR) based on local regression and variable parameter analysis was proposed by Fotheringham and other scholars (Fig. 1).

Fig. 1

Location map of the study area (Zarandi et al., 2021)

Geographically weighted regression (GWR) technique is a model building tool, which provides relatively more precise results if the relationship between variables shows spatial nonstationary characteristics (Fotheringham et al., 2002; Wang & Wang, 2021). Since the GWR technique calculates the relationship between variables with a different coefficient for each geographical unit, it is possible to map where the relations are weak and strong, significant, and insignificant. It is more explanatory than global regression since it can explain via maps the relationship between two variables for each observation in geography with a different LOG ODD (odds ratio, OR) coefficient (Torun et al., 2020).

The expression of the geographically weighted regression model is as follows:

In this equation, y is the dependent variable, xi is the independent variable, β0 and β1 are the estimated coefficients, є is the error component, ui and vi are the latitude and longitude of the point i, and βk (ui, vi) is the implementation of the factor examined on a continuous level (Alibakhshi et al., 2020; Chu, 2012; Fotheringham et al., 2015; Mahmoud Ahmadi et al., 2018a, 2018b).

Results

Distribution of infected patients’ condition to COVID-19

As shown in Fig. 2, COVID-19 patients are found in almost all areas of Tehran. Although the number of patients is somewhat higher in the eastern regions, the deaths do not follow a specific spatial pattern, but most cases are in the western or western-oriented areas. Density in some densely populated areas, such as regions 4, 8, and 13, as well as the presence of some disease clusters in neighborhoods such as Tehranpars, Ayat Streets, and Piroozi, are significant. On the contrary, the low frequency of patients and death in the western regions of Tehran, such as regions 9, 18, 21, and 22, is one of the findings of the study. However, in these areas, there are also clusters of patients in the eastern regions of the 18th region and the eastern regions of the 21st region.

Fig. 2

Distribution of infected patients condition to COVID-19

As shown in Fig. 3a, people with the COVID-19 virus are almost equal in number, but their distribution is disproportionate. As seen here in Fig. 2, IDW-based interpolation, according to gender, shows that in the western regions with fewer population distribution are industrial zones with mostly male, and that's why most patients in these areas are men. In general, about 60% of those patients and 69% of those who died were male.

Fig. 3

Distribution of infected patients condition to COVID-19 in terms of hospitalization discharge and homes area

Figure 4a shows the distribution of the patient's occasion with underlying diseases (hypertension, heart disease, respiratory and asthma, diabetes, cancer, liver, kidney, nervous), and also shown in these figures, a large number of patients had these underlying diseases, among which a large number had hypertension. The tendency of chronic diseases is concentrated in the eastern regions of Tehran, with a high rate of exposure, which is one of the effective factors in increasing hospitalization and death in these areas. Hypertension (47.3), cardiovascular disease (16.3), and diabetes (15.2) were the most common in deaths from COVID-19 in Tehran. Furthermore, respiratory, neurological, and renal diseases appear to be present in border areas 13 and 14.

Fig. 4

Land-use of Tehran metropolis

Figure 5a shows the patient's condition in terms of discharge, hospitalization, and death that normalized by various underlying chronic diseases (hypertension, heart disease, diabetes, neurological disease, and cancer). From the following figures, it can be concluded that the COVID-19 patients had various underlying diseases of hypertension, heart disease, diabetes, neurological diseases, and cancer, respectively; which have the same consequence in terms of death. The geographical distribution of patients and deaths from hypertension, heart disease, and diabetes in different districts has been almost proportional, with the disproportionate central regions also infected with COVID-19 suffering from neurological and cancer diseases.

Fig. 5

Distribution of infected patients condition to COVID-19 in terms of hospitalization discharge and density of different land-use in urban areas

Distribution of infected patients’ condition to COVID-19 in terms of hospitalization discharge and Population

Fig. 6 The total population of the 22 regions is currently about 9 million people with an average age of about 33 years, and a median age of 32 years. More than a third of this population lives in 4-person households (the dimension of household size in Tehran is 3). The total population density of Tehran is about 6700 people per square kilometer, and the population in areas 4 and 5 is higher than in other areas. Figure 6 includes simple population distribution, population density, and population hotspots (hot and cold spots according to their surroundings, which have been studied by Getis-Ord-Gi). As a result, it can be concluded that most patients live in densely populated areas, which has increased the probability that they will come into exposure with another infected person

Fig. 6

Distribution of infected patients condition to COVID-19 in terms of hospitalization discharge and Population

Distribution of infected patients’ condition to COVID-19 in terms of hospitalization discharge and homes area

The 2016 census of Tehran shows that about 41% of Tehran's households are tenants, about 90% of whom live in 50 m’ homes. About 80% of Tehran's population resides in apartments. The average per capita housing in Tehran is about 18 square meters. As shown in Fig. 3, most people infected with COVID-19 reside in a low home area, which can lead to more contact between family members due to lack of space and, consequently, an increased risk of transmitting the virus to each other. The figure shows that we have the lowest area of houses in the eastern and central areas, which could lead to an increase in the number of patients for these areas in the coming days.

Distribution of infected patients’ condition to COVID-19 in terms of hospitalization discharge and density of different land-use in urban areas

Figure 5 shows the distribution of patients in terms of discharge status and density of different land uses in urban areas. In these figures, areas with a high density of residential lands, commercial, and administrative lands, have a high number of patients, but there is a less number of patients in the area with a high density of industries. This indicates the number of infected people is low in industrial areas, and contrariwise, more people live or move in densely populated areas of residential, commercial and administrative area, that increases probability of exposure with the COVID-19.

Discussion

The United Nations has expressed COVID-19 as a human, social, and economic crisis (Esri, 2020; Mahmoud Ahmadi et al., 2018b). The socioeconomic impacts and disease burden are especially evident in developing countries (Esri, 2020; Mahmoud Ahmadi et al., 2018b). It is predicted that the annual global gross domestic product will decrease by 24%, meaning that it is projected to diminish by 2% each month (Chu, 2012; Esri, 2020). The predictions also estimate a 13% to 32% decline in global trade (Chu, 2012; Esri, 2020). Recent studies over the world and Iran have shown that various factors like humidity, temperature, age, sex, smoking, and air pollution may add to the cruelty and spread rate of COVID-19. Taghizadeh-Hesary et al. (2020) shows negative impacts of smoking on the health outcomes of COVID-19 patients. Wang et al. (2020) indicated that humidity and temperature could influence the transmission of COVID-19 when compared to other respiratory viruses, suggesting a decline in disease spread. Also, out of other studies done about prevalent comorbidity and COVID-19, the majority showed a positive relationship between patients with COVID-19 and other diseases, including hypertension, diabetes, cardiovascular diseases (Guan et al., 2020; Li et al., 2020; Taghizadeh-Hesary & Akbari, 2020; Wang et al., 2020; Xu et al., 2020; Zheng et al., 2020). In the meantime, studies were done that related to origin, transmission, and epidemiological characteristics of COVID-19 (Guo et al., 2020; Surveillances, 2020).

The results of patients' studies in Tehran showed that the majority of patients and deaths were men, but the death rate was higher in women than in men. 59% of definite patients and 63% of deceased patients are men. Of course, higher mortality rates in men than women can be justified. Because the percentage of cases in men is higher, the death rate is consequently higher in men. According to the results, it can be argued that the presence of more men in the outdoor environment is the cause of more exposure and an increase in the rate of infection and death (in Iran, the presence of men outside the home is more likely than women). It is not possible to say with certainty why women have a higher ratio of death to infection than men. Given that the results of the present study have been performed only on visitors to the hospital, it is not possible to draw a correct conclusion from this section.

It cannot be stated with certainty that the incidence of the disease is higher in men, or that the severity of the disease is more severe in men unless the higher rate of exposure to men for job reasons or men outside the home is the cause of the difference. On the other hand, for biological reasons, women may have a milder form and less frequent visits, and it requires comprehensive serological studies in the epidemiology of the disease. It seems that more studies are needed in women to determine if women are more likely to have a subclinical condition and mild symptoms of the disease and not be diagnosed because of the low severity of the disease.

The results show a disproportionate distribution of patients in Tehran, although in the eastern regions the number of infected people is higher than in other districts. The eastern areas have a high population density as well as residential land use. These results also indicate a high relationship between population density in residential districts and administrative-commercial in all regions. The eastern areas with a high density of residential land-use, central and northern regions with a high density of administrative-commercial land-use, respectively, have the highest number of COVD-19 cases. As a result, they have more deaths, but in the western high-density industrial districts, the number of patients is lower than in other regions (Figs. 5,  7 and Table 1). There is a direct relationship between the area of the houses and the infected rate to COVID-19; most of the infected live in houses with an area of fewer than 100 m. The results of this study showed that most of the patients in their family did not have an infected person and often were family members of 3 people. Figure 7 shows Local R2 Between Population Density and infected patients to COVID-19. The range of local R2 between population density and infected patients to COVID-19 was 0.74 to 0.98, which is high in all areas and shows a direct relationship between population density and the number of patients.

Fig. 7

Local R2 Between Population and infected patients to COVID-19

Table 1

Some relevant data of Tehran (2017).

Source: Statistical Center of Iran

Regions of Tehran	Area of regions (M2)	Total of green area (M2)	(Total of green area/Area of regions) × 100	Total of Sport Area (M2)	(Total of sport area/Area of regions) × 100	Population (person)	Population of males (person)	Population of females (person)	Population density	Estimation of active population in 2017		Mortality (person)	Percent of mortality (person)
										Male working population	Female working population
1	46,612,194	16,361,578	35.101	349,500	0.749	493,889	241,805	252,084	0.010596	11,257	37,687	2682	0.543037
2	47,005,370	14,477,667	30.800	380,650	0.809	692,579	338,136	354,443	0.014734	145,261	62,735	2854	0.412083
3	29,216,831	5,198,517	17.792	268,490	0.918	330,004	158,054	171,950	0.011295	73,106	28,772	1489	0.451207
4	61,554,522	20,009,565	32.507	577,015	0.937	917,261	456,394	460,867	0.014902	217,865	51,440	3489	0.380372
5	53,161,331	12,536,590	23.582	455,907	0.857	856,565	420,431	436,134	0.016113	190,558	69,991	3580	0.417948
6	21,367,428	3,111,889	14.563	238,460	1.115	250,753	122,804	127,949	0.011735	52,697	25,323	1645	0.656024
7	15,335,212	1383,122	9.019	169,100	1.102	312,002	151,882	160,120	0.020345	73,574	25,820	1318	0.422433
8	13,156,442	1704,307	12.954	219,530	1.668	425,044	208,284	216,760	0.032307	100,577	29,983	1226	0.288441
9	19,746,502	2832,648	14.345	136,162	0.689	174,115	88,092	86,023	0.008818	42,052	8725	736	0.422709
10	8,185,467	824,570	10.0735	166,655	2.035	326,885	162,035	164,850	0.039935	81,596	21,105	1229	0.375973
11	12,031,459	1,536,301	12.769	556,206	4.622	308,176	154,516	153,660	0.025614	77,733	19,971	1041	0.337794
12	16,007,106	1,442,505	9.011	251,224	1.569	240,909	122,121	118,788	0.01505	64,188	9061	2935	1.218302
13	12,862,735	3,045,542	23.677	290,030	2.254	253,054	125,617	127,437	0.019673	57,493	12,519	718	0.283734
14	14,552,708	3,935,550	27.043	260,550	1.790	489,101	244,500	244,601	0.033609	125,727	21,399	3174	0.648946
15	27,740,625	10,675,530	38.483	769,212	2.772	659,468	335,314	324,154	0.023773	177,048	18,833	2045	0.310098
16	16,515,696	3,117,697	18.877	538,750	3.262	267,678	134,250	133,428	0.016207	66,040	8600	1545	0.577186
17	8,251,825	914,657	11.084	237,282	2.875	278,354	140,131	138,223	0.033732	71,534	8770	1012	0.363566
18	37,869,105	6,329,703	16.714	561,325	1.482	419,249	213,518	205,731	0.011071	109,585	12,011	1352	0.322481
19	20,341,427	6,854,017	33.694	986,574	4.850	255,533	130,203	125,330	0.012562	68,650	5703	740	0.289591
20	23,580,300	6,635,274	28.139	437,291	1.854	367,600	184,224	183,376	0.015589	91,744	14,488	2531	0.68852
21	51,525,318	7,000,510	13.586	156,750	0.304	186,319	93,739	92,580	0.003616	42,805	10,152	456	0.244742
22	59,003,253	10,321,647	17.493	345,660	0.585	175,398	89,146	86,252	0.002973	42,231	9589	521	0.297039
Sum of total	615,622,854	140,249,386	22.781	8,352,323	1.356	8,679,936	4,315,196	4,364,740	0.014099	2,084,221	512,679	38,318	0.441455

According to Fig. 6 and Fig. 10, it seems that the southern and southwestern regions had a lower number of cases as well as lower infected patients to COVID-19. Conversely, the more we go to the northern regions, the infected patients to COVID-19 are increasing, which can be directly related to the economic situation of these areas. Non-referral to medical centers due to unfavorable economic conditions has been the reason for the lower number of patients in the southern and southwestern regions compared to the northern areas. Inversely, the number of patients in the northern regions was high due to better economic ability and more visits to medical centers. In Iran, medical expenses are very high compared to income. In general, the average annual net cost and income for an urban household was 479.946 and 515.35 million Rials in Tehran (2017), respectively. Also, the price of the dollar varied between 30,000 and 50,000 Rials during 2017, and the same price for the dollar in 2020 reached over 300,000 Rials in Iran. It has increased almost ten times compared to the previous three years, which in turn has reduced people's economic ability. However, people's incomes haven't risen over the past three years, which affected families with weaker economies. It seems that if the economic situation of the people is better, the number of patients would be proportionally higher, which is due to the number of visits to medical centers. If we look at the data in Table 1, we find that the highest mortality rate was in the southern and central regions of Tehran, which have a weaker economy than the northern parts of the city. Also, these areas are facing environmental pollution, including air pollution, noise pollution, and traffic, more than other regions. While the northern areas of the city have a better climate, as well as more green space than other regions, it is worth mentioning that the people of the northern regions of Tehran have a much better economy than the southern areas (Table 2).

Fig. 10

Std Residual and Local R2 Between Total of Chronic Diseases and infected patients to COVID-19

Table 2

Infected patients condition to COVID-19 by regions.

Source: Ministry of Health and Medical Education of Iran

Regions of Tehran	The percentage of COVID-19 patients	The mortality rate of patients	The hospitalization rate of patients	The discharge rate of patients
1	6.4	17.04	21.59	61.36
2	8	16.36	28.18	55.45
3	5.67	24.35	25.64	50
4	9.47	15.38	26.92	57.69
5	8.16	4.46	31.25	64.28
6	2.41	18.18	15.15	66.67
7	4.51	11.29	22.58	66.13
8	5.74	10.12	32.91	56.96
9	1.53	9.52	14.28	76.19
10	2.77	10.52	23.68	65.79
11	4.22	10.34	27.58	62.07
12	4.95	16.17	17.64	66.17
13	4.08	16.07	26.78	57.14
14	7.27	16	19	65
15	5.31	16.43	28.76	54.79
16	4	16.36	18.18	65.45
17	2.19	23.33	13.33	63.33
18	2.47	11.76	23.52	64.70
19	3.43	10.63	14.89	74.46
20	3.64	14	18	68
21	1.53	9.52	71.42	19.04
22	2.25	9.67	29.03	61.29
Sum of total	100

The results of this study showed that a significant number of patients with chronic disease had a higher share of hypertension, heart disease, and diabetes than other underlying disorders.

The results of local R2 were interesting among patients and underlying disorders. Local R2 between hypertension and neurological diseases was 0.91 and 0.79, respectively (Figs. 7a,  10a), which was higher than other disorders. The highest rates of local R2 for diabetes and heart disease were 0.67 and 0.55, respectively (Figs. 8a,  9a). However, as shown in Fig. 10, local R2 was higher than 0.91 for all of the patients with comorbidity in all areas. These results should be evaluated clinically to determine the exact cause of underlying patients' vulnerability to COVID-19 relative to healthy individuals.

Fig. 8

Infected patients condition to COID-19 by regions

Fig. 9

Local Moran's I for Total infected patients to COVID-19

We used Std Residual and Local Moran's I statistics to validate the results of local R2. The results of the std residual show that the GWR model is generally a suitable model for all the studied regions. Based on the obtained results, all areas matched with the GWR model. Likewise, most of these regions were within the Std range of the − 2.5 to + 2.5.

The spatial autocorrelation was also observed when the analyses were conducted according to distribution of patient’s number (Moran’s I index 0.261858, z-score 2.050243, P = 0.040341). When the analyses were conducted according to number of infected patients to COVID-19, the correlation was found in the 18 region, with one cluster of high values (high-high). No significant spatial patterns were observed in further areas.

Conclusion

It can be concluded that there is a strong relationship between patients and their location of residence in terms of population and home area, type of land-use around the living home of patients, and underlying diseases. In general, the study results show that in areas with high population density and districts with commercial and residential land-use, the number of patients is high, and there is a direct relationship between them. This study also found that patients with comorbidity have a higher risk of death and infection to COVID-19 than healthy individuals. The people with chronic diseases, especially those with hypertension and neurological disorders, were more vulnerable to COVID-19 than healthy people, which is noteworthy for patients with hypertension, and neurological disorders; they need to be given more attention than other people. It emphasizes the importance of maintaining social distance and limiting commute in the city, especially in densely populated districts and commercial-administrative regions. It should be noted that it shouldn't limit the traffic passage to only sick people. Rather, the restrictions should be considered for all people. With a large number of Tehran citizens starting their working days in Tehran, there are probably many concerns about the rise in the number of people carrying COVID-19 and their exposure to different surfaces of the city, resulting in a very high increase in the coming days ahead.