Disclosure of Invention
In order to solve the problems in the prior art, the application provides a method, a device and a storage medium for obtaining the land area, which can restore the actual range of the linear land object, the actual area of the land type pattern spot and the actual deducted land type area on the map surface, and improve the accuracy of the area analysis.
In order to achieve the above object, a first aspect of the present application provides a method for acquiring a land area, including:
superposing a spatial analysis layer of a target area and a map type spot layer corresponding to the target area to obtain an area superposition result;
acquiring a corresponding ground type map spot superposition result and a linear ground object superposition result according to the area superposition result;
acquiring an updated deduction relation according to the linear ground object superposition result;
and calculating and acquiring the actual area of the target area according to the land type pattern spot superposition result, the linear ground object superposition result and the updated deduction relation.
Optionally, the obtaining a linear ground object superposition result according to the area superposition result includes:
intersecting the area superposition result with the linear ground object to obtain an intersecting result;
and acquiring a linear ground object superposition result according to the spatial relationship of the linear ground objects in the intersection result.
Optionally, the calculating and obtaining the actual area of the target region according to the land type map spot superposition result, the linear ground object superposition result, and the updated deduction relationship includes:
obtaining deduction times of the linear ground object according to the linear ground object superposition result;
recalculating the actual area of the linear ground object in the target area according to the deduction times; calculating and acquiring the actual area of the land type pattern spots in the target area according to the overlapping result of the land type pattern spots and a preset calculation formula;
and calculating and acquiring the actual area of the target area according to the actual area of the linear ground object in the target area and the actual area of the ground class pattern spot in the target area.
Optionally, the method further comprises:
backfilling the actual total area of the linear ground object into the region superposition result to obtain a backfilling result;
and acquiring the actual deducted land area in the target area according to the actual area of the land pattern spot, the backfill result and the known deducted land coefficient.
Optionally, the calculating and obtaining the actual area of the target region according to the land type map spot superposition result, the linear ground object superposition result, and the updated deduction relationship includes:
And calculating and obtaining the actual area of the target area pattern spot land according to the land pattern spot superposition result, the linear land feature superposition result, the deduction land area and the updated deduction relation.
This application second aspect provides a land area acquisition device, includes: the system comprises a superposition unit, an acquisition unit and a calculation unit, wherein the superposition unit is used for superposing a spatial analysis map layer of a target area and a geographical map speckle map layer corresponding to the target area to acquire an area superposition result;
the acquisition unit is used for acquiring a corresponding ground type image spot superposition result and a linear ground object superposition result according to the area superposition result; acquiring an updated deduction relation according to the linear ground object superposition result;
and the calculating unit is used for calculating and acquiring the actual area of the target area according to the land type map spot superposition result, the linear ground object superposition result and the updated deduction relation.
Optionally, the obtaining unit is configured to intersect the area superposition result with the linear ground object to obtain an intersection result; and acquiring a linear ground object superposition result according to the spatial relationship of the linear ground objects in the intersection result.
Optionally, the calculating unit is configured to obtain a deduction number of the linear feature according to the linear feature superposition result; recalculating the actual area of the linear ground object in the target area according to the deduction times; calculating and acquiring the actual area of the geographical map spots in the target area according to the superposition result of the geographical map spots and a preset calculation formula;
and calculating and acquiring the actual area of the target area according to the actual area of the linear ground object in the target area and the actual area of the ground class pattern spot in the target area.
Optionally, the apparatus further comprises: the backfilling unit is used for backfilling the actual total area of the linear ground object into the region superposition result to obtain a backfilling result;
and the acquisition unit is used for acquiring the actual deducted land area in the target area according to the actual area of the land pattern spot, the backfill result and the known deducted land coefficient.
Optionally, the calculating unit is configured to calculate and obtain an actual area of the target region map spot according to the map spot superposition result, the linear ground object superposition result, the deducted ground area, and the updated deduction relationship.
The third aspect of the present application provides a land area acquisition apparatus, comprising: a processor, a storage medium and a bus, the storage medium storing machine-readable instructions executable by the processor, the processor and the storage medium communicating via the bus when the apparatus is operated, the processor executing the machine-readable instructions to perform the steps of the method of the first aspect of the application.
A fourth aspect of the present application provides a storage medium having stored thereon a computer program for executing the steps of the method of the first aspect of the present application when the computer program is run by a processor.
According to the method, the device and the storage medium for acquiring the land area, the space analysis layer of the target area is overlapped with the map layer of the land type corresponding to the target area, and an area overlapping result is acquired; acquiring a corresponding ground type map spot superposition result and a linear ground object superposition result according to the area superposition result; acquiring an updated deduction relation according to the linear ground object superposition result; and calculating and acquiring the actual area of the target area according to the land type map spot superposition result, the linear ground object superposition result and the updated deduction relation. By means of the superposition operation among the graph layers and the mode of updating the linear ground object deduction relation, the actual range of the linear ground object and the actual area of the ground type graph spot can be restored on the graph surface, and the accuracy of area analysis is improved.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are some embodiments of the present application, but not all embodiments. The components of the embodiments of the present application, as generally described and illustrated in the figures herein, may be arranged and designed in a wide variety of different configurations.
Thus, the following detailed description of the embodiments of the present application, as presented in the figures, is not intended to limit the scope of the claimed application, but is merely representative of selected embodiments of the application. All other embodiments obtained by a person of ordinary skill in the art without any inventive work based on the embodiments in the present application are within the scope of protection of the present application.
It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures.
The embodiment of the application provides a land area obtaining method aiming at the problem that the error between the land area obtained by adopting the existing deduction relation and the actual area is large.
Fig. 1 is a schematic flowchart of a method for acquiring a land area according to an embodiment of the present application, where the method may be implemented by a computer, a processor, and other devices with a computing function, and is not limited herein.
As shown in fig. 1, the method includes:
s101, overlapping the spatial analysis layer of the target area with the geographical map speckle pattern layer corresponding to the target area to obtain an area overlapping result.
In the embodiment of the application, the land type map layer corresponding to the target area is obtained according to the "land utilization database standard (TD/T1016-:
TABLE 1-map layer attribute structure description of geographical type (DLTB for short)
TABLE 2-description table of attribute structure of linear ground object layer (abbreviate-XZDW)
The map layer of land type map (DLTB) represents a kind of vector data, which may be, for example, second national land survey data (binary data). The spatial analysis layer (ClipFC for short) is also used as vector data, specifically, the spatial analysis layer can be third-time national land survey data (third-tone data), and is different from the land type map layer, and the spatial analysis layer refers to survey data to be analyzed.
It should be noted that the layer overlay performed in the embodiment of the present application may be implemented based on a scalable and comprehensive Geographic Information System (GIS), for example, the layer overlay may be an ArcGis platform.
In the embodiment of the application, according to the standard of land utilization database (TD/T1016-2007), the relevant areas in the map layer of the geological map and the linear surface feature need to be ellipsoid areas. When the ellipsoid areas of the divided image spots are calculated, in order to improve the calculation speed and avoid repeated calculation, the ellipsoid areas after division are calculated by adopting a method of dividing the proportion of the images before and after division. Before the space analysis, namely the space analysis map layer and the land type map layer corresponding to the target area are superposed, relevant fields are added in the land type map layer and the linear ground object layer for storing the image area and the image length before the division.
Illustratively, a field is added in the map layer of the geographical type map spot for saving the graphic area of the geographical type map spot before analysis, wherein:
area, shape area (TXMJ);
adding a field for storing the graph length of the linear ground object before analysis on the linear ground object layer;
length, length.
And performing superposition operation on a map layer (DLTB) of the land type map and a map layer (Clip FC) of the spatial analysis range, and outputting a region superposition result (DLTB _ Clip for short), wherein the graphic data in the DLTB _ Clip map layer is a graphic divided by the DLTB according to the Clip FC map layer, and the attribute data in the DLTB _ Clip map layer is a union of the DLTB and the Clip FC.
And S102, acquiring a corresponding ground type map spot superposition result and a linear ground object superposition result according to the area superposition result.
In this embodiment of the present application, the obtaining of the corresponding overlay result of the geographical map spots may specifically include:
and adding a spot area field UNION-TBMJ in the DLTB _ Clip for recalculating the spot area of the DLTB _ Clip after the combination.
Acquiring a corresponding linear ground object superposition result may include:
and (3) performing overlapping intersection operation on the linear ground object (XZDW), the spatial analysis layer and the ground map spot layer, namely (DLTB _ Clip), wherein the output layer is a result (XZDW _ Clip) of overlapping intersection of the linear ground object, the graphic data in the XZDW _ Clip layer is a graphic divided by the DLTB _ Clip layer, and the attribute data in the XZDW _ Clip layer is a union of the XZDW and the DLTB _ Clip.
And S103, acquiring an updated deduction relation according to the linear ground object superposition result.
Optionally, in this embodiment of the present application, the updated linear feature length (CD) is the linear feature length before update (CD) × shape. Length may be the updated linear ground object graph length.
And adding a linear ground object area field INTERSECT _ XZDWMJ in the XZDW _ Clip layer for recalculating the linear ground object area of the XZDW _ Clip after intersection.
The area of the line feature (INTERSECT _ XZDWMJ) is the updated length (CD) of the line feature and the width (KD).
The repetition number N of each linear feature element is counted according to the attribute BSM (identification code) and shape in the XZDW _ Clip layer.
And S104, calculating and acquiring the actual area of the target area according to the ground type pattern spot superposition result, the linear ground object superposition result and the updated deduction relation.
In this embodiment of the present application, the obtaining the actual area of the target region specifically includes: and acquiring the actual area of the linear ground object in the target area and the actual area of the ground pattern spot. Wherein, obtaining the actual area of the linear ground object in the target area comprises:
The recalculated linear ground object area (INTERSECT _ XZDWMJ) is equal to the linear ground object area (INTERSECT _ XZDWMJ)/N.
In the land area obtaining method provided by this embodiment, a spatial analysis layer of a target area is superimposed on a land type map patch layer corresponding to the target area, so as to obtain an area superimposition result; acquiring a corresponding ground type map spot superposition result and a linear ground object superposition result according to the region superposition result; acquiring an updated deduction relation according to the linear ground object superposition result; and calculating and acquiring the actual area of the target area according to the terrain map spot superposition result, the linear ground object superposition result and the updated deduction relation. By means of the superposition operation among the graph layers and the mode of updating the linear ground object deduction relation, the actual range of the linear ground object and the actual area of the ground type graph spot can be restored on the graph surface, and the accuracy of area analysis is improved.
It should be noted that, in the embodiment of the present application, the step S102 obtaining the linear ground object superposition result according to the area superposition result specifically includes: intersecting the area superposition result with the linear ground object to obtain an intersection result; and acquiring a linear ground object superposition result according to the spatial relationship of the linear ground objects in the intersection result.
The intersection operation generates two results, i.e. within the patch or by-pass (on the edge of the patch), the subtraction factor is 1 if the linear feature is within the patch, and 0.5 if the linear feature is on the edge of the patch.
As shown in fig. 2, in some embodiments of the present application, step S104 may specifically include: steps S1041 to S1044 are as follows:
and S1041, obtaining the deduction times of the linear ground object according to the superposition result of the linear ground object.
In the embodiment of the present application, first, a linear feature (XZDW) and a superposition result DLTB _ Clip of a spatial analysis range map layer and a feature map are used to perform a superposition intersection operation, and a result XZDW _ Clip after the superposition intersection is output, and then, the number N of times of repetition of each linear feature element is counted according to an attribute BSM (identification code) and a shape.
And S1042, recalculating the actual area of the linear ground object in the target area according to the deduction times.
The area of the linear ground object (INTERSECT _ XZDWMJ) is calculated by the length (CD) and the width (KD), and the area of the linear ground object (INTERSECT _ XZDWMJ) is calculated by the width (KD).
And S1043, calculating and obtaining the actual area of the geographical map spots in the target area according to the geographical map spot superposition result and a preset calculation formula.
And adding a spot area field UNION-TBMJ in the land type spot superposition result (DLTB _ Clip) for recalculating the spot area of the DLTB _ Clip after the UNION.
Optionally, the updated ground class patch area (UNION _ TBMJ) ═ the pre-update ground class patch area (TBMJ) · shape.
Area represents an updated graphic area in the embodiment of the present application.
And S1044, calculating and acquiring the actual area of the target area according to the actual area of the linear ground object in the target area and the actual area of the ground pattern spot in the target area.
In this embodiment, the actual area of the target region is the actual area of the ground-like pattern patch — the actual area of the linear ground object.
In the present embodiment, the actual area of the linear feature refers to the actual total area of the linear features in the target region.
Optionally, in another embodiment, in order to more accurately obtain the actual area of the target region, the actual subtracted geo-area in the target region may also be taken into account. That is, on the basis of the above embodiment, the method further includes: backfilling the actual total area of the linear ground object into the region superposition result to obtain a backfilling result; and acquiring the actual deducted geo-area in the target area according to the actual area of the geo-map spot, the backfill result and the known deducted geo-coefficient.
It should be noted that there may be two cases of subtracting the land-like coefficient: 1 or more or less than 1. The subtracted terrain area (TKMJ) may be calculated according to the following formula:
when the subtracted terrain class coefficient (TKXS) > 1: expressed as a percentage of denominator 100;
deducted land area (TKMJ) ((land pattern spot area (UNION _ TBMJ) -line land area (UNION _ XZDWMJ))) deducted land coefficient (TKXS)/100.
When the subtracted terrain class coefficient (TKXS) < 1: expressed as a percentage of denominator 1;
deducted land area (TKMJ) ((land pattern spot area (UNION _ TBMJ) -line land area (UNION _ XZDWMJ))) deducted land coefficient (TKXS).
Alternatively, the deduction relationship is reconstructed according to the conditions that the attribute fields of the result (XZDW _ Clip) after the line ground object intersection deduct the patch number 1(KCTBBH1), deduct the patch unit code 1(KCTBDWDM1), deduct the patch number 2(KCTBBH2), deduct the patch unit code 2(KCTBDWDM2), the patch number (TBBH) and the seating unit code (ZLDWDM).
1) The deduction relationship of the attribute field is as follows:
note: considering that the subtracted spot number and the spot number are not unique, the subtracted spot number and the subtracted spot unit code are analyzed together as a unique value.
A1 is in (a,0), that is, a1 is equal to a, and at this time, the linear ground object area is given to the deduction patch number + deduction patch unit code corresponding to a;
b1 is not in (A,0), namely B1 is not equal to A, and the linear ground object area is deducted to the deduction pattern spot number corresponding to A + the largest pattern spot in deduction pattern unit codes;
a1 is in (a, B), that is, a1 ═ a or a1 ═ B, at this time, the linear feature area is deducted to the deduction pattern number + deduction pattern unit code corresponding to a 1;
and fourthly, B1 is in (a, B), that is, B1 is equal to a or B1 is equal to B, and at this time, the linear ground object area is deducted to the deduction pattern number + deduction pattern unit code corresponding to B1.
2) And adding a UNION _ XZDWMJ field in the land type map spot superposition result (DLTB _ Clip).
3) And (2) correspondingly backfilling the linear ground object area (INTERSECT _ XZDWMJ) in the result (XZDW _ Clip) after the linear ground objects are intersected into the UNION _ XZDWMJ field in the result (DLTB _ Clip) after the combination of the ground type map spots according to the deduction relation of the attribute field in the step 1).
4) Updating the deduction land area:
when the subtracted terrain class coefficient (TKXS) > 1: expressed as a percentage of denominator 100;
deducted land area (TKMJ) ((land pattern spot area (UNION _ TBMJ) -line land area (UNION _ XZDWMJ))) deducted land coefficient (TKXS)/100.
When the subtraction ratio land class coefficient (TKXSKCBL) < 1: expressed as a percentage of denominator 1;
deducted land area (TKMJ) ((land pattern spot area (UNION _ TBMJ) -line land area (UNION _ XZDWMJ))) deducted land coefficient (TKXS).
It should be noted that, in this embodiment, the subtracted ground area is taken into consideration, where when the subtraction type of the subtracted ground area is "TK", the subtracted ground area represents the subtracted sill area; when the subtraction type is not "TK", the subtracted area represents a proportionally subtracted hashed area of other places, and in the embodiment of the present application, the subtracted area may be specifically a field area.
Correspondingly, calculating and acquiring the actual area of the target area according to the superposition result of the ground type map spots, the superposition result of the linear ground object and the updated deduction relationship, wherein the calculation comprises the following steps: and calculating and obtaining the actual area of the target area pattern spot land according to the land pattern spot superposition result, the linear land feature superposition result, the deduction land area and the updated deduction relation.
Wherein, the area of the land type image spot can be calculated according to the superposition result of the land type image spot; and calculating the area of the linear ground object according to the linear ground object superposition result and the updated deduction relation.
Alternatively, the actual area of the target area spot land (TBDLMJ) ═ land area of the land (UNION _ TBMJ) -linear land area (UNION _ XZDWMJ) -deducted land area (TKMJ) fig. 3 is a land area obtaining device provided by an embodiment of the present application, which may be a computer, a server, or other computing device, and corresponds to the land area obtaining method, and as shown in fig. 3, the device includes: a superposition unit 301, an acquisition unit 302, and a calculation unit 303, wherein;
the superimposing unit 301 is configured to superimpose the spatial analysis layer of the target area and the map-based spot layer corresponding to the target area, so as to obtain an area superimposing result.
An obtaining unit 302, configured to obtain a corresponding ground type pattern spot superposition result and a linear ground object superposition result according to the area superposition result; and acquiring an updated deduction relation according to the linear ground object superposition result.
And the calculating unit 303 is configured to calculate and obtain an actual area of the target region according to the land-like pattern spot superposition result, the linear ground object superposition result, and the updated deduction relationship.
Optionally, the obtaining unit 302 is configured to intersect the region overlapping result with the linear ground object to obtain an intersection result; and acquiring a linear ground object superposition result according to the spatial relationship of the linear ground objects in the intersection result.
Optionally, the calculating unit 303 is configured to obtain the deduction times of the linear feature according to the linear feature superposition result; recalculating the actual area of the linear ground object in the target area according to the deduction times; calculating and acquiring the actual area of the ground type image spots in the target area according to the superposition result of the ground type image spots and a preset calculation formula; and calculating and acquiring the actual area of the target area according to the actual area of the linear ground object in the target area and the actual area of the ground class pattern spot in the target area.
Fig. 4 is a land area acquiring device according to another embodiment of the present application; as shown in fig. 4, the apparatus further includes: a backfill unit 304, wherein:
and the backfilling unit 304 is configured to backfill the actual total area of the linear feature into the region overlapping result to obtain a backfilling result.
An obtaining unit 302, configured to obtain an actual deducted land area in the target region according to the actual area of the land pattern, the backfill result, and the known deducted land coefficient.
Optionally, the calculating unit 303 is configured to calculate and obtain an actual area of the map spot land of the target area according to the land pattern spot superposition result, the linear land feature superposition result, the deducted land pattern area, and the updated deduction relationship.
Fig. 5 is a schematic structural diagram of an apparatus for acquiring an earth area according to an embodiment of the present application, which corresponds to an earth area acquiring method; the land area obtaining apparatus may include: a processor 510, a storage medium 520, and a bus 530, wherein the storage medium 520 stores machine-readable instructions executable by the processor 510, and when the above-mentioned apparatus is operated, the processor 510 communicates with the storage medium 520 via the bus 530, and the processor 510 executes the machine-readable instructions to perform the steps of the above-mentioned method embodiment. The specific implementation and technical effects are similar, and are not described in detail herein.
The embodiment of the application provides a storage medium, wherein a computer program is stored on the storage medium, and the computer program is executed by a processor to execute the method.
In the several embodiments provided in the present application, it should be understood that the disclosed apparatus and method may be implemented in other ways. For example, the above-described apparatus embodiments are merely illustrative, and for example, the division of the units is only one logical division, and other divisions may be realized in practice, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may also be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present application may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit can be realized in a form of hardware, and can also be realized in a form of hardware and a software functional unit.
It is noted that relational terms such as "first" and "second," and the like, may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. And, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
The above description is only a preferred embodiment of the present application and is not intended to limit the present application, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the present application shall be included in the protection scope of the present application. It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures. The above description is only a preferred embodiment of the present application and is not intended to limit the present application, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the present application shall be included in the protection scope of the present application.