Weight-Aware Sensor Deployment in Wireless Sensor Networks for Smart Cities

During the construction of wireless sensor networks (WSNs) for smart cities, a preliminary survey of the relative criticalness within the monitored area can be performed. It is a challenge for deterministic sensor deployment to balance the tradeoff of sensing reliability and cost. In this paper, based on the sensing accuracy of the sensor, we establish a reliability model of the sensing area which is divided into sensing grids, and different weights are allocated to those grids. We employ a practical evaluation criterion using seesaw mapping for determining the weights of sensing grids. We further formulate and solve an optimization problem for maximizing the trust degree of the WSNs. With our proposed method, the efficient deployment of sensors can be realized. Simulation results show that our proposed deployment strategy can achieve higher trust degree with reduced sensor deployment cost and lower number of sensors at a certain miss probability threshold.


Introduction
There are many wireless sensor networks (WSNs) applications in smart cities, such as residential security, factory pollution monitoring, elevator monitoring, lamp control, power grid monitoring, fire monitoring, medical care, building monitoring, and air monitoring.During the construction of wireless sensor networks for smart cities, we can survey the terrain situations in advance.In the monitored region, data in some areas is more critical.For instance, in traffic accident monitoring, some sections, such as curves and ramps, are accident-prone areas.Since the accuracy of the sensors cannot be the same in each monitored area, the data collection capability of the critical areas should be considered in the deployment of sensors in these regions.
At present, there are many research works about sensor deployment.They focus on various aspects.Some studies are about routing of senses with the consideration of the energy problem [1].The authors have optimized the sensor placement for indoor localization [2].In the process of sensor deployment, the weight of the sensing region should be considered.In [3], Zou and Chakrabarty considered the most vulnerable regions in the sensor field when designing algorithms.In [4], the authors think it is necessary to locate the sensor nodes in the field of interest which refers to the preferred field of coverage.It is of practical significance that the sensor deployment should consider the criticalness of each subregion of the sensing area.It is not necessary to deploy sensors in every corner.We need to weigh the criticalness of each subregion in the monitored region in order to ensure the reliability of the data collected.
In our work, we present an optimization framework for sensor deployment.The proposed approach aims to improve the comprehensive performance of the sensor network and ensure the effective collection of data.
The contributions of this paper are listed as follows.
(1) In order to collect data effectively and reliably, considering the actual engineering requirements, this paper proposes a weight-aware sensor deployment method based on the weight of each sensing grid in the monitored region, which is the meshed region according to the sensor's sensing region.
(2) In this paper, the weight of the sensing grid is estimated from three aspects: the probability of abnormal events, the consequence of abnormal events, and the tolerance time of abnormal events.A seesaw mapping is proposed to calculate the weights by taking into account the combined effect of the influencing factors, the cost of damage resulting from an abnormal event, and the tolerance time of an abnormal event.
(3) In this paper, the sensing accuracy and coverage degree are synthesized and the sensing trust degree model of the whole WSNs is proposed by taking the weight of sensing grid into account.
(4) Considering the fact that sensing areas are often subject to obstacles and sensing accuracy is limited, the corresponding algorithms are proposed, which makes the sensing trust degree calculation model more widely useful.
The remainder of this paper is organized as follows.In Section 2, we review prior research on the WSNs coverage and deployment.The model for the optimized sensor deployment is proposed in Section 3, and the estimation for the weight is in Section 4. The computing for the sensing trust degree of whole WSNs is also given in detail in Section 5.In Section 6, we describe the weight-aware sensor deployment algorithms in detail.In Section 7, we present the simulation results to evaluate the network performance for the proposed algorithms.Finally, we summarize the main results and give concluding remarks in Section 8.

Related Work
When sensor nodes are deployed in the monitored area, it is critical to calculate the weight of the sensing region.We divide the sensing region into many grids.An irregular sensor field is modeled as a collection of grids.The sensing matrix is formed by the weights of the sensing grids.Much research has considered dividing the sensing region into many grids.Xiao and Zhang proposed a divide-and-conquerbased surveillance framework, in which a large surveillance area is divided into small areas by critical points and critical lines [5].In [6], the authors consider an optimization problem of how to partition the given field into multiple parcels and to deploy sensor nodes in some selected parcels such that the field information covered by the deployed sensor nodes meets the requirement.In [7], the performance analysis of different grid types (e.g., triangle, square, pentagon, hexagon, heptagon, and octagon) is discussed by computing their coverage area and efficient coverage area ratios.
In the smart cities, sensor deployment must take into account the nature of the terrain, for example, obstacles such as buildings and trees in the line of vision for IR sensors, uneven surfaces and elevations for hilly terrains, redundancy due to the likelihood of sensor failures, and the power needed to transmit information between deployed sensors and between a deployed sensor and the cluster head [8].
Sensor deployment methods consider obstacles.In [9], the authors develop a series of mechanisms to handle the obstacles in the 3D environment and propose a greedy heuristic.Chang et al. [10] presented the efficient obstacle-resistant robot deployment (ORRD) algorithm, which involves the design of a node localization policy, a serpentine movement policy, obstacle-handling rules, and boundary rules.
To ensure the efficient collection of data, the minimum reliability of wireless sensor networks needs to be considered.
Long et al. [11] present an efficient data gathering scheme that guarantees the quality of service and optimizes the following network performance metrics as well as the end-to-end reliability in WSNs.In [12], Deif and Gadallah proposed an ant colony optimization (ACO) algorithm to solve the minimum cost reliability constrained sensor node deployment problem.
In summary, there has been a lot of work on the sensor deployment of the WSNs.Some papers are to solve the coverage problem [13][14][15][16].However, most of the researches did not discuss the weight of the sensing field in detail.In fact, the weight of the sensing field has a great influence on the reliability and accuracy of data sampling.The relationship between factors which are involved in the estimation of weights is complicated.In order to solve this problem, we develop a novel seesaw mapping.The major issue of deploying sensors is equivalent to the problem of finding the location combination of the sensors.Considering the weights, we proposed a series of weight-aware sensor deployment algorithms in three cases which are the given total number of sensors, minimum sensing reliability, and some obstacles in the sensing region.

