Terrain Aided Navigation(TAN)technology has become increasingly important due to its effectiveness in environments where Global Positioning System(GPS)is unavailable.In recent years,TAN systems have been extensively r...Terrain Aided Navigation(TAN)technology has become increasingly important due to its effectiveness in environments where Global Positioning System(GPS)is unavailable.In recent years,TAN systems have been extensively researched for both aerial and underwater navigation applications.However,many TAN systems that rely on recursive Unmanned Aerial Vehicle(UAV)position estimation methods,such as Extended Kalman Filters(EKF),often face challenges with divergence and instability,particularly in highly non-linear systems.To address these issues,this paper proposes and investigates a hybrid two-stage TAN positioning system for UAVs that utilizes Particle Filter.To enhance the system’s robustness against uncertainties caused by noise and to estimate additional system states,a Fuzzy Particle Filter(FPF)is employed in the first stage.This approach introduces a novel terrain composite feature that enables a fuzzy expert system to analyze terrain non-linearities and dynamically adjust the number of particles in real-time.This design allows the UAV to be efficiently localized in GPS-denied environments while also reducing the computational complexity of the particle filter in real-time applications.In the second stage,an Error State Kalman Filter(ESKF)is implemented to estimate the UAV’s altitude.The ESKF is chosen over the conventional EKF method because it is more suitable for non-linear systems.Simulation results demonstrate that the proposed fuzzy-based terrain composite method achieves high positional accuracy while reducing computational time and memory usage.展开更多
Laser scanning technology has been widely used in landslide aspects.However,the existing deformation analysis based on terrain laser scanners can only provide limited information,which is insufficient for understandin...Laser scanning technology has been widely used in landslide aspects.However,the existing deformation analysis based on terrain laser scanners can only provide limited information,which is insufficient for understanding landslide kinematics and failure mechanisms.To overcome this limitation,this paper proposes an automated method for processing point clouds collected in landslide physical modeling.This method allows the acquisition of quantitative three-dimensional(3D)deformation field information.The results show the organized and spatially related point cloud segmentation in terms of spherical targets.The segmented point clouds can be fitted to determine the locations of all preset targets and their corresponding location changes.The proposed method has been validated based on theoretical analysis and numerical and physical tests,which indicates that it can batch-process massive data sets with high computational efficiency and good noise resistance.Compared to existing methods,this method shows a significant potential for understanding landslide kinematics and failure mechanisms and advancing the application of 3D laser scanning in geotechnical modeling.展开更多
Numerical weather prediction(NWP)models have always presented large forecasting errors of surface wind speeds over regions with complex terrain.In this study,surface wind forecasts from an operational NWP model,the SM...Numerical weather prediction(NWP)models have always presented large forecasting errors of surface wind speeds over regions with complex terrain.In this study,surface wind forecasts from an operational NWP model,the SMS-WARR(Shanghai Meteorological Service-WRF ADAS Rapid Refresh System),are analyzed to quantitatively reveal the relationships between the forecasted surface wind speed errors and terrain features,with the intent of providing clues to better apply the NWP model to complex terrain regions.The terrain features are described by three parameters:the standard deviation of the model grid-scale orography,terrain height error of the model,and slope angle.The results show that the forecast bias has a unimodal distribution with a change in the standard deviation of orography.The minimum ME(the mean value of bias)is 1.2 m s^(-1) when the standard deviation is between 60 and 70 m.A positive correlation exists between bias and terrain height error,with the ME increasing by 10%−30%for every 200 m increase in terrain height error.The ME decreases by 65.6%when slope angle increases from(0.5°−1.5°)to larger than 3.5°for uphill winds but increases by 35.4%when the absolute value of slope angle increases from(0.5°−1.5°)to(2.5°−3.5°)for downhill winds.Several sensitivity experiments are carried out with a model output statistical(MOS)calibration model for surface wind speeds and ME(RMSE)has been reduced by 90%(30%)by introducing terrain parameters,demonstrating the value of this study.展开更多
This article proposes a VGG network with histogram of oriented gradient(HOG) feature fusion(HOG-VGG) for polarization synthetic aperture radar(PolSAR) image terrain classification.VGG-Net has a strong ability of deep ...This article proposes a VGG network with histogram of oriented gradient(HOG) feature fusion(HOG-VGG) for polarization synthetic aperture radar(PolSAR) image terrain classification.VGG-Net has a strong ability of deep feature extraction,which can fully extract the global deep features of different terrains in PolSAR images,so it is widely used in PolSAR terrain classification.However,VGG-Net ignores the local edge & shape features,resulting in incomplete feature representation of the PolSAR terrains,as a consequence,the terrain classification accuracy is not promising.In fact,edge and shape features play an important role in PolSAR terrain classification.To solve this problem,a new VGG network with HOG feature fusion was specifically proposed for high-precision PolSAR terrain classification.HOG-VGG extracts both the global deep semantic features and the local edge & shape features of the PolSAR terrains,so the terrain feature representation completeness is greatly elevated.Moreover,HOG-VGG optimally fuses the global deep features and the local edge & shape features to achieve the best classification results.The superiority of HOG-VGG is verified on the Flevoland,San Francisco and Oberpfaffenhofen datasets.Experiments show that the proposed HOG-VGG achieves much better PolSAR terrain classification performance,with overall accuracies of 97.54%,94.63%,and 96.07%,respectively.展开更多
A rainstorm caused by mesoscale convective system (MCS) in Guizhou Province in June 25-26 in 2005 was simulated with the MM5 model. Based on the good simulated results of the MCS developing and the clouds physics proc...A rainstorm caused by mesoscale convective system (MCS) in Guizhou Province in June 25-26 in 2005 was simulated with the MM5 model. Based on the good simulated results of the MCS developing and the clouds physics process, and by means of reducing the height of Yunnan-Guizhou Plateau and cutting off the middle-east of the Yunnan-Guizhou Plateau on the simulated tests, the question as how the ladder terrain on the west of Yunnan-Guizhou Plateau impact on the rainstorm of Guizhou was studied. The analysis results showed that the second ladder terrain of Yunnan-Guizhou Plateau only affected the development of convective clouds on its backward position,and hardly affected the rain on its upward. The whole terrain of the Yunnan-Guizhou Plateau had a distinct impact not only on the windward slope rainfall of the west of the plateau, but also on the rainfall distribution, intensity and continuing time of the convective clouds on the middle-east of the plateau.展开更多
To improve the navigation accuracy of an autonomous underwater vehicle (AUV), a novel terrain passive integrated navigation system (TPINS) is presented. According to the characteristics of the underwater environme...To improve the navigation accuracy of an autonomous underwater vehicle (AUV), a novel terrain passive integrated navigation system (TPINS) is presented. According to the characteristics of the underwater environment and AUV navigation requirements of low cost and high accuracy, a novel TPINS is designed with a configuration of the strapdown inertial navigation system (SINS), the terrain reference navigation system (TRNS), the Doppler velocity sonar (DVS), the magnetic compass and the navigation computer utilizing the unscented Kalman filter (UKF) to fuse the navigation information from various navigation sensors. Linear filter equations for the extended Kalman filter (EKF), nonlinear filter equations for the UKF and measurement equations of navigation sensors are addressed. It is indicated from the comparable simulation experiments of the EKF and the UKF that AUV navigation precision is improved substantially with the proposed sensors and the UKF when compared to the EKF. The TPINS designed with the proposed sensors and the UKF is effective in reducing AUV navigation position errors and improving the stability and precision of the AUV underwater integrated navigation.展开更多
In order to slove a realistic test problem of TF/TA algorithm, it is necessary to make a discrete terrain model. The model has adjustable roughness parameters so as to test the optimization procedure for different typ...In order to slove a realistic test problem of TF/TA algorithm, it is necessary to make a discrete terrain model. The model has adjustable roughness parameters so as to test the optimization procedure for different types of terrain. In this paper, an algorithm to generate random terrain data is given.展开更多
Most of the current computing methods used to determine the magnetic field of a uniformly magnetized cuboid assume that the observation point is located in the upper half space without a source. However, such methods ...Most of the current computing methods used to determine the magnetic field of a uniformly magnetized cuboid assume that the observation point is located in the upper half space without a source. However, such methods may generate analytical singularities for conditions of undulating terrain. Based on basic geomagnetic field theories, in this study an improved magnetic field expression is derived using an integration method of variable substitution, and all singularity problems for the entire space without a source are discussed and solved. This integration process is simpler than that of previous methods, and final integral results with a more uniform form. AT at all points in the source-flee space can be calculated without requiring coordinate transformation; thus forward modeling is also simplified. Corresponding model tests indicate that the new magnetic field expression is more correct because there is no analytical singularity and can be used with undulating terrain.展开更多
[目的]分析雅砻河流域生态系统服务价值的变化,为流域生态文明建设提供科学依据。[方法]基于高精度土地利用数据,运用地形梯度分级、标准差椭圆和地理探测器,探究了近20年流域生态系统服务价值(Ecosystem service value,ESV)的时空演变...[目的]分析雅砻河流域生态系统服务价值的变化,为流域生态文明建设提供科学依据。[方法]基于高精度土地利用数据,运用地形梯度分级、标准差椭圆和地理探测器,探究了近20年流域生态系统服务价值(Ecosystem service value,ESV)的时空演变和地形梯度分异特征。[结果](1)草地是流域主要的土地利用类型,占比达92.45%。研究期内林地、灌木和草地面积有所减少,其他类型有所增加;(2)流域ESV在研究期间增长了0.13%(0.16亿元),ESV在空间上呈现出东南部高西北部低的特征。随着海拔、坡度和地形起伏度的抬升,ESV呈先增后降的分布规律,坡向梯度呈阴坡高于阳坡的特征。(3)流域ESV的标准差椭圆和重心在研究期内逐渐向西南方向移动,并趋于聚集。ESV的空间分异主要受自然和经济因子的影响,其中年均地温是主导因子(q=0.24),任意因子交互均增强了其分异性。[结论]依据雅砻河流域ESV的时空和地形分布特征,采取因地制宜的生态保护措施,促进流域生态环境的可持续发展。展开更多
The Hengshan complex forms part of the central zone of the North China Craton and consists predominantly of ductilely-deformed late Archaean to Palaeoproterozoic high-grade, partly migmatitic, granitoid orthogneisses,...The Hengshan complex forms part of the central zone of the North China Craton and consists predominantly of ductilely-deformed late Archaean to Palaeoproterozoic high-grade, partly migmatitic, granitoid orthogneisses, intruded by mafic dykes of gabbroic composition. Many highly strained rocks were previously misinterpreted as supracrustal sequences and represent mylonitized granitoids and sheared dykes. Our single zircon dating documents magmatic granitoid emplacement ages between 2.52 Ga and 2.48 Ga, with rare occurrences of 2.7 Ga gneisses, possibly reflecting an older basement. A few granitic gneisses have emplacement ages between 2.35 and 2.1 Ga and show the same structural features as the older rocks, indicating that the main deformation occurred after -2.1 Ga. Intrusion of gabbroic dykes occurred at -1920 Ma, and all Hengshan rocks underwent granulite-facies metamorphism at 1.88-1.85 Ga, followed by retrogression, sheafing and uplift. We interpret the Hengshan and adjacent Fuping granitoid gneisses as the lower, plutonic, part of a late Archaean to early Palaeoproterozoic Japan-type magmatic arc, with the upper, volcanic part represented by the nearby Wutai complex. Components of this arc may have evolved at a continental margin as indicated by the 2.7 Ga zircons. Major deformation and HP metamorphism occurred in the late Palaeoproterozoic during the Luliang orogeny when the Eastern and Western blocks of the North China Craton collided to form the Trans-North China orogen. Shear zones in the Hengshan are interpreted as major lower crustal discontinuities post-dating the peak of HP metamorphism, and we suggest that they formed during orogenic collapse and uplift of the Hengshan complex in the late Palaeoproterozoic (〈1.85 Ga).展开更多
The Loess positive and negative terrains (P-N terrains), which are widely distributed on the Loess Plateau, are discussed for the first time by introducing its characteristic, demarcation as well as extraction metho...The Loess positive and negative terrains (P-N terrains), which are widely distributed on the Loess Plateau, are discussed for the first time by introducing its characteristic, demarcation as well as extraction method from high-resolution Digital Elevation Models. Using 5 m-resolution DEMs as original test data, P-N terrains of 48 geomorphological units in different parts of Shaanxi Loess Plateau are extracted accurately. Then six indicators for depicting the geomorphologic landscape and spatial configuration characteristic of P-N terrains are proposed. The spatial distribution rules of these indicators and the relationship between the P-N terrains and Loess relief are discussed for further understanding of Loess landforms. Finally, with the integration of P-N terrains and traditional terrain indices, a series of un-supervised classification methods are applied to make a proper landform classification in northern Shaanxi. Results show that P-N terrains are an effect clue to reveal energy and substance distribution rules on the Loess Plateau. A continuous change of P-N terrains from south to north in Shaanxi Loess Plateau shows an obvious spatial difference of Loess land-forms and the positive terrain area only accounted for 60.5% in this region. The P-N terrains participant landform classification method increases validity of the result, especially in the Loess tableland, Loess tableland-ridge and the Loess low-hill area. This research is significant on the study of Loess landforms with the Digital Terrains Analysis methods.展开更多
文摘Terrain Aided Navigation(TAN)technology has become increasingly important due to its effectiveness in environments where Global Positioning System(GPS)is unavailable.In recent years,TAN systems have been extensively researched for both aerial and underwater navigation applications.However,many TAN systems that rely on recursive Unmanned Aerial Vehicle(UAV)position estimation methods,such as Extended Kalman Filters(EKF),often face challenges with divergence and instability,particularly in highly non-linear systems.To address these issues,this paper proposes and investigates a hybrid two-stage TAN positioning system for UAVs that utilizes Particle Filter.To enhance the system’s robustness against uncertainties caused by noise and to estimate additional system states,a Fuzzy Particle Filter(FPF)is employed in the first stage.This approach introduces a novel terrain composite feature that enables a fuzzy expert system to analyze terrain non-linearities and dynamically adjust the number of particles in real-time.This design allows the UAV to be efficiently localized in GPS-denied environments while also reducing the computational complexity of the particle filter in real-time applications.In the second stage,an Error State Kalman Filter(ESKF)is implemented to estimate the UAV’s altitude.The ESKF is chosen over the conventional EKF method because it is more suitable for non-linear systems.Simulation results demonstrate that the proposed fuzzy-based terrain composite method achieves high positional accuracy while reducing computational time and memory usage.
基金the National Natural Science Foundation of China(Grant No.42020104006).
文摘Laser scanning technology has been widely used in landslide aspects.However,the existing deformation analysis based on terrain laser scanners can only provide limited information,which is insufficient for understanding landslide kinematics and failure mechanisms.To overcome this limitation,this paper proposes an automated method for processing point clouds collected in landslide physical modeling.This method allows the acquisition of quantitative three-dimensional(3D)deformation field information.The results show the organized and spatially related point cloud segmentation in terms of spherical targets.The segmented point clouds can be fitted to determine the locations of all preset targets and their corresponding location changes.The proposed method has been validated based on theoretical analysis and numerical and physical tests,which indicates that it can batch-process massive data sets with high computational efficiency and good noise resistance.Compared to existing methods,this method shows a significant potential for understanding landslide kinematics and failure mechanisms and advancing the application of 3D laser scanning in geotechnical modeling.
基金supported by the National Natural Science Foundation of China(No.U2142206).
文摘Numerical weather prediction(NWP)models have always presented large forecasting errors of surface wind speeds over regions with complex terrain.In this study,surface wind forecasts from an operational NWP model,the SMS-WARR(Shanghai Meteorological Service-WRF ADAS Rapid Refresh System),are analyzed to quantitatively reveal the relationships between the forecasted surface wind speed errors and terrain features,with the intent of providing clues to better apply the NWP model to complex terrain regions.The terrain features are described by three parameters:the standard deviation of the model grid-scale orography,terrain height error of the model,and slope angle.The results show that the forecast bias has a unimodal distribution with a change in the standard deviation of orography.The minimum ME(the mean value of bias)is 1.2 m s^(-1) when the standard deviation is between 60 and 70 m.A positive correlation exists between bias and terrain height error,with the ME increasing by 10%−30%for every 200 m increase in terrain height error.The ME decreases by 65.6%when slope angle increases from(0.5°−1.5°)to larger than 3.5°for uphill winds but increases by 35.4%when the absolute value of slope angle increases from(0.5°−1.5°)to(2.5°−3.5°)for downhill winds.Several sensitivity experiments are carried out with a model output statistical(MOS)calibration model for surface wind speeds and ME(RMSE)has been reduced by 90%(30%)by introducing terrain parameters,demonstrating the value of this study.
基金Sponsored by the Fundamental Research Funds for the Central Universities of China(Grant No.PA2023IISL0098)the Hefei Municipal Natural Science Foundation(Grant No.202201)+1 种基金the National Natural Science Foundation of China(Grant No.62071164)the Open Fund of Information Materials and Intelligent Sensing Laboratory of Anhui Province(Anhui University)(Grant No.IMIS202214 and IMIS202102)。
文摘This article proposes a VGG network with histogram of oriented gradient(HOG) feature fusion(HOG-VGG) for polarization synthetic aperture radar(PolSAR) image terrain classification.VGG-Net has a strong ability of deep feature extraction,which can fully extract the global deep features of different terrains in PolSAR images,so it is widely used in PolSAR terrain classification.However,VGG-Net ignores the local edge & shape features,resulting in incomplete feature representation of the PolSAR terrains,as a consequence,the terrain classification accuracy is not promising.In fact,edge and shape features play an important role in PolSAR terrain classification.To solve this problem,a new VGG network with HOG feature fusion was specifically proposed for high-precision PolSAR terrain classification.HOG-VGG extracts both the global deep semantic features and the local edge & shape features of the PolSAR terrains,so the terrain feature representation completeness is greatly elevated.Moreover,HOG-VGG optimally fuses the global deep features and the local edge & shape features to achieve the best classification results.The superiority of HOG-VGG is verified on the Flevoland,San Francisco and Oberpfaffenhofen datasets.Experiments show that the proposed HOG-VGG achieves much better PolSAR terrain classification performance,with overall accuracies of 97.54%,94.63%,and 96.07%,respectively.
文摘A rainstorm caused by mesoscale convective system (MCS) in Guizhou Province in June 25-26 in 2005 was simulated with the MM5 model. Based on the good simulated results of the MCS developing and the clouds physics process, and by means of reducing the height of Yunnan-Guizhou Plateau and cutting off the middle-east of the Yunnan-Guizhou Plateau on the simulated tests, the question as how the ladder terrain on the west of Yunnan-Guizhou Plateau impact on the rainstorm of Guizhou was studied. The analysis results showed that the second ladder terrain of Yunnan-Guizhou Plateau only affected the development of convective clouds on its backward position,and hardly affected the rain on its upward. The whole terrain of the Yunnan-Guizhou Plateau had a distinct impact not only on the windward slope rainfall of the west of the plateau, but also on the rainfall distribution, intensity and continuing time of the convective clouds on the middle-east of the plateau.
基金Pre-Research Program of General Armament Department during the11th Five-Year Plan Period (No51309020503)the National Defense Basic Research Program of China (973Program)(No973-61334)+1 种基金the National Natural Science Foundation of China(No50575042)Specialized Research Fund for the Doctoral Program of Higher Education (No20050286026)
文摘To improve the navigation accuracy of an autonomous underwater vehicle (AUV), a novel terrain passive integrated navigation system (TPINS) is presented. According to the characteristics of the underwater environment and AUV navigation requirements of low cost and high accuracy, a novel TPINS is designed with a configuration of the strapdown inertial navigation system (SINS), the terrain reference navigation system (TRNS), the Doppler velocity sonar (DVS), the magnetic compass and the navigation computer utilizing the unscented Kalman filter (UKF) to fuse the navigation information from various navigation sensors. Linear filter equations for the extended Kalman filter (EKF), nonlinear filter equations for the UKF and measurement equations of navigation sensors are addressed. It is indicated from the comparable simulation experiments of the EKF and the UKF that AUV navigation precision is improved substantially with the proposed sensors and the UKF when compared to the EKF. The TPINS designed with the proposed sensors and the UKF is effective in reducing AUV navigation position errors and improving the stability and precision of the AUV underwater integrated navigation.
文摘In order to slove a realistic test problem of TF/TA algorithm, it is necessary to make a discrete terrain model. The model has adjustable roughness parameters so as to test the optimization procedure for different types of terrain. In this paper, an algorithm to generate random terrain data is given.
基金supported by China Geological Survey Northeastern Tarim Aeromagnetic and Aerogravity comprehensive survey project(No.12120115039401)
文摘Most of the current computing methods used to determine the magnetic field of a uniformly magnetized cuboid assume that the observation point is located in the upper half space without a source. However, such methods may generate analytical singularities for conditions of undulating terrain. Based on basic geomagnetic field theories, in this study an improved magnetic field expression is derived using an integration method of variable substitution, and all singularity problems for the entire space without a source are discussed and solved. This integration process is simpler than that of previous methods, and final integral results with a more uniform form. AT at all points in the source-flee space can be calculated without requiring coordinate transformation; thus forward modeling is also simplified. Corresponding model tests indicate that the new magnetic field expression is more correct because there is no analytical singularity and can be used with undulating terrain.
文摘[目的]分析雅砻河流域生态系统服务价值的变化,为流域生态文明建设提供科学依据。[方法]基于高精度土地利用数据,运用地形梯度分级、标准差椭圆和地理探测器,探究了近20年流域生态系统服务价值(Ecosystem service value,ESV)的时空演变和地形梯度分异特征。[结果](1)草地是流域主要的土地利用类型,占比达92.45%。研究期内林地、灌木和草地面积有所减少,其他类型有所增加;(2)流域ESV在研究期间增长了0.13%(0.16亿元),ESV在空间上呈现出东南部高西北部低的特征。随着海拔、坡度和地形起伏度的抬升,ESV呈先增后降的分布规律,坡向梯度呈阴坡高于阳坡的特征。(3)流域ESV的标准差椭圆和重心在研究期内逐渐向西南方向移动,并趋于聚集。ESV的空间分异主要受自然和经济因子的影响,其中年均地温是主导因子(q=0.24),任意因子交互均增强了其分异性。[结论]依据雅砻河流域ESV的时空和地形分布特征,采取因地制宜的生态保护措施,促进流域生态环境的可持续发展。
基金funded by the German Science Foundation(DFG,grant Kr 590/62)the National Natural Science Foundation of China(grant Nos.49832030,49772143 and 49572140 to Li Jianghai)
文摘The Hengshan complex forms part of the central zone of the North China Craton and consists predominantly of ductilely-deformed late Archaean to Palaeoproterozoic high-grade, partly migmatitic, granitoid orthogneisses, intruded by mafic dykes of gabbroic composition. Many highly strained rocks were previously misinterpreted as supracrustal sequences and represent mylonitized granitoids and sheared dykes. Our single zircon dating documents magmatic granitoid emplacement ages between 2.52 Ga and 2.48 Ga, with rare occurrences of 2.7 Ga gneisses, possibly reflecting an older basement. A few granitic gneisses have emplacement ages between 2.35 and 2.1 Ga and show the same structural features as the older rocks, indicating that the main deformation occurred after -2.1 Ga. Intrusion of gabbroic dykes occurred at -1920 Ma, and all Hengshan rocks underwent granulite-facies metamorphism at 1.88-1.85 Ga, followed by retrogression, sheafing and uplift. We interpret the Hengshan and adjacent Fuping granitoid gneisses as the lower, plutonic, part of a late Archaean to early Palaeoproterozoic Japan-type magmatic arc, with the upper, volcanic part represented by the nearby Wutai complex. Components of this arc may have evolved at a continental margin as indicated by the 2.7 Ga zircons. Major deformation and HP metamorphism occurred in the late Palaeoproterozoic during the Luliang orogeny when the Eastern and Western blocks of the North China Craton collided to form the Trans-North China orogen. Shear zones in the Hengshan are interpreted as major lower crustal discontinuities post-dating the peak of HP metamorphism, and we suggest that they formed during orogenic collapse and uplift of the Hengshan complex in the late Palaeoproterozoic (〈1.85 Ga).
基金Key Project of National Natural Science Foundation of China, No.40930531 National Youth Science Foundation of China, No.40801148 Anhui Provincial Natural Science Foundation. No. 090412062
文摘The Loess positive and negative terrains (P-N terrains), which are widely distributed on the Loess Plateau, are discussed for the first time by introducing its characteristic, demarcation as well as extraction method from high-resolution Digital Elevation Models. Using 5 m-resolution DEMs as original test data, P-N terrains of 48 geomorphological units in different parts of Shaanxi Loess Plateau are extracted accurately. Then six indicators for depicting the geomorphologic landscape and spatial configuration characteristic of P-N terrains are proposed. The spatial distribution rules of these indicators and the relationship between the P-N terrains and Loess relief are discussed for further understanding of Loess landforms. Finally, with the integration of P-N terrains and traditional terrain indices, a series of un-supervised classification methods are applied to make a proper landform classification in northern Shaanxi. Results show that P-N terrains are an effect clue to reveal energy and substance distribution rules on the Loess Plateau. A continuous change of P-N terrains from south to north in Shaanxi Loess Plateau shows an obvious spatial difference of Loess land-forms and the positive terrain area only accounted for 60.5% in this region. The P-N terrains participant landform classification method increases validity of the result, especially in the Loess tableland, Loess tableland-ridge and the Loess low-hill area. This research is significant on the study of Loess landforms with the Digital Terrains Analysis methods.