In China′s Loess Plateau area, gully head is the most active zone of a drainage system in gully areas. The differentiation of loess gully head follows geospatial patterns and reflects the process of the loess landfor...In China′s Loess Plateau area, gully head is the most active zone of a drainage system in gully areas. The differentiation of loess gully head follows geospatial patterns and reflects the process of the loess landform development and evolution of its drainage system to some extent. In this study, the geomorphic meaning, basic characteristics, morphological structure and the basic types of loess gully heads were systematically analysed. Then, the loess gully head′s conceptual model was established, and an extraction method based on Digital Elevation Model(DEM) for loess gully head features and elements was proposed. Through analysing the achieved statistics of loess gully head features, loess gully heads have apparently similar and different characteristics depending on the different loess landforms where they are found. The loess head characteristics reflect their growth period and evolution tendency to a certain degree, and they indirectly represent evolutionary mechanisms. In addition, the loess gully developmental stages and the evolutionary processes can be deduced by using loess gully head characteristics. This study is of great significance for development and improvement of the theoretical system for describing loess gully landforms.展开更多
The properties and feasibility of L-band differential InSAR for detecting and monitoring mining-induced subsidence were systematically analyzed and demonstrated. The largest monitored subsidence gradient of 7.9×1...The properties and feasibility of L-band differential InSAR for detecting and monitoring mining-induced subsidence were systematically analyzed and demonstrated. The largest monitored subsidence gradient of 7.9×10-3 and magnitude of 91 cm were firstly derived by theoretical derivation. Then, the stronger phase maintaining capacity and weaker sensitivity to minor land subsidence compared with C-band DInSAR were illustrated by phase simulation of the actual mine subsidence. Finally, the data processing procedure of two-pass DInSAR was further refined to accurately observe subsidence of a coalfield of Jining in Northern China using 7 ALOS PALSAR images. The largest monitored subsidence magnitude of 39.22 cm and other properties were better investigated by testing results interpretation and subsidence analysis, and the absolute difference varying from 0.5 mm to 17.9 mm was obtained by comparison with leveling-measured subsidence. All of results show that L-band DInSAR technique can investigate the location, amount, area and other detailed subsidence information with relatively higher accuracy.展开更多
Mesoscale eddies are an important component of oceanic features.How to automatically identify these mesoscale eddies from available data has become an important research topic.Through careful examination of existing m...Mesoscale eddies are an important component of oceanic features.How to automatically identify these mesoscale eddies from available data has become an important research topic.Through careful examination of existing methods,we propose an improved,SSH-based automatic identification method.Using the inclusion relation of enclosed SSH contours,the mesoscale eddy boundary and core(s) can be automatically identified.The time evolution of eddies can be examined by a threshold search algorithm and a tracking algorithm based on similarity.Sea-surface height(SSH) data from Naval Research Laboratory Layered Ocean Model(NLOM) and sea-level anomaly(SLA) data from altimeter are used in the many experiments,in which different automatic identification methods are compared.Our results indicate that the improved method is able to extract the mesoscale eddy boundary more precisely,retaining the multiple-core structure.In combination with the tracking algorithm,this method can capture complete mesoscale eddy processes.It can thus provide reliable information for further study of reconstructing eddy dynamics,merging,splitting,and evolution of a multi-core structure.展开更多
Associations between autumn Arctic sea ice concentration(SIC) and early winter precipitation in China are studied using singular value decomposition analysis. The results show that a reduced SIC almost everywhere in t...Associations between autumn Arctic sea ice concentration(SIC) and early winter precipitation in China are studied using singular value decomposition analysis. The results show that a reduced SIC almost everywhere in the Arctic Ocean, except the northern Greenland Sea and Canadian Basin, are accompanied by dry conditions over central China, extending northeast from the Tibetan Plateau toward the Japan Sea, the Bohai Sea and the Yellow Sea, and wet conditions over South China and North China. Atmospheric circulation anomalies associated with SIC variability show two wave-train structures, which are persistent from autumn to winter, leading to the identified relationship between autumn Arctic SIC and early winter precipitation in China. Given that the decline in autumn SIC in the Arctic Ocean is expected to continue as the climate warms, this relationship provides a possible long-term outlook for early winter precipitation in China.展开更多
Vertical errors often present in multibeam swath bathymetric data. They are mainly sourced by sound refraction, internal wave disturbance, imperfect tide correction, transducer mounting, long period heave, static draf...Vertical errors often present in multibeam swath bathymetric data. They are mainly sourced by sound refraction, internal wave disturbance, imperfect tide correction, transducer mounting, long period heave, static draft change, dynamic squat and dynamic motion residuals, etc. Although they can be partly removed or reduced by specific algorithms, the synthesized depth biases are unavoidable and sometimes have an important influence on high precise utilization of the final bathymetric data. In order to confidently identify the decimeter-level changes in seabed morphology by MBES, we must remove or weaken depth biases and improve the precision of multibeam bathymetry further. The fixed-interval profiles that are perpendicular to the vessel track are generated to adjust depth biases between swaths. We present a kind of postprocessing method to minimize the depth biases by the histogram of cumulative depth biases. The datum line in each profile can be obtained by the maximum value of histogram. The corrections of depth biases can be calculated according to the datum line. And then the quality of final bathymetry can be improved by the corrections. The method is verified by a field test.展开更多
We developed a three-step classification approach for forest road extraction utilizing LiDAR data. The first step employed the IDW method to interpolate LiDAR point data(first and last pulses) to achieve DSM, DTM and ...We developed a three-step classification approach for forest road extraction utilizing LiDAR data. The first step employed the IDW method to interpolate LiDAR point data(first and last pulses) to achieve DSM, DTM and DNTM layers(at 1 m resolution). For this interpolation RMSE was 0.19 m. In the second step, the Support Vector Machine(SVM) was employed to classify the LiDAR data into two classes, road and non-road. For this classification, SVM indicated the merged distance layer with intensity data and yielded better identification of the road position. Assessments of the obtained results showed 63% correctness, 75% completeness and 52% quality of classification. In the next step, road edges were defined in the LiDAR-extracted layers, enabling accurate digitizing of the centerline location. More than 95% of the LiDAR-derived road was digitized within 1.3 m to the field surveyed normal. The proposed approach can provide thorough and accurate road inventory data to support forest management.展开更多
Associations between the autumn Arctic sea ice concentration(SIC) and Asian winter temperature are discussed using the singular value decomposition analysis.Results show that in recent 33 years reduced autumn Arctic s...Associations between the autumn Arctic sea ice concentration(SIC) and Asian winter temperature are discussed using the singular value decomposition analysis.Results show that in recent 33 years reduced autumn Arctic sea ice is accompanied by Asian winter temperature decrease except in the Tibetan plateau and the Arctic Ocean and the North Pacific Ocean coast.The autumn SIC reduction excites two geopotential height centers in Eurasia and the north Arctic Ocean,which are persistent from autumn to winter.The negative center is in Barents Sea/Kara Sea.The positive center is located in Mongolia.The anomalous winds are associated with geopotential height centers,providing favorable clod air for the Asian winter temperature decreasing in recent 33 years.This relationship indicates a potential long-term outlook for the Asian winter temperature decrease as the decline of the autumn sea ice in the Arctic Ocean is expected to continue as climate warms.展开更多
The de-coherence phenomena such as Low-SNR radar signal, shadows and layover caused by topography, etc., causing phase data discontinuity, makes the result of unwrapping phase inaccuracy or even completely wrong. Base...The de-coherence phenomena such as Low-SNR radar signal, shadows and layover caused by topography, etc., causing phase data discontinuity, makes the result of unwrapping phase inaccuracy or even completely wrong. Based on the analysis of influencing factors to weight choice, this thesis develops a new method to choose the weights based on the measure of the confidence in the frequency domain. Experiments show that it could overcome the defect of sub-estimate to the slope of least squares method very well, which has a better rationale, stability and performance.展开更多
The decorelation phenomena such as Low-SNR radar signal,shadows and layover caused by topography etc,causes phase data discontinuous and makes the result of unwrapping phase inaccurate or completely wrong.Based on the...The decorelation phenomena such as Low-SNR radar signal,shadows and layover caused by topography etc,causes phase data discontinuous and makes the result of unwrapping phase inaccurate or completely wrong.Based on the analysis of influencing factors to the weight selection,this paper develops a new method to choose the weights based on the measurement of confidence in frequency domain.Results show that it is more precise and robust than other methods,and can make up for the defect of sub-estimate to the slope of least squares method.展开更多
With the development of mining industry,people have obtained profits from it,but they are facing environmental damages.In order to monitor these environmental changes,a spectral library is set up for the spectrum data...With the development of mining industry,people have obtained profits from it,but they are facing environmental damages.In order to monitor these environmental changes,a spectral library is set up for the spectrum data organization and management of mine typical objects.Most of the spectrum data come from the long-term field measuring in mining area and other spectral libraries.For the data quality control and error detection in the measuring data,an inner precision calculation method is presented and a series of interactive graphical controls are developed for the spectrum visualization and analysis.Through extracting and saving spectrum characters for the mine typical objects,realizs spectrum matching and classification for new measured spectrum samples are realized by using Euclidean distance,Aitchison distance,Pearson correlation coefficient and vector angular cosine methods.Based on the matching result,this work is able to gather dynamically physicochemical environment parameters from the library and gives an early warning for the mine environmental changes.展开更多
This study aimed to design a new ecological model to be used in ecological land evaluation studies, using careful planning and sustainable resource in the region to develop a sustainable and comprehensive guide. The r...This study aimed to design a new ecological model to be used in ecological land evaluation studies, using careful planning and sustainable resource in the region to develop a sustainable and comprehensive guide. The results show that the study area has grade 1 and 2 urban sustainable planning. Capability Grade 1 has area 149,724.25 hectares equal 22.2% that is stretched from East to West Range. Capability Grade2 has area about 18533.28 hectares equal 2.7% that located in the South, East and South-East. The results show that the greater extent of this range is the unsustainable site, which extends over 506,664.81 hectares, equivalent to 75.1% that areas compliance by the Environmental Protection area. Another important result of this study is unsuitable area for urban development and the range of services already seating workshops and industrial units. This is due to the arrival of adverse environmental effects. Therefore, this study has been focusing on reform of the status of land of use and the optimal use of the ecological potential.展开更多
The sensitivities of the initial value and the sampling information to the accuracy of a high accuracy surface modeling(HASM) are investigated and the implementations of this new modeling method are modified and enhan...The sensitivities of the initial value and the sampling information to the accuracy of a high accuracy surface modeling(HASM) are investigated and the implementations of this new modeling method are modified and enhanced. Based on the fundamental theorem of surface theory, HASM is developed to correct the error produced in geographical information system and ecological modeling process. However, the earlier version of HASM is theoretically incomplete and its initial value must be produced by other surface modeling methods, such as spline, which limit its promotion. In other words, we must use other interpolators to drive HASM. According to the fundamental theorem of surface theory, we modify HASM, namely HASM.MOD, by adding another important nonlinear equation to make it independent of other methods and, at the same time, have a complete and solid theory foundation. Two mathematic surfaces and monthly mean temperature of 1951–2010 are used to validate the effectiveness of the new method. Experiments show that the modified version of HASM is insensitive to the selection of initial value which is particular important for HASM. We analyze the sensitivities of sampling error and sampling ratio to the simulation accuracy of HASM.MOD. It is found that sampling information plays an important role in the simulation accuracy of HASM.MOD. Another feature of the modified version of HASM is that it is theoretically perfect as it considers the third equation of the surface theory which reflects the local warping of the surface. The modified HASM may be useful with a wide range of spatial interpolation as it would no longer rely on other interpolation methods.展开更多
Surface modeling with very large data sets is challenging. An efficient method for modeling massive data sets using the high accuracy surface modeling method(HASM) is proposed, and HASM_Big is developed to handle very...Surface modeling with very large data sets is challenging. An efficient method for modeling massive data sets using the high accuracy surface modeling method(HASM) is proposed, and HASM_Big is developed to handle very large data sets. A large data set is defined here as a large spatial domain with high resolution leading to a linear equation with matrix dimensions of hundreds of thousands. An augmented system approach is employed to solve the equality-constrained least squares problem(LSE) produced in HASM_Big, and a block row action method is applied to solve the corresponding very large matrix equations.A matrix partitioning method is used to avoid information redundancy among each block and thereby accelerate the model.Experiments including numerical tests and real-world applications are used to compare the performances of HASM_Big with its previous version, HASM. Results show that the memory storage and computing speed of HASM_Big are better than those of HASM. It is found that the computational cost of HASM_Big is linearly scalable, even with massive data sets. In conclusion,HASM_Big provides a powerful tool for surface modeling, especially when there are millions or more computing grid cells.展开更多
基金Under the auspices of National Youth Science Foundation of China(No.41001294)Key Project of National Natural Science Foundation of China(No.40930531)Research Fund of State Key Laboratory Resources and Environment Information System(No.2010KF0002SA)
文摘In China′s Loess Plateau area, gully head is the most active zone of a drainage system in gully areas. The differentiation of loess gully head follows geospatial patterns and reflects the process of the loess landform development and evolution of its drainage system to some extent. In this study, the geomorphic meaning, basic characteristics, morphological structure and the basic types of loess gully heads were systematically analysed. Then, the loess gully head′s conceptual model was established, and an extraction method based on Digital Elevation Model(DEM) for loess gully head features and elements was proposed. Through analysing the achieved statistics of loess gully head features, loess gully heads have apparently similar and different characteristics depending on the different loess landforms where they are found. The loess head characteristics reflect their growth period and evolution tendency to a certain degree, and they indirectly represent evolutionary mechanisms. In addition, the loess gully developmental stages and the evolutionary processes can be deduced by using loess gully head characteristics. This study is of great significance for development and improvement of the theoretical system for describing loess gully landforms.
基金Projects(41274007,40874001)supported by the National Natural Science Foundations of ChinaProjects(ZR2012DM001,ZR2010DQ020)supported by Shandong Province Natural Science Foundation,China+2 种基金Project(2011B04)supported by the Key Laboratory of Surveying and Mapping Technology on Island and Reef,National Administration of Surveying,Mapping and Geoinformation,ChinaProject(2011KYTD103)supported by SDUST Research Fund,ChinaProject(BS2013F013)supported by Shangdong Province Outstanding Youth Scientist Foundation,China
文摘The properties and feasibility of L-band differential InSAR for detecting and monitoring mining-induced subsidence were systematically analyzed and demonstrated. The largest monitored subsidence gradient of 7.9×10-3 and magnitude of 91 cm were firstly derived by theoretical derivation. Then, the stronger phase maintaining capacity and weaker sensitivity to minor land subsidence compared with C-band DInSAR were illustrated by phase simulation of the actual mine subsidence. Finally, the data processing procedure of two-pass DInSAR was further refined to accurately observe subsidence of a coalfield of Jining in Northern China using 7 ALOS PALSAR images. The largest monitored subsidence magnitude of 39.22 cm and other properties were better investigated by testing results interpretation and subsidence analysis, and the absolute difference varying from 0.5 mm to 17.9 mm was obtained by comparison with leveling-measured subsidence. All of results show that L-band DInSAR technique can investigate the location, amount, area and other detailed subsidence information with relatively higher accuracy.
基金jointly supported by a grant from the National Natural Science Foundation of China(General Program)(41071250)Innovation Program of State Key Laboratory of Resources and Environmental Information System,Institute of Geographic SciencesNatural Resources Research,Chinese Academy of Sciences(088RA500KA)
文摘Mesoscale eddies are an important component of oceanic features.How to automatically identify these mesoscale eddies from available data has become an important research topic.Through careful examination of existing methods,we propose an improved,SSH-based automatic identification method.Using the inclusion relation of enclosed SSH contours,the mesoscale eddy boundary and core(s) can be automatically identified.The time evolution of eddies can be examined by a threshold search algorithm and a tracking algorithm based on similarity.Sea-surface height(SSH) data from Naval Research Laboratory Layered Ocean Model(NLOM) and sea-level anomaly(SLA) data from altimeter are used in the many experiments,in which different automatic identification methods are compared.Our results indicate that the improved method is able to extract the mesoscale eddy boundary more precisely,retaining the multiple-core structure.In combination with the tracking algorithm,this method can capture complete mesoscale eddy processes.It can thus provide reliable information for further study of reconstructing eddy dynamics,merging,splitting,and evolution of a multi-core structure.
基金The Chinese Polar Environment Comprehensive Investigation and Assessment Programmes,State Oceanic Administration under contact Nos CHINARE2014-03-01 and CHINARE2014-04-03the Public Science and Technology Research Funds Projects of Ocean under contact No.201205007the Basic Research Operating Funds of the First Institute of Oceanography,State Oceanic Administration under contact Nos 2014T02 and 2014G02
文摘Associations between autumn Arctic sea ice concentration(SIC) and early winter precipitation in China are studied using singular value decomposition analysis. The results show that a reduced SIC almost everywhere in the Arctic Ocean, except the northern Greenland Sea and Canadian Basin, are accompanied by dry conditions over central China, extending northeast from the Tibetan Plateau toward the Japan Sea, the Bohai Sea and the Yellow Sea, and wet conditions over South China and North China. Atmospheric circulation anomalies associated with SIC variability show two wave-train structures, which are persistent from autumn to winter, leading to the identified relationship between autumn Arctic SIC and early winter precipitation in China. Given that the decline in autumn SIC in the Arctic Ocean is expected to continue as the climate warms, this relationship provides a possible long-term outlook for early winter precipitation in China.
基金supported by the project of Shandong Provincial National Science Foundation for Distinguished Young Scholars(Grant No.JQ201113)SDUST's National Science Foundation for Distinguished Young Scholars(Grant No.2010KYJQ102)
文摘Vertical errors often present in multibeam swath bathymetric data. They are mainly sourced by sound refraction, internal wave disturbance, imperfect tide correction, transducer mounting, long period heave, static draft change, dynamic squat and dynamic motion residuals, etc. Although they can be partly removed or reduced by specific algorithms, the synthesized depth biases are unavoidable and sometimes have an important influence on high precise utilization of the final bathymetric data. In order to confidently identify the decimeter-level changes in seabed morphology by MBES, we must remove or weaken depth biases and improve the precision of multibeam bathymetry further. The fixed-interval profiles that are perpendicular to the vessel track are generated to adjust depth biases between swaths. We present a kind of postprocessing method to minimize the depth biases by the histogram of cumulative depth biases. The datum line in each profile can be obtained by the maximum value of histogram. The corrections of depth biases can be calculated according to the datum line. And then the quality of final bathymetry can be improved by the corrections. The method is verified by a field test.
基金supported by Tarbiat Modares University(TMU)of Iran
文摘We developed a three-step classification approach for forest road extraction utilizing LiDAR data. The first step employed the IDW method to interpolate LiDAR point data(first and last pulses) to achieve DSM, DTM and DNTM layers(at 1 m resolution). For this interpolation RMSE was 0.19 m. In the second step, the Support Vector Machine(SVM) was employed to classify the LiDAR data into two classes, road and non-road. For this classification, SVM indicated the merged distance layer with intensity data and yielded better identification of the road position. Assessments of the obtained results showed 63% correctness, 75% completeness and 52% quality of classification. In the next step, road edges were defined in the LiDAR-extracted layers, enabling accurate digitizing of the centerline location. More than 95% of the LiDAR-derived road was digitized within 1.3 m to the field surveyed normal. The proposed approach can provide thorough and accurate road inventory data to support forest management.
基金The Basic Research Operating Funds of First Institute of Oceanography,State Oceanic Administration under contract No.2014T02 and 2014G02the Chinese Polar Environment Comprehensive Investigation and Assessment Programmes,State Oceanic Administration under contract Nos CHINARE2015-03-01 and CHINARE2015-04-03the Public science and Technology Research Funds Projects of Ocean under contract No.201205007
文摘Associations between the autumn Arctic sea ice concentration(SIC) and Asian winter temperature are discussed using the singular value decomposition analysis.Results show that in recent 33 years reduced autumn Arctic sea ice is accompanied by Asian winter temperature decrease except in the Tibetan plateau and the Arctic Ocean and the North Pacific Ocean coast.The autumn SIC reduction excites two geopotential height centers in Eurasia and the north Arctic Ocean,which are persistent from autumn to winter.The negative center is in Barents Sea/Kara Sea.The positive center is located in Mongolia.The anomalous winds are associated with geopotential height centers,providing favorable clod air for the Asian winter temperature decreasing in recent 33 years.This relationship indicates a potential long-term outlook for the Asian winter temperature decrease as the decline of the autumn sea ice in the Arctic Ocean is expected to continue as climate warms.
基金supported by the National Natural Science Fundation of China(40874001)Key Laboratory of Surveying and Mapping Technology on Island and Reef,National Administration of Surveying,Mapping and Geoinformation(2010A01)
文摘The de-coherence phenomena such as Low-SNR radar signal, shadows and layover caused by topography, etc., causing phase data discontinuity, makes the result of unwrapping phase inaccuracy or even completely wrong. Based on the analysis of influencing factors to weight choice, this thesis develops a new method to choose the weights based on the measure of the confidence in the frequency domain. Experiments show that it could overcome the defect of sub-estimate to the slope of least squares method very well, which has a better rationale, stability and performance.
基金supported by the National Natural Science Foundation(40874001)Key Laboratory of Surveying and Mapping Technology on Island and Reef,State Bureau of Surveying and Mapping(2010A01)
文摘The decorelation phenomena such as Low-SNR radar signal,shadows and layover caused by topography etc,causes phase data discontinuous and makes the result of unwrapping phase inaccurate or completely wrong.Based on the analysis of influencing factors to the weight selection,this paper develops a new method to choose the weights based on the measurement of confidence in frequency domain.Results show that it is more precise and robust than other methods,and can make up for the defect of sub-estimate to the slope of least squares method.
基金Supported by the National Key Technology R&D Program of China(No.2012BAH27B04,2011BAC03B03)the National Natural Science Foundation of China(No.41471330)+1 种基金Research Fund for the Doctoral Program of Higher Education of China(20113718110001)SDUST Research Fund(2011KYTD103)
文摘With the development of mining industry,people have obtained profits from it,but they are facing environmental damages.In order to monitor these environmental changes,a spectral library is set up for the spectrum data organization and management of mine typical objects.Most of the spectrum data come from the long-term field measuring in mining area and other spectral libraries.For the data quality control and error detection in the measuring data,an inner precision calculation method is presented and a series of interactive graphical controls are developed for the spectrum visualization and analysis.Through extracting and saving spectrum characters for the mine typical objects,realizs spectrum matching and classification for new measured spectrum samples are realized by using Euclidean distance,Aitchison distance,Pearson correlation coefficient and vector angular cosine methods.Based on the matching result,this work is able to gather dynamically physicochemical environment parameters from the library and gives an early warning for the mine environmental changes.
基金Acknowledgments The research is supported by the National Science Foundation of China (40874001) and National 863 plans projects of China (2009AA12Z147). The authors would like to express thanks to ESA (European Space Agency) for providing ENVISAT satellite data.
文摘This study aimed to design a new ecological model to be used in ecological land evaluation studies, using careful planning and sustainable resource in the region to develop a sustainable and comprehensive guide. The results show that the study area has grade 1 and 2 urban sustainable planning. Capability Grade 1 has area 149,724.25 hectares equal 22.2% that is stretched from East to West Range. Capability Grade2 has area about 18533.28 hectares equal 2.7% that located in the South, East and South-East. The results show that the greater extent of this range is the unsustainable site, which extends over 506,664.81 hectares, equivalent to 75.1% that areas compliance by the Environmental Protection area. Another important result of this study is unsuitable area for urban development and the range of services already seating workshops and industrial units. This is due to the arrival of adverse environmental effects. Therefore, this study has been focusing on reform of the status of land of use and the optimal use of the ecological potential.
基金supported by the Program of National Natural Science Foundation of China(Grant No.91325204)National Basic Research Priorities Program of Ministry of Science and Technology of the People’s Republic of China(Grant No.2010CB950904)+1 种基金National High-tech R&D Program of the Ministry of Science and Technology of the People’s Republic of China(Grant No.2013AA122003)the Key Program of National Natural Science of China(Grant No.41023010)
文摘The sensitivities of the initial value and the sampling information to the accuracy of a high accuracy surface modeling(HASM) are investigated and the implementations of this new modeling method are modified and enhanced. Based on the fundamental theorem of surface theory, HASM is developed to correct the error produced in geographical information system and ecological modeling process. However, the earlier version of HASM is theoretically incomplete and its initial value must be produced by other surface modeling methods, such as spline, which limit its promotion. In other words, we must use other interpolators to drive HASM. According to the fundamental theorem of surface theory, we modify HASM, namely HASM.MOD, by adding another important nonlinear equation to make it independent of other methods and, at the same time, have a complete and solid theory foundation. Two mathematic surfaces and monthly mean temperature of 1951–2010 are used to validate the effectiveness of the new method. Experiments show that the modified version of HASM is insensitive to the selection of initial value which is particular important for HASM. We analyze the sensitivities of sampling error and sampling ratio to the simulation accuracy of HASM.MOD. It is found that sampling information plays an important role in the simulation accuracy of HASM.MOD. Another feature of the modified version of HASM is that it is theoretically perfect as it considers the third equation of the surface theory which reflects the local warping of the surface. The modified HASM may be useful with a wide range of spatial interpolation as it would no longer rely on other interpolation methods.
基金supported by the National Natural Science Foundation of China (Grant Nos. 41541010, 41701456, 41421001, 41590840 & 91425304)the Key Programs of the Chinese Academy of Sciences (Grant No. QYZDY-SSW-DQC007)the Cultivate Project of Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences (Grant No. TSYJS03)
文摘Surface modeling with very large data sets is challenging. An efficient method for modeling massive data sets using the high accuracy surface modeling method(HASM) is proposed, and HASM_Big is developed to handle very large data sets. A large data set is defined here as a large spatial domain with high resolution leading to a linear equation with matrix dimensions of hundreds of thousands. An augmented system approach is employed to solve the equality-constrained least squares problem(LSE) produced in HASM_Big, and a block row action method is applied to solve the corresponding very large matrix equations.A matrix partitioning method is used to avoid information redundancy among each block and thereby accelerate the model.Experiments including numerical tests and real-world applications are used to compare the performances of HASM_Big with its previous version, HASM. Results show that the memory storage and computing speed of HASM_Big are better than those of HASM. It is found that the computational cost of HASM_Big is linearly scalable, even with massive data sets. In conclusion,HASM_Big provides a powerful tool for surface modeling, especially when there are millions or more computing grid cells.
基金partially supported by the National Key Research and Development Program of China(2016YFA0600104)supported by donations made by Delos Living LLC,and the Cyrus Tang Foundation+2 种基金supported by the National Natural Science Foundation of China(41471419)Beijing Institute of Urban Planningsupported by the Fundamental Research Funds for the Central Universities(CCNU19TD002).