Method and device for determining saddle region, electronic equipment and storage medium
By identifying saddle areas in the digital elevation model of a geographic region and using the curve parameters of the initial grid and topographic profile lines, the problem of inaccurate saddle area identification was solved, enabling accurate positioning and severity assessment of saddle areas, thus improving the efficiency and safety of fire fighting.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- TSINGHUA UNIVERSITY
- Filing Date
- 2023-12-25
- Publication Date
- 2026-06-19
AI Technical Summary
Existing technologies make it difficult to accurately determine saddle areas, leading to inaccurate identification of saddle areas during forest fire fighting and increasing the risk of casualties.
By acquiring a digital elevation model of a geographic region, using each initial grid as the central grid, N initial grids are determined as target grid regions. Initial terrain profile lines are determined along the four axes of symmetry of the target grid regions. The curve parameters of each profile line are calculated to determine whether a local geographic region is a saddle region and to determine the degree of saddleness.
Accurate identification of saddle areas improves the accuracy of saddle area determination, helps predict the difficulty of fire fighting, and improves the efficiency and safety of forest fire fighting.
Smart Images

Figure CN117830376B_ABST
Abstract
Description
Technical Field
[0001] This disclosure relates to the field of geographic information spatial analysis technology, specifically to a method, apparatus, electronic device, and storage medium for determining a saddle region. Background Technology
[0002] The saddle area refers to the middle part of a connected mountain range. The saddle area has a relatively high altitude and steep sides. When the wind crosses the saddle area of the ridge, it will form horizontal and vertical vortices, known as saddle eddies. Affected by the eddies, the forest fire situation in the saddle area is extremely complex and can easily cause casualties. It is an area that must be given high priority in the process of fighting forest fires.
[0003] Therefore, it is necessary to accurately determine the saddle area. Summary of the Invention
[0004] This disclosure aims to at least partially address one of the technical problems in the related art.
[0005] Therefore, the purpose of this disclosure is to provide a method, apparatus, electronic device and storage medium for determining saddle regions, which can accurately determine whether a geographical area is a saddle region and effectively improve the determination effect of saddle regions.
[0006] To achieve the above objectives, the method for determining the saddle region proposed in the first aspect of this disclosure includes:
[0007] Obtain a digital elevation model of the geographic region, wherein the digital elevation model includes: multiple initial grids;
[0008] Determine N using each initial grid as the center grid. N initial grids are used as the target grid region, where N is an odd number greater than 1;
[0009] Four initial terrain profile lines are determined along the four axes of symmetry of the target grid region;
[0010] Determine the curve parameters for each initial terrain profile line, where the curve parameters describe the slope change of the initial terrain profile line as it passes through the central grid.
[0011] Based on the curve parameters of each initial topographic profile line, determine whether the local geographic region corresponding to the initial grid is a saddle region.
[0012] In some embodiments of this disclosure, determining the curve parameters for each initial topographic profile line includes:
[0013] If there are no extreme points of the initial curve in the initial terrain profile, then the curve parameter corresponding to the initial terrain profile is determined to be 0.
[0014] In some embodiments of this disclosure, determining the curve parameters of each initial topographic profile line further includes:
[0015] If there are initial curve extrema in the initial terrain profile, then determine the target curve extrema that is closest to the center point of the central grid from at least one initial curve extrema.
[0016] Extract monotonically changing portions of the initial terrain profile from both sides of the extreme points of the target curve from the initial terrain profile.
[0017] The initial terrain profile lines that change monotonically on both sides of the extreme point of the target curve are used together as the target terrain profile line.
[0018] Based on the target terrain profile, determine the curve parameters of the initial terrain profile.
[0019] In some embodiments of this disclosure, curve parameters are determined based on the target terrain profile, including:
[0020] If the extreme point of the target curve in the target terrain profile is a minimum point, then the curve parameter is determined to be negative.
[0021] If the extreme point of the target curve in the target terrain profile is a maximum point, then the curve parameter is determined to be a positive number.
[0022] In some embodiments of this disclosure, determining whether a local geographic region corresponding to an initial grid is a saddle region based on the curve parameters of each initial terrain profile line includes:
[0023] Determine the parametric product between the curve parameters of the initial terrain profile corresponding to each pair of mutually perpendicular axes of symmetry;
[0024] If the product of at least one parameter is less than 0, then the local geographic region corresponding to the initial grid is determined to be a saddle region.
[0025] If the product of the two parameters is greater than 0, or the product of the two parameters is equal to 0, then the local geographic region corresponding to the initial grid is determined not to be a saddle region.
[0026] In some embodiments of this disclosure, it also includes:
[0027] Determine the saddle degree information corresponding to each saddle region, whereby the saddle degree information is used to describe the topographic relief of the saddle region.
[0028] In some embodiments of this disclosure, determining saddle degree information corresponding to each saddle region includes:
[0029] If the parameter product is less than 0, then a pair of target terrain profiles corresponding to the parameter product are determined;
[0030] Based on a pair of target terrain profiles, determine the target saddle degree value corresponding to each saddle region;
[0031] The target saddle depth value is used as the saddle depth information.
[0032] In some embodiments of this disclosure, determining the target saddle degree value corresponding to each saddle region based on a pair of target terrain profiles includes:
[0033] Determine the first average slope of each target terrain profile line in a pair of target terrain profile lines, the first elevation value of the extreme point of the target curve in each target terrain profile line, and the horizontal distance between the extreme point of the target curve and the center point in each target terrain profile line;
[0034] Determine the maximum and minimum elevation values of a pair of target terrain profiles, and the second average slope of the first terrain profile in each target terrain profile, wherein the first terrain profile is a portion of the target terrain profile between the extreme point and the center point of the target curve in the target terrain profile.
[0035] Determine the second elevation value of the central grid;
[0036] The candidate saddle depth value is determined by multiplying the first average slope, the second average slope, the maximum elevation value, the minimum elevation value, the horizontal distance, the first elevation value, the second elevation value, and the parameters corresponding to a pair of target terrain profile lines.
[0037] Based on the candidate saddle severity values, determine the target saddle severity value corresponding to each saddle region.
[0038] In some embodiments of this disclosure, the candidate saddle depth value is determined based on the product of a first average slope, a second average slope, a maximum elevation value, a minimum elevation value, a horizontal distance, a first elevation value, a second elevation value, and a pair of parameters corresponding to target terrain profile lines, including:
[0039] The candidate saddle point severity value is calculated using the following formula:
[0040] ;
[0041] Where D is the candidate saddle point severity value, and T is the parameter product. , This represents the horizontal distance between a pair of target terrain profile lines. This is the second elevation value. , These are the first elevation values corresponding to a pair of target terrain profile lines. and Let be the first average slope corresponding to a pair of target terrain profile lines. For the maximum elevation value, The minimum elevation value, = + , and These are the second average slopes corresponding to each of the first terrain profiles in a pair of target terrain profiles.
[0042] In some embodiments of this disclosure, determining a target saddle severity value corresponding to each saddle region based on candidate saddle severity values includes:
[0043] If the number of candidate saddle severity values is 1, then the candidate saddle severity value is used as the target saddle severity value.
[0044] If the number of candidate saddle severity values is 2, then the candidate saddle severity value with the larger value among the two candidate saddle severity values is taken as the target saddle severity value.
[0045] To achieve the above objectives, a device for determining a saddle region is provided in a second aspect embodiment of this disclosure. The device includes:
[0046] The acquisition module is used to acquire the digital elevation model of a geographic area, wherein the digital elevation model includes: multiple initial grids;
[0047] The first determining module is used to determine N, with each initial grid as the center grid. N initial grids are used as the target grid region, where N is an odd number greater than 1;
[0048] The second determining module is used to determine four initial terrain profile lines along the four axes of symmetry of the target grid region, respectively.
[0049] The third determination module is used to determine the curve parameters of each initial terrain profile line. The curve parameters are used to describe the slope change of the initial terrain profile line when it passes through the central grid.
[0050] The fourth determination module is used to determine whether the local geographic region corresponding to the initial grid is a saddle region based on the curve parameters of each initial terrain profile line.
[0051] In some embodiments of this disclosure, the third determining module is further configured to:
[0052] If there are no extreme points of the initial curve in the initial terrain profile, then the curve parameter corresponding to the initial terrain profile is determined to be 0.
[0053] In some embodiments of this disclosure, the third determining module is further configured to:
[0054] If there are initial curve extrema in the initial terrain profile, then determine the target curve extrema that is closest to the center point of the central grid from at least one initial curve extrema.
[0055] Extract monotonically changing portions of the initial terrain profile from both sides of the extreme points of the target curve from the initial terrain profile.
[0056] The initial terrain profile lines that change monotonically on both sides of the extreme point of the target curve are used together as the target terrain profile line.
[0057] Determine the curve parameters of the initial terrain profile based on the target terrain profile.
[0058] In some embodiments of this disclosure, the third determining module is further configured to:
[0059] If the extreme point of the target curve in the target terrain profile is a minimum point, then the curve parameter is determined to be negative.
[0060] If the extreme point of the target curve in the target terrain profile is a maximum point, then the curve parameter is determined to be a positive number.
[0061] In some embodiments of this disclosure, the fourth determining module is further configured to:
[0062] Determine the parametric product between the curve parameters of the initial terrain profile corresponding to each pair of mutually perpendicular axes of symmetry;
[0063] If the product of at least one parameter is less than 0, then the local geographic region corresponding to the initial grid is determined to be a saddle region.
[0064] If the product of the two parameters is greater than 0, or the product of the two parameters is equal to 0, then the local geographic region corresponding to the initial grid is determined not to be a saddle region.
[0065] In some embodiments of this disclosure, the fourth determining module is further configured to:
[0066] Determine the saddle degree information corresponding to each saddle region, whereby the saddle degree information is used to describe the topographic relief of the saddle region.
[0067] In some embodiments of this disclosure, the fourth determining module is further configured to:
[0068] If the parameter product is less than 0, then a pair of target terrain profiles corresponding to the parameter product are determined;
[0069] Based on a pair of target terrain profiles, determine the target saddle degree value corresponding to each saddle region;
[0070] The target saddle depth value is used as the saddle depth information.
[0071] In some embodiments of this disclosure, the fourth determining module is further configured to:
[0072] Determine the first average slope of each target terrain profile line in a pair of target terrain profile lines, the first elevation value of the extreme point of the target curve in each target terrain profile line, and the horizontal distance between the extreme point of the target curve and the center point in each target terrain profile line;
[0073] Determine the maximum and minimum elevation values of a pair of target terrain profiles, and the second average slope of the first terrain profile in each target terrain profile, wherein the first terrain profile is a portion of the target terrain profile between the extreme point and the center point of the target curve in the target terrain profile.
[0074] Determine the second elevation value of the central grid;
[0075] The candidate saddle depth value is determined by multiplying the first average slope, the second average slope, the maximum elevation value, the minimum elevation value, the horizontal distance, the first elevation value, the second elevation value, and the parameters corresponding to a pair of target terrain profile lines.
[0076] Based on the candidate saddle severity values, determine the target saddle severity value corresponding to each saddle region.
[0077] In some embodiments of this disclosure, the fourth determining module is further configured to:
[0078] The candidate saddle point severity value is calculated using the following formula:
[0079] ;
[0080] Where D is the candidate saddle point severity value, and T is the parameter product. , This represents the horizontal distance between a pair of target terrain profile lines. This is the second elevation value. , These are the first elevation values corresponding to a pair of target terrain profile lines. and Let be the first average slope corresponding to a pair of target terrain profile lines. For the maximum elevation value, The minimum elevation value, = + , and These are the second average slopes corresponding to each of the first terrain profiles in a pair of target terrain profiles.
[0081] In some embodiments of this disclosure, the fourth determining module is further configured to:
[0082] If the number of candidate saddle severity values is 1, then the candidate saddle severity value is used as the target saddle severity value.
[0083] If the number of candidate saddle severity values is 2, then the candidate saddle severity value with the larger value among the two candidate saddle severity values is taken as the target saddle severity value.
[0084] The electronic device proposed in the third aspect of this disclosure includes: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the method for determining the saddle region as proposed in the first aspect of this disclosure.
[0085] A fourth aspect of this disclosure provides a non-transitory computer-readable storage medium having a computer program stored thereon that, when executed by a processor, implements the method for determining the saddle region as described in the first aspect of this disclosure.
[0086] A fifth aspect of this disclosure provides a computer program product in which, when instructions in the computer program product are executed by a processor, the method for determining a saddle region as described in a first aspect of this disclosure is performed.
[0087] The method, apparatus, electronic device, and storage medium for determining the saddle region provided in this disclosure have at least the following beneficial effects: They acquire a digital elevation model of the geographic region, wherein the digital elevation model includes: multiple initial grids, and then, with each initial grid as the center grid, determine N... N initial grids are used as the target grid region, where N is an odd number greater than 1. Four initial terrain profile lines are then determined along the four axes of symmetry of the target grid region. The curve parameters of each initial terrain profile line are then determined, which describe the slope change of the initial terrain profile line when it passes through the central grid. Based on the curve parameters of each initial terrain profile line, it is determined whether the local geographic area corresponding to the initial grid is a saddle region. This allows for accurate determination of whether a geographic area is a saddle region, effectively improving the determination of saddle regions.
[0088] Additional aspects and advantages of this disclosure will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of this disclosure. Attached Figure Description
[0089] The above and / or additional aspects and advantages of this disclosure will become apparent and readily understood from the following description of the embodiments taken in conjunction with the accompanying drawings, in which:
[0090] Figure 1 This is a flowchart illustrating a method for determining the saddle region according to an embodiment of this disclosure;
[0091] Figure 2A This is a schematic diagram of the saddle region proposed in an embodiment of this disclosure;
[0092] Figure 2B This is a schematic diagram of a target grid region proposed in an embodiment of this disclosure;
[0093] Figure 3 This is a flowchart illustrating a method for determining the saddle region according to another embodiment of this disclosure;
[0094] Figure 4 This is a schematic diagram of an initial terrain profile line proposed in an embodiment of this disclosure;
[0095] Figure 5 This is a schematic diagram of the structure of a saddle region determination device according to an embodiment of the present disclosure;
[0096] Figure 6 A block diagram of an exemplary electronic device suitable for implementing embodiments of the present disclosure is shown. Detailed Implementation
[0097] Embodiments of this disclosure are described in detail below, with examples of embodiments illustrated in the accompanying drawings, wherein the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the accompanying drawings are exemplary and are used only to explain this disclosure, and should not be construed as limiting this disclosure. Rather, embodiments of this disclosure include all variations, modifications, and equivalents falling within the spirit and scope of the appended claims.
[0098] Figure 1 This is a schematic flowchart of a method for determining the saddle region according to an embodiment of this disclosure.
[0099] It should be noted that the execution subject of the method for determining the saddle region in this embodiment is the device for determining the saddle region. This device can be implemented by software and / or hardware. This device can be configured in an electronic device, which may include, but is not limited to, a terminal, a server, etc. For example, the terminal may be a mobile phone, a PDA, etc.
[0100] like Figure 1 As shown, the method for determining the saddle region includes:
[0101] S101: Obtain the digital elevation model of the geographic region, wherein the digital elevation model includes: multiple initial grids.
[0102] In this embodiment of the disclosure, firstly, in conjunction with Figure 2A For a detailed explanation of the saddle area, see [link to relevant documentation]. Figure 2A , Figure 2AThis is a schematic diagram of a saddle region proposed in an embodiment of the present disclosure. The saddle region is the middle part of two connected mountain ranges. The relative height of the saddle region is relatively high and the sides are relatively steep. When the wind passes over the saddle of the mountain ridge, it will form a saddle-shaped vortex.
[0103] The geographical region can be a pre-defined geographical region, or it can be any geographical region that needs to be determined to be a saddle region; there are no restrictions on this.
[0104] Among them, the Digital Elevation Model (DEM) of a geographic region is a commonly used data model in geographic information systems. It is used to describe the elevation information of the earth's surface, that is, to realize the digital simulation of the ground terrain (i.e., the digital expression of the terrain surface morphology) through limited terrain elevation data. The digital elevation model can be either a rectangular grid or an irregular triangular grid. Correspondingly, the grids that make up the rectangular grid or the triangular grid are the initial grids.
[0105] In this embodiment of the disclosure, a digital elevation model will be used as a rectangular grid to illustrate the subsequent embodiments of the disclosure, but this is not intended to be limiting.
[0106] In this embodiment of the disclosure, the digital elevation model of a geographic area can be obtained by measuring the surface elevation data of the geographic area through various methods such as radar and laser scanning. Alternatively, it can be collected from existing topographic maps, such as grid reading method, manual tracking with a digitizer, and semi-automatic collection by a scanner, and then generating the digital elevation model of the geographic area through interpolation. There are no limitations on this.
[0107] S102: Determine N using each initial grid as the center grid. N initial grids are used as the target grid region, where N is an odd number greater than 1.
[0108] In this embodiment of the disclosure, after obtaining a digital elevation model of a geographic area, the digital elevation model includes: multiple initial grids, and N can be determined with each initial grid as the center grid. N initial grids are used as the target grid region, where N is an odd number greater than 1.
[0109] In this embodiment of the disclosure, the initial grid can be, for example, a square grid. Accordingly, N can be determined with each initial grid as the center grid. N initial grids form a square target grid region.
[0110] For example, see Figure 2B , Figure 2BThis is a schematic diagram of a target mesh region proposed in an embodiment of this disclosure, which can be determined by taking the initial mesh as the center mesh and defining 7 The 7 initial grids consist of, as follows Figure 2B The target grid area is shown. Then, the elevation data of the target grid area can be combined to accurately determine whether the local geographic area corresponding to the central grid is a saddle area. For details, please refer to the subsequent embodiments, which will not be repeated here.
[0111] S103: Determine the four initial terrain profile lines along the four axes of symmetry of the target grid region.
[0112] A topographic profile line is a line segment or curve used to show changes in surface elevation, which moves from one location to another along a selected path or specific route.
[0113] In this embodiment of the disclosure, see the above. Figure 2B This can be done by determining four initial terrain profile lines along the four axes of symmetry of the target grid region, which means determining the elevation data of each point on the horizontal plane along the four axes of symmetry of the target grid region, and generating four initial terrain profile lines based on the elevation data of each point on the four axes of symmetry.
[0114] S104: Determine the curve parameters for each initial terrain profile line, where the curve parameters describe the slope change of the initial terrain profile line as it passes through the central grid.
[0115] The curve parameter describes the slope change of the initial terrain profile line as it passes through the central grid.
[0116] In some embodiments, determining the curve parameters of each initial terrain profile line can be achieved by determining the curve equation corresponding to each initial terrain profile line, then differentiating the curve equation to determine the derivative of the curve equation as the curve parameter of the initial terrain profile line. Alternatively, considering that the initial terrain profile line may have many undulations, determining the curve parameters of each initial terrain profile line can also be achieved by considering the changing trend of the curves closer to the center point of the central grid. That is, it can be achieved by separately determining the slope of the local initial terrain profile lines corresponding to the two grids on the left and right sides of the central grid, and then using the slopes of the local initial terrain profile lines corresponding to the two grids together as the curve parameters of the initial terrain profile line. There are no restrictions on this.
[0117] In some embodiments, the curve parameters of each initial terrain profile line can be determined as follows: when the slope of the initial terrain profile line changes from positive to negative as it passes through the central grid, the curve parameter is determined to be positive; when the slope of the initial terrain profile line changes from negative to positive as it passes through the central grid, the curve parameter is determined to be negative; and when there is no change in the slope of the initial terrain profile line as it passes through the central grid, the curve parameter is determined to be 0.
[0118] For example, the curve parameter of the initial terrain profile can be set to 1 when the slope of the initial terrain profile changes from positive to negative as it passes through the central grid; -1 when the slope of the initial terrain profile changes from negative to positive as it passes through the central grid; and 0 when the slope of the initial terrain profile does not change from positive to negative as it passes through the central grid. There are no restrictions on this.
[0119] Optionally, in some embodiments, determining the curve parameters of each initial terrain profile line can be done by setting the curve parameters of the initial terrain profile line to 0 when there are no extreme points of the initial curve in the initial terrain profile line.
[0120] S105: Based on the curve parameters of each initial topographic profile line, determine whether the local geographic region corresponding to the initial grid is a saddle region.
[0121] In this embodiment of the disclosure, after determining the curve parameters of each initial terrain profile line, it is possible to determine whether the local geographic region corresponding to the initial grid is a saddle region based on the curve parameters of each initial terrain profile line.
[0122] In some embodiments, determining whether the local geographic region corresponding to the initial grid is a saddle region based on the curve parameters of each initial topographic profile line can be achieved by inputting the curve parameters of each initial topographic profile line into a pre-trained saddle region determination model (wherein, the saddle region determination model has been pre-learned the mapping relationship between the curve parameters of each initial topographic profile line of different initial grids and the saddle region determination results) to obtain the saddle region determination result output by the saddle region determination model (wherein, the saddle region determination result is used to indicate whether the local geographic region corresponding to the initial grid is a saddle region), and there are no restrictions on this.
[0123] In this embodiment of the disclosure, a digital elevation model of a geographical region is obtained. The digital elevation model includes multiple initial grids, and N is determined using each initial grid as the center grid. N initial grids are used as the target grid region, where N is an odd number greater than 1. Four initial terrain profile lines are then determined along the four axes of symmetry of the target grid region. The curve parameters of each initial terrain profile line are then determined, which describe the slope change of the initial terrain profile line when it passes through the central grid. Based on the curve parameters of each initial terrain profile line, it is determined whether the local geographic area corresponding to the initial grid is a saddle region. This allows for accurate determination of whether a geographic area is a saddle region, effectively improving the determination of saddle regions.
[0124] Figure 3 This is a flowchart illustrating a method for determining the saddle region according to another embodiment of this disclosure.
[0125] like Figure 3 As shown, the method for determining the saddle region includes:
[0126] S301: Obtain the digital elevation model of the geographic region, wherein the digital elevation model includes: multiple initial grids.
[0127] S302: Determine N using each initial grid as the center grid. N initial grids are used as the target grid region, where N is an odd number greater than 1.
[0128] S303: Determine four initial terrain profile lines along the four axes of symmetry of the target grid region.
[0129] For detailed descriptions of S301-S303, please refer to the above embodiments, which will not be repeated here.
[0130] S304: If there are initial curve extrema in the initial terrain profile, then determine the target curve extrema that is closest to the center point of the central grid from at least one initial curve extrema.
[0131] In this embodiment of the disclosure, see Figure 4 , Figure 4 This is a schematic diagram of an initial terrain profile line proposed in an embodiment of the present disclosure. Specifically, when it is determined that there are initial curve extreme points (maximum points and / or minimum points) in the initial terrain profile line, the target curve extreme point (point b) that is closest to the center point (point a) of the central grid is determined from at least one initial curve extreme point.
[0132] In this embodiment of the disclosure, determining the target curve extreme point closest to the center point of the central grid from at least one initial curve extreme point can be done by determining the curve extreme point closest to the center point in Euclidean distance as the target curve extreme point. When there are two or more curve extreme points with the same Euclidean distance to the center point, the target curve extreme point can be further determined from the curve extreme points with the same Euclidean distance to the center point in the horizontal direction as the target curve extreme point. When there are two or more curve extreme points with the same horizontal distance to the center point in the horizontal direction, the target curve extreme point can be further selected from the curve extreme points with the same horizontal distance to the center point in the horizontal direction as the curve extreme point with the same vertical distance to the center point in the vertical direction as the target curve extreme point. There are no restrictions on this.
[0133] S305: Extract monotonically changing portions of the initial terrain profile from the initial terrain profile lines on both sides of the extreme points of the target curve.
[0134] In this embodiment of the disclosure, after determining the target curve extreme point that is closest to the center point of the central grid from at least one initial curve extreme point, it may be possible to extract monotonically changing portions of the initial terrain profile line on both sides of the target curve extreme point from the initial terrain profile line.
[0135] See above Figure 4 One approach is to determine an extreme point of the curve (point c and point d) on both sides of point b in the initial terrain profile, and then extract the portion of the initial terrain profile between point c and point b, and the portion of the initial terrain profile between point d and point b.
[0136] S306: The initial terrain profile lines that change monotonically on both sides of the extreme point of the target curve are used together as the target terrain profile line.
[0137] In this embodiment of the disclosure, see the above. Figure 4 The initial terrain profile can be extracted from the initial terrain profile, and the monotonically changing portions of the initial terrain profile on both sides of the extreme point of the target curve can be taken as the portion of the initial terrain profile between point c and point d.
[0138] S307: Determine the curve parameters of the initial terrain profile line based on the target terrain profile line.
[0139] In this embodiment of the disclosure, after determining the target terrain profile, the curve parameters of the initial terrain profile can be determined based on the target terrain profile.
[0140] Optionally, in some embodiments, the curve parameters of the initial terrain profile are determined based on the target terrain profile. This can be done by determining the curve parameters as negative if the extreme point of the target curve in the target terrain profile is a minimum point, and as positive if the extreme point of the target curve in the target terrain profile is a maximum point.
[0141] For example, in the embodiments of this disclosure, see the above. Figure 4 If the extreme point (point b) of the target curve in the target terrain profile is a minimum point, then the curve parameter is determined to be a negative number (e.g., -1). If the extreme point (point b) of the target curve in the target terrain profile is a maximum point, then the curve parameter is determined to be a positive number (e.g., 1). There are no restrictions on this.
[0142] S308: Based on the curve parameters of each initial topographic profile line, determine whether the local geographic region corresponding to the initial grid is a saddle region.
[0143] Optionally, in some embodiments, determining whether the local geographic region corresponding to the initial grid is a saddle region based on the curve parameters of each initial topographic profile line can be achieved by determining the product of the curve parameters of the initial topographic profile lines corresponding to each pair of mutually perpendicular axes of symmetry, and determining that the local geographic region corresponding to the initial grid is a saddle region when at least one parameter product is less than 0, or determining that the local geographic region corresponding to the initial grid is not a saddle region when both parameter products are greater than 0 or both parameter products are equal to 0.
[0144] In other words, in this embodiment of the present disclosure, the parameter product between the curve parameters of the initial terrain profile line corresponding to each pair of mutually perpendicular axes of symmetry can be determined. When it is determined that at least one parameter product is less than 0, the local geographic region corresponding to the initial grid is determined to be a saddle region. When it is determined that two parameter products are both greater than 0, or two parameter products are both equal to 0, the local geographic region corresponding to the initial grid is determined not to be a saddle region.
[0145] S309: Determine the saddle degree information corresponding to each saddle region, wherein the saddle degree information is used to describe the topographic relief of the saddle region.
[0146] Specifically, when the terrain undulations in a saddle area are significant, the changes in the saddle-shaped eddies are complex, leading to highly complex forest fire patterns and making firefighting in that area more difficult. Conversely, when the terrain undulations in a saddle area are relatively small, the changes in the saddle-shaped eddies are clearer, making firefighting in that area relatively easier. Therefore, pre-determining the saddle degree information corresponding to each saddle area allows for advance assessment of the work difficulty in forest fire fighting or other operations, enabling the preparation of corresponding work plans in advance and facilitating the implementation of forest fire fighting and other related tasks.
[0147] Optionally, in some embodiments, determining the saddle degree information corresponding to each saddle region may involve determining a pair of target terrain profiles corresponding to the parameter product when the parameter product is less than 0, then determining the target saddle degree value corresponding to each saddle region based on the pair of target terrain profiles, and then using the target saddle degree value as the saddle degree information. Thus, the terrain undulation amplitude of the saddle region can be quantitatively described based on the target saddle degree value, thereby effectively improving the reference value of the saddle degree information.
[0148] Among them, the target saddle degree value quantifies the topographic relief of the saddle area. The larger the target saddle degree value, the greater the topographic relief of the saddle area, and vice versa.
[0149] In other words, in this embodiment of the present disclosure, when the parameter product is less than 0, a pair of target terrain profile lines corresponding to the parameter product are determined, and then the target saddle degree value corresponding to each saddle region is determined based on the pair of target terrain profile lines.
[0150] Optionally, in some embodiments, determining the target saddle degree value corresponding to each saddle region based on a pair of target terrain profiles may involve determining the first average slope of each target terrain profile in the pair, the first elevation value of the extreme point of the target curve in each target terrain profile, and the horizontal distance between the extreme point of the target curve and the center point in each target terrain profile; determining the maximum elevation value and the minimum elevation value of the pair of target terrain profiles; determining the second average slope of the first terrain profile in each target terrain profile; determining the second elevation value of the central grid; determining candidate saddle degree values based on the product of the first average slope, the second average slope, the maximum elevation value, the minimum elevation value, the horizontal distance, the first elevation value, the second elevation value, and the parameters corresponding to the pair of target terrain profiles; and then determining the target saddle degree value corresponding to each saddle region based on the candidate saddle degree values.
[0151] See above. Figure 4The horizontal distances between the extreme points and the center point of the target curve corresponding to each of the two target terrain profiles are d1 and d2, respectively. The first elevation values of the extreme points of the target curve corresponding to each of the two target terrain profiles are respectively... , .
[0152] Among them, the second elevation value of the central grid The second elevation value can be the elevation value at the center point of the central grid, or it can be the average elevation value of the central grid, without any restrictions.
[0153] Among them, the maximum elevation value and minimum elevation value It is the unique maximum elevation value and the unique minimum elevation value determined by the joint determination of a target terrain profile.
[0154] Among them, the average slope of each target terrain profile in a pair of target terrain profiles is respectively and .
[0155] Among them, the second average slope corresponding to the first terrain profile line in each target terrain profile line is respectively and The first terrain profile is the portion of the target terrain profile between the extreme points and the center point of the target curve.
[0156] Optionally, in some embodiments, the candidate saddle depth value is determined based on the product of a first average slope, a second average slope, a maximum elevation value, a minimum elevation value, a horizontal distance, a first elevation value, a second elevation value, and the parameters corresponding to a pair of target terrain profile lines. The candidate saddle depth value can be calculated using the following formula:
[0157] ;
[0158] Where D is the candidate saddle point severity value, and T is the parameter product. , This represents the horizontal distance between a pair of target terrain profile lines. This is the second elevation value. , These are the first elevation values corresponding to a pair of target terrain profile lines. and Let be the first average slope corresponding to a pair of target terrain profile lines. For the maximum elevation value, The minimum elevation value, = + , and These are the second average slopes corresponding to each of the first terrain profiles in a pair of target terrain profiles.
[0159] Optionally, in some embodiments, the target saddle degree value corresponding to each saddle region is determined based on the candidate saddle degree value. This can be done by using the candidate saddle degree value as the target saddle degree value when the number of candidate saddle degree values is 1, or by using the candidate saddle degree value with the largest value among the two candidate saddle degree values as the target saddle degree value when the number of candidate saddle degree values is 2.
[0160] In this embodiment of the disclosure, after calculating the candidate saddle degree value, when the number of candidate saddle degree values is 1, the candidate saddle degree value can be used as the target saddle degree value, or when the number of candidate saddle degree values is 2, the candidate saddle degree value with the largest value among the two candidate saddle degree values can be used as the target saddle degree value. Thus, the unique target saddle degree value corresponding to each saddle region can be accurately determined.
[0161] In this embodiment of the disclosure, a digital elevation model of a geographical region is obtained. The digital elevation model includes multiple initial grids, and N is determined using each initial grid as the center grid. N initial grids are used as the target grid region, where N is an odd number greater than 1. Four initial terrain profile lines are then determined along the four axes of symmetry of the target grid region. If an initial curve extremum exists in one of the initial terrain profile lines, the target curve extremum point closest to the center point of the central grid is determined from at least one of the initial curve extremum points. Monotonically changing portions of the initial terrain profile lines on both sides of the target curve extremum point are then extracted from the initial terrain profile lines. These monotonically changing portions of the initial terrain profile lines on both sides of the target curve extremum point are combined to form the target terrain profile line. Finally, based on the target... The initial topographic profile is determined by defining the curve parameters of each initial topographic profile. Based on these parameters, it is then determined whether the local geographic area corresponding to the initial grid is a saddle region. Finally, the saddle degree information corresponding to each saddle region is determined. This saddle degree information describes the topographic relief of the saddle region. By pre-determining the saddle degree information for each saddle region, the difficulty of working in that saddle region can be known in advance during forest fire fighting or other operations. This allows for the preparation of corresponding work plans in advance, thereby facilitating the implementation of forest fire fighting and other related work.
[0162] Figure 5 This is a schematic diagram of the structure of a saddle region determination device according to an embodiment of this disclosure.
[0163] like Figure 5 As shown, the device 50 for determining the saddle region includes:
[0164] The acquisition module 501 is used to acquire the digital elevation model of a geographic area, wherein the digital elevation model includes: multiple initial grids;
[0165] The first determining module 502 is used to determine N, with each initial grid as the center grid. N initial grids are used as the target grid region, where N is an odd number greater than 1;
[0166] The second determining module 503 is used to determine four initial terrain profile lines along the four axes of symmetry of the target grid region, respectively.
[0167] The third determining module 504 is used to determine the curve parameters of each initial terrain profile line, wherein the curve parameters are used to describe the slope change of the initial terrain profile line when it passes through the central grid.
[0168] The fourth determination module 505 is used to determine whether the local geographic region corresponding to the initial grid is a saddle region based on the curve parameters of each initial terrain profile line.
[0169] In some embodiments of this disclosure, the third determining module 504 is further configured to:
[0170] If there are no extreme points of the initial curve in the initial terrain profile, then the curve parameter corresponding to the initial terrain profile is determined to be 0.
[0171] In some embodiments of this disclosure, the third determining module 504 is further configured to:
[0172] If there are initial curve extrema in the initial terrain profile, then determine the target curve extrema that is closest to the center point of the central grid from at least one initial curve extrema.
[0173] Extract monotonically changing portions of the initial terrain profile from both sides of the extreme points of the target curve from the initial terrain profile.
[0174] The initial terrain profile lines that change monotonically on both sides of the extreme point of the target curve are used together as the target terrain profile line.
[0175] Determine the curve parameters of the initial terrain profile based on the target terrain profile.
[0176] In some embodiments of this disclosure, the third determining module 504 is further configured to:
[0177] If the extreme point of the target curve in the target terrain profile is a minimum point, then the curve parameter is determined to be negative.
[0178] If the extreme point of the target curve in the target terrain profile is a maximum point, then the curve parameter is determined to be a positive number.
[0179] In some embodiments of this disclosure, the fourth determining module 505 is further configured to:
[0180] Determine the parametric product between the curve parameters of the initial terrain profile corresponding to each pair of mutually perpendicular axes of symmetry;
[0181] If the product of at least one parameter is less than 0, then the local geographic region corresponding to the initial grid is determined to be a saddle region.
[0182] If the product of the two parameters is greater than 0, or the product of the two parameters is equal to 0, then the local geographic region corresponding to the initial grid is determined not to be a saddle region.
[0183] In some embodiments of this disclosure, the fourth determining module 505 is further configured to:
[0184] Determine the saddle degree information corresponding to each saddle region, whereby the saddle degree information is used to describe the topographic relief of the saddle region.
[0185] In some embodiments of this disclosure, the fourth determining module 505 is further configured to:
[0186] If the parameter product is less than 0, then a pair of target terrain profiles corresponding to the parameter product are determined;
[0187] Based on a pair of target terrain profiles, determine the target saddle degree value corresponding to each saddle region;
[0188] The target saddle depth value is used as the saddle depth information.
[0189] In some embodiments of this disclosure, the fourth determining module 505 is further configured to:
[0190] Determine the first average slope of each target terrain profile line in a pair of target terrain profile lines, the first elevation value of the extreme point of the target curve in each target terrain profile line, and the horizontal distance between the extreme point of the target curve and the center point in each target terrain profile line;
[0191] Determine the maximum and minimum elevation values of a pair of target terrain profiles, and the second average slope of the first terrain profile in each target terrain profile, wherein the first terrain profile is a portion of the target terrain profile between the extreme point and the center point of the target curve in the target terrain profile.
[0192] Determine the second elevation value of the central grid;
[0193] The candidate saddle depth value is determined by multiplying the first average slope, the second average slope, the maximum elevation value, the minimum elevation value, the horizontal distance, the first elevation value, the second elevation value, and the parameters corresponding to a pair of target terrain profile lines.
[0194] Based on the candidate saddle severity values, determine the target saddle severity value corresponding to each saddle region.
[0195] In some embodiments of this disclosure, the fourth determining module 505 is further configured to:
[0196] The candidate saddle point severity value is calculated using the following formula:
[0197] ;
[0198] Where D is the candidate saddle point severity value, and T is the parameter product. , This represents the horizontal distance between a pair of target terrain profile lines. This is the second elevation value. , These are the first elevation values corresponding to a pair of target terrain profile lines. and Let be the first average slope corresponding to a pair of target terrain profile lines. For the maximum elevation value, The minimum elevation value, = + , and These are the second average slopes corresponding to each of the first terrain profiles in a pair of target terrain profiles.
[0199] In some embodiments of this disclosure, the fourth determining module 505 is further configured to:
[0200] If the number of candidate saddle severity values is 1, then the candidate saddle severity value is used as the target saddle severity value.
[0201] If the number of candidate saddle severity values is 2, then the candidate saddle severity value with the larger value among the two candidate saddle severity values is taken as the target saddle severity value.
[0202] It should be noted that the aforementioned explanation of the method for determining the saddle region also applies to the device for determining the saddle region in this embodiment, and will not be repeated here.
[0203] In this embodiment of the disclosure, a digital elevation model of a geographical region is obtained. The digital elevation model includes multiple initial grids, and N is determined using each initial grid as the center grid. N initial grids are used as the target grid region, where N is an odd number greater than 1. Four initial terrain profile lines are then determined along the four axes of symmetry of the target grid region. The curve parameters of each initial terrain profile line are then determined, which describe the slope change of the initial terrain profile line when it passes through the central grid. Based on the curve parameters of each initial terrain profile line, it is determined whether the local geographic area corresponding to the initial grid is a saddle region. This allows for accurate determination of whether a geographic area is a saddle region, effectively improving the determination of saddle regions.
[0204] Figure 6 A block diagram of an exemplary electronic device suitable for implementing embodiments of the present disclosure is shown. Figure 6 The electronic device 12 shown is merely an example and should not impose any limitation on the functionality and scope of use of the embodiments disclosed herein.
[0205] like Figure 6 As shown, the electronic device 12 is represented in the form of a general-purpose computing device. The components of the electronic device 12 may include, but are not limited to: one or more processors or processing units 16, system memory 28, and bus 18 connecting different system components (including system memory 28 and processing unit 16).
[0206] Bus 18 represents one or more of several bus architectures, including a memory bus or memory controller, a peripheral bus, a graphics acceleration port, a processor, or a local bus using any of the various bus architectures. Examples of these architectures include, but are not limited to, the Industry Standard Architecture (ISA) bus, the Micro Channel Architecture (MAC) bus, the Enhanced ISA bus, the Video Electronics Standards Association (VESA) local bus, and the Peripheral Component Interconnect (PCI) bus.
[0207] Electronic device 12 typically includes a variety of computer system readable media. These media can be any available media that can be accessed by electronic device 12, including volatile and non-volatile media, removable and non-removable media.
[0208] Memory 28 may include computer system readable media in the form of volatile memory, such as Random Access Memory (RAM) 30 and / or cache memory 32. Electronic device 12 may further include other removable / non-removable, volatile / non-volatile computer system storage media. By way of example only, storage system 34 may be used to read and write non-removable, non-volatile magnetic media (… Figure 6 Not shown; usually referred to as a "hard drive".
[0209] although Figure 6 As not shown, a disk drive for reading and writing to a removable non-volatile disk (e.g., a "floppy disk") and an optical disc drive for reading and writing to a removable non-volatile optical disc (e.g., a compact disc read-only memory (CD-ROM), a digital video disc read-only memory (DVD-ROM), or other optical media) may be provided. In these cases, each drive may be connected to bus 18 via one or more data media interfaces. Memory 28 may include at least one program product having a set (e.g., at least one) of program modules configured to perform the functions of the embodiments of this disclosure.
[0210] A program / utility 40 having a set (at least one) of program modules 42 may be stored, for example, in memory 28. Such program modules 42 include, but are not limited to, an operating system, one or more application programs, other program modules, and program data. Each or some combination of these examples may include an implementation of a network environment. Program modules 42 typically perform the functions and / or methods described in the embodiments of this disclosure.
[0211] Electronic device 12 can also communicate with one or more external devices 14 (e.g., keyboard, pointing device, display 24, etc.), and with one or more devices that enable human interaction with electronic device 12, and / or with any device that enables electronic device 12 to communicate with one or more other computing devices (e.g., network card, modem, etc.). This communication can be performed via input / output (I / O) interface 22. Furthermore, electronic device 12 can also communicate with one or more networks (e.g., local area network (LAN), wide area network (WAN), and / or public networks, such as the Internet) via network adapter 20. As shown, network adapter 20 communicates with other modules of electronic device 12 via bus 18. It should be understood that, although not shown in the figures, other hardware and / or software modules can be used in conjunction with electronic device 12, including but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, and data backup storage systems.
[0212] The processing unit 16 executes various functional applications and data processing by running programs stored in the system memory 28, such as implementing the method for determining the saddle region mentioned in the foregoing embodiments.
[0213] To implement the above embodiments, this disclosure also proposes a non-transitory computer-readable storage medium storing a computer program that, when executed by a processor, implements the method for determining the saddle region as proposed in the foregoing embodiments of this disclosure.
[0214] To implement the above embodiments, this disclosure also proposes a computer program product that, when the instruction processor in the computer program product is executed, performs the method for determining the saddle region as proposed in the foregoing embodiments of this disclosure.
[0215] Other embodiments of this disclosure will readily occur to those skilled in the art upon consideration of the specification and practice of the invention disclosed herein. This disclosure is intended to cover any variations, uses, or adaptations of this disclosure that follow the general principles of this disclosure and include common knowledge or customary techniques in the art not disclosed herein. The specification and examples are to be considered exemplary only, and the true scope and spirit of this disclosure are indicated by the following claims.
[0216] It should be understood that this disclosure is not limited to the precise structures described above and shown in the accompanying drawings, and various modifications and changes can be made without departing from its scope. The scope of this disclosure is limited only by the appended claims.
[0217] It should be noted that in the description of this disclosure, the terms "first," "second," etc., are used for descriptive purposes only and should not be construed as indicating or implying relative importance. Furthermore, in the description of this disclosure, unless otherwise stated, "a plurality of" means two or more.
[0218] Any process or method description in the flowchart or otherwise herein can be understood as representing a module, segment, or portion of code comprising one or more executable instructions for implementing a particular logical function or process, and the scope of preferred embodiments of this disclosure includes additional implementations in which functions may be performed not in the order shown or discussed, including substantially simultaneously or in reverse order depending on the function involved, as will be understood by those skilled in the art to which embodiments of this disclosure pertain.
[0219] It should be understood that various parts of this disclosure can be implemented using hardware, software, firmware, or a combination thereof. In the above embodiments, multiple steps or methods can be implemented using software or firmware stored in memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, it can be implemented using any one or a combination of the following techniques known in the art: discrete logic circuits having logic gates for implementing logical functions on data signals, application-specific integrated circuits (ASICs) having suitable combinational logic gates, programmable gate arrays (PGAs), field-programmable gate arrays (FPGAs), etc.
[0220] Those skilled in the art will understand that all or part of the steps of the methods in the above embodiments can be implemented by a program instructing related hardware. The program can be stored in a computer-readable storage medium, and when executed, the program includes one or a combination of the steps of the method embodiments.
[0221] Furthermore, the functional units in the various embodiments of this disclosure can be integrated into a processing module, or each unit can exist physically separately, or two or more units can be integrated into a module. The integrated module can be implemented in hardware or as a software functional module. If the integrated module is implemented as a software functional module and sold or used as an independent product, it can also be stored in a computer-readable storage medium.
[0222] The storage media mentioned above can be read-only memory, disk, or optical disk, etc.
[0223] In the description of this specification, references to terms such as "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., indicate that a specific feature, structure, material, or characteristic described in connection with that embodiment or example is included in at least one embodiment or example of this disclosure. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples.
[0224] Although embodiments of the present disclosure have been shown and described above, it is to be understood that the above embodiments are exemplary and should not be construed as limiting the present disclosure. Those skilled in the art can make changes, modifications, substitutions and variations to the above embodiments within the scope of the present disclosure.
Claims
1. A method of determining a saddle region, characterized in that, The method includes: Obtain a digital elevation model of a geographic region, wherein the digital elevation model includes: multiple initial grids; Determine N using each of the initial grids as the center grid. N initial grids are used as the target grid region, where N is an odd number greater than 1; Four initial terrain profile lines are determined along the four axes of symmetry of the target grid region; Determine the curve parameters for each of the initial terrain profile lines, wherein the curve parameters are used to describe the slope change of the initial terrain profile line as it passes through the central grid; Based on the curve parameters of each initial terrain profile line, determine whether the local geographic region corresponding to the initial grid is a saddle region; Determine saddle degree information corresponding to each saddle region, wherein the saddle degree information is used to describe the topographic relief amplitude of the saddle region; The method for obtaining the saddle depth information includes: Determine the parametric product between the curve parameters of the initial terrain profile lines corresponding to each pair of mutually perpendicular axes of symmetry; If the product of the parameters is less than 0, then a pair of target terrain profiles corresponding to the product of the parameters are determined; Determine the first average slope of each of the target terrain profiles in the pair of target terrain profiles, the first elevation value of the extreme point of the target curve in each of the target terrain profiles, and the horizontal distance between the extreme point of the target curve and the center point in each of the target terrain profiles; Determine the maximum elevation value, minimum elevation value, and second average slope of the first terrain profile in each of the target terrain profiles, wherein the first terrain profile is a portion of the target terrain profile between the extreme point of the target curve and the center point in the target terrain profile; Determine the second elevation value of the central grid; The candidate saddle point severity value is calculated using the following formula: ; Where D is the candidate saddle point degree value, and T is the product of the parameters. , The horizontal distances corresponding to a pair of target terrain profile lines are respectively. This is the second elevation value. , Let be the first elevation value corresponding to each pair of the target terrain profile lines. and Let be the first average slope corresponding to each pair of the target terrain profile lines. The maximum elevation value, The minimum elevation value, = + , and These are the second average slopes corresponding to each first terrain profile line in a pair of target terrain profile lines; Based on the candidate saddle severity values, a target saddle severity value corresponding to each saddle region is determined; The target saddle degree value is used as the saddle degree information.
2. The method as described in claim 1, characterized in that, The determination of the curve parameters for each of the initial terrain profile lines includes: If there is no extreme point of the initial curve in the initial terrain profile, then the curve parameter corresponding to the initial terrain profile is determined to be 0.
3. The method as described in claim 1, characterized in that, The determination of the curve parameters for each of the initial terrain profile lines further includes: If there are initial curve extreme points in the initial terrain profile, then determine the target curve extreme point that is closest to the center point of the central grid from at least one of the initial curve extreme points; Extract monotonically changing portions of the initial terrain profile from both sides of the extreme point of the target curve from the initial terrain profile; The initial terrain profile lines that change monotonically on both sides of the extreme point of the target curve are used together as the target terrain profile line. Based on the target terrain profile, determine the curve parameters of the initial terrain profile.
4. The method as described in claim 3, characterized in that, The step of determining the curve parameters of the initial terrain profile based on the target terrain profile includes: If the extreme point of the target curve in the target terrain profile is a minimum point, then the curve parameter is determined to be negative. If the extreme point of the target curve in the target terrain profile is a maximum point, then the curve parameter is determined to be a positive number.
5. The method as described in claim 4, characterized in that, The step of determining whether the local geographic region corresponding to the initial grid is a saddle region based on the curve parameters of each initial terrain profile line includes: Determine the parametric product between the curve parameters of the initial terrain profile lines corresponding to each pair of mutually perpendicular axes of symmetry; If the product of at least one of the parameters is less than 0, then the local geographic region corresponding to the initial grid is determined to be the saddle region; If the product of the two parameters is greater than 0, or the product of the two parameters is equal to 0, then the local geographic region corresponding to the initial grid is determined not to be the saddle region.
6. The method as described in claim 1, characterized in that, The step of determining the target saddle severity value corresponding to each saddle region based on the candidate saddle severity values includes: If the number of candidate saddle severity values is 1, then the candidate saddle severity value is used as the target saddle severity value; If the number of candidate saddle severity values is 2, then the candidate saddle severity value with the largest value among the two candidate saddle severity values is taken as the target saddle severity value.
7. A device for determining a saddle region, characterized in that, The device includes: An acquisition module is used to acquire a digital elevation model of a geographic area, wherein the digital elevation model includes: multiple initial grids; The first determining module is used to determine N, with each of the initial grids as the center grid. N initial grids are used as the target grid region, where N is an odd number greater than 1; The second determining module is used to determine four initial terrain profile lines along the four axes of symmetry of the target grid region, respectively. The third determining module is used to determine the curve parameters of each of the initial terrain profile lines, wherein the curve parameters are used to describe the slope change of the initial terrain profile line when it passes through the central grid. The fourth determining module is used to determine whether the local geographic region corresponding to the initial grid is a saddle region based on the curve parameters of each initial terrain profile line. The fourth determining module is further configured to determine saddle degree information corresponding to each of the saddle regions, wherein the saddle degree information is used to describe the topographic relief of the saddle region; The method for obtaining the saddle depth information includes: Determine the parametric product between the curve parameters of the initial terrain profile lines corresponding to each pair of mutually perpendicular axes of symmetry; If the product of the parameters is less than 0, then a pair of target terrain profiles corresponding to the product of the parameters are determined; Determine the first average slope of each of the target terrain profiles in the pair of target terrain profiles, the first elevation value of the extreme point of the target curve in each of the target terrain profiles, and the horizontal distance between the extreme point of the target curve and the center point in each of the target terrain profiles; Determine the maximum elevation value, minimum elevation value, and second average slope of the first terrain profile in each of the target terrain profiles, wherein the first terrain profile is a portion of the target terrain profile between the extreme point of the target curve and the center point in the target terrain profile; Determine the second elevation value of the central grid; The candidate saddle point severity value is calculated using the following formula: ; Where D is the candidate saddle point degree value, and T is the product of the parameters. , The horizontal distances corresponding to a pair of target terrain profile lines are respectively. This is the second elevation value. , Let be the first elevation value corresponding to each pair of the target terrain profile lines. and Let be the first average slope corresponding to each pair of the target terrain profile lines. The maximum elevation value, The minimum elevation value, = + , and These are the second average slopes corresponding to each first terrain profile line in a pair of target terrain profile lines; Based on the candidate saddle severity values, a target saddle severity value corresponding to each saddle region is determined; The target saddle degree value is used as the saddle degree information.
8. An electronic device, characterized in that, include: At least one processor; as well as A memory communicatively connected to the at least one processor; wherein, The memory stores instructions that can be executed by the at least one processor to enable the at least one processor to perform the method of any one of claims 1-6.
9. A non-transitory computer-readable storage medium storing computer instructions, characterized in that, in, The computer instructions are used to cause the computer to perform the method according to any one of claims 1-6.
Citation Information
Patent Citations
Construction method for digital elevation model for discrete expression of landform on earth surface
CN101900546A
Complex geological curved surface feature extraction and reconstruction method
CN106446910A