Deployment Model
The deployment of the sensor network should ensure that all the monitored regions are covered.The number of sensors and the cost of the network are also considered.We can estimate the covered area and the position of the sensors using the method based on the grids.The monitored region is divided into several grids in which the sensors are deployed [17].
The probability of abnormal events is very high in certain fields of the monitored region.In monitoring, it is necessary to obtain data from these areas more precisely.In order to obtain accurate monitoring data of an interesting area, several sensors are usually deployed.Sensor deployment requires accurate and computationally feasible sensor detection models.When the sensor is deployed in the critical area, the minimum sensing radius of the sensor   that does not obstruct the monitoring is considered as the effective sensing threshold of the sensor node of the whole network.The monitored area is divided into  ×  sensing grids.We assume that each sensing grid is a regular inscribed quadrilateral of a circle whose diameter is 2  .The index of the sensing grid is (, ), where 1 ≤  ≤ , 1 ≤  ≤ .The model is illustrated in Figure 1.The sensing data matrix of the monitored area is set as  ∈ R × .
In this work, we determine a sensing grid to deploy only one sensor.If an area needs more sensors, we can take measures of selecting smaller   or increase the sensor layout density.
Here, we adopt a mask operator Γ(⋅) to represent the procedure of deploying the sensors: where   represents the matrix of the sensing grids' weights and ⊗ represents the elementwise product of two matrices. is an  ×  matrix defined by the following form: 1 A sensor has been deployed in the grid which is the (, ) th position of the  0 otherwise, where (, ) denotes the index of the sensing grid in .The sensors localization set  is obtained by combining all the coordinates of the sensing grid that satisfy (, ) = 1: Now, we can use  to represent the total number of the elements of .It is obvious that  is the total number of sensors.
The index of elements of  is numbered by line priority.The index of the numbers is from 1 to  × .The essence of weight-aware sensor deterministic deployment optimization is a combinatorial optimization issue for the number of sensing grids in which the sensors are located.The combination scheme for the number of sensor positions with the best monitoring effect is selected from many combinations.The set of various combinations of sensors is where   stands for the th combination of  numbers which is from 1 to  × . is the index of various combinations.
There are many combinations of sensor localization.It becomes a critical issue for the weight-aware sensor deployment to choose the best one.What criterion is the best combination of the sensor localization?To solve this problem, we need to consider the weight of the sensing grid in the matrix  and ensure the effective and reliable data collection.We introduce the concept of the sensing trust degree of the whole WSN to evaluate deployment.It is related to the sensing reliability degree of the sensing grid, the sensing coverage rate of the whole WSN, and the weight of the sensing grid.In the following, we introduce these concepts.
Definition 1 (the sensing reliability degree of the sensing grid).This is defined as the probability that the data in the sensing grid can be effectively and successfully collected and transmitted to the sink node.The sensing reliability degree of the (, )th sensing grid is represented by (, ).
(, ) is the dual concept of data miss probability in sensing grid.The more the data points collected in the sensing region, the smaller the possibility of data loss and the greater the value of (, ).(, ) is the result of the combined action of all sensors that can sense the (, )th sensing grid.It is related to the total number of sensors that can cover the sensing grid and the sensors' sensing accuracies.
The more the sensor nodes that collect the same data are, the higher the reliability of the data acquisition (, ) is.When one piece of data of the grid is collected by multiple sensors, if one of the sensors fails, it can be acquired by other sensors.For the critical areas of the monitored region, the reliability of data acquisition must be guaranteed.
Definition 2 (the sensing coverage degree of the sensing grid).The more powerful the sensor, the greater the area that it extends.The coverage ability of sensors is reflected by their sensing accuracy.The sensing coverage degree of the sensing grid is related to the sensing accuracy of sensors covered in this grid.The effective coverage degree produced by the most powerful sensor in this grid is used as the sensing coverage degree.The most powerful sensor's sensing accuracy is maximum.The sensing coverage degree of the (, )th sensing grid is represented by (, ).(, ) × (, ) is defined as the effective coverage area of the (, )th sensing grid, where (, ) denotes the area of the (, )th grid.
The sensing coverage rate of the whole WSNs is denoted by .
where  denotes the area of the whole sensing region. denotes the effective coverage area of the whole WSN.
The area where sensors overlap increases as the density of sensor deployment increases.As a result, the sensing coverage of the whole network is becoming lower. reflects the density degree of sensor deployment.
Definition 3 (the sensing trust degree of the whole WSNs).It is the sum of all sensing reliabilities of sensing grids in a certain sensing coverage rate of the whole WSN.It is denoted by (⋅).The higher the value of (⋅) is, the higher the reliability of the data collected is, that is, the better the performance of the monitoring network is.(i, j) (p, q): The coordinate of the sensing grid (i, j): The sensor location Finally, in the case that the number of sensors is given, we can get the best sensor deployment scheme by solving the following optimization problem: In actual engineering, to ensure the effective collection of data from the monitored region, we also encounter the minimum requirement for grid sensing.It is discussed in the following.
In practice, the minimum sensor reliability  min of sensing grid and the data loss probability threshold  min in [8] are a pair of dual concepts.The larger the value of  min is, the smaller the value of  min is, and vice versa.In the case that the minimum sensor reliability is defined, the problem of weight-aware sensor deterministic deployment problem becomes

Estimation of the Weights
Weight estimation is a complex work which is investigated with many methods.In [18], Fujimoto et al. proposed a model reestimation method that involves the extraction of reliable characteristics using Gaussian pruning with weight normalization.In this paper, we hold the main factors that affect the weights.However, the relationship between the weight and the factors is not simply cumulative.
The weight of each sensing grid is affected by many factors.In this paper, we have found the three main factors: PL, the probability of an abnormal event occurrence; , the cost of damage resulting from an abnormal event, that is, the extent of the impact; and , the tolerance time of an abnormal event, that is, the period time over which the monitored object exhibits an abnormal event to cause great losses.These three quantities of monitoring objects can be estimated by fuzzy comprehensive evaluation method or by Delphi method.
Among them, the relationship between the cost of damage resulting from an abnormal event and the tolerance time of an abnormal event is very special.The two variables are independent and unrelated to each other, and they have the opposite direction of development.The problem of calculating their effects on computation of the weight is getting more and more attention.Since the two variables are independent and unrelated, their coordinates are perpendicular to each other.It is not appropriate to calculate the common effects of two independent variables in a summation method.In this paper, the characteristics of two independent and unrelated variables are analyzed, and a novel seesaw mapping is introduced on the basis of the seesaw model in [19,20].

