A method, device, terminal device and storage medium for determining terrain complexity

Abstract

The present application belongs to the field of meteorology, and in particular relates to a method, apparatus and terminal device for determining terrain complexity, the method comprising: obtaining surface elevation data of a pre-study area; determining height data of a target terrain based on the surface elevation data and a preset reference height; determining the terrain complexity of the target terrain within a preset radius based on the height data of the target terrain and a preset radius, wherein the preset radius takes the center position of the typhoon as the sphere center, that is, the present application can determine the height data of the target terrain based on the surface elevation data and the preset reference height, and further determine the terrain complexity of the target terrain within the preset radius, and use the terrain complexity of the target terrain within the preset radius as a quantitative analysis indicator of the impact of the underlying surface terrain complexity on the wind field before and after the typhoon lands, thereby improving the fitting degree of the regression model of the wind circle radius before and after the typhoon lands.

Description

Technical Field

[0001] The present application relates to the field of meteorology, and in particular to a method, apparatus, terminal device and storage medium for determining terrain complexity. Background Art

[0002] A typhoon is a strong low-pressure vortex that generates and develops on tropical or subtropical oceans. It often brings severe weather such as strong winds, heavy rains and storm surges, causing serious casualties and huge property losses to the ocean (sea) and land it passes through. It is one of the most common meteorological disasters in the world.

[0003] The major disasters caused by typhoons are often caused before and after they make landfall. This is the period before and after a typhoon makes landfall when the typhoon's environmental field and underlying topography change most dramatically, and when the typhoon's structure shifts most strongly. The underlying land and sea surface and topography can cause significant changes in the typhoon's wind field distribution during landfall.

[0004] Currently, my country's typhoon monitoring and early warning systems still fail to meet the needs of disaster prevention and mitigation, and research on the wind field distribution and changes caused by landfalling typhoons is scarce. The main factors influencing the distribution of landfalling typhoons include the environmental field, the complexity of the underlying terrain, and the vortex structure of the typhoon itself. The complexity of the underlying terrain has a significant impact on the radius of the typhoon's wind circle, but there are currently no quantitative indicators or parameters to analyze the impact of the underlying terrain complexity on the radius of the typhoon's wind circle. Summary of the Invention

[0005] The embodiments of the present application provide a method, apparatus, terminal device, and storage medium for determining terrain complexity, thereby improving the fit of a wind circle radius regression model before and after a typhoon makes landfall.

[0006] In a first aspect, an embodiment of the present application provides a method for determining terrain complexity, the method comprising:

[0007] Obtain surface elevation data of the pre-study area;

[0008] Determining height data of a target terrain according to the surface elevation data and a preset reference height, wherein the target terrain is a terrain corresponding to the preset reference height in the pre-study area;

[0009] The terrain complexity of the target terrain within the preset radius is determined according to the height data of the target terrain and a preset radius, wherein the preset radius takes the center position of the typhoon as the center of the sphere.

[0010] In a possible implementation of the first aspect, obtaining surface elevation data of the pre-study area includes:

[0011] Acquire a digital elevation data file of a first resolution of the pre-study area;

[0012] Performing raster processing and grayscale processing on the digital elevation data file to obtain grayscale raster data of a first resolution of the pre-study area;

[0013] The surface elevation data of the pre-study area is acquired according to the grayscale raster data of the first resolution of the pre-study area.

[0014] The step of determining the height data of the target terrain according to the surface elevation data and the preset reference height includes:

[0015] Converting the grayscale raster data of the first resolution into grayscale raster data of a second resolution, wherein the second resolution is greater than the first resolution;

[0016] The height data of the target terrain is determined according to the grayscale raster data of the second resolution and the preset reference height.

[0017] The step of determining the height data of the target terrain based on the grayscale raster data of the second resolution and the preset reference height includes:

[0018] Determining that the surface elevation data obtained from the grayscale raster data of the first resolution and the number of pixel points of the grayscale raster data of the first resolution included in each pixel of the grayscale raster data of the second resolution are a first number;

[0019] Determining that the number of pixel points in each pixel of the grayscale raster data of the second resolution, in which the surface elevation data obtained from the grayscale raster data of the first resolution is greater than the preset reference height, is a second number;

[0020] The ratio of the second number to the first number is determined as the height data of the target terrain.

[0021] The step of determining the height data of the target terrain according to the surface elevation data and the preset reference height includes:

[0022] The height data of the target terrain is determined according to the grayscale raster data of the first resolution and the preset reference height.

[0023] The step of determining the height data of the target terrain based on the grayscale raster data of the first resolution and the preset reference height includes:

[0024] If the surface elevation data corresponding to a pixel of the grayscale raster data of the first resolution is greater than the preset reference height, determining that the height data corresponding to the pixel is a first binary digit;

[0025] If the surface elevation data corresponding to the pixel of the grayscale raster data of the first resolution is less than or equal to the preset reference height, determining that the height data corresponding to the pixel is a second binary number;

[0026] The height data of the target terrain is determined according to the binary number corresponding to each pixel.

[0027] The step of determining the terrain complexity of the target terrain within the preset radius based on the height data of the target terrain and the preset radius includes:

[0028] Determining a first distance between each pixel in the grayscale raster data of the first resolution or the grayscale raster data of the second resolution and the center position of the typhoon;

[0029] Determining the number of target pixels as a third number, the target pixels being pixels whose first distance from the center of the typhoon is less than the preset radius;

[0030] Determining the terrain complexity of the target pixel according to the height data of the target pixel and the target terrain, wherein the terrain complexity of the target pixel is the sum of the height data of the target terrain of the target pixel;

[0031] The ratio of the terrain complexity of the target pixel to the third number of the target pixel is determined as the terrain complexity of the target terrain within the preset radius.

[0032] In a second aspect, an embodiment of the present application provides a device for determining terrain complexity, the device comprising:

[0033] Acquisition module, used to obtain surface elevation data of the pre-study area;

[0034] A first determining module is configured to determine height data of a target terrain based on the surface elevation data and a preset reference height, wherein the target terrain is a terrain corresponding to the preset reference height in the pre-study area;

[0035] The second determining module is used to determine the terrain complexity of the target terrain within the preset radius according to the height data of the target terrain and a preset radius, wherein the preset radius takes the center position of the typhoon as the center of the sphere.

[0036] In a third aspect, an embodiment of the present application provides a terminal device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the method for determining the terrain complexity as described in any one of the first aspects is implemented.

[0037] In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method for determining the terrain complexity as described in any one of the first aspects is implemented.

[0038] Compared with the prior art, the beneficial effects of the embodiments of the present application are as follows: in the technical solution of the present application, by obtaining the surface elevation data of the pre-study area; determining the height data of the target terrain based on the surface elevation data and the preset reference height; determining the terrain complexity of the target terrain within the preset radius based on the height data of the target terrain and the preset radius, the preset radius takes the center position of the typhoon as the center of the sphere, that is, the present application can determine the height data of the target terrain based on the surface elevation data and the preset reference height, and then determine the terrain complexity of the target terrain within the preset radius, and use the terrain complexity of the target terrain within the preset radius as a quantitative analysis indicator of the impact of the underlying surface terrain complexity on the wind field before and after the typhoon lands, thereby improving the fitting degree of the regression model of the wind circle radius before and after the typhoon lands. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following briefly introduces the drawings required for use in the embodiments or descriptions of the prior art. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0040] Figure 1 This is a schematic diagram of an application scenario of the method for determining terrain complexity provided in one embodiment of the present application;

[0041] Figure 2 is a schematic flow chart of a method for determining the terrain complexity of a target terrain within a preset radius provided by an embodiment of the present application;

[0042] Figure 3 is a schematic flow chart of a method for obtaining surface elevation data of a pre-study area provided in an embodiment of the present application;

[0043] Figure 4 is a schematic flow chart of a method for determining height data of a target terrain provided in an embodiment of the present application;

[0044] Figure 5a is a schematic flow chart of a specific method of S402 provided in an embodiment of the present application;

[0045] Figure 5b This is an example diagram of grayscale raster data of a first resolution and grayscale raster data of a second resolution provided by an embodiment of the present application;

[0046] Figure 6a is a schematic flow chart of another specific method for determining height data of a target terrain provided in an embodiment of the present application;

[0047] Figure 6b This is an example diagram of grayscale raster data of a first resolution provided in an embodiment of the present application;

[0048] Figure 7 is a schematic flow chart of a specific method of S203 provided in an embodiment of the present application;

[0049] Figure 8 is a structural diagram of a device for determining terrain complexity provided in an embodiment of the present application;

[0050] Figure 9 It is a structural diagram of a terminal device provided in one embodiment of the present application. DETAILED DESCRIPTION

[0051] In the following description, specific details such as specific system structures and technologies are provided for the purpose of illustration rather than limitation, so as to provide a thorough understanding of the embodiments of the present application. However, it should be clear to those skilled in the art that the present application can also be implemented in other embodiments without these specific details. In other cases, detailed descriptions of well-known systems, devices, circuits, and methods are omitted to avoid obstructing the description of the present application with unnecessary details. In other cases, the specific technical details of the various embodiments can be referenced to each other, and specific systems not described in one embodiment can be referenced to other embodiments.

[0052] It should be understood that when used in the present specification and the appended claims, the term "comprising" indicates the presence of described features, integers, steps, operations, elements and / or components, but does not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components and / or collections thereof.

[0053] It will also be understood that the term "and / or" used in this specification and the appended claims refers to and includes any and all possible combinations of one or more of the associated listed items.

[0054] References to "an embodiment of the present application" or "some embodiments" in the present specification mean that one or more embodiments of the present application include specific features, structures, or characteristics described in conjunction with the embodiment. Therefore, phrases such as "in other embodiments," "an embodiment of the present application," and "other embodiments of the present application" appearing in different places in this specification do not necessarily refer to the same embodiment, but rather mean "one or more but not all embodiments," unless otherwise specifically emphasized. The terms "including," "comprising," "having," and their variations all mean "including but not limited to," unless otherwise specifically emphasized.

[0055] In addition, in the description of the present application specification and the appended claims, the terms "first", "second", etc. are only used to distinguish the descriptions and cannot be understood as indicating or implying relative importance.

[0056] Currently, my country's typhoon monitoring and early warning systems still fall short of meeting the needs of disaster prevention and mitigation, and research on the distribution and changes in wind fields caused by landfalling typhoons is scarce. The main factors influencing the distribution of landfalling typhoons include the environmental field, underlying topography, and the vortex structure of the typhoon itself. The complexity of the underlying topography has a significant impact on the radius of the typhoon's wind circle, but there are currently no quantitative indicators or parameters available to analyze the impact of underlying topography complexity on the radius of the typhoon's wind circle.

[0057] In order to solve the above-mentioned defects, the invention of this application is as follows:

[0058] This application can determine the height data of the target terrain based on the surface elevation data and the preset benchmark height, and then determine the terrain complexity within the preset radius. The terrain complexity within the preset radius is used as a quantitative analysis indicator of the impact of the terrain complexity of the underlying surface on the wind field before and after the typhoon lands, thereby improving the fitting degree of the regression model of the wind circle radius before and after the typhoon lands.

[0059] In order to illustrate the technical solution of the present application, specific embodiments are provided below.

[0060] Please refer to Figure 1 , Figure 1 This is a schematic diagram of an application scenario of the method for determining terrain complexity provided by an embodiment of the present application. For the sake of convenience, only the part related to the present application is shown. The application scenario includes: a terminal device 100.

[0061] The terminal device 100 includes but is not limited to a mobile phone, a tablet computer, a wearable device, a vehicle-mounted device, a laptop computer, an ultra-mobile personal computer (UMPC), a netbook, a personal digital assistant (PDA), etc. The embodiment of the present application does not impose any restrictions on the specific type of the terminal device.

[0062] The terminal device 100 is used to obtain surface elevation data of a pre-study area, and determine the terrain complexity of a target terrain within a preset radius based on the surface elevation data.

[0063] In the embodiment of the present application, please refer to the specific method for determining the terrain complexity of the target terrain within the preset radius. Figure 2 , Figure 2 This is a schematic flowchart of a method for determining the terrain complexity of a target terrain within a preset radius provided in an embodiment of the present application. Figure 2 The execution subject of the method in can be Figure 1 The terminal device 100 in Figure 2 As shown, the method includes: S201 to S203.

[0064] S201. The terminal device obtains surface elevation data of the pre-study area.

[0065] In the embodiment of the present application, the pre-study area specifically refers to the area within the range of 100° east longitude to 125° east longitude and 15° north latitude to 30° north latitude, and the surface elevation data refers to the height of the surface of the pre-study area from the sea level.

[0066] In the embodiment of the present application, the terminal device can obtain the surface elevation data of the pre-study area from the specific method. Figure 3 , Figure 3 This is a schematic flowchart of a method for obtaining surface elevation data of a pre-study area provided in an embodiment of the present application. Figure 3 The execution subject of the method in can be Figure 1 The terminal device 100 in Figure 3 As shown, the method includes: S301 to S303.

[0067] S301. The terminal device obtains a digital elevation data file of a first resolution of a pre-study area.

[0068] Specifically, the terminal device can obtain digital elevation data files from elevation collection websites, which include: Geospatial Data Cloud, National Earth System Science Data Sharing Platform, and Resource and Environmental Science and Data Center, etc. The embodiment of this application does not limit the way in which the terminal device obtains digital elevation data files.

[0069] A digital elevation data file, also known as a digital elevation model (DEM), is a physical ground model that represents ground elevation in the form of an ordered array of numerical values. DEM data can represent elevation data in a variety of ways, such as using a rectangular grid or a triangulated irregular network. This embodiment of the present application uses a rectangular grid as an example.

[0070] DEM resolution is an important indicator of the accuracy of DEM terrain depiction. It refers to the length of the smallest DEM cell. Since DEM data is discrete, the (X, Y) coordinates are small squares. The length of each small square is the DEM resolution, and each small square is labeled with an elevation value.

[0071] The first resolution of the embodiment of the present application does not specifically refer to a specific numerical resolution; it is merely used to distinguish it from different resolutions that may appear in subsequent embodiments. For example, the first resolution DEM data of the embodiment of the present application may be 0.00028°*0.00028° or 0.01°*0.01°, where both 0.00028° and 0.01° represent angles.

[0072] S302: The terminal device performs raster processing and grayscale processing on the digital elevation data file to obtain grayscale raster data of the first resolution of the pre-study area.

[0073] In the embodiments of the present application, since the digital elevation data file represents elevation data in the form of a rectangular grid, in some embodiments, the digital elevation data file obtained by S301 is obtained by obtaining multiple rectangular grids, and the multiple rectangular grids are merged to obtain the digital elevation data file of the pre-study area.

[0074] In an embodiment of the present application, the terminal device merges multiple rectangular grids through grid processing.

[0075] Exemplarily, a data management tool is configured in the terminal device, wherein the data management tool has the function of mosaicking a rectangular grid into a new grid, and grid processing can be completed through the data management tool.

[0076] In some embodiments, after performing raster processing on the digital elevation data file, the terminal device may use a data management tool to crop the merged grid data to remove elevation data outside the pre-study area.

[0077] In other embodiments, after the terminal device crops the merged grid data, it can project the cropped grid data through a data management tool. The purpose of projection is to prevent the raster data from being distorted (before projection, the terrain in the raster data is flatter than the real terrain shape, and this problem can be solved by projection processing).

[0078] In some embodiments, the terminal device is configured with an attribute extraction tool, and the attribute extraction tool performs negative value processing on the raster data to obtain raster data that is not negative. The purpose of performing negative value processing on the raster data is to remove abnormal data.

