CN115511947B - Real-time accurate acre measurement and correction method for land parcels - Google Patents

Real-time accurate acre measurement and correction method for land parcels Download PDF

Info

Publication number
CN115511947B
CN115511947B CN202211229682.4A CN202211229682A CN115511947B CN 115511947 B CN115511947 B CN 115511947B CN 202211229682 A CN202211229682 A CN 202211229682A CN 115511947 B CN115511947 B CN 115511947B
Authority
CN
China
Prior art keywords
value
acre
convex polygon
virtual
area
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202211229682.4A
Other languages
Chinese (zh)
Other versions
CN115511947A (en
Inventor
董立柱
侯祥英
乌云塔娜
谢杨青
路伟
张长明
陶言民
曲盛林
刘健
孙杰
常艳
张培培
张晓鸥
张鸿川
王笑笑
杨丽华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong Kexiang Intelligent Technology Co ltd
Original Assignee
Shandong Kexiang Intelligent Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong Kexiang Intelligent Technology Co ltd filed Critical Shandong Kexiang Intelligent Technology Co ltd
Priority to CN202211229682.4A priority Critical patent/CN115511947B/en
Publication of CN115511947A publication Critical patent/CN115511947A/en
Application granted granted Critical
Publication of CN115511947B publication Critical patent/CN115511947B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/64—Analysis of geometric attributes of convexity or concavity
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30181—Earth observation
    • G06T2207/30188—Vegetation; Agriculture
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Geometry (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Evolutionary Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a land parcel real-time accurate mu measurement and correction method, and belongs to the technical field of measurement. The invention is obtained by
Figure 464657DEST_PATH_IMAGE001
Value and area
Figure 178535DEST_PATH_IMAGE002
The mapping relation between the two land parcel measuring models can be quickly solved based on the mapping relation
Figure 302480DEST_PATH_IMAGE002
Then by polygonal for each virtual space
Figure 166531DEST_PATH_IMAGE003
Corresponding area prediction value
Figure 839958DEST_PATH_IMAGE002
Summing and taking the sum value
Figure 131262DEST_PATH_IMAGE004
To verify the first land parcel survey pattern to obtain
Figure 438441DEST_PATH_IMAGE005
The value being obtained by solving independent variables
Figure 496396DEST_PATH_IMAGE006
The value is
Figure 712744DEST_PATH_IMAGE007
By the verification mode, the accuracy of mu measurement of the convex polygon area is greatly improved. Moreover, only the mu measurement result of the convex polygon area needs to be obtained when the mu measurement result is predicted
Figure 174950DEST_PATH_IMAGE001
Value sum
Figure 788334DEST_PATH_IMAGE005
The value is enough, and the calculation speed is greatly improved. And the expansion or contraction of the convex polygon area can be dependent through the calculation of the vertex deviation, so that the regularity of the shape of the expanded or contracted area is ensured.

Description

Real-time accurate acre measurement and correction method for land parcels
Technical Field
The invention relates to the technical field of measurement, in particular to a real-time accurate acre measurement and correction method for land parcels.
Background
In the scenes of farmland division, land measurement and the like, whether a certain area or a plurality of areas on the land reach a mu or not needs to be determined, in the existing method, a manual mu measuring mode is generally adopted, namely, personnel measure mu along the boundary of the areas in a tape pulling mode, the manual mode is time-consuming and labor-consuming to operate, and when the number of areas needing to be measured is large, the defect is more prominent. Moreover, the area of the area to be measured is usually in an irregular shape, which may be a rectangle, a pentagon, a hexagon, etc., after the length of each boundary of the area to be measured with an irregular shape is measured by a tape, when the area is calculated, the area to be measured is also required to be correspondingly divided, such as a combination of a regular rectangle and a triangle, or a combination of a trapezoid or a triangle, etc., then the area of each regular shape is calculated and summed to obtain the area of the area to be measured, but the artificial division of the area to be measured into the combination of regular shapes takes longer time, the artificial area calculation is performed on each regular shape after measuring the boundary length of each regular shape one by one, which is very complex and troublesome.
Moreover, if the measured area of the mu to be measured does not reach the area of one mu or exceeds the area of one mu to be measured and needs to be contracted, the area of the mu to be measured with irregular shape is usually expected to be as regular as possible after the area of one mu to be measured is expanded or contracted, and the influence on the other areas outside the area of mu to be measured is reduced as much as possible, but a measurer on site does not have a high-altitude visual angle, so that irregular points or boundaries in the area of mu to be measured cannot be rapidly judged, the point or boundaries from which to start to be expanded or contracted are not known, and whether the point or boundaries in the area of mu to be measured lack scientific judgment basis is not regular or not, so that the area of mu to be measured after the area of mu to be measured is artificially expanded or contracted by a tape is difficult to reach the expected value of ' regular ' shape ' of people.
Disclosure of Invention
The invention aims to simplify the acre measurement calculation and area correction process, improve acre measurement precision and acre measurement and area correction automation degree, and provides a land block real-time accurate acre measurement and correction method.
To achieve the purpose, the invention adopts the following technical scheme:
the method for accurately measuring acre and correcting land parcels in real time comprises the following steps:
S1, drawing a corresponding virtual convex polygon of a convex polygon area determined on a land block to be measured
Figure 805755DEST_PATH_IMAGE001
；
S2, predicting the virtual convex polygon
Figure 199827DEST_PATH_IMAGE001
Area of (2)
Figure 266878DEST_PATH_IMAGE002
；
S3, calculating
Figure 797217DEST_PATH_IMAGE002
Difference from one mu of area
Figure 409464DEST_PATH_IMAGE003
And calculates the difference
Figure 166198DEST_PATH_IMAGE003
Is expressed as the absolute value of
Figure 663038DEST_PATH_IMAGE004
；
S4, judging the difference value
Figure 172517DEST_PATH_IMAGE003
Whether or not it is negative，
If yes, the area expansion flow is switched to;
if not, turning to an area contraction flow;
the area expansion flow comprises the following steps:
m1, calculating the virtual convex polygon
Figure 831032DEST_PATH_IMAGE001
Each vertex on
Figure 308018DEST_PATH_IMAGE005
Degree of deviation of (2)
Figure 467604DEST_PATH_IMAGE006
The calculation method is expressed by the following formula (15):
Figure 238114DEST_PATH_IMAGE007
in the formula (15) of the present invention,
Figure 942896DEST_PATH_IMAGE008
representing vertices
Figure 798857DEST_PATH_IMAGE005
To the virtual convex polygon
Figure 762133DEST_PATH_IMAGE001
Other vertices on
Figure 855991DEST_PATH_IMAGE009
A linear distance therebetween;
Figure 495789DEST_PATH_IMAGE010
representing the virtual convex polygon
Figure 573466DEST_PATH_IMAGE001
The number of vertices thereon;
m2, select the degree of deviation
Figure 74855DEST_PATH_IMAGE006
The vertex with the smallest value is taken as the extension initial vertex
Figure 429744DEST_PATH_IMAGE011
And will be connected with
Figure 866542DEST_PATH_IMAGE011
The vertex with relatively small deviation degree of two adjacent vertexes is used as an extension initial vertex
Figure 290570DEST_PATH_IMAGE012
Then obtain
Figure 471015DEST_PATH_IMAGE011
、
Figure 919226DEST_PATH_IMAGE012
The coordinates of the vertices are respectively noted as
Figure 651559DEST_PATH_IMAGE013
、
Figure 703828DEST_PATH_IMAGE014
，
M3 according to the coordinates
Figure 297752DEST_PATH_IMAGE013
、
Figure 751867DEST_PATH_IMAGE014
Computing connections
Figure 655101DEST_PATH_IMAGE011
And
Figure 194667DEST_PATH_IMAGE012
the length of the first line of the vertex is recorded as
Figure 90816DEST_PATH_IMAGE015
；
M4, according to absolute value
Figure 399438DEST_PATH_IMAGE004
、
Figure 473573DEST_PATH_IMAGE015
Calculating the respective slave
Figure 969276DEST_PATH_IMAGE011
、
Figure 170582DEST_PATH_IMAGE012
The vertexes being connected at right angles to
Figure 333710DEST_PATH_IMAGE011
And
Figure 578746DEST_PATH_IMAGE012
the height of two extension lines extending in the direction of the first straight line between the vertexes
Figure 92904DEST_PATH_IMAGE016
So that
Figure 330857DEST_PATH_IMAGE011
Vertex, slave
Figure 348491DEST_PATH_IMAGE011
Termination point of vertex extension
Figure 764429DEST_PATH_IMAGE017
From the
Figure 375670DEST_PATH_IMAGE012
Termination point of vertex extension
Figure 43412DEST_PATH_IMAGE018
、
Figure 774607DEST_PATH_IMAGE012
The area of the rectangle formed by enclosing between 4 points of the vertex is equal to the absolute value solved in the step S3
Figure 502392DEST_PATH_IMAGE004
The method comprises the steps of carrying out a first treatment on the surface of the The area contraction flow comprises the following steps:
n1, calculating the virtual convex polygon
Figure 99464DEST_PATH_IMAGE001
Each vertex on
Figure 305318DEST_PATH_IMAGE005
Degree of deviation of (2)
Figure 156599DEST_PATH_IMAGE006
The calculation is expressed by the following formula (16):
Figure 930651DEST_PATH_IMAGE019
in the formula (16) of the present invention,
Figure 375539DEST_PATH_IMAGE008
representing vertices
Figure 509717DEST_PATH_IMAGE005
To the virtual convex polygon
Figure 622029DEST_PATH_IMAGE001
The apex on
Figure 528500DEST_PATH_IMAGE009
A linear distance therebetween;
Figure 726263DEST_PATH_IMAGE010
representing the virtual convex polygon
Figure 398553DEST_PATH_IMAGE001
The number of vertices thereon;
n2, select degree of deviation
Figure 365372DEST_PATH_IMAGE006
Maximum valueVertex point
Figure 215648DEST_PATH_IMAGE005
As a first initial vertex of area contraction, and selecting the virtual convex polygon
Figure 900707DEST_PATH_IMAGE001
Any one of a first adjacent vertex and a second adjacent vertex adjacent to the first initial vertex is taken as a second initial vertex with area shrinkage;
n3, obtaining the coordinates of the first initial vertex and the second initial vertex, calculating the length of a second straight line connected between the first initial vertex and the second initial vertex according to the obtained coordinates of the first initial vertex and the second initial vertex, and marking as
Figure 376688DEST_PATH_IMAGE020
The second straight line is used as the waist of an isosceles triangle to be contracted;
N4, at the absolute value
Figure 306335DEST_PATH_IMAGE004
As the area of the isosceles triangle to be contracted, and according to the length of the second straight line
Figure 717725DEST_PATH_IMAGE020
And the coordinates of the first initial vertex and the second initial vertex calculate the base length of the isosceles triangle to be contracted
Figure 749135DEST_PATH_IMAGE021
；
N5, shrinking the first initial vertex toward another adjacent vertex which is not used as the second initial vertex and is adjacent to the first initial vertex
Figure 373014DEST_PATH_IMAGE021
Distance, obtain the contraction point
Figure 189792DEST_PATH_IMAGE022
；
N6, the first initial vertex, the second initial vertex and the contraction point
Figure 506503DEST_PATH_IMAGE022
Removing the enclosed isosceles triangle; s5, expanding or contracting the area of the virtual convex polygon
Figure 25209DEST_PATH_IMAGE001
Mapping to the convex polygon area after area expansion or contraction under the physical space.
Preferably, in step S2, the virtual convex polygon is calculated by the following method steps
Figure 718359DEST_PATH_IMAGE001
Area of (2)
Figure 622599DEST_PATH_IMAGE002
：
L1, for the virtual convex polygon
Figure 500425DEST_PATH_IMAGE001
Performing similarity matching with each virtual convex polygon in the virtual convex polygon database,
if the matching is successful, the step L2 is carried out;
if the matching fails, terminating the mu measuring flow of the convex polygon area;
l2, the virtual convex polygon
Figure 381793DEST_PATH_IMAGE001
Discrete into a plurality of the virtual space polygons
Figure 488421DEST_PATH_IMAGE023
Then the virtual convex polygon is obtained
Figure 138845DEST_PATH_IMAGE001
Corresponding to
Figure 921993DEST_PATH_IMAGE024
Data pairs and acquisition of each of the virtual space polygons
Figure 290658DEST_PATH_IMAGE023
Corresponding to
Figure 699511DEST_PATH_IMAGE025
The data pair is used to determine the data pair,
Figure 63496DEST_PATH_IMAGE026
、
Figure 158491DEST_PATH_IMAGE010
respectively represent the virtual convex polygon
Figure 889818DEST_PATH_IMAGE001
Is the central site of (2)
Figure 462882DEST_PATH_IMAGE027
The distance average value and the number of the vertexes;
Figure 212532DEST_PATH_IMAGE028
representing each of the virtual space polygons separately
Figure 212849DEST_PATH_IMAGE023
Is the central site of (2)
Figure 670287DEST_PATH_IMAGE029
The distance average value and the number of the vertexes;
l3, will
Figure 312621DEST_PATH_IMAGE026
Input to
Figure 916778DEST_PATH_IMAGE010
In the corresponding first land block mu measuring model, the model outputs the corresponding
Figure 963362DEST_PATH_IMAGE030
The value is recorded as
Figure 528336DEST_PATH_IMAGE002
The method comprises the steps of carrying out a first treatment on the surface of the Will be
Figure 567836DEST_PATH_IMAGE031
Input to
Figure 167445DEST_PATH_IMAGE032
In the corresponding second land block acre measuring model, the model outputs the corresponding
Figure 149045DEST_PATH_IMAGE030
The value is recorded as
Figure 935735DEST_PATH_IMAGE033
；
L4, the pair is opposite to the virtual convex polygon
Figure 778926DEST_PATH_IMAGE001
Each of the virtual space polygons having a discrete relationship
Figure 233041DEST_PATH_IMAGE023
Corresponding area prediction value
Figure 887008DEST_PATH_IMAGE033
Summing to obtain a sum value
Figure 426574DEST_PATH_IMAGE034
；
L5, judging the sum value
Figure 73456DEST_PATH_IMAGE034
And (3) with
Figure 382077DEST_PATH_IMAGE002
Whether the absolute value of the difference is smaller than a preset first difference threshold,
if yes, will
Figure 439901DEST_PATH_IMAGE002
The value is used as a mu measurement result of the convex polygon area;
if not, form
Figure 325817DEST_PATH_IMAGE035
Data pair addition to
Figure 651756DEST_PATH_IMAGE010
Corresponding first fitting point set and form
Figure 690251DEST_PATH_IMAGE036
Data pair addition to
Figure 810653DEST_PATH_IMAGE032
And a corresponding second fitting point set.
Preferably, in step L2, the virtual convex polygon is obtained
Figure 183866DEST_PATH_IMAGE001
Corresponding said
Figure 47917DEST_PATH_IMAGE024
The method steps of the data pair comprise:
a1, respectively installing a distance sensor at each vertex of the convex polygon area, and recording that each installation position is at
Figure 439453DEST_PATH_IMAGE037
The coordinates in the axis coordinate system are respectively marked as
Figure 855390DEST_PATH_IMAGE038
,
Figure 935473DEST_PATH_IMAGE039
Representing the number of vertices of the convex polygon area, and then drawing the virtual convex polygon corresponding to the convex polygon area on a computer according to the recorded coordinates
Figure 337636DEST_PATH_IMAGE001
；
A2, according to the virtual convex polygon
Figure 334410DEST_PATH_IMAGE001
Calculating the coordinates of each vertex of the virtual convex polygon
Figure 62195DEST_PATH_IMAGE001
Is defined by the central site of (2)
Figure 145687DEST_PATH_IMAGE027
At the position of
Figure 741753DEST_PATH_IMAGE037
The coordinates in the axis coordinate system are noted as
Figure 812609DEST_PATH_IMAGE040
The central site is then calculated
Figure 304770DEST_PATH_IMAGE027
And the virtual convex polygon where the virtual convex polygon is located
Figure 15237DEST_PATH_IMAGE001
The average value of the distance between each vertex is recorded as
Figure 398683DEST_PATH_IMAGE026
Thereby obtaining the described
Figure 245416DEST_PATH_IMAGE024
And (3) data pairs.
Preferably, each of the virtual space polygons is acquired
Figure 439637DEST_PATH_IMAGE023
Corresponding said
Figure 716029DEST_PATH_IMAGE025
The method of the data pair further comprises the following steps on the basis of the steps A1-A2:
a3, for the virtual convex polygon
Figure 529264DEST_PATH_IMAGE001
Is equally divided and each equally divided point is calculated
Figure 355137DEST_PATH_IMAGE041
At the position of
Figure 330047DEST_PATH_IMAGE037
Coordinates in an axial coordinate system, and dividing each of the equal points
Figure 389007DEST_PATH_IMAGE041
And the central site
Figure 740354DEST_PATH_IMAGE027
After connecting the lines, the virtual convex polygon
Figure 420734DEST_PATH_IMAGE001
Discrete into a plurality of the virtual space polygons
Figure 441911DEST_PATH_IMAGE023
，
Figure 348687DEST_PATH_IMAGE042
；
A4, calculating each virtual space polygon
Figure 628359DEST_PATH_IMAGE023
Each vertex coordinate on the model is respectively marked as
Figure 146933DEST_PATH_IMAGE043
，
Figure 729225DEST_PATH_IMAGE032
Representing the virtual space polygon
Figure 247931DEST_PATH_IMAGE023
Is the number of vertices of (a);
a5, according to the virtual space polygon
Figure 675501DEST_PATH_IMAGE023
Coordinates of each vertex on the virtual space polygon are calculated
Figure 346785DEST_PATH_IMAGE023
Is defined by the central site of (2)
Figure 99977DEST_PATH_IMAGE029
At the position of
Figure 574821DEST_PATH_IMAGE037
The coordinates in the axis coordinate system are noted as
Figure 451422DEST_PATH_IMAGE044
The central site is then calculated
Figure 101846DEST_PATH_IMAGE029
And the virtual space polygon where it is
Figure 150574DEST_PATH_IMAGE023
The average value of the distance between each vertex is recorded as
Figure 253659DEST_PATH_IMAGE031
Thereby obtaining an association of each of the virtual space polygons
Figure 163977DEST_PATH_IMAGE023
Is described in (2)
Figure 668908DEST_PATH_IMAGE025
And (3) data pairs.
Preferably, in step L3, the method steps of constructing the first land parcel measuring model and the second land parcel measuring model include:
c1, obtaining virtual convex polygons which are respectively corresponding to a plurality of convex polygon areas with different shapes and are determined on a land block with measured acre
Figure 888537DEST_PATH_IMAGE001
Associated with
Figure 478918DEST_PATH_IMAGE045
A data pair, wherein,
Figure 691462DEST_PATH_IMAGE046
representing the virtual convex polygon
Figure 316479DEST_PATH_IMAGE001
Is a part of the area of (2);
c2, each virtual convex polygon
Figure 441429DEST_PATH_IMAGE001
Discrete into a plurality of the virtual space polygons
Figure 394473DEST_PATH_IMAGE023
Then, obtaining and each virtual convex polygon
Figure 505649DEST_PATH_IMAGE001
Each of the virtual space polygons having a discrete relationship
Figure 109805DEST_PATH_IMAGE023
Corresponding to
Figure 546603DEST_PATH_IMAGE047
The data pair is used to determine the data pair,
Figure 485478DEST_PATH_IMAGE048
representing the virtual space polygon
Figure 665923DEST_PATH_IMAGE023
Is a part of the area of (2);
C3，
Figure 390166DEST_PATH_IMAGE049
the first land block mu measuring model and the second land block mu measuring model respectively correspond to each otherLand block acre measurement model;
c4, in several forms
Figure 997865DEST_PATH_IMAGE050
In data pairs
Figure 394342DEST_PATH_IMAGE026
As an argument, in a data pair
Figure 503112DEST_PATH_IMAGE046
Solving the first land block acre measurement model for the dependent variable to obtain a first parameter value of a first acre measurement parameter; using several
Figure 691648DEST_PATH_IMAGE051
In data pairs
Figure 844150DEST_PATH_IMAGE031
As an argument, in a data pair
Figure 118136DEST_PATH_IMAGE048
Solving the second land block acre measurement model for the dependent variable to obtain a first parameter value of a second acre measurement parameter;
c5, substituting the first parameter values of the first acre measuring parameters and the second acre measuring parameters into the first land block acre measuring model and the second land block acre measuring model respectively, and substituting each virtual convex polygon into the first land block acre measuring model and the second land block acre measuring model
Figure 765018DEST_PATH_IMAGE001
Corresponding to
Figure 683427DEST_PATH_IMAGE026
The value is input into the first land parcel measuring model, and the corresponding model is output
Figure 632928DEST_PATH_IMAGE030
The value is recorded as
Figure 784424DEST_PATH_IMAGE002
And each of the virtual pairsQuasi-space polygon
Figure 844784DEST_PATH_IMAGE023
Corresponding to
Figure 641533DEST_PATH_IMAGE031
The value is input into the second land parcel measuring model, and the model outputs corresponding
Figure 496357DEST_PATH_IMAGE030
The value is recorded as
Figure 869569DEST_PATH_IMAGE033
；
C6, judging each
Figure 733620DEST_PATH_IMAGE002
The value corresponds to the true value
Figure 892200DEST_PATH_IMAGE046
Whether the absolute value of the difference is smaller than a preset second difference threshold,
if yes, the input and output data pair of the first land parcel measuring model in the step C5 is stored
Figure 42559DEST_PATH_IMAGE035
Adding the first fitting point set as a fitting point;
if not, discarding the input/output data pair
Figure 778433DEST_PATH_IMAGE035
；
Simultaneously judging each of the
Figure 820076DEST_PATH_IMAGE033
The value corresponds to the true value
Figure 957797DEST_PATH_IMAGE048
Whether the absolute value of the difference is smaller than a preset third difference threshold,
if so, save the step C5Input/output data pair of second land block acre measurement model
Figure 544636DEST_PATH_IMAGE036
Added as fitting points to the second set of fitting points,
if not, discarding the input/output data pair
Figure 767807DEST_PATH_IMAGE036
；
C7, opposite to the virtual convex polygon
Figure 114606DEST_PATH_IMAGE001
Each of the virtual space polygons having a discrete relationship
Figure 841253DEST_PATH_IMAGE023
Corresponding model predictive value
Figure 598994DEST_PATH_IMAGE033
Summing to obtain a sum value
Figure 309461DEST_PATH_IMAGE034
；
C8, judging the sum value
Figure 692906DEST_PATH_IMAGE034
With said virtual convex polygon having a discrete relationship therewith
Figure 398694DEST_PATH_IMAGE001
True value of area of (2)
Figure 468281DEST_PATH_IMAGE046
Whether the absolute value of the difference is smaller than a preset fourth difference threshold,
if yes, go to step C9,
if not, filtering the virtual convex polygon from the first fitting point set
Figure 275832DEST_PATH_IMAGE001
The input/output data pair with association relation
Figure 89067DEST_PATH_IMAGE035
And filtering the virtual convex polygon from the second fitting point set
Figure 914940DEST_PATH_IMAGE001
The input/output data pair with association relation
Figure 155429DEST_PATH_IMAGE036
；
C9, fitting each fitting point in the first fitting point set and the second fitting point set respectively through an interpolation method of a Lagrangian interpolation polynomial to obtain a first fitting curve corresponding to the first fitting point set and a second fitting curve corresponding to the second fitting point set;
C10, solving second parameter values of the first acre measuring parameters of the first land parcel acre measuring model according to the first fitting curve, and solving second parameter values of the second acre measuring parameters of the second land parcel acre measuring model according to the second fitting curve;
and C11, updating and correcting second parameter values of the first acre measurement parameters and the second acre measurement parameters, and substituting the updated second parameter values into the corresponding first land parcel acre measurement model or the second land parcel acre measurement model to complete the construction of the first land parcel acre measurement model and the second land parcel acre measurement model.
Preferably, the following steps are continued after step A5 to solve the virtual space polygon
Figure 948810DEST_PATH_IMAGE023
Area of (2)
Figure 34578DEST_PATH_IMAGE048
：
A6, for the virtual convex polygon
Figure 980537DEST_PATH_IMAGE001
Two of said bisecting points on adjacent sides of (a)
Figure 1714DEST_PATH_IMAGE041
After direct connection, each virtual space polygon is obtained
Figure 174069DEST_PATH_IMAGE023
Further discretizing into a plurality of virtual triangles, denoted as
Figure 188162DEST_PATH_IMAGE052
；
A7, according to each virtual triangle
Figure 863994DEST_PATH_IMAGE052
Calculating the area of the vertex coordinates of (a)
Figure 560466DEST_PATH_IMAGE053
And for each of said virtual space polygons
Figure 954538DEST_PATH_IMAGE023
Each of the virtual triangles having a discrete relationship
Figure 772322DEST_PATH_IMAGE052
Area of (2)
Figure 568239DEST_PATH_IMAGE053
Summing, the sum value is taken as the corresponding virtual space polygon
Figure 931219DEST_PATH_IMAGE023
Area of (2)
Figure 671642DEST_PATH_IMAGE048
。
Preferably, the following steps are continued after step A7 to solve the virtual convex polygon
Figure 168482DEST_PATH_IMAGE001
Area of (2)
Figure 192808DEST_PATH_IMAGE046
：
A8, for each virtual space polygon
Figure 116901DEST_PATH_IMAGE023
Area of (2)
Figure 344620DEST_PATH_IMAGE048
Summing, the sum value is taken as the virtual convex polygon
Figure 113993DEST_PATH_IMAGE001
Area of (2)
Figure 759869DEST_PATH_IMAGE046
。
Preferably, in step C3, when
Figure 854864DEST_PATH_IMAGE054
At the time of acquisition of
Figure 569879DEST_PATH_IMAGE010
The corresponding first land parcel acre measurement model is expressed by the following formula (1):
Figure 516845DEST_PATH_IMAGE055
when (when)
Figure 141861DEST_PATH_IMAGE056
At the time of acquisition of
Figure 532391DEST_PATH_IMAGE010
The corresponding first land parcel acre measurement model is expressed by the following formula (2):
Figure 610068DEST_PATH_IMAGE057
when (when)
Figure 596610DEST_PATH_IMAGE058
At the time of acquisition of
Figure 76133DEST_PATH_IMAGE010
The corresponding first land parcel acre measurement model is expressed by the following formula (3):
Figure 637564DEST_PATH_IMAGE059
in the formulas (1) to (3),
Figure 310860DEST_PATH_IMAGE060
、
Figure 491306DEST_PATH_IMAGE061
、
Figure 949969DEST_PATH_IMAGE062
、
Figure 557668DEST_PATH_IMAGE063
、
Figure 219724DEST_PATH_IMAGE064
、
Figure 203861DEST_PATH_IMAGE065
、
Figure 517030DEST_PATH_IMAGE066
、
Figure 295631DEST_PATH_IMAGE067
、
Figure 961097DEST_PATH_IMAGE068
、
Figure 342399DEST_PATH_IMAGE069
、
Figure 916600DEST_PATH_IMAGE070
、
Figure 475889DEST_PATH_IMAGE071
solving each first acre measuring parameter for the parameter value to be calculated;
Figure 502750DEST_PATH_IMAGE072
acquired, acquired
Figure 687744DEST_PATH_IMAGE032
The corresponding second land parcel acre measurement model is expressed by the following formula (4):
Figure 850872DEST_PATH_IMAGE073
in the formula (4) of the present invention,
Figure 345176DEST_PATH_IMAGE074
、
Figure 718389DEST_PATH_IMAGE075
、
Figure 582440DEST_PATH_IMAGE076
、
Figure 475440DEST_PATH_IMAGE077
and solving each second acre measuring parameter for the parameter value to be calculated.
Preferably, in step C9, the method of fitting the first fitted curve or the second fitted curve by using the interpolation method of the lagrangian interpolation polynomial is expressed by the following formula (5):
Figure 766744DEST_PATH_IMAGE078
in the formula (5) of the present invention,
Figure 673258DEST_PATH_IMAGE079
representing the first fitting point set or the second fitting point set
Figure 465634DEST_PATH_IMAGE079
Fitting points;
Figure 603354DEST_PATH_IMAGE080
representing the number of fitting points in the first fitting point set or the second fitting point set;
Figure 675346DEST_PATH_IMAGE081
representing the first land block acre measuring model or the second land block acre measuring model according to the input first land block acre measuring model
Figure 898517DEST_PATH_IMAGE079
Predictive output of individual fitting points
Figure 963425DEST_PATH_IMAGE030
A value;
Figure 329554DEST_PATH_IMAGE082
represents a lagrangian basis function expressed by the following expression (6):
Figure 228239DEST_PATH_IMAGE083
Figure 63340DEST_PATH_IMAGE084
、
Figure 541726DEST_PATH_IMAGE085
representing the first fitting point set or the second fitting point set, respectively
Figure 529405DEST_PATH_IMAGE079
Fitting pointAnd (d)
Figure 333413DEST_PATH_IMAGE086
Fitting points of
Figure 390230DEST_PATH_IMAGE087
The value is
Figure 317647DEST_PATH_IMAGE026
Value or
Figure 284466DEST_PATH_IMAGE031
Values.
Preferably, in step C11, the method for updating and correcting the second parameter value of each of the first acre measurement parameters includes the steps of:
d1, calculating parameter solving error
Figure 649588DEST_PATH_IMAGE088
The calculation method is expressed by the following formula (7):
Figure 803489DEST_PATH_IMAGE089
in the formula (7) of the present invention,
Figure 30202DEST_PATH_IMAGE090
c4, a first parameter value of the first acre measurement parameter obtained by solving in the step is represented;
Figure 851528DEST_PATH_IMAGE091
c10, obtaining a second parameter value of the same first mu measuring parameter by solving;
d2, predicted for step C5
Figure 387551DEST_PATH_IMAGE002
Calculating prediction error
Figure 294328DEST_PATH_IMAGE092
The calculation method is implemented byThe following formula (8) expresses:
Figure 557688DEST_PATH_IMAGE093
d3, judging
Figure 499099DEST_PATH_IMAGE088
Whether or not to follow
Figure 940444DEST_PATH_IMAGE092
Is increased by the increase of (a),
if yes, correcting a second parameter value of the first mu measuring parameter by the following formula (9):
Figure 209883DEST_PATH_IMAGE094
If not, correcting the second parameter value of the first acre measurement parameter by the following formula (10):
Figure 637453DEST_PATH_IMAGE095
in the formulas (9) to (10),
Figure 292425DEST_PATH_IMAGE096
a second parameter value representing the corrected first acre measurement parameter;
in step C11, the method for updating and correcting the second parameter value of each second acre measurement parameter includes the steps of:
f1, calculating parameter solving error
Figure 45618DEST_PATH_IMAGE097
The calculation method is expressed by the following formula (11):
Figure 300887DEST_PATH_IMAGE098
in the formula (11), the color of the sample is,
Figure 656782DEST_PATH_IMAGE099
c4, representing the first parameter value of the second mu measuring parameter obtained by solving in the step;
Figure 307207DEST_PATH_IMAGE100
c10, obtaining a second parameter value of the same second mu measuring parameter by solving;
f2, for each of the predictions calculated in step C5
Figure 841087DEST_PATH_IMAGE033
Calculating a prediction error mean
Figure 209752DEST_PATH_IMAGE101
The calculation method is expressed by the following formula (12):
Figure 369338DEST_PATH_IMAGE102
in the formula (12) of the present invention,
Figure 874268DEST_PATH_IMAGE103
representing the virtual convex polygon
Figure 77585DEST_PATH_IMAGE001
Middle (f)
Figure 933546DEST_PATH_IMAGE104
Each of the virtual space polygons
Figure 896823DEST_PATH_IMAGE023
Is a prediction error of the area of (2);
f3, judge
Figure 866047DEST_PATH_IMAGE097
Whether or not to follow
Figure 131943DEST_PATH_IMAGE101
Is increased by the increase of (a),
if yes, correcting a second parameter value of the second acre measurement parameter by the following formula (13):
Figure 334254DEST_PATH_IMAGE105
if not, correcting a second parameter value of the second acre measurement parameter by the following formula (14):
Figure 976588DEST_PATH_IMAGE106
in the formulas (13) to (14),
Figure 824153DEST_PATH_IMAGE107
A second parameter value representing the corrected second acre measurement parameter.
The invention has the following beneficial effects:
1. the corresponding virtual convex polygon is drawn for the convex polygon area under the physical space, so that the acre measurement work under the real scene is transferred to the computer space, and the automatic acre measurement is possible;
2. the virtual convex polygon is scattered into a plurality of virtual space polygons to serve as acre measuring units, the acre measuring units are thinned, the area of each acre measuring unit is calculated and summed, the obtained sum value serves as the area of the virtual convex polygon, the error of the acre measuring result is reduced, and the acre measuring precision is improved;
3. finding the center site of the virtual convex polygon by constructing a first land block acre measurement model
Figure 260951DEST_PATH_IMAGE027
Average value of distance from each vertex thereof
Figure 950558DEST_PATH_IMAGE026
Mu measurement result of the virtual convex polygon
Figure 865424DEST_PATH_IMAGE046
The mapping relation between the two areas is that the convex polygon area determined on the land block to be measured is obtained
Figure 74820DEST_PATH_IMAGE026
The value can be used for quickly solving the mu measurement result corresponding to the convex polygon area
Figure 682519DEST_PATH_IMAGE046
The method solves the problems that the method for manually dividing an irregular area into a plurality of areas with regular shapes and carrying out area calculation and summation on each regular area is too time-consuming and too complex and the mu measurement error is larger;
4. Finding out the central site of the polygon in the virtual space by constructing a second land block mu measuring model
Figure 593843DEST_PATH_IMAGE029
Average value of distance from each vertex thereof
Figure 686302DEST_PATH_IMAGE031
Mu measurement result of polygon in virtual space
Figure 140417DEST_PATH_IMAGE048
Mapping relation among the virtual space polygons is obtained only by obtaining the corresponding virtual space polygons
Figure 778072DEST_PATH_IMAGE031
The value can be used for quickly solving the mu measurement result corresponding to the polygon in the virtual space
Figure 317637DEST_PATH_IMAGE048
Then, the mu measurement result of each virtual space polygon related to the same virtual convex polygon
Figure 715252DEST_PATH_IMAGE048
Summing and according to the sum and the value
Figure 882928DEST_PATH_IMAGE046
The difference value of the first land block acre measurement model can be used for judging whether the acre measurement result output by the first land block acre measurement model is correct or not, and the second verification of the acre measurement result output by the first land block acre measurement model is realized through the second land block acre measurement model.
5. By setting the first difference threshold, in step S4, when determining the sum value
Figure 832429DEST_PATH_IMAGE034
And (3) with
Figure 967613DEST_PATH_IMAGE002
If the absolute value of the difference is greater than or equal to the preset first difference threshold value
Figure 293552DEST_PATH_IMAGE035
Data pair
Figure 315735DEST_PATH_IMAGE036
The data pairs are added into the corresponding first fitting point set and second fitting point set respectively, so that the number of fitting points is increased, a first fitting curve or a second fitting curve obtained by fitting by a Lagrangian interpolation polynomial interpolation method is smoother, and the obtained second parameter values of each first measured mu parameter or each second measured parameter are more accurate, so that the accuracy of a measured mu result is improved;
6. By setting a second difference threshold, a third difference threshold and a fourth difference threshold to find noise fitting points in the first fitting point set and the second fitting point set, the accuracy of second parameter values corresponding to each first mu measuring parameter and each second mu measuring parameter obtained by reversely pushing a first fitting curve and a second fitting curve obtained by fitting by a subsequent interpolation method using a Lagrange interpolation polynomial is further improved;
7. by searching for
Figure 436138DEST_PATH_IMAGE088
And (3) with
Figure 825662DEST_PATH_IMAGE092
、
Figure 689713DEST_PATH_IMAGE097
And
Figure 831981DEST_PATH_IMAGE101
the parameter correction relation between the two is used for correcting the second parameter value of each first acre measurement parameter or each second acre measurement parameter according to the parameter correction relation, so that the parameter correction accuracy is improved, and the acre measurement accuracy of the convex polygon area is further improved;
8. by calculating virtual convex polygons
Figure 388864DEST_PATH_IMAGE001
Degree of deviation of each vertex on the table
Figure 498641DEST_PATH_IMAGE006
Determining the virtual convex polygon according to the magnitude of the deviation value
Figure 166382DEST_PATH_IMAGE001
The initial vertexes of expansion or contraction are carried out, so that the shape of the convex polygon area in the land block to be measured after the area expansion or contraction is as regular as possible;
9. area predicted by land block acre measurement model
Figure 897578DEST_PATH_IMAGE002
Absolute value of difference from area of one mu
Figure 500729DEST_PATH_IMAGE004
For the area to be expanded or contracted, after the initial vertex of expansion or contraction is determined, the expansion is carried out by using a parallelogram, and the contraction is carried out by using an isosceles triangle, so that the calculation complexity of the expansion and contraction of the area is considered, the calculation speed is improved, and the virtual convex polygon after the expansion or contraction of the area is ensured as much as possible
Figure 723900DEST_PATH_IMAGE001
Is a degree of regularity of (2);
10. for virtual convex polygon in computer
Figure 788808DEST_PATH_IMAGE001
After the area expansion or contraction, the virtual convex polygon is formed
Figure 781034DEST_PATH_IMAGE001
The coordinates of (2) are converted from a virtual space coordinate system to a physical space coordinate system, and the extension termination point position can be obtained
Figure 59481DEST_PATH_IMAGE017
、
Figure 504369DEST_PATH_IMAGE108
Shrinkage point
Figure 372968DEST_PATH_IMAGE022
In the position coordinates in the real physical space, the measuring personnel can rapidly realize the expansion or contraction of the area of the convex polygon area by only finding the two extension termination points or the contraction points on site.
Drawings
In order to more clearly illustrate the technical solution of the embodiments of the present invention, the drawings that are required to be used in the embodiments of the present invention will be briefly described below. It is evident that the drawings described below are only some embodiments of the present invention and that other drawings may be obtained from these drawings without inventive effort for a person of ordinary skill in the art.
Fig. 1 is a step diagram of implementing a method for accurately measuring acre and correcting land parcels in real time according to an embodiment of the present invention;
fig. 2 is a schematic diagram after performing virtual space polygon dispersion on a virtual convex quadrilateral which is drawn correspondingly to a convex quadrilateral region determined on a land parcel to be measured and performing secondary dispersion on each virtual space polygon;
fig. 3 is a schematic diagram after performing virtual space multi-deformation dispersion on a virtual convex pentagon which is drawn corresponding to a convex pentagon region determined on a land area to be measured and performing secondary dispersion on each virtual space multi-deformation;
fig. 4 is a schematic diagram after performing virtual space polygon discretization on a virtual hexagon drawn corresponding to a convex hexagon area determined on a land area to be measured and performing secondary discretization on each virtual space polygon;
fig. 5 is a flowchart of the calculation of the virtual convex polygon area according to the present embodiment.
Detailed Description
The technical scheme of the invention is further described below by the specific embodiments with reference to the accompanying drawings.
Wherein the drawings are for illustrative purposes only and are shown in schematic, non-physical, and not intended to be limiting of the present patent; for the purpose of better illustrating embodiments of the invention, certain elements of the drawings may be omitted, enlarged or reduced and do not represent the size of the actual product; it will be appreciated by those skilled in the art that certain well-known structures in the drawings and descriptions thereof may be omitted.
The same or similar reference numbers in the drawings of embodiments of the invention correspond to the same or similar components; in the description of the present invention, it should be understood that, if the terms "upper", "lower", "left", "right", "inner", "outer", etc. indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings, only for convenience in describing the present invention and simplifying the description, rather than indicating or implying that the apparatus or elements being referred to must have a specific orientation, be constructed and operated in a specific orientation, so that the terms describing the positional relationships in the drawings are merely for exemplary illustration and should not be construed as limiting the present patent, and that the specific meaning of the terms described above may be understood by those of ordinary skill in the art according to specific circumstances.
In the description of the present invention, unless explicitly stated and limited otherwise, the term "coupled" or the like should be interpreted broadly, as it may be fixedly coupled, detachably coupled, or integrally formed, as indicating the relationship of components; can be mechanically or electrically connected; can be directly connected or indirectly connected through an intermediate medium, and can be communication between the two parts or interaction relationship between the two parts. The specific meaning of the above terms in the present invention will be understood in specific cases by those of ordinary skill in the art.
The method for accurately measuring acre and correcting land parcels in real time provided by the embodiment of the invention, as shown in fig. 1, comprises the following steps:
s1, drawing a corresponding virtual convex polygon of a convex polygon area determined on a land block to be measured
Figure 360646DEST_PATH_IMAGE001
；
S2, predicting virtual convex polygon
Figure 164654DEST_PATH_IMAGE001
Area of (2)
Figure 487051DEST_PATH_IMAGE002
；
S3, calculating
Figure 34707DEST_PATH_IMAGE002
Difference from one mu of area
Figure 109848DEST_PATH_IMAGE003
And calculate the difference
Figure 350337DEST_PATH_IMAGE003
Is expressed as the absolute value of
Figure 894451DEST_PATH_IMAGE004
；
S4, judging the difference value
Figure 245798DEST_PATH_IMAGE003
Whether it is negative or not,
if yes, the area expansion flow is switched to;
if not, turning to an area contraction flow;
s5, expanding or contracting the area of the virtual convex polygon
Figure 942489DEST_PATH_IMAGE001
Mapping to the area-expanded or contracted convex polygon region under physical space (by virtual convex polygon
Figure 353879DEST_PATH_IMAGE001
The coordinate conversion of each vertex in the virtual space is realized by converting the coordinate in the physical space, and the specific conversion process applies the existing coordinate conversion method, so that the specific description is not made.
In this embodiment, the key point of real-time accurate acre measurement and correction is to ensure
Figure 119710DEST_PATH_IMAGE002
The speed and accuracy of calculation are calculated, so before the specific technical scheme of the accurate acre measurement and correction method for land parcels is introduced, the method for representing the land parcels is constructed
Figure 383070DEST_PATH_IMAGE026
Value and value of
Figure 324481DEST_PATH_IMAGE046
（
Figure 500247DEST_PATH_IMAGE002
True values of (a) mapping relationship between land parcel survey models and constructing a first land parcel survey model for characterization
Figure 894320DEST_PATH_IMAGE031
Value and value of
Figure 462835DEST_PATH_IMAGE048
The method of the second land parcel acre measurement model of the mapping relation is explained.
In this embodiment, the method steps of constructing the first land parcel measuring model and the second land parcel measuring model include:
c1, obtain under testVirtual convex polygons respectively corresponding to a plurality of convex polygon areas with different shapes determined on each mu of land parcels
Figure 258753DEST_PATH_IMAGE001
Associated with
Figure 871000DEST_PATH_IMAGE045
A data pair, wherein,
Figure 126270DEST_PATH_IMAGE026
、
Figure 357531DEST_PATH_IMAGE109
respectively represent virtual convex polygons
Figure 132589DEST_PATH_IMAGE001
Is the central site of (2)
Figure 791103DEST_PATH_IMAGE027
The mean value of the distance from each vertex, the number of vertices and the area;
Figure 35134DEST_PATH_IMAGE024
the data pair acquisition method comprises the following steps:
a1, respectively installing a distance sensor at each vertex of the convex polygon area, and recording each installation position in
Figure 70086DEST_PATH_IMAGE037
The coordinates in the axis coordinate system are respectively marked as
Figure 699650DEST_PATH_IMAGE038
,
Figure 794645DEST_PATH_IMAGE039
Representing the number of vertexes of the convex polygon area, and then drawing the virtual convex polygon corresponding to the convex polygon area on a computer according to the recorded coordinates
Figure 776506DEST_PATH_IMAGE001
The method comprises the steps of carrying out a first treatment on the surface of the When (when)
Figure 474204DEST_PATH_IMAGE110
When the virtual convex polygon is drawn
Figure 833641DEST_PATH_IMAGE001
As shown in fig. 2, 3 and 4, respectively;
a2, according to the virtual convex polygon
Figure 974903DEST_PATH_IMAGE001
Calculating virtual convex polygons by coordinates of each vertex of (3)
Figure 52581DEST_PATH_IMAGE001
Is the central site of (2)
Figure 819549DEST_PATH_IMAGE027
At the position of
Figure 299071DEST_PATH_IMAGE037
The coordinates in the axis coordinate system are noted as
Figure 109770DEST_PATH_IMAGE040
Then calculate the central site
Figure 674744DEST_PATH_IMAGE027
And the virtual convex polygon where the same is located
Figure 448665DEST_PATH_IMAGE001
The average value of the distance between each vertex is recorded as
Figure 923640DEST_PATH_IMAGE026
Thereby obtaining
Figure 531339DEST_PATH_IMAGE024
And (3) data pairs. Here, the virtual convex polygon is described as
Figure 708242DEST_PATH_IMAGE001
After each vertex coordinate of (2) is determined, a conventional mathematical algorithm can be applied to calculate the central locus
Figure 331859DEST_PATH_IMAGE027
The specific calculation process is not described here.
Figure 599024DEST_PATH_IMAGE046
The calculation process of (2) will be described below and is not yet repeated here.
Obtaining each virtual convex polygon
Figure 502257DEST_PATH_IMAGE001
Corresponding to
Figure 776244DEST_PATH_IMAGE024
After the data pair, the steps are carried out:
c2, each virtual convex polygon
Figure 672394DEST_PATH_IMAGE001
Discretizing into a plurality of virtual space polygons
Figure 981015DEST_PATH_IMAGE023
Then obtain and each virtual convex polygon
Figure 55150DEST_PATH_IMAGE001
Each virtual space polygon having a discrete relationship
Figure 816433DEST_PATH_IMAGE023
Corresponding to
Figure 17738DEST_PATH_IMAGE047
The data pair is used to determine the data pair,
Figure 180866DEST_PATH_IMAGE111
the representation represents each virtual space polygon separately
Figure 160324DEST_PATH_IMAGE023
Is the central site of (2)
Figure 54242DEST_PATH_IMAGE029
The mean value of the distance from each vertex, the number of vertices and the area;
Figure 918293DEST_PATH_IMAGE025
the data pair acquisition method further comprises the following steps on the basis of the steps A1-A2:
a3, for the virtual convex polygon
Figure 60562DEST_PATH_IMAGE001
Is equally divided and each equally divided point is calculated
Figure 351866DEST_PATH_IMAGE041
At the position of
Figure 963107DEST_PATH_IMAGE037
Coordinates in the axial coordinate system (each bisecting point
Figure 365269DEST_PATH_IMAGE041
The coordinates of (2) can be calculated from the coordinates of the two vertices on the bisector segment by using conventional mathematical algorithms, not illustrated herein, and dividing each of the bisectors
Figure 627623DEST_PATH_IMAGE041
With a central site
Figure 355408DEST_PATH_IMAGE027
After connecting the lines, the virtual convex polygon
Figure 421322DEST_PATH_IMAGE001
Discretizing into a plurality of virtual space polygons
Figure 17388DEST_PATH_IMAGE023
，
Figure 9615DEST_PATH_IMAGE042
；
A4, calculating each virtual space polygon
Figure 518088DEST_PATH_IMAGE023
Each vertex coordinate on the model is respectively marked as
Figure 962976DEST_PATH_IMAGE043
(As such, these vertex coordinates may be based on the isocenter coordinates, the center locus
Figure 97154DEST_PATH_IMAGE027
Coordinates and virtual convex polygons
Figure 943887DEST_PATH_IMAGE001
The coordinates of each vertex are calculated by a conventional mathematical algorithm, so the detailed calculation process is not described),
Figure 387375DEST_PATH_IMAGE032
representing virtual space polygons
Figure 178614DEST_PATH_IMAGE023
Is the number of vertices of (a); virtual space polygons
Figure 991849DEST_PATH_IMAGE023
Refer to FIG. 2 for an example
A5, according to the virtual space polygon
Figure 834034DEST_PATH_IMAGE023
Computing virtual space polygons from coordinates of each vertex on the polygon
Figure 808944DEST_PATH_IMAGE023
Is the central site of (2)
Figure 353057DEST_PATH_IMAGE029
At the position of
Figure 704404DEST_PATH_IMAGE037
The coordinates in the axis coordinate system are noted as
Figure 899631DEST_PATH_IMAGE044
(based on the virtual space polygon
Figure 45442DEST_PATH_IMAGE023
The 4 vertex coordinates of (1) are calculated by a conventional mathematical algorithm, the specific calculation process is not described), and then the central position is calculated
Figure 811272DEST_PATH_IMAGE029
And the virtual space polygon where it is
Figure 841676DEST_PATH_IMAGE023
The average value of the distance between each vertex is recorded as
Figure 517508DEST_PATH_IMAGE031
Thereby obtaining the associated each virtual space polygon
Figure 224433DEST_PATH_IMAGE023
Is described in (2)
Figure 352926DEST_PATH_IMAGE025
And (3) data pairs.
Solving virtual space polygons
Figure 414118DEST_PATH_IMAGE023
Area of (2)
Figure 210036DEST_PATH_IMAGE048
The method of (a) is that after the step A5, the following steps are continuously executed:
A6, for the virtual convex polygon
Figure 87862DEST_PATH_IMAGE001
Two bisecting points on adjacent edges of (a)
Figure 313438DEST_PATH_IMAGE041
After direct connection (indicated by a broken line in the direct connection diagram 2), each virtual space polygon is formed
Figure 810278DEST_PATH_IMAGE023
Further discretizing into a plurality of virtual triangles, denoted as
Figure 585336DEST_PATH_IMAGE052
；
A7, according to each virtual triangle
Figure 509430DEST_PATH_IMAGE052
Calculating the area of the vertex coordinates of (a)
Figure 251996DEST_PATH_IMAGE053
And for each virtual space polygon
Figure 21369DEST_PATH_IMAGE023
Each virtual triangle with discrete relationship
Figure 650933DEST_PATH_IMAGE052
Area of (2)
Figure 745928DEST_PATH_IMAGE053
Summing, and taking the obtained sum value as a corresponding virtual space polygon
Figure 211675DEST_PATH_IMAGE023
Area of (2)
Figure 909373DEST_PATH_IMAGE048
。
Through the steps A1-A7, each virtual space polygon is obtained
Figure 534389DEST_PATH_IMAGE023
Corresponding to
Figure 908608DEST_PATH_IMAGE047
And (3) data pairs.
Solving for virtual convex polygonsShape of a Chinese character
Figure 986285DEST_PATH_IMAGE001
Area of (2)
Figure 487674DEST_PATH_IMAGE046
The method of (a) is that after step A7, the steps are continuously executed:
a8, for each virtual space polygon
Figure 967197DEST_PATH_IMAGE023
Area of (2)
Figure 279360DEST_PATH_IMAGE048
Summing, the sum is taken as a virtual convex polygon
Figure 703389DEST_PATH_IMAGE001
Area of (2)
Figure 883834DEST_PATH_IMAGE046
。
Through the steps A1-A8, the virtual convex polygon is obtained
Figure 857344DEST_PATH_IMAGE001
Corresponding to
Figure 465043DEST_PATH_IMAGE045
And (3) data pairs.
Acquisition of
Figure 110788DEST_PATH_IMAGE047
Data pair
Figure 94924DEST_PATH_IMAGE045
After the data pair, the method for constructing the first land parcel measuring model and the second land parcel measuring model is transferred to the steps:
c3, obtain
Figure 424406DEST_PATH_IMAGE112
The first land block mu measuring model and the second land block mu measuring model which correspond respectively, namely
Figure 62060DEST_PATH_IMAGE112
When the values of the land areas are different, the land area measurement models are respectively corresponding to the land area measurement models; as shown for example in figure 3 of the drawings,
Figure 336047DEST_PATH_IMAGE056
、
Figure 972477DEST_PATH_IMAGE072
The method comprises the steps of carrying out a first treatment on the surface of the And in FIG. 4
Figure 281098DEST_PATH_IMAGE054
、
Figure 355234DEST_PATH_IMAGE072
。
Figure 382095DEST_PATH_IMAGE010
The larger the number of the convex polygon area determined on the mu area to be measured is, the more complex the shape of the convex polygon area is, the more the number of virtual space polygons which need to be scattered and the more the number of virtual triangles which need to be further scattered for each virtual space polygon is, the more the number of the discrete virtual space polygons and the number of the virtual triangles is, the higher the complexity of area block calculation on the convex polygon area is represented, the more calculation errors are easy to occur, and in order to balance the area calculation errors and the calculation speed, the embodiment adopts different land block measurement mu models to calculate the areas of the corresponding areas respectively for the virtual convex polygon with different vertexes and the virtual space polygons. In order to ensure the accuracy of area solution, the embodiment correlates the order of the high-order equation with the number of vertices of the virtual convex polygon and the virtual space polygon to obtain the number of vertices of the virtual convex polygon
Figure 317822DEST_PATH_IMAGE010
The orders of the first higher-order square equation serving as the first land parcel measuring model are respectively corresponding to the numbers of the vertexes of the virtual space polygons
Figure 605583DEST_PATH_IMAGE032
The order of the second higher-order equation serving as the corresponding second land parcel measuring mu model is determined.
The expression of the first higher-order equation as the first plot acre model is as follows:
When (when)
Figure 725986DEST_PATH_IMAGE054
At the time of acquisition of
Figure 82887DEST_PATH_IMAGE010
The corresponding first land parcel acre measurement model is expressed by the following formula (1):
Figure 946938DEST_PATH_IMAGE055
when (when)
Figure 354785DEST_PATH_IMAGE056
At the time of acquisition of
Figure 646089DEST_PATH_IMAGE010
The corresponding first land parcel acre measurement model is expressed by the following formula (2):
Figure 991751DEST_PATH_IMAGE057
when (when)
Figure 784127DEST_PATH_IMAGE058
At the time of acquisition of
Figure 921847DEST_PATH_IMAGE010
The corresponding first land parcel acre measurement model is expressed by the following formula (3):
Figure 492375DEST_PATH_IMAGE059
in the formulas (1) to (3),
Figure 715545DEST_PATH_IMAGE060
、
Figure 311612DEST_PATH_IMAGE061
、
Figure 913626DEST_PATH_IMAGE062
、
Figure 546732DEST_PATH_IMAGE063
、
Figure 381833DEST_PATH_IMAGE064
、
Figure 125798DEST_PATH_IMAGE065
、
Figure 346433DEST_PATH_IMAGE066
、
Figure 416020DEST_PATH_IMAGE067
、
Figure 207259DEST_PATH_IMAGE068
、
Figure 895860DEST_PATH_IMAGE069
、
Figure 597100DEST_PATH_IMAGE070
、
Figure 962222DEST_PATH_IMAGE071
and solving each first mu measuring parameter for the parameter value to be calculated.
It should be noted here that,
Figure 116123DEST_PATH_IMAGE010
the values of (2) are not limited to 4, 5, 6, but may be other values, where 4, 5, 6 are defined for convenience of expression
Figure 858949DEST_PATH_IMAGE010
Corresponding first placeA block-measuring mu model is adopted,
Figure 680275DEST_PATH_IMAGE113
when the corresponding first land measurement mu model is provided in the expression provided in the formula (1)
Figure 950719DEST_PATH_IMAGE114
And then add
Figure 732861DEST_PATH_IMAGE115
，
Figure 887899DEST_PATH_IMAGE116
At the time of
Figure 688365DEST_PATH_IMAGE115
After adding
Figure 5077DEST_PATH_IMAGE117
And so on, no further description is given.
The expression of the second higher-order equation as the second plot acre model is as follows:
Figure 773050DEST_PATH_IMAGE072
acquired, acquired
Figure 466200DEST_PATH_IMAGE032
The corresponding second plot acre measurement model is expressed by the following formula (4):
Figure 121172DEST_PATH_IMAGE073
in the formula (4) of the present invention,
Figure 484151DEST_PATH_IMAGE074
、
Figure 99941DEST_PATH_IMAGE075
、
Figure 721415DEST_PATH_IMAGE076
、
Figure 371839DEST_PATH_IMAGE077
and solving each second mu measuring parameter for the parameter value to be calculated.
Likewise, the number of the cells to be processed,
Figure 404255DEST_PATH_IMAGE032
not limited to the case of =4,
Figure 631974DEST_PATH_IMAGE032
when=5, the expression is formula (2),
Figure 401347DEST_PATH_IMAGE032
when=6, the expression is formula (1), and so on, and will not be described again.
Obtained is obtained
Figure 47223DEST_PATH_IMAGE112
After the corresponding models are respectively constructed, the method steps of constructing a first land parcel measuring model and a second land parcel measuring model are transferred to the steps of:
c4, in several forms
Figure 142218DEST_PATH_IMAGE050
In data pairs
Figure 591654DEST_PATH_IMAGE026
As an argument, in a data pair
Figure 804198DEST_PATH_IMAGE046
Solving a first land block acre measurement model for the dependent variable to obtain a first parameter value of a first acre measurement parameter; using several
Figure 429214DEST_PATH_IMAGE051
In data pairs
Figure 819744DEST_PATH_IMAGE031
As an argument, in a data pair
Figure 631843DEST_PATH_IMAGE048
Solving a second land block acre measurement model for the dependent variable to obtain a first parameter value of a second acre measurement parameter;
for example, for the first land parcel measuring model expressed by formulas (1) - (3), as long as
Figure 883964DEST_PATH_IMAGE050
The number of the data pairs is enough, so that first parameter values corresponding to the first acre measurement parameters in the formula can be solved, and the solving process is a conventional mathematical algorithm, so that no exchange is carried out in the specific calculating process.
C5, substituting the first parameter values of the first acre measuring parameters and the second acre measuring parameters into the first land block acre measuring model and the second land block acre measuring model respectively, and substituting each virtual convex polygon
Figure 363486DEST_PATH_IMAGE001
Corresponding to
Figure 659339DEST_PATH_IMAGE026
The value is input into a first land block mu measuring model, and the model outputs corresponding data
Figure 604073DEST_PATH_IMAGE030
The value is recorded as
Figure 518939DEST_PATH_IMAGE002
And polygonal each virtual space
Figure 977602DEST_PATH_IMAGE023
Corresponding to
Figure 585301DEST_PATH_IMAGE031
The value is input into a second land block mu measuring model, and the model outputs corresponding data
Figure 512937DEST_PATH_IMAGE030
The value is recorded as
Figure 356128DEST_PATH_IMAGE033
；
C6, judging each
Figure 544664DEST_PATH_IMAGE002
The value corresponds to the true value
Figure 697166DEST_PATH_IMAGE046
Whether the absolute value of the difference is smaller than a preset second difference threshold,
if yes, the input and output data pair of the first land parcel measuring model in the step C5 is saved
Figure 236731DEST_PATH_IMAGE035
Adding the first fitting point set as a fitting point;
if not, discarding the input/output data pair
Figure 883613DEST_PATH_IMAGE035
；
Simultaneously judging each
Figure 926656DEST_PATH_IMAGE033
The value corresponds to the true value
Figure 751523DEST_PATH_IMAGE048
Whether the absolute value of the difference is smaller than a preset third difference threshold,
if yes, the input and output data pair of the second land parcel measuring model in the step C5 is stored
Figure 637440DEST_PATH_IMAGE036
Added as fitting points to the second set of fitting points,
if not, discarding the input/output data pair
Figure 963379DEST_PATH_IMAGE036
；
C7, pairing with virtual convex polygon
Figure 234829DEST_PATH_IMAGE001
Each virtual space polygon having discrete relationship
Figure 620811DEST_PATH_IMAGE023
Corresponding model predictive value
Figure 994024DEST_PATH_IMAGE033
Summing to obtain a sum value
Figure 592495DEST_PATH_IMAGE034
；
C8, judging the sum value
Figure 751075DEST_PATH_IMAGE034
With virtual convex polygons having discrete relationships therewith
Figure 42379DEST_PATH_IMAGE001
True value of area of (2)
Figure 902888DEST_PATH_IMAGE046
Whether the absolute value of the difference is smaller than a preset fourth difference threshold,
if yes, go to step C9,
if not, filtering out the virtual convex polygon from the first fitting point set
Figure 678952DEST_PATH_IMAGE001
Input/output data pair with association relation
Figure 816672DEST_PATH_IMAGE035
And filtering out the virtual convex polygon from the second fitting point set
Figure 403511DEST_PATH_IMAGE001
Input/output data pair with association relation
Figure 361103DEST_PATH_IMAGE036
；
C9, fitting each fitting point in the first fitting point set and the second fitting point set respectively through an interpolation method of a Lagrangian interpolation polynomial to obtain a first fitting curve corresponding to the first fitting point set and a second fitting curve corresponding to the second fitting point set;
specifically, the method for fitting the first fitted curve or the second fitted curve by using the interpolation method of the lagrangian interpolation polynomial is expressed by the following formula (5):
Figure 442322DEST_PATH_IMAGE078
in the formula (5) of the present invention,
Figure 434549DEST_PATH_IMAGE079
representing the first fitting point set or the second fitting point set
Figure 457869DEST_PATH_IMAGE079
Fitting points;
Figure 270799DEST_PATH_IMAGE080
representing the number of fitting points in the first fitting point set or the second fitting point set;
Figure 14764DEST_PATH_IMAGE081
representing the first land block acre measuring model or the second land block acre measuring model according to the input first land block acre measuring model
Figure 251710DEST_PATH_IMAGE079
Predictive output of individual fitting points
Figure 790139DEST_PATH_IMAGE030
A value;
Figure 863268DEST_PATH_IMAGE082
represents a lagrangian basis function expressed by the following expression (6):
Figure 269979DEST_PATH_IMAGE083
Figure 236798DEST_PATH_IMAGE084
、
Figure 851188DEST_PATH_IMAGE085
representing the first fitting point set or the second fitting point set, respectively
Figure 270668DEST_PATH_IMAGE079
Individual fitting points and the th
Figure 746648DEST_PATH_IMAGE086
Fitting points of
Figure 302395DEST_PATH_IMAGE087
The value is
Figure 589151DEST_PATH_IMAGE026
Value (if)
Figure 495927DEST_PATH_IMAGE079
、
Figure 510019DEST_PATH_IMAGE086
The fitting points in the first fitting point set are
Figure 559752DEST_PATH_IMAGE026
) Or (b)
Figure 142043DEST_PATH_IMAGE031
Value (if)
Figure 395170DEST_PATH_IMAGE079
、
Figure 88320DEST_PATH_IMAGE086
The fitting points in the second fitting point set are
Figure 494024DEST_PATH_IMAGE031
）。
It should be noted that, the fitting data in the first fitting point set or the second fitting point set is substituted into the expression (6) and the expression (5), so as to obtain a corresponding fitting curve, and specific data substitution process is not described here.
After the first fitting curve and the second fitting curve are obtained, the method for constructing the first land parcel measuring model and the second land parcel measuring model is transferred to the steps:
c10, reversely solving second parameter values of all the first acre measuring parameters of the first land parcel acre measuring model according to the first fitting curve, reversely solving second parameter values of all the second acre measuring parameters of the second land parcel acre measuring model according to the second fitting curve;
here, the method of solving the term coefficients (parameters) of the higher-order equation by the inverse of the fitted curve is a conventional mathematical operation method, for example, for a unitary primary equation given a value having the term coefficients, a straight line may be obtained in the XY axis coordinate system, and the term coefficients of the respective terms of the unitary primary equation may be obtained by the inverse of the given straight line in the XY axis coordinate system. Different from the above, in this embodiment, the higher-order equation is used as the first land parcel measuring model and the second land parcel measuring model, when the number of fitting points is insufficient, the curve is not smooth enough, the term coefficient (the second parameter) obtained by back-pushing has a larger error with the first parameter of the same parameter obtained by solving in the step C4, and when the number of fitting points is sufficient, the curve is smoother, and at this time, the term coefficient obtained by back-pushing may be more accurate than the first parameter of the same parameter obtained by solving in the step C4.
However, each parameter in the first land parcel measuring model and the second land parcel measuring model is based on the first parameter value calculated in the step C4 or based on the second parameter value reversely calculated in the step C10, and needs to be supported by a corresponding theoretical basis. Therefore, after step C10, a process for calibrating the model parameters is also required, namely the steps of:
and C11, updating and correcting second parameter values of the first acre measuring parameters and the second acre measuring parameters, and substituting the updated second parameter values into the corresponding first land parcel acre measuring model or second land parcel acre measuring model to complete the construction of the first land parcel acre measuring model and the second land parcel acre measuring model.
In the invention, we find through repeated experimental summary that, whether the first measured mu parameter in the first land parcel measured mu model or the second measured mu parameter in the second land parcel measured mu model, for the first parameter value calculated in the step C4 and the second parameter value reversely obtained in the step C10 of the same parameter, the parameter solving error of the first parameter value and the second parameter value has a corresponding relation on trend with the model prediction error, the corresponding relation on the trend is utilized to correct the second parameter value of each parameter, and finally the corrected model is predicted to output
Figure 247217DEST_PATH_IMAGE030
The value is more similar to the true value. The correction scheme of the present embodiment for each parameter is specifically as follows:
in step C11, the method for updating and correcting the second parameter value of each first mu measuring parameter includes the steps of:
d1, calculating parameter solving error
Figure 253219DEST_PATH_IMAGE088
The calculation method is expressed by the following formula (7):
Figure 750059DEST_PATH_IMAGE089
in the formula (7) of the present invention,
Figure 774385DEST_PATH_IMAGE090
c4, a first parameter value of a first acre measurement parameter obtained by solving is represented;
Figure 291954DEST_PATH_IMAGE091
representing the second parameter value of the same first mu measuring parameter obtained in the step C10 in a back-pushing way；
D2, predicted for step C5
Figure 660618DEST_PATH_IMAGE002
Calculating prediction error
Figure 570937DEST_PATH_IMAGE092
The calculation method is expressed by the following formula (8):
Figure 341447DEST_PATH_IMAGE093
d3, judging
Figure 29917DEST_PATH_IMAGE088
Whether or not to follow
Figure 885877DEST_PATH_IMAGE092
Is increased by the increase of (a),
if yes, correcting the second parameter value of the first mu measuring parameter by the following formula (9):
Figure 104281DEST_PATH_IMAGE094
if not, correcting the second parameter value of the first acre measurement parameter by the following formula (10):
Figure 322773DEST_PATH_IMAGE095
in the formulas (9) to (10),
Figure 588669DEST_PATH_IMAGE096
a second parameter value representing the corrected first acre measurement parameter.
The method for updating and correcting the second parameter value of each second acre measuring parameter comprises the following steps:
f1, calculating parameter solving error
Figure 541713DEST_PATH_IMAGE097
The calculation method is expressed by the following formula (11):
Figure 184047DEST_PATH_IMAGE098
in the formula (11), the color of the sample is,
Figure 522624DEST_PATH_IMAGE099
c4, a first parameter value of a second mu measuring parameter obtained by solving is represented;
Figure 959422DEST_PATH_IMAGE100
C10, reversely pushing the second parameter value of the same second mu measuring parameter;
f2, for each of the predictions calculated in step C5
Figure 898297DEST_PATH_IMAGE033
Calculating a prediction error mean
Figure 78742DEST_PATH_IMAGE101
The calculation method is expressed by the following formula (12):
Figure 271826DEST_PATH_IMAGE102
in the formula (12) of the present invention,
Figure 754891DEST_PATH_IMAGE103
representing a pair of virtual convex polygons
Figure 807161DEST_PATH_IMAGE001
Middle (f)
Figure 650352DEST_PATH_IMAGE104
Multiple virtual space polygons
Figure 104467DEST_PATH_IMAGE023
Is a face of (2)Integrating the prediction error;
f3, judge
Figure 991389DEST_PATH_IMAGE097
Whether or not to follow
Figure 530955DEST_PATH_IMAGE101
Is increased by the increase of (a),
if yes, correcting a second parameter value of the second mu measuring parameter by the following formula (13):
Figure 912258DEST_PATH_IMAGE105
if not, correcting the second parameter value of the second acre measurement parameter by the following formula (14):
Figure 96246DEST_PATH_IMAGE106
in the formulas (13) to (14),
Figure 780168DEST_PATH_IMAGE107
a second parameter value representing the corrected second acre measurement parameter.
By integrating the scheme, the invention completes the construction of the first land parcel measuring model and the second land parcel measuring model, and for the convex polygon area determined on the land parcel to be measured, the measuring result of the convex polygon area can be quickly and automatically solved by applying the two models. The method for accurately measuring acre in real time for the land according to the embodiment of the invention, as shown in fig. 5, comprises the following steps:
l1, for virtual convex polygon
Figure 666084DEST_PATH_IMAGE001
Performing similarity matching with each virtual convex polygon in the virtual convex polygon database,
If the matching is successful, the step L2 is carried out;
if the matching fails, terminating the mu measuring flow of the convex polygon area;
it should be noted that there are many existing methods for matching the regional similarity, so the specific matching process of the regional similarity matching method adopted in the present application is not described here.
L2, the virtual convex polygon
Figure 992023DEST_PATH_IMAGE001
Discrete into a plurality of the virtual space polygons
Figure 263474DEST_PATH_IMAGE023
Then the virtual convex polygon is obtained
Figure 508510DEST_PATH_IMAGE001
Corresponding to
Figure 22668DEST_PATH_IMAGE024
Data pairs and acquisition of each of the virtual space polygons
Figure 762085DEST_PATH_IMAGE023
Corresponding to
Figure 514141DEST_PATH_IMAGE025
The data pair is used to determine the data pair,
Figure 195658DEST_PATH_IMAGE026
、
Figure 931532DEST_PATH_IMAGE010
respectively represent the virtual convex polygon
Figure 725174DEST_PATH_IMAGE001
Is the central site of (2)
Figure 721949DEST_PATH_IMAGE027
Distance average value and vertex number of each vertex;
Figure 449734DEST_PATH_IMAGE028
representing each of said virtual representationsQuasi-space polygon
Figure 548271DEST_PATH_IMAGE023
Is the central site of (2)
Figure 488545DEST_PATH_IMAGE029
An average of the distances from each of its vertices;
it should be noted here how to obtain
Figure 605406DEST_PATH_IMAGE024
Data pair
Figure 504092DEST_PATH_IMAGE025
The data pairs have been described in detail in constructing the first land parcel measuring pattern and the second land parcel measuring pattern, and are not described in detail herein.
L3, will
Figure 322881DEST_PATH_IMAGE026
Input to
Figure 332425DEST_PATH_IMAGE010
In the corresponding first land block mu measuring model, the model outputs the corresponding
Figure 303792DEST_PATH_IMAGE030
The value is recorded as
Figure 983166DEST_PATH_IMAGE002
The method comprises the steps of carrying out a first treatment on the surface of the Will be
Figure 915350DEST_PATH_IMAGE031
Input to
Figure 587640DEST_PATH_IMAGE032
In the corresponding second land block acre measuring model, the model outputs the corresponding
Figure 554459DEST_PATH_IMAGE030
The value is recorded as
Figure 637690DEST_PATH_IMAGE033
；
L4, pairing with virtual convex polygon
Figure 447383DEST_PATH_IMAGE001
Each virtual space polygon having discrete relationship
Figure 798730DEST_PATH_IMAGE023
Corresponding area prediction value
Figure 495422DEST_PATH_IMAGE033
Summing to obtain a sum value
Figure 906812DEST_PATH_IMAGE034
；
L5, judge the sum value
Figure 672643DEST_PATH_IMAGE034
And (3) with
Figure 562101DEST_PATH_IMAGE002
Whether the absolute value of the difference is smaller than a preset first difference threshold,
if yes, will
Figure 611834DEST_PATH_IMAGE002
The value is used as the mu measurement result of the convex polygon area;
if not, form
Figure 53180DEST_PATH_IMAGE035
Data pair addition to
Figure 181673DEST_PATH_IMAGE010
Corresponding first fitting point set and form
Figure 750189DEST_PATH_IMAGE036
Data pair addition to
Figure 546107DEST_PATH_IMAGE032
And a corresponding second fitting point set.
In addition, the center site is found
Figure 423933DEST_PATH_IMAGE027
、
Figure 305301DEST_PATH_IMAGE029
There are many existing methods of (a) such as to make a certain point inside the convex polygon area as close to the distance of each vertex as possible, and this point can be determined as a central point.
By the scheme, the model rapidly calculates the area of the convex polygon area
Figure 650744DEST_PATH_IMAGE002
When step S4 determines the difference
Figure 425802DEST_PATH_IMAGE003
When negative, the description
Figure 208950DEST_PATH_IMAGE002
If a customer needs to expand the area to one mu and the shape of the area after expansion is expected to be as regular as possible, the embodiment provides the following area expansion scheme, and the method comprises the following steps:
m1, calculating a virtual convex polygon
Figure 951516DEST_PATH_IMAGE001
Each vertex on
Figure 845522DEST_PATH_IMAGE005
Degree of deviation of (2)
Figure 350453DEST_PATH_IMAGE006
The calculation method is expressed by the following formula (15):
Figure 320814DEST_PATH_IMAGE007
In the formula (15) of the present invention,
Figure 911196DEST_PATH_IMAGE008
representing vertices
Figure 608893DEST_PATH_IMAGE005
To a virtual convex polygon
Figure 342232DEST_PATH_IMAGE001
Other vertices on
Figure 608128DEST_PATH_IMAGE009
A linear distance therebetween;
Figure 810439DEST_PATH_IMAGE010
representing a virtual convex polygon
Figure 452773DEST_PATH_IMAGE001
The number of vertices thereon;
m2, select the degree of deviation
Figure 542083DEST_PATH_IMAGE006
The vertex with the smallest value is taken as the extension initial vertex
Figure 978881DEST_PATH_IMAGE011
And will be connected with
Figure 668488DEST_PATH_IMAGE011
The vertex with relatively small deviation degree of two adjacent vertexes is used as an extension initial vertex
Figure 685817DEST_PATH_IMAGE012
Then obtain
Figure 285426DEST_PATH_IMAGE011
、
Figure 17758DEST_PATH_IMAGE012
The coordinates of the vertices are respectively noted as
Figure 804449DEST_PATH_IMAGE013
、
Figure 398372DEST_PATH_IMAGE014
，
M3 according to the coordinates
Figure 852487DEST_PATH_IMAGE013
、
Figure 755721DEST_PATH_IMAGE014
Computing connections
Figure 403609DEST_PATH_IMAGE011
And
Figure 925857DEST_PATH_IMAGE012
the length of the first line of the vertex is recorded as
Figure 359113DEST_PATH_IMAGE015
；
M4, according to absolute value
Figure 308614DEST_PATH_IMAGE004
、
Figure 679684DEST_PATH_IMAGE015
Calculating the respective slave
Figure 5623DEST_PATH_IMAGE011
、
Figure 293385DEST_PATH_IMAGE012
The vertexes being connected at right angles to
Figure 413787DEST_PATH_IMAGE011
And
Figure 36268DEST_PATH_IMAGE012
high of two extension lines extending in the direction of the first straight line between the vertexes
Figure 24952DEST_PATH_IMAGE016
So that
Figure 42587DEST_PATH_IMAGE011
Vertex, slave
Figure 474836DEST_PATH_IMAGE011
Termination point of vertex extension
Figure 210711DEST_PATH_IMAGE017
From the
Figure 471928DEST_PATH_IMAGE012
Termination point of vertex extension
Figure 875228DEST_PATH_IMAGE018
、
Figure 711334DEST_PATH_IMAGE012
The area of the rectangle formed by enclosing between 4 points of the vertex is equal to the absolute value solved in the step S3
Figure 668926DEST_PATH_IMAGE004
。
It should be noted here that the area of the rectangle to be expanded is known, the length of the long side (i.e. the first straight line) is known, the height
Figure 999413DEST_PATH_IMAGE016
Namely, is
Figure 867006DEST_PATH_IMAGE004
And (3) with
Figure 500113DEST_PATH_IMAGE015
From which the slave is determined by dividing the value of (a)
Figure 69634DEST_PATH_IMAGE011
、
Figure 79179DEST_PATH_IMAGE012
The length of the extension line from which the apex extends.
And when step S4 determines the difference
Figure 305673DEST_PATH_IMAGE003
To be positive, explain
Figure 375260DEST_PATH_IMAGE002
If the area is required to be shrunk to one mu and the shape of the area after shrinkage is expected to be as regular as possible, the embodiment provides the following area shrinkage scheme, and the method comprises the following steps:
N1, calculating a virtual convex polygon
Figure 166498DEST_PATH_IMAGE001
Each vertex on
Figure 855100DEST_PATH_IMAGE005
Degree of deviation of (2)
Figure 821919DEST_PATH_IMAGE006
The calculation method is expressed by the following formula (15):
Figure 921462DEST_PATH_IMAGE019
in the formula (16) of the present invention,
Figure 340942DEST_PATH_IMAGE008
representing vertices
Figure 66190DEST_PATH_IMAGE005
To a virtual convex polygon
Figure 887516DEST_PATH_IMAGE001
The apex on
Figure 892381DEST_PATH_IMAGE009
A linear distance therebetween;
Figure 674523DEST_PATH_IMAGE010
representing a virtual convex polygon
Figure 829561DEST_PATH_IMAGE001
The number of vertices thereon;
n2, select degree of deviation
Figure 630027DEST_PATH_IMAGE006
Vertex with maximum value
Figure 212318DEST_PATH_IMAGE005
As a first initial vertex of area contraction, and selecting a virtual convex polygon
Figure 714712DEST_PATH_IMAGE001
Any one of a first adjacent vertex and a second adjacent vertex adjacent to the first initial vertex is taken as a second initial vertex with area shrinkage;
n3, obtaining the coordinates of the first initial vertex and the second initial vertex, calculating the length of a second straight line connected between the first initial vertex and the second initial vertex according to the obtained coordinates of the first initial vertex and the second initial vertex, and marking as
Figure 407862DEST_PATH_IMAGE020
The second straight line is used as the waist of the isosceles triangle to be contracted;
n4, in absolute value
Figure 328413DEST_PATH_IMAGE004
As an isosceles triangle to be contracted and according to the length of the second straight line
Figure 81606DEST_PATH_IMAGE020
And coordinates of the first initial vertex and the second initial vertex, calculating a base length of the isosceles triangle to be contracted
Figure 307182DEST_PATH_IMAGE021
；
Figure 928656DEST_PATH_IMAGE021
Is of the meter(s)The calculation is simple mathematical operation, and the specific calculation process is not described;
N5, shrinking the first initial vertex toward another adjacent vertex adjacent to the second initial vertex
Figure 579080DEST_PATH_IMAGE021
Distance, obtain the contraction point
Figure 877075DEST_PATH_IMAGE022
；
N6, the first initial vertex, the second initial vertex and the contraction point
Figure 980160DEST_PATH_IMAGE022
The enclosed isosceles triangle is removed.
In conclusion, the invention is achieved by
Figure 139746DEST_PATH_IMAGE031
Value and area
Figure 254464DEST_PATH_IMAGE033
The mapping relation between the two land parcel measuring models can be quickly solved based on the mapping relation
Figure 349459DEST_PATH_IMAGE033
Then by polygonal for each virtual space
Figure 64474DEST_PATH_IMAGE023
Corresponding area prediction value
Figure 903117DEST_PATH_IMAGE033
Summing and taking the sum value
Figure 732113DEST_PATH_IMAGE034
To verify the first land parcel survey pattern to obtain
Figure 732430DEST_PATH_IMAGE026
The value being obtained by solving independent variables
Figure 934741DEST_PATH_IMAGE030
The value is
Figure 311496DEST_PATH_IMAGE002
By the verification mode, the accuracy of mu measurement of the convex polygon area is greatly improved. Moreover, only the mu measurement result of the convex polygon area needs to be obtained when the mu measurement result is predicted
Figure 666385DEST_PATH_IMAGE031
Value sum
Figure 962237DEST_PATH_IMAGE026
And the value is achieved, the area division of the convex polygon areas is not needed, the area calculation and summation are respectively carried out on each divided area, the calculation speed is greatly improved, and the method is particularly suitable for scenes in which the prediction of measuring mu is carried out on a plurality of convex polygon areas with different shapes at the same time.
It should be understood that the above description is only illustrative of the preferred embodiments of the present invention and the technical principles employed. It will be apparent to those skilled in the art that various modifications, equivalents, variations, and the like can be made to the present invention. However, such modifications are intended to fall within the scope of the present invention without departing from the spirit of the present invention. In addition, some terms used in the specification and claims of the present application are not limiting, but are merely for convenience of description.

Claims (10)

1. The real-time accurate land area mu measurement and correction method is characterized by comprising the following steps:
s1, drawing a corresponding virtual convex polygon of a convex polygon area determined on a land block to be measured
Figure FDA0004218807980000011
S2, predicting the virtual convex polygon
Figure FDA0004218807980000012
Area of->
Figure FDA0004218807980000013
S3, calculating
Figure FDA0004218807980000014
Calculating the difference DV between the two areas of one mu, and calculating the absolute value of the difference DV, and recording the absolute value as AVD;
s4, judging whether the difference DV is negative,
if yes, the area expansion flow is switched to;
if not, turning to an area contraction flow;
the area expansion flow comprises the following steps:
m1, calculating the virtual convex polygon
Figure FDA0004218807980000015
Each vertex p on k Deviation of->
Figure FDA0004218807980000016
The calculation method is expressed by the following formula (15): />
Figure FDA0004218807980000017
In the formula (15), L k-r Representing vertex p k To the virtual convex polygon
Figure FDA0004218807980000018
Other vertices p on r A linear distance therebetween;
n represents the virtual convex polygon
Figure FDA0004218807980000019
The number of vertices thereon;
m2, select the degree of deviation
Figure FDA00042188079800000110
The vertex with the smallest value is taken as an extension initial vertex p k1 And will be in contact with p k1 The vertex with relatively small deviation degree of two adjacent vertexes serves as an extension initial vertex p k2 Then obtain p k1 、p k2 The coordinates of the vertices are marked +.>
Figure FDA00042188079800000111
M3 according to the coordinates
Figure FDA00042188079800000112
Calculating the connection p k1 And p k2 The length of the first line of the apex, denoted L k ；
M4, according to absolute values AVD, L k Calculating the values from p k1 、p k2 The apex being connected at p perpendicularly to k1 And p k2 The height h of two extension lines extending in the direction of the first straight line between the vertexes 1 So that p k1 Vertex, slave p k1 Termination point p of vertex extension k1 ' from p k2 Termination point p of vertex extension k2 '、p k2 The area of a rectangle formed by enclosing between 4 points of the vertex is equal to the absolute value AVD solved in the step S3; the area contraction flow comprises the following steps:
n1, calculating the virtual convex polygon
Figure FDA0004218807980000021
Each vertex p on k Deviation of->
Figure FDA0004218807980000022
The calculation is expressed by the following formula (16): />
Figure FDA0004218807980000023
In the formula (16), L k-r Representing vertex p k To the virtual convex polygon
Figure FDA0004218807980000024
The upper vertex p r A linear distance therebetween;
n represents the virtual convex polygon
Figure FDA0004218807980000025
The number of vertices thereon;
n2, select degree of deviation
Figure FDA0004218807980000026
The vertex p with the greatest value k As a first initial vertex of area contraction, and selecting said virtual convex polygon ++>
Figure FDA0004218807980000027
Any one of a first adjacent vertex and a second adjacent vertex adjacent to the first initial vertex is taken as a second initial vertex with area shrinkage;
n3, obtaining the coordinates of the first initial vertex and the second initial vertex, calculating the length of a second straight line connected between the first initial vertex and the second initial vertex according to the obtained coordinates of the first initial vertex and the second initial vertex, and marking as L R The second straight line is used as the waist of an isosceles triangle to be contracted;
n4, taking the absolute value AVD as the area of the isosceles triangle to be contracted, and according to the length L of the second straight line R The coordinates of the first initial vertex and the second initial vertex calculate the bottom length w of the isosceles triangle to be contracted;
n5, shrinking the first initial vertex by w distance to the direction of another adjacent vertex which is not used as the second initial vertex and is adjacent to the first initial vertex, and obtaining a shrinking point p Shrinkage of ；
N6, the first initial vertex, the second initial vertex and the contraction point p Shrinkage of Removing the enclosed isosceles triangle; s5, expanding the areaThe virtual convex polygon after expansion or contraction
Figure FDA0004218807980000028
Mapping to the convex polygon area after area expansion or contraction under the physical space.
2. The method for real-time accurate acre measurement and correction of land parcels according to claim 1, wherein in step S2, the virtual convex polygon is calculated by the following method steps
Figure FDA0004218807980000029
Area of->
Figure FDA00042188079800000210
L1, for the virtual convex polygon
Figure FDA00042188079800000211
Performing similarity matching with each virtual convex polygon in the virtual convex polygon database,
if the matching is successful, the step L2 is carried out;
If the matching fails, terminating the mu measuring flow of the convex polygon area;
l2, the virtual convex polygon
Figure FDA0004218807980000031
Discretized into a plurality of virtual space polygons>
Figure FDA0004218807980000032
Then the virtual convex polygon is obtained +.>
Figure FDA0004218807980000033
Corresponding L 1 -n data pairs and obtaining each of said virtual space polygons +.>
Figure FDA0004218807980000034
Corresponding L 2 -m data pairs, L 1 N respectively represent the virtual convex polygon +.>
Figure FDA0004218807980000035
Central site O of (2) 1 The distance average value and the number of the vertexes; l (L) 2 M represents each of said virtual space polygons +.>
Figure FDA0004218807980000036
Central site O of (2) 2 The distance average value and the number of the vertexes;
l3, L 1 Inputting the data into a first land block mu measuring model corresponding to n, outputting a corresponding y value by the model, and recording the y value as
Figure FDA0004218807980000037
Will L 2 Inputting the measured area model into a second land block corresponding to m, outputting a corresponding y value by the model, and marking the y value as +.>
Figure FDA0004218807980000038
L4, the pair is opposite to the virtual convex polygon
Figure FDA00042188079800000323
Each of said virtual space polygons having a discrete relationship +.>
Figure FDA0004218807980000039
Corresponding area prediction value +.>
Figure FDA00042188079800000310
Summing to obtain a sum +.>
Figure FDA00042188079800000311
L5, judging the sum value
Figure FDA00042188079800000312
And->
Figure FDA00042188079800000313
Whether the absolute value of the difference is smaller than a preset first difference threshold,
if yes, will
Figure FDA00042188079800000314
The value is used as a mu measurement result of the convex polygon area;
if not, form
Figure FDA00042188079800000315
The data pairs are added into the first fitting point set corresponding to n and form +.>
Figure FDA00042188079800000316
Figure FDA00042188079800000317
The data pairs are added to the second fitting point set corresponding to m.
3. The method for real-time accurate acre measurement and correction of land parcels according to claim 2, wherein in step L2, the virtual convex polygon is obtained
Figure FDA00042188079800000318
Corresponding to the L 1 The method steps of the n data pairs include:
a1, respectively installing a distance sensor at each vertex of the convex polygon area, and recording the coordinates of each installation position under an XY axis coordinate system, which are respectively marked as (x) 1 ,y 1 ),(x 2 ,y 2 ),…,(x n ,y n ) N represents the number of vertices of the convex polygon area, and then the corresponding points of the convex polygon area are drawn on a computer according to the recorded coordinatesThe virtual convex polygon
Figure FDA00042188079800000319
A2, according to the virtual convex polygon
Figure FDA00042188079800000320
Calculating the coordinates of each vertex of the virtual convex polygon ++>
Figure FDA00042188079800000321
Is defined by the central site O of (2) 1 The coordinates in the XY-axis coordinate system are denoted as (x o1 ,y o1 ) Then calculate the central site O 1 And the virtual convex polygon where it is located ≡>
Figure FDA00042188079800000322
The average value of the distance between each vertex is denoted as L 1 Thereby obtaining the L 1 -n data pairs.
4. The method for real-time accurate acre measurement and correction of land parcels according to claim 3, wherein each virtual space polygon is obtained
Figure FDA0004218807980000041
Corresponding to the L 2 The method of m data pairs further comprises, on the basis of steps A1-A2:
a3, for the virtual convex polygon
Figure FDA0004218807980000042
Equally dividing each side of (2) and calculating each equally dividing point +.>
Figure FDA0004218807980000043
Coordinates in the XY-axis coordinate system, and +.>
Figure FDA0004218807980000044
And the central site O 1 After connecting the lines, the virtual convex polygon
Figure FDA0004218807980000045
Discretized into a plurality of virtual space polygons>
Figure FDA0004218807980000046
Figure FDA0004218807980000047
A4, calculating each virtual space polygon
Figure FDA0004218807980000048
Each vertex coordinate on the model is respectively marked as
Figure FDA0004218807980000049
m represents the virtual space polygon +.>
Figure FDA00042188079800000410
Is the number of vertices of (a);
a5, according to the virtual space polygon
Figure FDA00042188079800000411
Coordinates of each vertex on the virtual space polygon are calculated
Figure FDA00042188079800000412
Is defined by the central site O of (2) 2 The coordinates in the XY-axis coordinate system are denoted as (x o2 ,y o2 ) Then calculate the central site O 2 And the virtual space polygon where it is located +.>
Figure FDA00042188079800000413
The average value of the distance between each vertex is denoted as L 2 Thereby obtaining the association of each of said virtual space polygons +.>
Figure FDA00042188079800000414
Is not less than the L 2 -m data pairs.
5. The method for real-time accurate land parcel measurement and correction according to claim 4, wherein in step L3, the method steps of constructing the first land parcel measurement model and the second land parcel measurement model include:
C1, obtaining virtual convex polygons which are respectively corresponding to a plurality of convex polygon areas with different shapes and are determined on a land block with measured acre
Figure FDA00042188079800000415
Associated +.>
Figure FDA00042188079800000416
Data pair, wherein->
Figure FDA00042188079800000417
Representing the virtual convex polygon
Figure FDA00042188079800000418
Is a real area of (2);
c2, each virtual convex polygon
Figure FDA00042188079800000419
Discretized into a plurality of virtual space polygons>
Figure FDA00042188079800000420
Then obtaining +/for each said virtual convex polygon>
Figure FDA00042188079800000421
Each of said virtual space polygons having a discrete relationship +.>
Figure FDA00042188079800000422
Corresponding to
Figure FDA00042188079800000423
Data pair (s)/(s)>
Figure FDA00042188079800000424
Representing the virtual space polygon ++>
Figure FDA00042188079800000425
Is a real area of (2);
c3, acquiring the first land block acre measurement model and the second land block acre measurement model which correspond to n and m respectively;
c4, in several forms
Figure FDA00042188079800000426
L in data pair 1 As an argument +.>
Figure FDA00042188079800000427
Solving the first land block acre measurement model for the dependent variable to obtain a first parameter value of a first acre measurement parameter; by several->
Figure FDA0004218807980000051
Figure FDA0004218807980000052
L in data pair 2 As an argument +.>
Figure FDA0004218807980000053
Solving the second land block acre measurement model for the dependent variable to obtain a first parameter value of a second acre measurement parameter;
c5, substituting the first parameter values of the first acre measuring parameters and the second acre measuring parameters into the first land block acre measuring parameters respectivelyIn the model and the second land parcel measuring model, each virtual convex polygon is obtained
Figure FDA0004218807980000054
Corresponding L 1 The value is input into the first land parcel measuring model, the model outputs a corresponding y value which is recorded as +.>
Figure FDA0004218807980000055
And each of said virtual space polygons +.>
Figure FDA0004218807980000056
Corresponding L 2 The value is input into the second land parcel acre measuring model, the model outputs a corresponding y value which is recorded as
Figure FDA0004218807980000057
/>
C6, judging each
Figure FDA0004218807980000058
The value corresponds to the true value +.>
Figure FDA0004218807980000059
Whether the absolute value of the difference is smaller than a preset second difference threshold,
if yes, the input-output data pair L of the first land parcel measuring model in the step C5 is stored 1 -
Figure FDA00042188079800000510
Adding the first fitting point set as a fitting point;
if not, discarding the input/output data pair
Figure FDA00042188079800000511
Simultaneously judging each of the
Figure FDA00042188079800000512
The value corresponds to the true value +.>
Figure FDA00042188079800000513
Whether the absolute value of the difference is smaller than a preset third difference threshold,
if yes, the input and output data pair of the second land parcel measuring model in the step C5 is stored
Figure FDA00042188079800000514
Figure FDA00042188079800000515
Added as fitting points to the second set of fitting points,
if not, discarding the input/output data pair
Figure FDA00042188079800000516
C7, opposite to the virtual convex polygon
Figure FDA00042188079800000517
Each of said virtual space polygons having a discrete relationship +.>
Figure FDA00042188079800000518
Corresponding model predictive value +.>
Figure FDA00042188079800000519
Summing to obtain a sum +.>
Figure FDA00042188079800000520
C8, judging the sum value
Figure FDA00042188079800000521
Said virtual convex polygon having a discrete relationship therewith ++ >
Figure FDA00042188079800000522
Area true value +.>
Figure FDA00042188079800000523
Whether the absolute value of the difference is smaller than a preset fourth difference threshold,
if yes, go to step C9,
if not, filtering the virtual convex polygon from the first fitting point set
Figure FDA00042188079800000524
Said input/output data pair having an association relationship +.>
Figure FDA0004218807980000061
And filtering the virtual convex polygon from the second fitting point set
Figure FDA0004218807980000062
Said input/output data pair having an association relationship +.>
Figure FDA0004218807980000063
C9, fitting each fitting point in the first fitting point set and the second fitting point set respectively through an interpolation method of a Lagrangian interpolation polynomial to obtain a first fitting curve corresponding to the first fitting point set and a second fitting curve corresponding to the second fitting point set;
c10, solving second parameter values of the first acre measuring parameters of the first land parcel acre measuring model according to the first fitting curve, and solving second parameter values of the second acre measuring parameters of the second land parcel acre measuring model according to the second fitting curve;
and C11, updating and correcting second parameter values of the first acre measurement parameters and the second acre measurement parameters, and substituting the updated second parameter values into the corresponding first land parcel acre measurement model or the second land parcel acre measurement model to complete the construction of the first land parcel acre measurement model and the second land parcel acre measurement model.
6. The method for real-time accurate acre measurement and correction of land parcels according to claim 5, wherein the following steps are continuously performed after the step A5 to solve the virtual space polygon
Figure FDA0004218807980000064
Area of->
Figure FDA0004218807980000065
A6, for the virtual convex polygon
Figure FDA0004218807980000066
Two of said bisection points +.>
Figure FDA0004218807980000067
After direct connection, each virtual space polygon is +.>
Figure FDA0004218807980000068
Further discretizing into several virtual triangles, denoted +.>
Figure FDA0004218807980000069
/>
A7, according to each virtual triangle
Figure FDA00042188079800000610
Calculating its area +.>
Figure FDA00042188079800000611
And for each of said virtual space polygons +.>
Figure FDA00042188079800000612
Each of said virtual triangles having a discrete relationship +.>
Figure FDA00042188079800000613
Area of->
Figure FDA00042188079800000614
Summing, the sum value obtained is used as the corresponding virtual space polygon +>
Figure FDA00042188079800000615
Area of->
Figure FDA00042188079800000616
7. The method for real-time accurate acre measurement and correction of land parcels of claim 6, wherein the following steps are continued after step A7 to solve for the virtual convex polygon
Figure FDA00042188079800000617
Area of->
Figure FDA00042188079800000618
A8, for each virtual space polygon
Figure FDA00042188079800000619
Area of->
Figure FDA00042188079800000620
Summing the resulting sum as said virtual convex polygon +.>
Figure FDA00042188079800000621
Area of->
Figure FDA00042188079800000622
8. The method for real-time accurate measurement and correction of land parcels according to claim 5, wherein in step C3, when n=6, the obtained first land parcels measurement model corresponding to n is expressed by the following formula (1):
y＝a 1 x n +b 1 x n-1 +c 1 x n-2 +d 1 x n-3 +e 1 x n-4 +f 1 x n-5 +K 1 Formula (1)
When n=5, the obtained first land parcel measuring model corresponding to n is expressed by the following formula (2):
y＝a 2 x n +b 2 x n-1 +c 2 x n-2 +d 2 x n-3 +e 2 x n-4 +K 2 formula (2)
When n=4, the obtained first land parcel measuring model corresponding to n is expressed by the following formula (3):
y＝a 3 x n +b 3 x n-1 +c 3 x n-2 +d 3 x n-3 +K 3 formula (3)
In the formulas (1) - (3), a 1 、b 1 、c 1 、d 1 、e 1 、f 1 、K 1 、a 2 、b 2 、c 2 、d 2 、e 2 、K 2 、a 3 、b 3 、c 3 、d 3 、K 3 Solving each first acre measuring parameter for the parameter value to be calculated;
m=4, and the obtained second land parcel acre measurement model corresponding to m is expressed by the following formula (4):
y＝a 4 x m +b 4 x m-1 +c 4 x m-2 +d 4 x m-3 +K 4 formula (4)
In the formula (4), a 4 、b 4 、c 4 、d 4 、K 4 Solving for the parameter values to be madeAnd the second acre measurement parameter.
9. The method for real-time accurate land parcel measurement and correction according to claim 5, wherein in step C9, the method for fitting the first fitting curve or the second fitting curve by using the interpolation method of the lagrangian interpolation polynomial is expressed by the following formula (5):
Figure FDA0004218807980000071
in formula (5), q represents the q-th fitting point in the first fitting point set or the second fitting point set;
q represents the number of fitting points in the first fitting point set or the second fitting point set;
y q the y value which is predicted and output by the first land parcel measuring model or the second land parcel measuring model according to the input q fitting point is represented;
l q (x represents a lagrangian basis function expressed by the following expression (6):
Figure FDA0004218807980000072
x q 、x p the L value of the q-th fitting point and the p-th fitting point in the first fitting point set or the second fitting point set is L 1 Value or L 2 Values.
10. The method for real-time accurate measurement and correction of land parcels according to claim 5, wherein in step C11, the method for updating and correcting the second parameter value of each of the first measurement parameters comprises the steps of:
d1, calculating parameter solving error E 1 The calculation method is expressed by the following formula (7):
Figure FDA0004218807980000081
in the formula (7), par 1 C4, a first parameter value of the first acre measurement parameter obtained by solving in the step is represented;
Par 2 c10, obtaining a second parameter value of the same first mu measuring parameter by solving;
d2, predicted for step C5
Figure FDA0004218807980000082
Calculating a prediction error E 2 The calculation method is expressed by the following formula (8):
Figure FDA0004218807980000083
d3, judge E 1 Whether or not to follow E 2 Is increased by the increase of (a),
if yes, correcting a second parameter value of the first mu measuring parameter by the following formula (9):
Par 2 '＝(1-E 1 )Par 2 formula (9)
If not, correcting the second parameter value of the first acre measurement parameter by the following formula (10):
Par 2 '＝(1+E 1 )Par 2 formula (10)
In formulas (9) - (10), par 2 ' a second parameter value representing the corrected first acre measurement parameter;
In step C11, the method for updating and correcting the second parameter value of each second acre measurement parameter includes the steps of:
f1, calculating parameter solving error E 3 The calculation method is expressed by the following formula (11):
Figure FDA0004218807980000084
in the formula (11), par 3 C4, representing the first parameter value of the second mu measuring parameter obtained by solving in the step;
Par 4 c10, obtaining a second parameter value of the same second mu measuring parameter by solving;
f2, for each of the predictions calculated in step C5
Figure FDA0004218807980000085
Calculating a prediction error mean mv, wherein the calculation method is expressed by the following formula (12):
Figure FDA0004218807980000091
in the formula (12), E i Representing the virtual convex polygon
Figure FDA0004218807980000092
I th said virtual space polygon +.>
Figure FDA0004218807980000093
Is a prediction error of the area of (2);
f3, judge E 3 Whether or not to increase with increasing mv,
if yes, correcting a second parameter value of the second acre measurement parameter by the following formula (13):
Par 4 '＝(1-E 3 )Par 4 formula (13)
If not, correcting a second parameter value of the second acre measurement parameter by the following formula (14):
Par 4 '＝(1+E 3 )Par 4 formula (14)
Par in formulas (13) - (14) 4 ' represents a second parameter value of the corrected second acre measurement parameter.
CN202211229682.4A 2022-10-09 2022-10-09 Real-time accurate acre measurement and correction method for land parcels Active CN115511947B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211229682.4A CN115511947B (en) 2022-10-09 2022-10-09 Real-time accurate acre measurement and correction method for land parcels

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211229682.4A CN115511947B (en) 2022-10-09 2022-10-09 Real-time accurate acre measurement and correction method for land parcels

Publications (2)

Publication Number Publication Date
CN115511947A CN115511947A (en) 2022-12-23
CN115511947B true CN115511947B (en) 2023-06-09

Family

ID=84508809

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211229682.4A Active CN115511947B (en) 2022-10-09 2022-10-09 Real-time accurate acre measurement and correction method for land parcels

Country Status (1)

Country Link
CN (1) CN115511947B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109883320A (en) * 2019-01-22 2019-06-14 湖南土流信息有限公司 A kind of Land area measure method and system

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002298146A (en) * 2001-04-02 2002-10-11 Fujitsu General Ltd Polygonal area calculating method
JP4501479B2 (en) * 2004-03-18 2010-07-14 カシオ計算機株式会社 Polygon model simplification method, image processing apparatus, image drawing apparatus, and program
US8508533B2 (en) * 2011-09-28 2013-08-13 Palantir Technologies, Inc. Simplifying a polygon
CN106127820A (en) * 2016-06-30 2016-11-16 北京合众思壮科技股份有限公司 The area of a kind of irregular terrain profiles determines method and device
CN106568412A (en) * 2016-11-09 2017-04-19 无锡汉和航空技术有限公司 Area measurement method for irregularly-shaped plot
CN108332652B (en) * 2018-01-15 2019-11-05 中国农业大学 A kind of farmland working area metering method and device

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109883320A (en) * 2019-01-22 2019-06-14 湖南土流信息有限公司 A kind of Land area measure method and system

Also Published As

Publication number Publication date
CN115511947A (en) 2022-12-23

Similar Documents

Publication Publication Date Title
CN106204557A (en) A kind of extracting method of the non-complete data symmetrical feature estimated with M based on extension Gaussian sphere
CN106709943B (en) A kind of point cloud registration method based on optimal transmission
CN110838145B (en) Visual positioning and mapping method for indoor dynamic scene
CN111369607A (en) Prefabricated part assembling and matching method based on picture analysis
CN113777931B (en) Icing wing type pneumatic model construction method, device, equipment and medium
Rouhani et al. Implicit polynomial representation through a fast fitting error estimation
CN113468706A (en) Laser point cloud power transmission line lead fitting method for distribution network live working robot
CN112033338A (en) Blade curved surface contact type scanning measurement probe radius surface compensation method
CN113570722A (en) Surrounding rock crack information extraction and integrity coefficient rapid determination method
CN115511947B (en) Real-time accurate acre measurement and correction method for land parcels
CN114323536A (en) Interpolation method for improving measurement precision of five-hole probe
CN109446614A (en) A kind of discrimination method and device of static reactive model parameter
CN114937184A (en) Training method and device for cardiac coronary vessel naming model and electronic equipment
CN112153564B (en) Efficient multi-hop positioning method based on combination of centralized and distributed computing
CN115235439B (en) Land parcel mu measuring method based on grid division and data correction
CN110046412B (en) Circular grating angle measurement error correction method based on optimized BP neural network
CN106695796A (en) Portable laser scanning measuring arm hand-eye calibration method based on RANSAC
CN112991300B (en) Single wood skeleton extraction and visualization method based on neighborhood characteristics
Al-Anssari et al. Three-dimensional laplacian spatial filter of a field of vectors for geometrical edges magnitude and direction detection in point cloud surfaces
CN113393507A (en) Unmanned aerial vehicle point cloud and ground three-dimensional laser scanner point cloud registration method
CN111765889A (en) Pose positioning method of mobile robot in production workshop based on multi-cell-ellipsoid dual-filtering
CN113029049A (en) Complex workpiece optical measurement method based on weighted positive and negative margin variance minimization algorithm
CN111210877A (en) Method and device for deducing physical property parameters
CN116309641B (en) Image area acquisition system
CN106570864A (en) Quadratic curve fitting method based on geometric error optimization in image

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant