Dynamic Inversion of Inland Aquaculture Water Quality Based on UAVs-WSN Spectral Analysis

: The inland aquaculture environment is an artiﬁcial ecosystem, where the water quality is a key factor which is closely related to the economic beneﬁts of inland aquaculture and the quality of aquatic products. Compared with marine aquaculture, inland aquaculture is normally smaller and susceptible to pollution, with poor self-puriﬁcation capacity. Considering its low cost and large-scale monitoring ability, many researches have developed spectrum sensor on-board satellite platforms to allow remote monitoring of inland water surface. However, there remain many problems, such as low image resolution, poor ﬂexible data acquisition, and anti-interference. Apart from that, the conventional forecasting model is of weak generalization ability and low accuracy. In our study, we combine unmanned aerial vehicles system (UAVs) with the wireless sensor network (WSN) to design a new ground water quality parameter and drone spectrum information acquisition approach, and to propose a novel dynamic network surgery-deep neural networks (DNS-DNNs) model based on multi-source feature fusion to forecast the distribution of dissolved oxygen (DO) and turbidity (TUB) in inland aquaculture areas. The result of using fused features, including characteristic spectrum, Gray-level co-occurrence matrix (GLCM) texture feature, and convolutional neural network (CNN) texture feature to build a model is that the characteristic spectrum + CNN texture fusion features were the best input items for DNS-DNNs when forecasting DO, with the determination coe ﬃ cient R 2 of the vertical set arriving at 0.8741, while the characteristic spectrum + GLCM texture + CNN texture fusion features were the best for TUB, with the R 2 reaching 0.8531. Compared with a variety of conventional models, our model had a better performance in the inversion of DO and TUB, and there was a strong correlation between predicted and real values: R 2 reached 0.8042 and 0.8346, whereas the root mean square error (RMSE) were only 0.1907 and 0.1794, separately. Our study provides a new insight about using remote sensing to rapidly monitor water quality in inland aquaculture regions.


Introduction
Inland aquaculture refers to the use of freshwater on the land surface to engage in freshwater fishery production activities such as fishing and breeding of aquatic products with economic value. It is one of the important components of the entire aquaculture industry [1]. The good development of aquatic products has strict requirements on the quality of water. In inland aquaculture, sustainable by Ai Yeshuang et al. [14], an adaptive sampling consistency limit learning machine algorithm was applied to invert and analyze the concentration of suspended solids in a water body, and thereby stable effects were obtained. The researchers developed the remote sensing monitoring of water quality from simple qualitative analysis into quantitative inversion, from an empirical model with spatial and temporal limitations to a widely applicable biological optical model, which continuously broadened the application prospects of remote sensing water quality monitoring [15,16]. In addition, with the emergence of various remote sensing data sources, especially hyperspectral and multispectral remote sensing data, and the deep understanding of the spectra of various water quality parameters, the parameters of remote sensing water quality monitoring continue to increase, and the inversion accuracy also constantly improves [17][18][19]. However, the past remote sensing monitoring was mainly based on satellite remote sensing, and its data source only had a spatial resolution of 10-30 m and a spectral resolution of 10-20 nm, so it isn't suitable for the mapping of inland small and medium inland water, especially the typical aquaculture area in Guangdong Province. In addition, it is necessary to eliminate the influence of atmospheric scattering, so it cannot be widely used in the monitoring of various water bodies. In fact, there hasn't been a perfect satellite remote sensing data source for inland water bodies up to now [20].
In recent years, the application field of unmanned aerial vehicles system (UAVs) has been continuously expanded. It has become a new trend to observe water quality characteristics through ground experiments and spectral radiometers by using UAVs equipped with spectral cameras, which has been favored by researchers [21][22][23], as they are rapid, efficient, and flexible spectrum information acquisition systems. Peter et al. [24] used UAVs remote sensing to obtain data, analyzed the impact of land leveling in catchment areas on linear soil erosion in ditches, and carried out an effective assessment of water and soil loss in water areas. Su and Chou [25] used the multispectral sensors installed on UAVs to map the nutritional status of small reservoirs, and the designed remote sensing system had a relatively high spatial resolution (0.1 m). After the ground measurement, UAVs spectral imaging was implemented as soon as possible to effectively shorten the time interval between remote sensing imaging and ground measurement. Taking the east lake on the campus of Zhejiang A & F University as the research object, Liu et al. [26] obtained spectral reflectance data by means of multispectral sensor mounted on UAVs, constructed the inversion models of total phosphorus (TP), suspended substance concentration (SS) and turbidity (TUB), and studied the spatial distribution of their concentrations, which provided a new technical means for the monitoring of water quality parameters. Compared with common satellite remote sensing, UAV remote sensing can more clearly obtain the spectrum information of water quality, because its sensors are closer to the surface of water, it can minimize the impact of the atmosphere in the radiation transmission and has higher spatial resolution and spectral data accuracy, so it is more suitable for dynamically monitoring the water quality factors of factory intensive aquaculture in inland water environment. However, UAVs spectral remote sensing combined with ground WSN monitoring data and the inversion water quality distribution approaches have rarely been reported.
Water quality remote sensing inversion algorithm can be roughly divided into empirical algorithm and model-based algorithm. The empirical algorithm mostly adopts a ratio method, such as hyperbolic exponential algorithm [27], binary quadratic polynomial algorithm [28], SeaWiFS algorithm [29] as well as the semi-empirical Carder model [30]. However, with the development of the bio-optical model, the quantitative inversion algorithm based on model has become the mainstream of retrieving the correlation between the water body components and the spectral radiation characteristics of water body. Quantitative inversion mostly adopts linear regression model. For example, Olmanson et al. [31] used unitary linear equations to invert the water quality evaluation index of the Mississippi River. However, it isn't enough for the linear model to describe the complex relationship between remote sensing observation indexes and water quality parameters, which may easily lead to information loss and distortion, and its practicability is limited. Therefore, nonlinear models such as artificial neural network (ANN), extreme learning machine (ELM) and so on are gradually adopted. Zhu Yunfang et al. [32] established the ANN model by using GF-1 WFV4 image and real-time ground sampling data, which indicated that the accuracy of establishing inversion model by the ANN was usually higher than that of the previous ratio method. Sun Deyong et al. [33] made an attempt on remote sensing inversion of the concentration of suspended substances in Taihu Lake by using ANN. Although ANN can approach arbitrarily complex nonlinear relations, there are still some problems, such as slow learning convergence speed, easiness to fall into local extremum, and difficulty to determine the network structure. As a new single-hidden layer forward neural network, ELM has the advantages that ANN describes nonlinear relations and its own network structure is easy to determine, and at the same time, it has fast learning speed, strong generalization ability [34]. However, due to numerous and complex uncertain factors in the ground synchronous measurement of water quality elements (such as measurement methods, weather conditions, and the laboratory technicians' proficiency in work standards), large data size, much noise, as well as complex relationship between mixed pixels, so that the conventional linear regression, ANN and ELM algorithm cannot adapt to these features, resulting in a low inversion precision. Therefore, a new information extraction method and inversion model are needed.
In this paper, combining with the ground WSN and UAVs spectral remote sensing technology, an acquisition scheme of water quality spectral elements suitable for the complex water bodies of aquaculture was designed. Based on the general packet radio service (GPRS) protocol, the ground WSN network built a real water quality parameter acquisition system on the ground, and the parameters include dissolved oxygen (DO) and turbidity (TUB) values, which are sensitive to aquaculture organisms, large changes in them may directly lead to the death of aquaculture organisms. The UAVs is equipped with a spectral imager to obtain the spectral data of water quality elements. Based on the measured data of UAVs spectral data and ground sensor monitoring, a dynamic network surgery-deep neural networks (DNS-DNNs) water quality spectral inversion model based on multi-source feature fusion was proposed to predict the distribution of dissolved oxygen and turbidity in aquaculture area.
The main innovations and contributions of our research work are as follows: (1) The UAV was used as the remote sensing platform to make the data acquired by multispectral sensor easier to process, which could avoid the difficult atmospheric correction process due to low reflectance response. UAVs and ground WSN were combined into a complete acquisition system, to realize automatic data acquisition and improve acquisition efficiency. (2) The texture features of characteristic spectrum images were extracted by GLCM and CNN methods, respectively, and different fusion features were formed by combining with the characteristic spectrum to explore the optimal fusion features suitable for different water quality parameters. (3) A DNS-DNNs water quality spectrum inversion model was proposed. Due to the high computational complexity of the DNNs model, the dynamic network surgery method was adopted to induce the regular term by updating the weight matrix, so as to make it more sparse, thus effectively compressing the network and reducing the complexity. (4) In order to visually display the spatial variation law of water quality parameters, the DNS-DNNs optimal estimation model based on multi-source feature fusion was used to calculate and generate the prediction distribution diagram of DO and TUB, and the three-dimensional visualization of water quality parameters was realized.

Materials and Methods
The method that we are proposing can be structured according to the following steps: UAVs overflight with multispectral cameral and processing of spectral data; Ground WSN water quality data collection and processing; Characteristic spectrum acquisition based on correlation analysis; training and vertification of the DNS_DNNs; and inversion of DO and TUB. The flowchart of the proposed method is depicted in Figure 1 and detailed in the following items.