The Seesaw Mapping.
In some application scenarios, two independent variables or entities affect another variable or entity together, but the direction of development of these two entities is opposite, and two variables have opposite orientations relative to the third variable or entity.It is necessary to determine the extent of the combined effect.
The condition of seesaw model in [19] is as follows: (1) two independent variables or entities are generally independent of each other; (2) two independent variables or entities affect another variable or entity together.For the two variables or entities that satisfy the condition of the seesaw model, a seesaw mapping can be applied.The seesaw model condition ensures that the two variables or entities can be mapped to Euclidean space.The values of the two coordinates in the Euclidean space reflect the values of the two variables at some certain state points and the values of two individual developments.It is necessary to map the common effects of two quantities in the Euclidean space to a single coordinate.We compare their values by the front and back relationships in a coordinate.This comparison is intuitive and clear.This series of processes is called the seesaw mapping.
We assume the variables  and  satisfied the condition of seesaw model.
The set of the seesaw pairs of  and  is set as  = {(  ,   ) |  ∈ }, where (  ,   ) is a seesaw pair.One   corresponds to one   ; namely, there is a one-to-one correspondence between   and   . is the distribution sequence obeyed by  and . ⊂ R. The distribution sequence can be either spatial or temporal or any other distribution.For example, we can use the seesaw mapping to study the changes of  and  over time.We can also use it to study the difference between  and  when they are at different points in space. is the index of (, ) in . is commonly {0, 1, 2, . ..} or [, ], where  can take 0 or −∞ and  can be +∞.If  is a countable set, we call  a discrete sequence.  denotes the value of variable  at point  in .  denotes the value of variable  at point  in .
Definition 4 (the state of the combined effect).This is the quantitative representation of the combined effect of  and  at the  point in .It is represented by   (  ,   ), and   (  ,   ) ∈ R.
Definition 5 (the state space).This is the set of   (  ,   ).One calls the state space of the combined effect of  and .It is denoted by   (, ).  (, ) ⊂ R.
The seesaw mapping is represented as a mapping: where ( * , * ) is mapping operator.( * , * ) :  ×  → R is a two-variable single valued function defined on  × .
The seesaw mapping is as follows.
where the function F( * ) is to normalize the variables.
If two variables are measured in different units, various normalization methods are used to map the data of each variable to a range of 0∼1.The dimensional expressions are transformed into dimensionless expressions, that is, scalar quantity.In many scenarios, the dimensions of the two variables are different, and it is difficult to compare the size directly, so the two variables need to be converted into scalar quantities.
Step 2 (perform coordinate transform).The developing vector of a variable is taken by a series of values in which the variable varies with .
Definition 6 (positive/negative vectors).For the variable that plays a positive role in the combined effect, one sets it as a positive vector.The other vector is a negative vector.
Definition 7 (equilibrium factor).When the ratio of the normalized variables reaches a constant, the effect of two variables on a third variable is always the same; that is, the two variables are in equilibrium.This constant is defined as the equilibrium factor.It is denoted by .Obviously, if Definition 8 (seesaw axis).It is a coordinate axis measuring the degree of the seesaw effect in two variables.The equilibrium state is the origin.
In Figure 2, a rectangular coordinate system is made up of the developing vectors of two variables.Draw a ray with an angle of tan −1  in the counterclockwise direction where the coordinates of the positive vector are located.We can take The diagram of the proof for saw transform.one point at the rays and make a vertical line.The pedal can be used as the origin of axis.The direction of this vertical line, which is near the positive vector, is the positive direction of the axis.The other is the negative direction.This vertical line is used as the seesaw coordinate of the two variables.The number of negative directions to the origin is negative, and the number from the origin to the positive is a positive number.

Theorem 9. The values of two variables A and B at a
point  in T are, respectively, normalized.Let   be the positive vector and   be a negative vector.The value of Proof.Assume that two variables  and , which are normalized, can be mapped into an Euclidean space.Let  norm be horizontal ordinate and  norm be vertical ordinate.It is shown in Figure 2 that the th state of  and    (  ,   ) is mapped as the point  in the Euclidean space.∠ = tan −1 .Take a spot   at the ray .Draw a vertical line     from the   pedal.The direction close to the coordinate  is chosen as the positive direction.Finally,     can be used as the seesaw axis.The projection of   (  ,   ) on the axis     is    .Then, the value At the maximum or minimum point, the seesaw mapping is reduced to a size comparison of the values of a variable.In this paper, it is assumed that the contribution of the common effect of  and  to the weight is not 0 in the process of weight calculation.