[0079] In the embodiment of the present application, the terminal device performs grayscale processing on the raster data to obtain grayscale raster data of the first resolution of the pre-study area. Grayscale processing methods include component method, maximum method, and average method, etc. The embodiment of the present application does not limit the grayscale processing method.

[0080] S303: The terminal device obtains the surface elevation data of the pre-study area based on the grayscale raster data of the first resolution of the pre-study area.

[0081] Specifically, the terminal device can obtain the elevation value corresponding to each pixel by reading the grayscale value of each pixel in the grayscale raster data of the first resolution of the pre-study area, and then obtain the surface elevation data of the South China area.

[0082] S202: The terminal device determines the height data of the target terrain based on the surface elevation data and the preset reference height.

[0083] Specifically, the target terrain is the terrain corresponding to the preset reference height in the pre-study area.

[0084] In the embodiments of the present application, the preset reference altitude refers to an altitude value from sea level pre-stored in the terminal device. For example, the preset reference altitude can be 0m, 20m, 50m, 100m, 200m, 500m, and 1000m, etc., where the specific value is the altitude value from sea level. The embodiments of the present application do not limit the number or specific values of the preset reference altitudes.

[0085] The preset reference height is configured in the embodiment of the present application to study the effect of terrain above the preset reference height on the radius of the typhoon wind circle. For example, the preset reference height is 20m to study the effect of terrain above 20m above sea level on the radius of the typhoon wind circle.

[0086] In some embodiments, the target terrain is terrain above a preset reference altitude. For example, if the preset reference altitude is 20 meters, the target terrain is terrain above 20 meters above sea level. For another example, if the preset reference altitude is 50 meters, the target terrain is terrain above 50 meters above sea level.

[0087] In the embodiment of the present application, for the specific method of determining the height data of the target terrain based on the surface elevation data and the preset reference height, please refer to Figure 4 , Figure 4 This is a schematic flowchart of a method for determining height data of a target terrain provided in an embodiment of the present application. Figure 4 The execution subject of the method in can be Figure 1 The terminal device 100 in Figure 4 As shown, the method includes: S401 to S402.

[0088] S401: The terminal device converts grayscale raster data of a first resolution into grayscale raster data of a second resolution.

[0089] Specifically, the second resolution is greater than the first resolution.

[0090] In an embodiment of the present application, if the first resolution of the grayscale raster data obtained by the terminal device is too small, the amount of grayscale raster data will be large, which in turn will result in a large amount of computation required for subsequent calculations using the grayscale raster data. In an embodiment of the present application, to reduce the amount of computation, the terminal device converts the grayscale raster data of the first resolution into grayscale raster data of the second resolution.

[0091] For example, the terminal device may convert grayscale raster data with a resolution of 0.00028°*0.00028° into grayscale raster data with a resolution of 0.01°*0.01°.

[0092] In the embodiment of the present application, the terminal device can convert the resolution of the grayscale raster data through the configured data management tool. The embodiment of the present application does not limit the resolution conversion method.

[0093] S402: The terminal device determines the height data of the target terrain according to the grayscale raster data of the second resolution and the preset reference height.

[0094] Specifically, the terminal device determines the height data of the target terrain based on the grayscale raster data of the second resolution (larger resolution) and the preset reference height, which can reduce the amount of calculation when determining the height data of the target terrain.

[0095] In the embodiment of this application, please refer to the specific method of S402. Figure 5a , Figure 5a It is a schematic flow chart of a specific method of S402 provided in an embodiment of the present application. Figure 5a The execution subject of the method in can be Figure 1 The terminal device 100 in Figure 5a As shown, the method includes: S501 to S503.

[0096] S501. The terminal device determines that the surface elevation data obtained from the grayscale raster data of the first resolution is included in each pixel of the grayscale raster data of the second resolution, and the number of pixel points of the grayscale raster data of the first resolution is a first number.

[0097] Please refer to Figure 5b , Figure 5b This is an example diagram of grayscale raster data of a first resolution and grayscale raster data of a second resolution provided in an embodiment of the present application.

[0098] exist Figure 5b In the figure, Figure A represents grayscale raster data with a primary resolution of 4*4, and Figure B represents a grayscale raster image with a secondary resolution of 2*2. For each pixel in Figure A, the 16 surface elevation data points obtained are: 23m, 19m, 13m, 12m, 17m, 16m, 39m, 57m, 57m, 87m, 14m, 11m, 66m, 96m, 10m, and 17m.

[0099] As can be seen from Figure B, one pixel of the grayscale raster data of the second resolution includes four pixels of the grayscale raster data of the first resolution.

[0100] In some embodiments, determining that the surface elevation data obtained from the grayscale raster data of the first resolution and included in each pixel of the grayscale raster data of the second resolution refers to the four elevation values contained in each pixel in the B image.

[0101] In some embodiments, the first number determined by the grayscale raster data of the first resolution included in each pixel of the grayscale raster data of the second resolution refers to the number of elevation values included in each pixel in the B image.

[0102] S502: The terminal device determines that the number of pixel points whose surface elevation data obtained from the grayscale raster data of the first resolution included in each pixel in the grayscale raster data of the second resolution is greater than a preset reference height is a second number.

[0103] Specifically, the second number refers to the number of the four elevation values contained in each pixel in the B image that is greater than the preset reference height.

[0104] The preset reference height in the embodiment of the present application is 20m. Of course, the method in the embodiment of the present application is also applicable to preset reference heights of other values.