Study Area
The experimental area is an aquaculture area that is located in Zengcheng Teaching and Research Base of South China Agricultural University, Guangzhou, China (113°38'14.5" E, 23°14 '36.3" N), as shown in the red area in Figure 2c. According to actual surveys, the area of aquaculture is 10,000 m 2 and the average depth between 2 and 2.5 m. About 6000 Mandarin fish are cultured in this area. Mandarin fish are sensitive to DO and TUB of water, The DO concentration value should be greater than 3.0 ppm (1 ppm = 1mg / L) and the TUB concentration value should be lower than 25.0 ppm in water [35]. This area meets the composition criteria for a small inland aquaculture area.

Wireless Sensor Network
The ground WSN acquisition node (AN) consists of DO sensor, TUB sensor, solar panel and system control box internally equipped with GPRS communication module, as shown in Figure 3, and it can monitor the water quality of aquafarm in real time. The system control box has a built-in STM32F407 controller, GPRS communication module and RS485 interface accessories, which were

Study Area
The experimental area is an aquaculture area that is located in Zengcheng Teaching and Research Base of South China Agricultural University, Guangzhou, China (113 • 38 14.5 E, 23 • 14 36.3 N), as shown in the red area in Figure 2c. According to actual surveys, the area of aquaculture is 10,000 m 2 and the average depth between 2 and 2.5 m. About 6000 Mandarin fish are cultured in this area. Mandarin fish are sensitive to DO and TUB of water, The DO concentration value should be greater than 3.0 ppm (1 ppm = 1mg / L) and the TUB concentration value should be lower than 25.0 ppm in water [35]. This area meets the composition criteria for a small inland aquaculture area.

Study Area
The experimental area is an aquaculture area that is located in Zengcheng Teaching and Research Base of South China Agricultural University, Guangzhou, China (113°38'14.5" E, 23°14 '36.3" N), as shown in the red area in Figure 2c. According to actual surveys, the area of aquaculture is 10,000 m 2 and the average depth between 2 and 2.5 m. About 6000 Mandarin fish are cultured in this area. Mandarin fish are sensitive to DO and TUB of water, The DO concentration value should be greater than 3.0 ppm (1 ppm = 1mg / L) and the TUB concentration value should be lower than 25.0 ppm in water [35]. This area meets the composition criteria for a small inland aquaculture area.

Wireless Sensor Network
The ground WSN acquisition node (AN) consists of DO sensor, TUB sensor, solar panel and system control box internally equipped with GPRS communication module, as shown in Figure 3, and it can monitor the water quality of aquafarm in real time. The system control box has a built-in STM32F407 controller, GPRS communication module and RS485 interface accessories, which were

Wireless Sensor Network
The ground WSN acquisition node (AN) consists of DO sensor, TUB sensor, solar panel and system control box internally equipped with GPRS communication module, as shown in Figure 3, and it can monitor the water quality of aquafarm in real time. The system control box has a built-in STM32F407 controller, GPRS communication module and RS485 interface accessories, which were burnt with the sensor driver and communication algorithm program. In order to continue field operation, the acquisition controller used solar panels made of monocrystalline silicon as the power supply. The acquisition interval is 3 min/time. If the battery is fully charged, its time of endurance is more than two days. The collected data are uploaded to the cloud for analysis and processing, and sent to the mobile terminal software for display. The measured value of water quality monitoring can be exported from the mobile terminal. burnt with the sensor driver and communication algorithm program. In order to continue field operation, the acquisition controller used solar panels made of monocrystalline silicon as the power supply. The acquisition interval is 3 min/time. If the battery is fully charged, its time of endurance is more than two days. The collected data are uploaded to the cloud for analysis and processing, and sent to the mobile terminal software for display. The measured value of water quality monitoring can be exported from the mobile terminal.

Unmanned Aerial Remote Sensing System
The spectral reflectance data were obtained by means of the multispectral sensor mounted on the UAVs, and the inversion models of DO and TUB were constructed to study the spatial distribution of their concentrations. The UAVs spectrum platform adopted is composed of DJI M600 pro UAVs, manufactured by DJI Innovations, K4 multi-spectrometer and ground control station, as shown in Figure 4. The wavelength range of K4 is 395~1000 nm, including 21 filters, and the spatial resolution is 0.155 m.

Data Acquisition and Processing
Under the condition of clear sky, the calm water surface is captured by the spectral sensor due to the flare generated by clouds and solar scintillation, which results in the changes in the spectral information of water quality and seriously affects the modeling accuracy [36]. The flare intensity decreases with the increase of the azimuth difference between the sun and the sensor, and rapidly decreases with the increase of the solar altitude. After surveying for many times, when the sun's

Unmanned Aerial Remote Sensing System
The spectral reflectance data were obtained by means of the multispectral sensor mounted on the UAVs, and the inversion models of DO and TUB were constructed to study the spatial distribution of their concentrations. The UAVs spectrum platform adopted is composed of DJI M600 pro UAVs, manufactured by DJI Innovations, K4 multi-spectrometer and ground control station, as shown in Figure 4. The wavelength range of K4 is 395~1000 nm, including 21 filters, and the spatial resolution is 0.155 m.
Remote Sens. 2019, 11, x FOR PEER REVIEW 6 of 18 burnt with the sensor driver and communication algorithm program. In order to continue field operation, the acquisition controller used solar panels made of monocrystalline silicon as the power supply. The acquisition interval is 3 min/time. If the battery is fully charged, its time of endurance is more than two days. The collected data are uploaded to the cloud for analysis and processing, and sent to the mobile terminal software for display. The measured value of water quality monitoring can be exported from the mobile terminal.

Unmanned Aerial Remote Sensing System
The spectral reflectance data were obtained by means of the multispectral sensor mounted on the UAVs, and the inversion models of DO and TUB were constructed to study the spatial distribution of their concentrations. The UAVs spectrum platform adopted is composed of DJI M600 pro UAVs, manufactured by DJI Innovations, K4 multi-spectrometer and ground control station, as shown in Figure 4. The wavelength range of K4 is 395~1000 nm, including 21 filters, and the spatial resolution is 0.155 m.

Data Acquisition and Processing
Under the condition of clear sky, the calm water surface is captured by the spectral sensor due to the flare generated by clouds and solar scintillation, which results in the changes in the spectral information of water quality and seriously affects the modeling accuracy [36]. The flare intensity decreases with the increase of the azimuth difference between the sun and the sensor, and rapidly decreases with the increase of the solar altitude. After surveying for many times, when the sun's

Data Acquisition and Processing
Under the condition of clear sky, the calm water surface is captured by the spectral sensor due to the flare generated by clouds and solar scintillation, which results in the changes in the spectral information of water quality and seriously affects the modeling accuracy [36]. The flare intensity decreases with the increase of the azimuth difference between the sun and the sensor, and rapidly decreases with the increase of the solar altitude. After surveying for many times, when the sun's altitude angle increases to 30 • and the azimuth angle increases to 70 • , the solar flare almost disappears. Therefore, in order to eliminate the influence of solar flares on the extraction of water quality information, we selected four time points (09:30 on 3 May 2019, 15:30 on 3 May 2015, 09:30 on 5 June 2019, and 15:30 on 5 June 2019) for sample measurement. To reduce the accumulated error, the calibrated reference was measured prior to each sample measurement. We acquired 20 spectral images with each filters, respectively, and used the empirical line method (ELM) [37,38] to calibrate them to surface reflectance. As a result, there are a total of 420 spectral images after radiometric calibration.
There were an aeration zone and a feeding zone in the study area. In order to reflect the diversity and representativeness of the sampling points, a total of 11 sampling points were set, including 3 offshore sampling points (SP9-SP11) and eight shore side sampling points (SP1-SP8). Combining with the ground measured data obtained by the developed ground WSN equipment, the test site and ground data sampling points are shown in Figure 5. The statistical distribution of the real values of DO and TUB during the four sampling times were shown in Figure 6, where 1 ppm = 1 mg/L. To reduce the accumulated error, the calibrated reference was measured prior to each sample measurement. We acquired 20 spectral images with each filters, respectively, and used the empirical line method (ELM) [37,38] to calibrate them to surface reflectance. As a result, there are a total of 420 spectral images after radiometric calibration.
There were an aeration zone and a feeding zone in the study area. In order to reflect the diversity and representativeness of the sampling points, a total of 11 sampling points were set, including 3 offshore sampling points (SP9-SP11) and eight shore side sampling points (SP1-SP8). Combining with the ground measured data obtained by the developed ground WSN equipment, the test site and ground data sampling points are shown in Figure 5. The statistical distribution of the real values of DO and TUB during the four sampling times were shown in Figure 6, where 1 ppm = 1 mg/L.  In the same four time periods, a UAVs with multi-spectrometers was used to obtain multispectral images 200 m above the study area. The UAVs was controlled to be in a hovering state, and the shooting angle of the spectrometer was vertically downward. To reduce the accumulated error, the calibrated reference was measured prior to each sample measurement. We acquired 20 spectral images with each filters, respectively, and used the empirical line method (ELM) [37,38] to calibrate them to surface reflectance. As a result, there are a total of 420 spectral images after radiometric calibration.
There were an aeration zone and a feeding zone in the study area. In order to reflect the diversity and representativeness of the sampling points, a total of 11 sampling points were set, including 3 offshore sampling points (SP9-SP11) and eight shore side sampling points (SP1-SP8). Combining with the ground measured data obtained by the developed ground WSN equipment, the test site and ground data sampling points are shown in Figure 5. The statistical distribution of the real values of DO and TUB during the four sampling times were shown in Figure 6, where 1 ppm = 1 mg/L.  In the same four time periods, a UAVs with multi-spectrometers was used to obtain multispectral images 200 m above the study area. The UAVs was controlled to be in a hovering state, and the shooting angle of the spectrometer was vertically downward. In the same four time periods, a UAVs with multi-spectrometers was used to obtain multispectral images 200 m above the study area. The UAVs was controlled to be in a hovering state, and the shooting angle of the spectrometer was vertically downward.
In order to obtain sensitive bands, 420 spectral images were selected visually. We totally selected one spectral image in each of the four time periods for every band, so there are 84 spectral images, and that cover all sampling point. A large amount of data taken was selected visually. The image directly above the sampling point was selected and imported into ENVI5.3 software. Then, 11 water surface sampling points were found according to according to the latitude and longitude coordinates, and a 16 × 16 (ppi) matrix centering on the sampling point was constructed as the region of interest (ROI). Later, taking the average spectral reflectance of all points in the region as the spectral reflectance data of this point, a total of 924 sets of spectral reflectance data of water quality corresponding to ground WSN monitoring data were obtained.

Image Denoise and Feature Spectrum Acquisition
The spectral data acquisition process is often influenced by instrument noise and other factors, and there are nonlinear absorbance, baseline variation and additional scattering variation [39], which inevitably produces errors. Therefore, in order to reduce the interference of noise with the model, Savitzky-Golay filtering method based on time-domain local polynomial least square fitting was adopted for preprocessing [40].
To establish an inversion model with high accuracy, the Pearson correlation coefficient was carried out to the measured DO and TUB and the normalized (N) and first-order differential (FoD) water spectral reflectance to obtain the sensitive bands of DO and TUB. The correlation coefficient is shown in Figure 7. In order to obtain sensitive bands, 420 spectral images were selected visually. We totally selected one spectral image in each of the four time periods for every band, so there are 84 spectral images, and that cover all sampling point. A large amount of data taken was selected visually. The image directly above the sampling point was selected and imported into ENVI5.3 software. Then, 11 water surface sampling points were found according to according to the latitude and longitude coordinates, and a 16 × 16 (ppi) matrix centering on the sampling point was constructed as the region of interest (ROI). Later, taking the average spectral reflectance of all points in the region as the spectral reflectance data of this point, a total of 924 sets of spectral reflectance data of water quality corresponding to ground WSN monitoring data were obtained.

Image Denoise and Feature Spectrum Acquisition
The spectral data acquisition process is often influenced by instrument noise and other factors, and there are nonlinear absorbance, baseline variation and additional scattering variation [39], which inevitably produces errors. Therefore, in order to reduce the interference of noise with the model, Savitzky-Golay filtering method based on time-domain local polynomial least square fitting was adopted for preprocessing [40].
To establish an inversion model with high accuracy, the Pearson correlation coefficient was carried out to the measured DO and TUB and the normalized (N) and first-order differential (FoD) water spectral reflectance to obtain the sensitive bands of DO and TUB. The correlation coefficient is shown in Figure 7. By comparison, it was found that there was a strong correlation(| | > 0.6)between DO and the first-order differential water spectral reflectance near the wave band of 550 nm, 590 nm, 615 nm and 632 nm, and there was a significant correlation between DO and the normalized spectral reflectance at 550 nm (| | = 0.6) and 850 nm (| | = 0.88), indicating that the variation of DO had a significant impact on the spectral reflectance of water body at 550 nm, 590 nm, 615 nm, 632 nm, and 850 nm. Meanwhile, there was a strong correlation between TUB and the first-order differential water spectral reflectance at 808 nm (| | = 0.80) and 850 nm (| | = 0.81), and there was a significant correlation between TUB and normalized spectral reflectance at 880 nm (| | = 0.66). Therefore, in this paper, the sensitive bands of DO were selected as 550 nm, 590 nm, 615 nm, 632 nm, 850 nm, and the sensitive bands of TUB were selected as 808 nm, 850 nm, and 880 nm.
In order to reduce the interference of background information and extract more effective spectral information, several popular band combination empirical indices, including average red green (ARG), red green ratio (RGR), near-infrared red ratio (NRR), red green(RG), sum red green (SRG), senescence reflectance (SR), anthocyanin reflectance (AR), vogelmann (VOG) and average nearinfrared (ANIR), were extracted as presented in Table 1. Those combination formulae were applied to the sensitive bands to obtain the characteristic spectral information. By comparison, it was found that there was a strong correlation(|r| > 0.6) between DO and the first-order differential water spectral reflectance near the wave band of 550 nm, 590 nm, 615 nm and 632 nm, and there was a significant correlation between DO and the normalized spectral reflectance at 550 nm (|r| = 0.6) and 850 nm (|r| = 0.88), indicating that the variation of DO had a significant impact on the spectral reflectance of water body at 550 nm, 590 nm, 615 nm, 632 nm, and 850 nm. Meanwhile, there was a strong correlation between TUB and the first-order differential water spectral reflectance at 808 nm (|r| = 0.80) and 850 nm (|r| = 0.81), and there was a significant correlation between TUB and normalized spectral reflectance at 880 nm (|r| = 0.66). Therefore, in this paper, the sensitive bands of DO were selected as 550 nm, 590 nm, 615 nm, 632 nm, 850 nm, and the sensitive bands of TUB were selected as 808 nm, 850 nm, and 880 nm.
In order to reduce the interference of background information and extract more effective spectral information, several popular band combination empirical indices, including average red green (ARG), red green ratio (RGR), near-infrared red ratio (NRR), red green(RG), sum red green (SRG), senescence reflectance (SR), anthocyanin reflectance (AR), vogelmann (VOG) and average near-infrared (ANIR), were extracted as presented in Table 1. Those combination formulae were applied to the sensitive bands to obtain the characteristic spectral information.

Image Texture Information Extraction
The decomposition of suspended organic matters containing chlorophyll a (Chl.a) in freshwater aquaculture water leads to a huge consumption of oxygen and produces a large amount of carbon dioxide, which eliminates the effect of photosynthesis of algae on DO, thus resulting in a negative correlation between the change of chlorophyll a and DO. However, the change of chlorophyll a content is reflected in the change of water texture characteristics, so the DO value change can be represented by the texture characteristics of chlorophyll a in water body [48]. Turbidity refers to the darkening of the light beam passing through the suspended substances in water, and there is a positive correlation with suspended substances, so the change of turbidity value can also be represented by the texture characteristics of suspended substances in water [49].
In this paper, the Gray level co-concurrence matrix (GLCM) [50] was used to extract the texture features of characteristic spectrum, and the convolutional neural network was applied to extract the CNN texture. Based on the GLCM texture feature, the extraction distance parameter was set as 1 pixel, and the direction was 0 • , 45 • , 90 • , and 135 • , respectively. 4 texture features, including energy, contrast, correlation and homogeneity, were extracted in the ROI (61px × 61px). Based on the pre-trained AlexNet [51] neural network, the texture features of spectral images under sensitive bands and band combinations were extracted.

DNS-DNNs Algorithm
DNNs is a deep learning model, which is one of the most popular research directions in the field of artificial intelligence in recent years. Its main principle is to automatically extract features in data, layer by layer, through deep neural network and feedback propagation. Set the network weight matrix as V, the bias vector as b, the total number of samples as m, the total number of network layers as n l−1 . Assuming the function is h, the network layer number is l (l = 1,2, . . . ,n l−1 ), the total number of neurons in the lth layer is s l , the weighted attenuation factor is λ, and the minimum loss function J is where V l ji is the weight matrix between the jth neuron in the layer l + 1 network and the ith neuron in the layer 1 network; P i is the ith input sample vector; y i is the ith output vector; h v,b is a hypothesis function with parameters V and b.
Due to the complex structure and large calculated amount of DNNs, the model's power consumption is high and regression speed is slow. In order to compress the network, dynamic network surgery [52] was used in the training process to carry out the induction of regular items for the updating of the weight matrix v to make it more sparse, so that most of the weights were 0. Dynamic network surgery includes two processes: pruning and splicing, in which pruning cuts the weights which are considered to be unimportant, while splicing verifies the importance of weights and recovers the important weights that have been cut. Updating v as the following scheme. That is, in which β indicates a positive learning rate. We update not only the important parameters, but also the ones corresponding to zero entries of b, which are considered unimportant and ineffective to decrease the network loss. The DNS-DNNs algorithm in this paper as Table 2.  Update v l by formula (2) and the current loss function gradient ∆L end for Update: iter ← iter + 1 and β ← f (α, iter) until iter reaches its desired maximum Two indexes, R 2 and RMSE, were used in the modeling evaluation. The mean value of the performance indexes was taken as the experimental result after running for 20 times.

Spectral Data Correlation Analysis and Texture Feature Extraction
Pearson correlation coefficient was carried out to each calculation formula and the DO content and TUB values actually measured at 11 sampling points, as shown in Table 3. Bands and band combination parameters with a strong correlation (r > 0.8) and significance level p < 0.05 were selected as the characteristic spectra for modeling analysis.
As shown in Table 3, the characteristic spectra conforming to the requirements of DO are S1, S3, S7, S10, S12 and S14, respectively, with 264 sets of data. The characteristic spectra conforming to the requirement of TUB are S1, S2, S3 and S5, respectively, with 174 sets of data.
Characteristic spectral images were extracted from the DO and TUB data samples conforming to the requirements, and 264 DO characteristic spectral images and 174 TUB characteristic spectral images were obtained. Four kinds of texture features were adopted for each image, and the collection directions of each feature were 0 • , 45 • , 90 • , and 135 • , respectively, and 16-dimensional GLCM texture was obtained. Therefore, the data volumes of GLCM texture features of DO and TUB are 4224 and 2784, respectively.
AlexNet extracted 160-dimensional CNN texture, and there were a total of 176-dimensional texture features. Therefore, the data volumes of CNN texture features of DO and TUB are 42,240 and 27,840 respectively. Further, 80% of characteristic parameters, including characteristic spectral, GLCM texture features, and CNN texture features, were randomly selected as the correction set for modeling, and the remaining 20% was used as the verification set for evaluating the model's performance.

Comparative Analysis of DNS-DNNs Modeling Results Based on Multi-Source Feature Fusion
The DNS-DNNs model was established for the two water quality parameters. After experiments, when the weight attenuation factor was adjusted to 0.9, the number of DNNs layers was set as 6, the number of hidden neurons was 20-10-10-20, the learning rate was α = 0.03, and the dropout intensity of the optimization objective function was set as 0.1, the best parameter performance was obtained. In the training process, the best iteration times were selected by using the early stop method, as shown in Figure 8.
As shown in Table 3, the characteristic spectra conforming to the requirements of DO are S1, S3, S7, S10, S12 and S14, respectively, with 264 sets of data. The characteristic spectra conforming to the requirement of TUB are S1, S2, S3 and S5, respectively, with 174 sets of data.
Characteristic spectral images were extracted from the DO and TUB data samples conforming to the requirements, and 264 DO characteristic spectral images and 174 TUB characteristic spectral images were obtained. Four kinds of texture features were adopted for each image, and the collection directions of each feature were 0°, 45°, 90°, and 135°, respectively, and 16-dimensional GLCM texture was obtained. Therefore, the data volumes of GLCM texture features of DO and TUB are 4224 and 2784, respectively. AlexNet extracted 160-dimensional CNN texture, and there were a total of 176dimensional texture features. Therefore, the data volumes of CNN texture features of DO and TUB are 42,240 and 27,840 respectively. Further, 80% of characteristic parameters, including characteristic spectral, GLCM texture features, and CNN texture features, were randomly selected as the correction set for modeling, and the remaining 20% was used as the verification set for evaluating the model's performance.

Comparative Analysis of DNS-DNNs Modeling Results Based on Multi-Source Feature Fusion
The DNS-DNNs model was established for the two water quality parameters. After experiments, when the weight attenuation factor was adjusted to 0.9, the number of DNNs layers was set as 6, the number of hidden neurons was 20-10-10-20, the learning rate wasα = 0.03, and the dropout intensity of the optimization objective function was set as 0.1, the best parameter performance was obtained. In the training process, the best iteration times were selected by using the early stop method, as shown in Figure 8. The characteristic spectrum was optimized by fusing with GLCM texture feature and CNN texture feature, the DNS-DNNs model based on characteristic spectrum +GLCM texture feature fusion, characteristic spectrum +CNN texture feature fusion, characteristic spectrum +GLCM texture feature +CNN texture feature fusion, and single characteristic spectrum was established, respectively, as shown in Table 4. The characteristic spectrum was optimized by fusing with GLCM texture feature and CNN texture feature, the DNS-DNNs model based on characteristic spectrum +GLCM texture feature fusion, characteristic spectrum +CNN texture feature fusion, characteristic spectrum +GLCM texture feature +CNN texture feature fusion, and single characteristic spectrum was established, respectively, as shown in Table 4.

Comparative Analysis of Water Quality Inversion Model Based on Feature Fusion
In order to compare with conventional inversion model performance, multivariable linear regression model (MLRM), ANN and ELM inversion models were established, and then the calibration set of the characteristic spectrum of DO + CNN texture feature, and the calibration set of the characteristic spectrum of TUB + GLCM texture + CNN texture fusion feature was imported into the above model, respectively. The determination coefficient R 2 and root mean square error (RMSE) of the inversion model corresponding to each parameter are shown in Table 5. The data of verification set was imported into each conventional model and DNS-DNNs model to predict the concentration values of DO and TUB, which were linearly fitted with the measured values. The test results are shown in Figure 9.

Comparative Analysis of Water Quality Inversion Model Based on Feature Fusion
In order to compare with conventional inversion model performance, multivariable linear regression model (MLRM), ANN and ELM inversion models were established, and then the calibration set of the characteristic spectrum of DO + CNN texture feature, and the calibration set of the characteristic spectrum of TUB + GLCM texture + CNN texture fusion feature was imported into the above model, respectively. The determination coefficient and root mean square error (RMSE) of the inversion model corresponding to each parameter are shown in Table 5.

DO and TUB Content Verification and Distribution Dynamic Inversion
Based on the characteristic spectrum image, obtained in Section 3.1, the background and shadows are removed by automatic binarization, and the entire water area is extracted as a new ROI. At this time, the size of the characteristic spectrum image is uniformly 610px × 610px. Based on the above conclusions, the DNS-DNNs model based on multi-source feature fusion was established to predict the distribution of DO and TUB in water body. The spectral curve corresponding to each pixel was all imported into the DNS-DNNs model, and the spectral information of each pixel and its related points could predict the concentration values of DO and TUB at that point, so as to further obtain the distribution of DO and TUB in the whole region.
We select 4 ROI which are from characteristic spectrum images at 09:30 and 15:30 on May 3, and use the trained DNS-DNNs model to forecast the DO and TUB concentration distributions in ROI, and present them in a three-dimensional space. The distribution results are shown in Figure 10.

DO and TUB Content Verification and Distribution Dynamic Inversion
Based on the characteristic spectrum image, obtained in Section 3.1, the background and shadows are removed by automatic binarization, and the entire water area is extracted as a new ROI. At this time, the size of the characteristic spectrum image is uniformly 610px × 610px. Based on the above conclusions, the DNS-DNNs model based on multi-source feature fusion was established to predict the distribution of DO and TUB in water body. The spectral curve corresponding to each pixel was all imported into the DNS-DNNs model, and the spectral information of each pixel and its related points could predict the concentration values of DO and TUB at that point, so as to further obtain the distribution of DO and TUB in the whole region.
We select 4 ROI which are from characteristic spectrum images at 09:30 and 15:30 on May 3, and use the trained DNS-DNNs model to forecast the DO and TUB concentration distributions in ROI, and present them in a three-dimensional space. The distribution results are shown in Figure 10.

DO and TUB Content Verification and Distribution Dynamic Inversion
Based on the characteristic spectrum image, obtained in Section 3.1, the background and shadows are removed by automatic binarization, and the entire water area is extracted as a new ROI. At this time, the size of the characteristic spectrum image is uniformly 610px × 610px. Based on the above conclusions, the DNS-DNNs model based on multi-source feature fusion was established to predict the distribution of DO and TUB in water body. The spectral curve corresponding to each pixel was all imported into the DNS-DNNs model, and the spectral information of each pixel and its related points could predict the concentration values of DO and TUB at that point, so as to further obtain the distribution of DO and TUB in the whole region.
We select 4 ROI which are from characteristic spectrum images at 09:30 and 15:30 on May 3, and use the trained DNS-DNNs model to forecast the DO and TUB concentration distributions in ROI, and present them in a three-dimensional space. The distribution results are shown in Figure 10.

DNS-DNNs Performance Analysis Based on Different Multi-Source Feature Fusion
As provided in Figure 8, when the number of iterations was selected as 3819 by using the early stop method, the model determination coefficient of the verification set reached the maximum, the root-mean-square error reached the minimum, and the model's comprehensive performance was the best.
The correction set and verification set of DNS-DNNs model based on characteristic spectrum for the two water quality parameters was less than 0.8 in Table 4. The model's performance didn't reach the ideal state, and the generalization ability was weak. The characteristic spectrum + CNN texture fusion feature had the best effect in the prediction of DO, and at this time, the correction set and verification set was 0.8741 and 0.8042, respectively, and RMSE was 0.1938 and 0.1907,

DNS-DNNs Performance Analysis Based on Different Multi-Source Feature Fusion
As provided in Figure 8, when the number of iterations was selected as 3819 by using the early stop method, the model determination coefficient of the verification set reached the maximum, the root-mean-square error reached the minimum, and the model's comprehensive performance was the best.
The correction set and verification set R 2 of DNS-DNNs model based on characteristic spectrum for the two water quality parameters was less than 0.8 in Table 4. The model's performance didn't reach the ideal state, and the generalization ability was weak. The characteristic spectrum + CNN texture fusion feature had the best effect in the prediction of DO, and at this time, the correction set and verification set R 2 was 0.8741 and 0.8042, respectively, and RMSE was 0.1938 and 0.1907, respectively. In the prediction of TUB, the effect of the characteristic spectrum + GLCM texture feature + CNN texture feature was the best. Compared with the modeling using characteristic spectrum alone, the stability of the model was improved and the generalization ability was enhanced to a certain extent under the fusion of various data including characteristic spectrum and image texture features.

Different Inversion Model Performance Analysis
After comparing Table 4 with Table 5, it could be seen that, among the training results of DO fusion features, ELM model had the worst performance, and its R 2 was only 0.6493, and the ANN model had the best performance, and its R 2 was only 0.7853, which was lower than that of DNS-DNNs (0.8741). In the comparative analysis on the training results of Do fusion features, it was found that, among the two water quality parameters: DO and TUB, for the MLRM model with the best performance, R 2 was 0.7727, which was also lower than that of DNS-DNNs model (0.8531).
It could be known that the conventional models all had different degrees of predicted deviation data points in Figure 9, mainly because the modeling involved in the model's bias-variance contradiction, the model's complexity wasn't enough, and the extracted effective features were not sufficient or they were not completely related to the original data, thus affecting the model's robustness and generalization ability. The predicted values of the two water quality parameters by the DNS-DNNs model were close to the trend line of the measured values, and R 2 was 0.8042 and 0.8346, respectively, and there were bias data points. Basically, it could be seen that the estimation level of the model was stable and close to the real value.

Effects of Dynamic Inversion Distribution
As shown in Figure 10a, at 09:00, the concentration of DO in the oxygenated zone was significantly higher than that in other regions, and gradually diffused to the surrounding area and decreased. In the feeding area, the average DO concentration was the lowest because there were more fish and oxygen consumption was fast. In Figure 10b, compared with the distribution of DO concentration at 09:00, the distribution at 15:00 was more uniform in the entire aquaculture area, while the overall average concentration was slightly lower than the former. The reason is that although the enhanced photosynthesis of the plants increased the natural DO on the water surface, the farmer reduced the power of the aerator which resulted into dramatic reduction in artificial DO. According to Figure 10c, the high turbidity value in this area was caused by the frequent fish movement and fast metabolism in the feeding area and oxygenated area, and the agitated water body drove the floating of impurities in water bottom. As the surrounding area of the study area was the shallow water area, and the spectral reflectance of bottom impurities was high, there was a higher average turbidity value. Compared with the TUB distribution in Figure 10c, that in Figure 10d does not change much, and the high TUB values are mainly concentrated in the aeration and feeding zones.

Conclusions
In the research, in view of the deficiencies existing in conventional water quality testing methods and water quality inversion models, and combined with internet of things and UAVs technology, a water quality parameter acquisition approach which is suitable for fresh water aquaculture was designed, and a prediction model of water quality parameter content distribution of deep learning based on multi-source feature fusion of spectral image and convolutional neural network was proposed. Both theoretical research and experimental verification proved that the application of UAVs remote sensing technology combined with ground WSN automatic monitoring technology was feasible in water quality factor inversion. UAVs remote sensing is a research carried out under the guidance of the inversion idea of water quality factors of satellite remote sensing, and it is proved that the data fusion conducted by the data extracted from multispectral sensors with different features can meet the inversion requirements of some water quality factors. The DNS-DNNs model constructed on the basis of deep learning has a better adaptability to multi-source feature fusion spectral information.
Compared with the conventional inversion algorithm, the model's stability was improved, and its generalization ability was also enhanced, so it can successfully predict the distribution of DO and TUB in inland aquaculture area.
At present, there are few applications in the field of inland aquaculture water monitoring using UAVs hyperspectral data. It is also important to choose the right algorithm for inland aquaculture water quality monitoring research. We hope that our research can provide technical support for the monitoring of water quality in inland water, especially in inland aquaculture. However, there is still some room for improvement. For example, on the one hand, we will explore new wireless transmission mode to transmit spectral data and real values of ground WSN to the server directly, and process water quality parameters in real time by combining with edge calculation. On the other hand, we will increase the number of spectral bands, and expand the amount of sensitive feature data for inversion accuracy improvement. Additionally, we also believe that extending our study to incorporate multi-UAVs collaborative operation to realize large-scale inversion of large areas of water, which will be interesting.