Computation of the Weight Based on the Seesaw Mapping.
The probability of the occurrence of an abnormal event and the effect of an abnormal event in each sensing grid are always different, so the importance of each sensing grid is different.Let the weight of the (, )th sensing grid be denoted by  , .The matrix of the sensing grids' weights is   .
The weight  , of each sensing grid is affected by many factors.In this paper, we list the three main factors: the probability of an abnormal event occurrence  , , the cost of damage resulting from an abnormal event, that is, the extent of the impact  , , and the tolerance time of an abnormal event, that is, the period time over which the monitored object exhibits an abnormal event to cause great losses  , .These three quantities of monitoring objects can be estimated by fuzzy comprehensive evaluation method or by Delphi method.The larger the value of  , of a sensing grid is, the smaller the value of  , is and the more important the sensing grid is. , plays a positive role in  , . , plays a negative role in  , . , , which is the premise of the combined effect of  , and  , , plays a decisive role in  , .Therefore,  , is directly proportional to  , and is proportional to the combined effect of  , and  , .Now, we need to study the extent of the combined effect of  , and  , on  , in each sensing grid.However, the units of measurement for  , and  , are different.The changes in  , and  , are independent of each other.These satisfy the definition of the seesaw model.We can use the seesaw mapping to compute the weights.
Step 1.Using the maximum range formula,  , is normalized into   , : Step 2. Using the maximum range formula,  , is normalized into   , : Step 3. Compute the weight  , of each sensing grid.The role of  , is to calculate the criticalness of a sensing grid relative to other grids, so we compute the weight using the normalized abnormal event tolerance time   , and the normalized extent of the impact of the abnormal events   , .If   , goes up in increments and   , goes toward decreasing direction, it is beneficial to  , .  , and   , are in a state of seesaw.The extent of the combined effect of   , and   , can be computed using the seesaw mapping.
The smaller the value of   , is, the greater the value of   , is and the more important the sensing grid is.  , can be as the positive vector.In this work,  is valued at 1; that is,  and  have the same criticalness.We regard that  and  have the same criticalness.
According to the seesaw mapping (Theorem 11), we know the extent of the combined effect of   , =  , ×  , . ( Step 4. Normalize  , for   , : The weights matrix   of the monitored region is composed of many weights   , of sensing grids.

The Sensing Trust Degree of WSN
As mentioned early, (, ) is related to the total number and sensing accuracy of sensors.This paper assumes that the sensor is arranged at the center of each sensing grid.
According to the sensing precision formula of single sensor node in [21], it is known that the sensing accuracy is attenuated with the increase of the distance.Assuming the position (, ) of the grid is arranged with a sensor, the sensing accuracy formula of position (, ) of the sensing grid is derived as where  (,)→(,) represents the distance between the (, )th and the (, )th sensing grid.The sensor is located in the (, )th sensing grid.The constants  and  are devicedependent parameters reflecting the physical features of a sensor.Generally,  ranges from 1 to 4 [17], and  is used as an adjustment parameter [16].Sensors collect data from the sensing grid.The precision of the data in the sensing grid becomes smaller as the distance between the sensing data grid and the sensor grid increases.When a sensing grid is equipped with a sensor, its sensing accuracy corresponding to the sensor is 1.Since  = 0, then  = 1.This denotes that the probability of sensor false positives is very small.For example, the closer you get to the camera, the clearer you see and the less likely you are to misjudge.A sensor is deployed in a sensing grid, and the sensing accuracy of the sensor grid for this sensor is 1.In this paper, we assume that sensors are deployed at the center of the sensing grid.
In the practical application of sensor networks, obstacles inevitably exist in the monitored region.These obstacles cause signal blocking and reduce sensor sensing efficiency.Here, we consider the obstruction to  (,)→(,) , assuming the known topography within the sensing region.If the obstacles appear in the line of sight from the grid point (, ) to the grid point (, ), namely, if the coordinates (, ) of the obstacle satisfy the equations of connecting (, ) and (, ) line, that is, ( − )/( − ) = ( − )/( − ), then  (,)→(,) = 0.
Local blocking can also be modeled by setting the sensing probability to very small nonzero values [21].Proof.In actual operations, when (, ) and (, ) are on the same line or the same column, there is − == 0 or − == 0. It is not allowed for the denominator to be zero.
Theorem 13.In the model, the distance between the (, )th sensing grid and the (, )th sensor grid is where (, ) are the coordinates of the sensor grid in which the sensor is located.
Proof.The sensor is located at the center of the sensing grid.The sensing grid is the inscribed quadrilateral of a circle with a radius of   .It is shown in Figure 4 that the sensor is placed at the point , and  (,)→(,) is the distance between the grid (, ) and the point .
∵  (,)→(,) = If a plurality of sensors are deployed in the sensing region, the sensing reliability degree of the (, )th sensing grid is Sensing reliability is related to whether data can be sampled by the sensor smoothly.The data in the many sensors' sensing accuracy overlapping field is read by several sensors.The reliability of the overlap field is higher because if one sensor fails, the data can also be sampled through other sensors.So, the sensing reliability is the accumulation of all sensors' sensing accuracies.Proof.
∵  ,→, = 0 ∴ 1 The sensing accuracy of sensors reflects the ability of sensors to collect data.The stronger the sensor, the greater the area that it covers.The coverage area of the whole WSN is the union of all sensors' coverage fields.In the overlap region, the coverage degree takes the maximum value of all sensor accuracies.The sensing coverage degree of the (, )th sensing grid is The total sensing coverage degree of the whole WSN is ∵ According to Definition 2, the sensing coverage rate of the whole WSN is It is shown in Figure 4 that The sensing accuracy of multiple sensors overlaps in some sensing grids, and the sensing reliability of the sensing grids increases, while the effective coverage area of the whole WSN is reduced (i.e., the sensing coverage rate becomes smaller).This is because the coverage degree takes the maximum of all sensors' sensing accuracy in the overlapping areas.
Finally, the sensing trust degree of the whole WSN is defined as We can use the algorithms in the next section to achieve the maximum of (Γ( ⊗   )).

Weight-Aware Sensor Deployment Algorithms
It is necessary to select the sensor location combination scheme whose sensing trust degree of the whole monitored region is the maximum, from the combination scheme of many sensor locations.If the same criterion is followed by each location scheme, then these combinations can be compared.This paper discusses algorithms to address the optimal sensor deployment in three cases: (1) the total number of sensors deployed in the monitored region is given; it is very hard to increase the total number of sensors; (2) the minimum sensing reliability of the sensing grid is set; (3) there are some obstacles in the sensing region.

Algorithm Notations.
To describe our proposed algorithms, the following notations and message types are used.
: the number of rows in the sensing matrix, : the number of columns in the sensing matrix, : the total number of sensors deployed in the monitored region, : adjustment parameter, : device-dependent parameter, : a list that stores the combinations of  numbers from 1 to  × , : the index of combination sequences in ,   : the th combination, : a list that stores a series of row coordinates and column coordinates of each sensor position corresponding to the th combination   ,  * : a list that stores a series of row coordinates and column coordinates of each sensor position corresponding to the combination whose sensing trust degree of the whole WSN is maximum, : a list that stores the sensing trust degree of the whole WSN gained from each sensor deployment, : weight matrix, : sensing accuracy degree of the whole WSN, : sensing coverage degree of the whole WSN, (, ): the sensing reliability degree of the (, )th sensing grid, (, ): the sensing coverage degree of the (, )th sensing grid, (, ): the coordinate of the sensing grid, (, ): the coordinate of the sensing grid deployed with sensor.

Algorithm Description
Case 1.The total number of sensors deployed in the monitored region is given (see Algorithm 1).
Step 1. Generate combinatorial  numbers from 1 to × and store these combinations into .
Step 2. The sensing accuracy of each sensor is added as the sensing reliability degree (, ) of this sensing grid.
Step 3. The maximum sensing accuracy of each sensor to this sensing grid is used as the sensing coverage degree (, ) of this sensing grid.
Step 4. The sensing coverage degree of each sensing grid is accumulated to obtain .
Step 5.The sensing reliability degree (, ) of each sensing grid is multiplied by the weight of the sensing grid and then accumulated to obtain .
Step 6. Use formula (30) () = /( × ) ×  to compute the maximum sensing trust degree of the whole WSN of the th layout.
Step 7. Go through all the combinations in .Put the sensing trust degree of each combination into .When the  value is the maximum, the corresponding sensor position combination  result is required.
Step 8. Using the function ( ), the sensor position combination  result is transformed into the specific coordinates of the sensor in the sensing data matrix.
Case 2. Minimum sensing reliability  min of the sensing grid is set.
In practical engineering, in order to ensure that the data of each grid is collected effectively and reliably, the minimum sensing reliability degree  min of each sensing grid is often set (see Algorithm 2).
Step 1.It is necessary to determine whether (, ) is less than  min each time after (, ) is computed.If (, ) is less than  min , the current selected combination is not appropriate.
There is a need to select the next combination and delete the currently inappropriate combination sequence   .
Step 2. If the total number of sensors is small, this will cause the sensing accuracy to be less than  min in the calculation process of each sensor position combination.After going through all the combinations, the sensing accuracy is still less than  min , such that  is empty.There is a need to increase the number of sensors.The new combinations are generated and stored in .
Case 3.There are some obstacles in the sensing region.
Before (, ) is computed, the coordinates of all the obstacles are judged whether they satisfy the line equations determined by two points at (, ) and (, ).If they satisfy the equation, then the value of (, ) is zero.Otherwise, the value of (, ) can be calculated by Algorithm 1.

Performance Validation
In the section, we first validate our proposed algorithm.The sensing trust degree of various sensor deployments caused by different matrix weights is compared.Then, our proposed algorithm is compared with the widely used sensor deployment algorithm MAX AVG COV [8].
7.1.Simulation Settings.The simulation experiment is completed on MATLAB.This paper assumes that the sensing grid is a regular square.It is assumed that each sensor has the same sensing accuracy; that is,  and  of each sensor are the same.Table 1 enumerates the common parameters of our simulation experiments.

Simulation Results.
In the first simulation, we study the weight's impact on the data collection effects of the WSN.The sensor deployment schemes with the same number of sensors have different sensing trust degrees of the whole WSN in the case of the different weight matrices.We choose the sensor deployment with the highest trust degree of the whole WSN as the best sensor deployment.The weight of the sensing grid in the sensing region is related to the loss cost, occurrence probability, and tolerance time of monitoring abnormal events.The weight of all the grids constitutes the weight matrix of WSN.In order to compare the different sensing trust degrees of the whole WSN caused by the input: , ,  output:  * Generate combinatorial  numbers from 1 to  × , and store these combination into .
for all   ∈  do  ← 0,  ← 0 for  = 1 to  do for  = 1 to  do for all (, ) ∈  do (, ) = (, ) + 1 End different weight matrices, we randomly set the estimated parameters in Tables 2 and 3 for the weight matrices  1 and  2 , respectively.
The number of sensors is set to 5. The number of abscissa grids is  = 3, and the number of ordinate grids is  = 3.We can see that the different sensing trust degrees are obtained by different weight matrices in Figure 5.When the number of sensors is 5, there are 126 sensor deployment methods, that is, 126 combinations.In Figure 5(a), the sensing trust degree of the WSN which is caused by the sensor deployment corresponding to the 118th combination is highest, up to 1.424.The sensing trust degree of the WSN which is caused by the sensor deployment corresponding to the 23rd combination is minimum, only 0.9707.In Figure 5(b), the sensing trust degree of the WSN which is caused by the sensor deployment corresponding to the 71st combination is highest, up to 1.504.The sensing trust degree of the WSN which is caused by the sensor deployment corresponding to the 71st combination is minimum, only 0.8902.The sensing trust degree of the best combination is about two times that of the worst one.Because of the weight difference, the effect of data collection on the same number of sensors is different.Finally, it is necessary to study the influence of grid weight on sensor deployment.
In order to further analyze the impact of grid weights on data collection, we compare the best deployed sensor locations with the worst locations in the case of the two different weight matrices.In Figure 6, these e denote the sensor locations with the highest sensing degree of the whole WSN, and these △ denote the sensor locations with minimal sensing degree.In Figure 6(a), the sensor deployment coordinates with the highest sensing degree of the whole As is shown in Figure 6, we can see that the sensors are deployed on the edge of the sensing area and in the grids with lower weights by the sensor deployment with the minimum trust degree of the whole WSN.However, the sensors are almost deployed on the center of the sensing area and in the grids with larger weights by the sensor deployment with the maximum trust degree.
In the second simulation, we compare our proposed algorithm with MAX AVG COV [8].We assume that the number of rows of the sensing matrix is 4 and the number of columns is 4. We measure the minimum number of sensors required, when miss probability threshold  min (corresponding to the minimum reliability of each grid  min ) varies from 0.1 to 0.7 with an increment of 0.1.
As we can see from Figure 7, the minimum number of sensors required decreases as the minimum reliability increases.When  min is 0.1, the minimum number of sensors required by the proposed algorithm is half the number in the MAX AVG COV algorithm.When  min = 0.7, the required number of sensors in our algorithm is less than 2 in the MAX AVG COV algorithm.Obviously, the algorithm proposed in this paper is superior to MAX AVG COV algorithm.
We compare the minimum number of sensors required by the two algorithms, as shown in Figure 10, when there are obstacles in the sensing region.
As we can see from Figure 10, the minimum number of sensors required varies when there are obstacles in the sensing area.The variations of the algorithms proposed in this paper are smaller.At the miss probability threshold 0.1, the minimum number of sensors required by the MAX AVG COV algorithm increases to 15, while that of our proposed algorithm is still 6, with a difference of 9.With the increase of miss probability threshold, the restriction of sensor deployment is getting smaller and smaller, and the minimum number of sensors required is also reduced.At the miss probability threshold 0.7, the minimum number of sensors required by our algorithm is 3, while the minimum number of sensors required by the MAX AVG COV algorithm is 5, with a difference of two.Since MAX AVG COV algorithm is a local optimum selection method, the localization of obstacles Obstacle X (i,j) : The sensing grid  has a great influence on the sensor deployment.The proposed algorithm is a global optimal selection method; as long as the computation time is sufficient, we can always choose an optimal layout scheme, even a way of sensor deployment that  can be chosen to avoid the effects of obstacles.From this, we can see that the algorithm proposed in this paper is obviously superior to MAX AVG COV algorithm.

Conclusion
In practical applications, especially for smart city wireless sensor network deployment, the inconsistency of data weights owing to relative criticalness in different sensing grids needs to be taken into account.In order to ensure the data collection of wireless sensor networks effectively, we propose a computing model of the sensing trust degree of the whole WSN.In this paper, the sensing region is divided into meshed grids.The weights of the sensing grids are estimated by seesaw mapping to reflect the combined effect of the variables related to abnormal events, including the probability of occurrence PC (, ), loss cost, that is, the degree of impact (, ), and the tolerance time (, ).The computing model based on sensing trust degree takes into account the sensor location, the impact of the weight of sensing grid, and the obstacles in the sensing region.We study the impact of different sensor deployments on the sensing trust degree of the whole WSN.Our model can obtain as much information as possible when the number of sensors is given.Based on our proposed computing model and formulated optimization problem, we propose the sensor deployment algorithms.Our algorithms are designed for three cases: (1) the number of sensors is limited, (2) the minimum sensing precision is defined, and (3) there are some obstacles in the sensing area.In the simulation experiments, we find out that the sensing affections of different sensor deployment schemes with different weight matrices are different under the condition of the same number of sensors.When the sensing accuracy is limited, the minimum number of sensors needed by our sensor deployment algorithm is less than that of the conventional MAX AVG COV algorithm.In this work, the sensing grid is square.In our following research, we will study the sensor deployment scheme for irregular grids and reduce the computational complexity.

Figure 1 :
Figure 1: The sensing grids of the monitored region.
The extended seesaw axis       in Figure2is shown in Figure3.The intersection of the seesaw axis and the horizontal ordinate  norm is   .The intersection of the seesaw axis and the horizontal ordinate  norm is   .At the point   , the value of axis reached the maximum, and then  norm  =  norm and  norm  = 0, such that  norm  ( norm  ,  norm  ) = 1.At the point   , the value of axis is minimum, and then  norm  =  norm and  norm  = 0, such that  norm  ( norm  ,  norm  ) = −1:

Figure 3 :
Figure 3: The diagram of the proof for saw transform property.

9 . 11 .
and there is  norm  ( norm  ,  norm  ) = 0.It requires the further evolution of Theorem Theorem Positive values of the common contribution    (   ,    ) of two variables  and  to the third variable are Proof.∵ The domain of the exponential function () =   is , the range is () > 0, and in the domain () increases with the increase of .∵  norm  ( norm  ,  norm  ) reflects the comparative relationship between the common effects of two variables at different points in .Since the exponential function () =   is an increasing function, it does not affect the comparative relationship between them.At the same time, it can ensure that  norm  ( norm  ,  norm  ) is more than 0. ∴ The domain of the exponential function () =   can be used as a quantized function of the common effect of two variables  and .

Figure 7 :
Figure 7: The comparison for minimum reliability.

Figure 8 :
Figure 8: The layout model of the obstacle of our proposed algorithm.

Figure 9 :
Figure 9: The layout model of the obstacle of MAX AVG COV algorithm.

Figure 10 :
Figure 10: The comparison for the existence of obstacles.
Algorithm 1: The weight-aware sensor deployment algorithm with the given number of sensors.

Table 2 :
The estimated parameters of  1 .

Table 3 :
The estimated parameters of  2 .