[0105] For example, in Figure B, among the four elevation values of the first pixel from the left in the upper row, the number of pixels with elevation values greater than 20 m is 1. This embodiment of the present application refers to this number as the second number.

[0106] It should be noted that the method for determining the second number of the surface elevation data in the second resolution is the same as the method in the above embodiment and will not be described in detail here.

[0107] S503: The terminal device determines the ratio of the second number to the first number as the height data of the target terrain.

[0108] For example, please refer to Figure 5b In Figure C, the terminal device calculates the ratio of the second number 1 to the first number 4 in the first pixel from the left in the top row to be 1 / 4. This 1 / 4 is then determined as the elevation of a pixel in the second-resolution surface elevation data, i.e., the target terrain's height. This method determines the elevation of each pixel in the second-resolution surface elevation data. The elevation of each pixel in the second-resolution surface elevation data is referred to as the target terrain's height.

[0109] In the above embodiment, the height data of the target terrain is determined when the first resolution of the grayscale raster data acquired by the terminal device is too small, resulting in a large amount of grayscale raster data. In other embodiments, the height data of the target terrain can also be determined when the first resolution of the grayscale raster data acquired by the terminal device is appropriate (for example, the first resolution is 0.01°*0.01°). Another method for determining the height data of the target terrain includes:

[0110] The height data of the target terrain is determined according to the grayscale raster data of the first resolution and a preset reference height.

[0111] For another method to determine the height data of the target terrain, please refer to Figure 6a , Figure 6a This is a schematic flowchart of another specific method for determining the height data of the target terrain provided in an embodiment of the present application. Figure 6a The execution subject of the method in can be Figure 1 The terminal device 100 in Figure 6a As shown, the method includes: S601 to S603.

[0112] S601: If the surface elevation data corresponding to a pixel of the grayscale raster data of the first resolution is greater than a preset reference height, determine that the height data corresponding to the pixel is a first binary digit.

[0113] For details, please refer to Figure 6b , Figure 6b This is an example diagram of grayscale raster data of a first resolution provided in an embodiment of the present application. Figure 6b In the figure, D represents the grayscale raster data of the first resolution with appropriate resolution.

[0114] The preset reference height of the embodiment of the present application is 20m. The surface elevation data corresponding to the first pixel from the left in the first row in Figure D is 23m, which is greater than the preset reference height of 20m. Therefore, the height data corresponding to the first pixel from the left in the first row is determined to be the first binary digit. For example, the first binary digit is 1. Of course, in other embodiments, the first binary digit can be 0. In this way, the surface elevation data of each pixel in the grayscale raster data of the first resolution that is greater than the preset reference height can be converted into the first binary digit. For an example of conversion, refer to Figure 6b Figure E in .

[0115] S602: If the surface elevation data corresponding to the pixel of the grayscale raster data of the first resolution is less than or equal to a preset reference height, determine that the height data corresponding to the pixel is a second binary digit.

[0116] Specifically, the surface elevation data corresponding to the second pixel from the left in the first row in Figure D is 19m, which is less than the preset reference height of 20m, and the height data corresponding to the second pixel from the left in the first row is determined to be the second binary digit. The surface elevation data corresponding to the third pixel from the left in the first row in Figure D is 20m, which is equal to the preset reference height of 20m, and the height data corresponding to the third pixel from the left in the first row is determined to be the second binary digit. For example: the second binary digit is 0. Of course, in other embodiments, the second binary digit may be 1. The second binary digit and the first binary digit in the embodiment of the present application are different binary digits. In this way, the surface elevation data in each pixel in the grayscale raster data of the first resolution that is less than or equal to the preset reference height can be converted into the second binary digit. For example, see . Figure 6b Figure E in .

[0117] S603: The terminal device determines the height data of the target terrain according to the binary number corresponding to each pixel.

[0118] Specifically, the first binary digit corresponding to each pixel is the height data of the target terrain.

[0119] S203: The terminal device determines the terrain complexity of the target terrain within the preset radius according to the height data of the target terrain and the preset radius.

[0120] Specifically, the preset radius takes the center of the typhoon as the center of the sphere. For example, the preset radius can be 400 km to 600 km, for example, the preset radius is 500 km.

[0121] In the embodiment of this application, please refer to the specific method of S203. Figure 7 , Figure 7 This is a schematic flow chart of a specific method of S203 provided in an embodiment of the present application. Figure 7The execution subject of the method in can be Figure 1 The terminal device 100 in Figure 7 As shown, the method includes: S701 to S704.

[0122] S701: The terminal device determines a first distance between each pixel in the grayscale raster data of the first resolution or the grayscale raster data of the second resolution and the center position of the typhoon.

[0123] Specifically, the longitude and latitude of each pixel can be obtained according to the grayscale raster data of the first resolution or the grayscale raster data of the second resolution.

[0124] For example, a pixel in the first-resolution grayscale raster data is called point A, and the center of the typhoon is called point B. The longitude of point A is α1 and the latitude is β1; the longitude of point B is α2 and the latitude is β2, where α1, α2, β1, and β2 are all radians, and the radius of the earth is R (km). It is also agreed that east longitude is positive and west longitude is negative; north latitude is positive and south latitude is negative. The distance from point A to point B can be calculated using the following formula:

[0125] dist=R*arcos(sin(β1)*sin(β2)+cos(β1)*cos(β2)*cos(α1-α2)).

[0126] Wherein, dist represents the distance from point A to point B, that is, the first distance of this application.

[0127] According to the above formula, the first distance between each pixel in the grayscale raster data of the first resolution or the grayscale raster data of the second resolution and the center position of the typhoon can be determined.

[0128] S702: The terminal device determines that the number of target pixels is a third number.

[0129] Specifically, the target pixel refers to a pixel whose first distance from the pixel point to the center position of the typhoon is less than a preset radius.

[0130] In this embodiment of the present application, a first distance between each pixel and the typhoon center is determined in step S701, and the first distance is compared with a preset radius (e.g., 500 km). Pixels corresponding to pixels whose first distance is less than the preset radius are target pixels. In this embodiment of the present application, the number of target pixels is referred to as the third number.

[0131] S703: The terminal device determines the terrain complexity of the target pixel based on the height data of the target pixel and the target terrain.

[0132] Specifically, the terrain complexity of the target pixel is the sum of the height data of the target terrain of the target pixel.

[0133] In some embodiments, if the resolution of the grayscale raster data acquired by the terminal device is appropriate, the terrain complexity of the target terrain determined based on the height data of the target terrain can be represented by the following example matrix:

[0134]

[0135] Through this matrix, the terrain complexity of the target terrain is determined to be: 1+0+0+0+0+0+1+1+1+1+0+0+1+1+0+0=7.

[0136] If the target pixel is Figure 6b In the second and third pixels from the left in the second row and the second and third pixels from the left in the third row of the D image, the terrain complexity of the target pixel can be determined based on the terrain complexity of the target terrain and can be represented by the following example matrix:

[0137]

[0138] Through this matrix, the terrain complexity of the target terrain is determined to be 0+1+1+0=2.

[0139] In other embodiments, if the resolution of the grayscale raster data acquired by the terminal device is relatively small, the terrain complexity of the target terrain determined based on the height data of the target terrain can be represented by the following example matrix:

[0140]

[0141] Through this matrix, the terrain complexity of the target terrain is determined to be 1 / 4+2 / 4+4 / 4+0 / 4=7 / 4.

[0142] If the target pixel is Figure 5b For the first and second pixels from the left in the first row of Figure B, the terrain complexity of the target pixel can be determined based on the terrain complexity of the target terrain and can be represented by the following example matrix:

[0143] [1 / 4 2 / 4].

[0144] Through this matrix, the terrain complexity of the target terrain is determined to be 1 / 4+2 / 4=3 / 4.

[0145] S704: The terminal device determines the ratio of the terrain complexity of the target pixel to the third number of the target pixel as the terrain complexity of the target terrain within the preset radius.

[0146] For example, if the terrain complexity of the target pixel is the sum of the data in the following matrix:

[0147]

[0148] The third number of the target pixel is 4, and the terrain complexity of the target terrain within the preset radius is (0+1+1+0) / 4=1 / 2.

[0149] If the terrain complexity of the target pixel is the sum of the data in the following matrix:

[0150] [1 / 4 2 / 4].

[0151] The third number of the target pixel is 2, and the terrain complexity of the target terrain within the preset radius is (1 / 4+2 / 4) / 2=3 / 8.

[0152] In an embodiment of the present application, after determining the terrain complexity of the target terrain within a preset radius, the terminal device can input it as a prediction factor into the multivariate regression model of the landing typhoon wind circle radius, thereby significantly improving the fitting degree of the multivariate regression model of the landing typhoon wind circle radius.

[0153] In summary, the technical solution of the present application obtains surface elevation data of the pre-study area; determines the height data of the target terrain based on the surface elevation data and the preset benchmark height; determines the terrain complexity of the target terrain within the preset radius based on the height data of the target terrain and the preset radius, and the preset radius takes the center position of the typhoon as the center of the sphere, that is, the present application can determine the height data of the target terrain based on the surface elevation data and the preset benchmark height, and then determine the terrain complexity of the target terrain within the preset radius, and use the terrain complexity of the target terrain within the preset radius as a quantitative analysis indicator of the impact of the underlying surface terrain complexity on the wind field before and after the typhoon lands, thereby improving the fitting degree of the regression model of the wind circle radius before and after the typhoon lands.

[0154] It should be understood that the size of the serial numbers of the steps in the above embodiments does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of this application.

[0155] Please refer to Figure 8 , Figure 8 : is a schematic structural diagram of a device for determining terrain complexity provided in an embodiment of the present application, the device comprising:

[0156] The acquisition module 81 is used to acquire the surface elevation data of the pre-study area.

[0157] The first determining module 82 is configured to determine height data of a target terrain according to the surface elevation data and a preset reference height. The target terrain is a terrain corresponding to the preset reference height in the pre-study area.

[0158] The second determining module 83 is configured to determine the terrain complexity of the target terrain within a preset radius according to the height data of the target terrain and a preset radius, wherein the preset radius takes the center of the typhoon as the center of the sphere.

[0159] The acquisition module 81 is further configured to acquire a digital elevation data file of a first resolution in the pre-study area;

[0160] Perform raster processing and grayscale processing on the digital elevation data file to obtain the grayscale raster data of the first resolution of the pre-study area;

[0161] The surface elevation data of the pre-study area is obtained based on the grayscale raster data of the first resolution of the pre-study area.

[0162] The first determining module 82 is further configured to convert the grayscale raster data of the first resolution into grayscale raster data of a second resolution, where the second resolution is greater than the first resolution.

[0163] The height data of the target terrain is determined according to the grayscale raster data of the second resolution and a preset reference height.

[0164] The first determining module 82 is further configured to determine that the surface elevation data obtained from the grayscale raster data of the first resolution is included in each pixel of the grayscale raster data of the second resolution, and the number of pixels of the grayscale raster data of the first resolution is a first number;

[0165] Determine the number of pixel points whose surface elevation data obtained from the grayscale raster data of the first resolution is greater than a preset reference height and is included in each pixel of the grayscale raster data of the second resolution as a second number;

[0166] The ratio of the second number to the first number is determined as the height data of the target terrain.

[0167] The first determining module 82 is further configured to determine the height data of the target terrain according to the grayscale raster data of the first resolution and a preset reference height.

[0168] The first determining module 82 is further configured to determine that the height data corresponding to a pixel of the grayscale raster data of the first resolution is a first binary digit if the surface elevation data corresponding to the pixel is greater than a preset reference height.

[0169] If the surface elevation data corresponding to the pixel of the grayscale raster data of the first resolution is less than or equal to the preset reference height, determining the height data corresponding to the pixel to be a second binary number;

[0170] The height data of the target terrain is determined based on the binary number corresponding to each pixel.

[0171] The second determining module 83 is further configured to determine a first distance between each pixel in the grayscale raster data of the first resolution or the grayscale raster data of the second resolution and the center position of the typhoon;

[0172] Determining the number of target pixels as the third number, where the target pixels refer to pixels whose first distance from the center of the typhoon is less than a preset radius;

[0173] Determining the terrain complexity of the target pixel according to the height data of the target pixel and the target terrain, wherein the terrain complexity of the target terrain is the sum of the height data of the target pixel and the target terrain;

[0174] The ratio of the terrain complexity of the target pixel to the third number of the target pixel is determined as the terrain complexity of the target terrain within the preset radius.

[0175] In summary, the technical solution of the present application can determine the terrain complexity of the target terrain within a preset radius, and use the terrain complexity of the target terrain within the preset radius as a quantitative analysis indicator of the impact of the underlying surface terrain complexity on the typhoon wind field, thereby improving the fitting degree of the regression model of the wind circle radius before and after the typhoon lands.

[0176] The technical solution of the present application can convert grayscale raster data of a first resolution into grayscale raster data of a second resolution, thereby reducing the amount of calculation while ensuring the integrity of the grayscale raster data of the first resolution.

[0177] The technical solution of the present application can determine the ratio of the second number of surface elevation data obtained from the grayscale raster data of the first resolution included in each pixel of the grayscale raster data of the second resolution that is greater than the preset reference height to the first number of the surface elevation data obtained from the grayscale raster data of the first resolution as the height data of the target terrain, and determine the height data of the target terrain based on the binary number corresponding to each pixel. The form of the determined height data of the target terrain is novel, which is convenient for calculating the terrain complexity of the target pixel within the preset radius and the terrain complexity of the target terrain.

[0178] Those skilled in the art can clearly understand that, for the convenience and brevity of description, only the division of the above-mentioned functional units and modules is used as an example for illustration. In actual applications, the above-mentioned functions can be distributed and completed by different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above. The functional units and modules in the embodiment can be integrated into one processing unit, or each unit can exist physically alone, or two or more units can be integrated into one unit. The above-mentioned integrated unit can be implemented in the form of hardware or in the form of software functional units. In addition, the specific names of the functional units and modules are only for the convenience of distinguishing each other, and are not used to limit the scope of protection of this application. The specific working process of the units and modules in the above-mentioned system can refer to the corresponding process in the aforementioned method embodiment, and will not be repeated here.

[0179] like Figure 9 As shown, an embodiment of the present application also provides a terminal device 200, including a memory 21, a processor 22, and a computer program 23 stored in the memory 21 and executable on the processor 22. When the processor 22 executes the computer program 23, the method for determining the terrain complexity of the above-mentioned embodiments is implemented.

[0180] The processor 22 may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor.

[0181] The memory 21 may be an internal storage unit of the terminal device 200. The memory 21 may also be an external storage device of the terminal device 200, such as a plug-in hard disk, a SmartMedia Card (SMC), a Secure Digital (SD) card, a flash memory card, etc. equipped on the terminal device 200. Furthermore, the memory 21 may include both an internal storage unit of the terminal device 200 and an external storage device. The memory 21 is used to store computer programs and other programs and data required by the terminal device 200. The memory 21 may also be used to temporarily store data that has been output or is about to be output.

[0182] An embodiment of the present application further provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the method for determining the terrain complexity of the above-mentioned embodiments is implemented.

[0183] An embodiment of the present application provides a computer program product. When the computer program product is run on a mobile terminal, the mobile terminal implements the method for determining the terrain complexity of the above-mentioned embodiments.

[0184] If the integrated unit is implemented as a software functional unit and sold or used as a standalone product, it can be stored in a computer-readable storage medium. Based on this understanding, the present application implements all or part of the process steps in the above-mentioned method embodiments by instructing the relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium. When executed by a processor, the computer program can implement the steps of each of the above-mentioned method embodiments. The computer program includes computer program code, which can be in source code form, object code form, executable file, or some intermediate form. The computer-readable storage medium may include at least: any entity or device capable of carrying computer program code to the camera / terminal device, recording medium, computer memory, read-only memory (ROM), random access memory (RAM), electric carrier signal, telecommunication signal, and software distribution medium. Examples include USB flash drives, removable hard drives, magnetic disks, or optical disks. In some jurisdictions, based on legislation and patent practice, computer-readable storage media cannot be electric carrier signals or telecommunication signals.

[0185] In the above embodiments, the description of each embodiment has its own focus. For parts that are not described or recorded in detail in a certain embodiment, reference can be made to the relevant description of other embodiments.

[0186] Those skilled in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0187] The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed across multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the embodiments of the present application.

[0188] The above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. These modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the various embodiments of the present application, and should all be included in the scope of protection of the present application.

Claims

1. A method for determining terrain complexity, characterized in that: The determination method includes: Obtain the first-resolution digital elevation data file of the pre-study area; Performing raster processing and grayscale processing on the digital elevation data file to obtain grayscale raster data of a first resolution of the pre-study area; Acquiring surface elevation data of the pre-study area based on the grayscale raster data of the first resolution of the pre-study area; If the surface elevation data corresponding to a pixel of the grayscale raster data of the first resolution is greater than a preset reference height, determining that the height data corresponding to the pixel is a first binary digit; If the surface elevation data corresponding to the pixel of the grayscale raster data of the first resolution is less than or equal to the preset reference height, determining that the height data corresponding to the pixel is a second binary number; Determining height data of a target terrain according to a binary number corresponding to each pixel, wherein the target terrain is a terrain corresponding to the preset reference height of the pre-study area; Determining a first distance between the grayscale raster data of the first resolution and the center position of the typhoon; Determining the number of target pixels as a third number, wherein the target pixels refer to pixels whose first distance from the center of the typhoon is less than a preset radius; Determining the terrain complexity of the target pixel according to the height data of the target pixel and the target terrain, wherein the terrain complexity of the target pixel is the sum of the height data of the target terrain of the target pixel; The ratio of the terrain complexity of the target pixel to the third number of the target pixel is determined as the terrain complexity of the target terrain within the preset radius, and the preset radius takes the center position of the typhoon as the center of the sphere.

2. The determination method according to claim 1, characterized in that The method further comprises: Converting the grayscale raster data of the first resolution into grayscale raster data of a second resolution, wherein the second resolution is greater than the first resolution; The height data of the target terrain is determined according to the grayscale raster data of the second resolution and the preset reference height.

3. The determination method according to claim 2, characterized in that: The step of determining the height data of the target terrain according to the grayscale grid data of the second resolution and the preset reference height includes: Determining that the surface elevation data obtained from the grayscale raster data of the first resolution and the number of pixel points of the grayscale raster data of the first resolution included in each pixel of the grayscale raster data of the second resolution are a first number; Determining that the number of pixel points of the surface elevation data obtained from the grayscale raster data of the first resolution included in each pixel of the grayscale raster data of the second resolution that is greater than the preset reference height is a second number; The ratio of the second number to the first number is determined as the height data of the target terrain.

4. The determination method according to claim 3, characterized in that: The method further comprises: Determining a first distance between each pixel in the grayscale raster data of the second resolution and the typhoon center position; Determining the number of target pixels as a third number, wherein the target pixels refer to pixels whose first distance from the center of the typhoon is less than the preset radius; Determining the terrain complexity of the target pixel according to the height data of the target pixel and the target terrain, wherein the terrain complexity of the target pixel is the sum of the height data of the target terrain of the target pixel; The ratio of the terrain complexity of the target pixel to the third number of the target pixel is determined as the terrain complexity of the target terrain within the preset radius.

5. A device for determining terrain complexity, characterized in that: The determining device comprises: an acquisition module, configured to acquire a digital elevation data file of a first resolution of a pre-study area; perform raster processing and grayscale processing on the digital elevation data file to obtain grayscale raster data of the first resolution of the pre-study area; and acquire surface elevation data of the pre-study area based on the grayscale raster data of the first resolution of the pre-study area; a first determining module configured to determine, if the surface elevation data corresponding to a pixel of the grayscale raster data of the first resolution is greater than a preset reference height, the elevation data corresponding to the pixel as a first binary digit; and, if the surface elevation data corresponding to a pixel of the grayscale raster data of the first resolution is less than or equal to the preset reference height, the elevation data corresponding to the pixel as a second binary digit; and determine, based on the binary digit corresponding to each pixel, elevation data of a target terrain, the target terrain being the terrain corresponding to the preset reference height in the pre-study area; The second determination module is used to determine the first distance between the first resolution grayscale raster data and the center position of the typhoon; determine the number of target pixels as a third number, the target pixel refers to a pixel whose first distance from the center position of the typhoon is less than a preset radius; determine the terrain complexity of the target pixel based on the height data of the target pixel and the target terrain, the terrain complexity of the target pixel being the sum of the height data of the target terrain of the target pixel; and determine the ratio of the terrain complexity of the target pixel to the third number of the target pixel as the terrain complexity of the target terrain within the preset radius, where the preset radius takes the center position of the typhoon as the sphere center.

6. A terminal device, characterized in that: The method comprises a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the method for determining the terrain complexity according to any one of claims 1 to 4 is implemented.

7. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the method for determining terrain complexity according to any one of claims 1 to 4 is implemented.

Images