CN106989664A - Soil subsoiling qualification rate acquisition methods with GPS track interpolation are scanned based on cell - Google Patents
Soil subsoiling qualification rate acquisition methods with GPS track interpolation are scanned based on cell Download PDFInfo
- Publication number
- CN106989664A CN106989664A CN201710379117.9A CN201710379117A CN106989664A CN 106989664 A CN106989664 A CN 106989664A CN 201710379117 A CN201710379117 A CN 201710379117A CN 106989664 A CN106989664 A CN 106989664A
- Authority
- CN
- China
- Prior art keywords
- point
- coordinate
- coordinates
- qualification rate
- longitude
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
- G01B7/26—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring depth
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Soil subsoiling qualification rate acquisition methods with GPS track interpolation are scanned based on cell, the present invention relates to soil subsoiling qualification rate acquisition methods.It is low the invention aims to solve existing method calculating soil subsoiling qualification rate precision, and there is the problem of practising fraud.Process is：Step one：Data acquisition：Agricultural machinery latitude and longitude coordinates and Subsoiler depth are gathered every the set time by the GPS positioning device on agricultural equipment and subsoiling sensor；Step 2：Coordinate Conversion：The agricultural machinery latitude and longitude coordinates of collection are converted to the XY coordinates be converted under plane right-angle coordinate according to gauss projection；Step 3：Piecewise interpolation；Step 4：Build circumscribed rectangle；Step 5：Scanning, obtains soil subsoiling qualification rate M ÷ N × 100%.The present invention is used for soil subsoiling qualification rate field.
Description
Technical field
The present invention relates to soil subsoiling qualification rate acquisition methods.
Background technology
Soil subsoiling refers to by tractor-drawn subsoiling equipment, chesson, breaks plough sole, improves topsoil structure,
Strengthen a farming technique of soil water storage soil moisture conservation and drought resisting water drainage ability.Carrying out subsoiling earth working is conducive to crops to give birth to
It is long, be improve crop yield important means it.The main calculation methods of current soil subsoiling qualification rate are constant durations
Subsoiling depth is recorded, if the total data points of record are W, wherein M data meet subsoiling requirement, and soil subsoiling is qualified
Rate is M ÷ W × 100%.Because agricultural machinery working speed is not at the uniform velocity, it is impossible to ensure that test point is evenly distributed on operation soil
On, it is low to calculate soil subsoiling qualification rate precision in this way, and there is " cheating " method.When practising fraud as producer's Subsoiler
Utilize existing method leak, the artificial means for improving soil subsoiling qualification rate.
The content of the invention
It is low the invention aims to solve existing method calculating soil subsoiling qualification rate precision, and there is asking for cheating
Topic, and propose to scan the soil subsoiling qualification rate acquisition methods with GPS track interpolation based on cell.
Scanned based on cell and the soil subsoiling qualification rate acquisition methods detailed process of GPS track interpolation is：
Step one：Data acquisition：By the GPS positioning device on agricultural equipment and subsoiling sensor every fixation
Time gathers agricultural machinery latitude and longitude coordinates and Subsoiler depth；
Set time is artificial setting；
Step 2：Coordinate Conversion：The agricultural machinery latitude and longitude coordinates of collection are converted into flat square according to gauss projection conversion
XY coordinates under coordinate system；
Step 3：Piecewise interpolation：Segment processing is carried out to each adjacent coordinates in coordinate points, k sections of coordinate point set C are obtainedi,
I=1 ..., k, k value are positive integer；Again by obtained k section coordinate point sets CiDistance is according to decile between middle adjacent coordinates point
Mode inserts new coordinate points, makes the distance between neighbor interpolation point take agricultural equipment cultivation width long less than or equal to set-point L, L
Degree；
Step 4：Build circumscribed rectangle：Judge the azimuthal point of each section of internal coordinate point, minimax is found according to azimuthal point
Four coordinate X of abscissa and minimax ordinatemin, Xmax, Ymin, Ymax, it is X to expand four coordinatesmin-L,Xmax+L,
Ymin-L,YmaxThis four coordinates are formed four boundary points by+L by combining, and make the flat of XY axles on the basis of this four boundary points
Line, builds circumscribed rectangle Ri；
Step 5：Scanning：By the length of side of the cultivation width L of agricultural equipment in every section of coordinate point set CiCircumscribed rectangle RiIn
Square shaped cells lattice are built, circumscribed rectangle R is paved with square shaped cells latticei, i=1 ..., k, successively to each square shaped cells
Lattice are traveled through, and are judged to whether there is the coordinate points after interpolation in this square shaped cells lattice, are added up if carrying out counting in the presence of if, most
After obtain every section of coordinate point set CiTotal number A containing square shaped cells latticei, i=1 ..., k, to each section of AiAdded up
Obtain overall operation square shaped cells lattice number in the periodN values are positive integer；One is taken in each square shaped cells lattice
The maximum data point of individual Subsoiler depth, judges whether the maximum data point of depth in square shaped cells lattice is qualified, if conjunction
The square shaped cells lattice number of lattice is M, then obtains soil subsoiling qualification rate, and soil subsoiling qualification rate is M ÷ N × 100%.
Beneficial effects of the present invention are：
It is of the invention effectively using the GIS-Geographic Information System in modern information technologies, the method scanned by cell is by subsoiling
Operation soil is divided into several junior unit lattice, and ensures there is at least one subsoiling data test point in each cell, i.e., real
Show subsoiling data test point being evenly distributed in operation soil, improve the precision for calculating soil subsoiling qualification rate, and avoid
The problem of cheating, thus the present invention be it is a kind of efficiently, the high solution of convenient, accuracy.
The soil subsoiling qualification rate that existing method is obtained is used in embodiment for 80%, and actual soil subsoiling qualification rate is only
Have 17%, and the present invention is scanned based on cell and the soil subsoiling qualification rate computational methods of GPS track interpolation calculate soil depth
Loose qualification rate is 16.6%, it is to avoid the error in existing method calculating process, obtained soil subsoiling qualification rate and actual soil
Ground subsoiling qualification rate is close, improves calculating soil subsoiling qualification rate precision.
Brief description of the drawings
Fig. 1 is the present invention based on cell scanning and GPS track interpolation soil subsoiling qualification rate acquisition methods flow chart；
Fig. 2 is to determine that cell scans schematic diagram；
Fig. 3 is that traditional work depth changes with time
Fig. 4 is that traditional work area changes with time.
Embodiment
Embodiment one：Illustrate present embodiment with reference to Fig. 1, present embodiment based on cell scanning and GPS
The soil subsoiling qualification rate acquisition methods detailed process of track interpolation is：
Step one：Data acquisition：By the GPS positioning device on agricultural equipment and subsoiling sensor every fixation
Time gathers agricultural machinery latitude and longitude coordinates and Subsoiler depth；
Set time is artificial setting；
Step 2：Coordinate Conversion：The agricultural machinery latitude and longitude coordinates of collection are converted into flat square according to gauss projection conversion
XY coordinates under coordinate system；
Step 3：Piecewise interpolation：Segment processing is carried out to each adjacent coordinates in coordinate points, k sections of coordinate point set C are obtainedi,
I=1 ..., k, k value are positive integer；Again by obtained k section coordinate point sets CiDistance is according to decile between middle adjacent coordinates point
Mode inserts new coordinate points, makes the distance between neighbor interpolation point take agricultural equipment cultivation width long less than or equal to set-point L, L
Degree；
So the dense degree of interpolation point can be determined according to the length of side of permissible accuracy and cell selection is calculated；Reason
Interpolation is more intensive within the specific limits by upper, and the precision of calculating is higher；
Step 4：Build circumscribed rectangle：Judge the azimuthal point of each section of internal coordinate point, minimax is found according to azimuthal point
Four coordinate X of abscissa and minimax ordinatemin, Xmax, Ymin, Ymax, in order that the coordinate points on border are also contained in list
It is X that four coordinates are expanded in first latticemin-L,Xmax+L,Ymin-L,YmaxThis four coordinates are formed four borders by+L by combining
Point, makees the parallel lines of XY axles on the basis of this four boundary points, builds circumscribed rectangle Ri；
Step 5：Scanning：By the length of side of the cultivation width L of agricultural equipment in every section of coordinate point set CiCircumscribed rectangle RiIn
Square shaped cells lattice are built, circumscribed rectangle R is paved with square shaped cells latticei, i=1 ..., k, successively to each square shaped cells
Lattice are traveled through, and are judged to whether there is the coordinate points after interpolation in this square shaped cells lattice, are added up if carrying out counting in the presence of if, most
After obtain every section of coordinate point set CiTotal number A containing square shaped cells latticei, i=1 ..., k, to each section of AiAdded up
Obtain overall operation square shaped cells lattice number in the periodN values are positive integer；One is taken in each square shaped cells lattice
The maximum data point of individual Subsoiler depth, judges whether the maximum data point of depth in square shaped cells lattice is qualified, if conjunction
The square shaped cells lattice number of lattice is M, then obtains soil subsoiling qualification rate, and soil subsoiling qualification rate is M ÷ N × 100%.Such as Fig. 2.
Embodiment two：Present embodiment from unlike embodiment one：Coordinate turns in the step 2
Change：The latitude and longitude coordinates of coordinate acquisition are converted to the XY coordinates be converted under plane right-angle coordinate according to gauss projection；Specifically
Process is：
The latitude and longitude coordinates of agricultural machinery are uploaded every the set time by the GPS positioning device on agricultural equipment, and
The latitude and longitude coordinates of upload are converted to the XY coordinates be converted under plane right-angle coordinate according to gauss projection；
Gauss projection full name converts for Gauss Kru＆4＆ger projection, is a kind of transverse elliptic cylindrical equalangle projection conversion, by becoming
The calculation process that coordinate after changing can just enter on row distance, and then calculate area.
Process is：
If the longitude and latitude that biography is engraved during GPS device is respectively longitude and latitude, unit is degree；
It is iPI=π/180.0 to calculate the radian corresponding to 1 °；
In formula, iPI is converted into the coefficient of radian for degree；
Using point method with projection, make the deformation of projecting edge will not be excessive：ZoneWide=6 is 6 ° of bands；
In formula, ZoneWide is through difference；
Coordinate system parameters were from 54 years Beijing coordinate system parameters：
Spheroid major semiaxis length：A=6378137.0
Spheroid semi-minor axis length：B=6356752.3142
Longitude and latitude is respectively converted into radian longitude1 and latitude1：
Longitude1=longitude*iPI
Latitude1=latitude*iPI
Search reel number：ProjNo=longitude/ZoneWide
Longitude0=ProjNo*ZoneWide+ZoneWide/2
In formula, ProjNo is longitude with being numbered through difference than obtained projection zone；Longitude0 is the longitude of central meridian；
Obtain central meridian：Longitude2=longitude0*iPI
In formula, longitude2 is the corresponding radian of longitude of central meridian；
The eccentricity of ellipsoid first：
The eccentricity of ellipsoid second：
Calculate radius of curvature in prime vertical：
Calculate Meridian arc length：
S=a* ((1-e2/4-3*e2 2/64-5*e2 3/256)*latitude1-(3*e2/8+3*e2 2/32+45*e2 3/
1024)*sin(2*latitude1)+(15*e2 2/256+45*e2 3/1024)*sin(4*latitude1)-(35*e2 3/3072)*
sin(6*latitude1))
In formula, * is multiplication sign；
With interior geodetic coordinates：
Xval=R* (A+ (1-T+C) * A3/6+(5-18*T+T2+72*C-58*e)*A5/120)
Yval=S+R*tan (latitude1) * (A2/2+(5-T+9*C+4*C2)*A4/24+(61-58*T+T2+600*C-
330*e)*A6/720)
In formula, xval is with interior the earth x-axis coordinate；Yval is with interior the earth y-axis coordinate；
T=tan2(latitude1)
C=e1*cos2(latitude1)
A=(longitude1-longitude0) * cos (latitude1)
T, C, A are intermediate variable；
Calculate the offset on abscissa direction：X0=1000000* (ProjNo+1)+500000；
Calculating obtains rectangular co-ordinate：X=xval+X0, Y=yval (unit is rice).
Other steps and parameter are identical with embodiment one.
Embodiment three：Present embodiment from unlike embodiment one or two：Divide in the step 3
Section interpolation：Segment processing is carried out to each adjacent coordinates in coordinate points, k sections of coordinate point set C are obtainedi, i=1 ..., k, k value
For positive integer；Again by obtained k section coordinate point sets CiDistance inserts new coordinate in the way of decile between middle adjacent coordinates point
Point, makes the distance between neighbor interpolation point take agricultural equipment cultivation width length less than or equal to set-point L, L；Detailed process is：
Piecewise interpolation：The distance in continuous coordinate points between each adjacent coordinates point, which is respectively less than, in temporal sequence is equal to M meters,
It is a continuous active section then to think continuous coordinate points, is divided into one section；If the distance between adjacent coordinates point is more than M meters,
Not think it is an active section, adjacent coordinates point is divided into two sections；K sections of coordinate point set C are obtained after segmented processingi, i=
1 ..., k, k value be positive integer；Again by obtained k section coordinate point sets CiDistance is in the way of decile between middle adjacent coordinates point
The new coordinate points of insertion, make the distance between neighbor interpolation point take agricultural equipment cultivation width length less than or equal to set-point L, L.
Other steps and parameter are identical with embodiment one or two.
Embodiment four：Unlike one of present embodiment and embodiment one to three：The step 5
Whether the maximum data point of the middle depth judged in square shaped cells lattice is qualified, and detailed process is：
The depth of the maximum data point of depth and Ministry of Agriculture regulation Subsoiler depth are compared, when depth maximum
When the depth of data point is more than or equal to Ministry of Agriculture's regulation Subsoiler depth, the data point of the depth maximum in square shaped cells lattice
It is otherwise as unqualified to be qualified.
For example：Subsoiler depth standards are 30cm as defined in Shanxi Agricultural portion；Subsoiling as defined in the Ministry of Agriculture of the Three Eastern Provinces is made
Industry depth is more than 25 centimetres, no more than 40 centimetres.
Other steps and parameter are identical with one of embodiment one to three.
Beneficial effects of the present invention are verified using following examples：
Embodiment one：
The present embodiment based on cell scan and GPS track interpolation soil subsoiling qualification rate acquisition methods be specifically according to
Prepared by following steps：
Fig. 3 and Fig. 4 list a kind of extreme case, it is intended to illustrate asking for current soil subsoiling qualification rate computational methods presence
Topic.Worked in the first two ten minutes of operation agricultural machinery with fast speed, subsoiling mean depth is 10cm, operation completes 83% area
Soil；Worked afterwards with slower speed within 80 minutes, complete the soil of remaining 17% area, Subsoiler mean depth is 35cm.Root
It is 25cm according to Shanxi Agricultural portion regulation Subsoiler depth, because subsoiling data are that constant duration is taken, therefore first 20 minutes
Subsoiling data point is number of non-compliances strong point, and the subsoiling data point of latter 80 minutes is qualified data point, is closed according to current soil subsoiling
Lattice rate acquisition methods, qualification rate isActual soil subsoiling qualification rate only has 17%, i.e.,
When obtaining this using soil subsoiling qualification rate acquisition methods of the present invention based on cell scanning and GPS track interpolation
Overall operation square shaped cells lattice number is 60 in section；The maximum data of a Subsoiler depth are taken in each square shaped cells lattice
Point, provides that Subsoiler depth is 25cm, qualified square shaped cells lattice number is 10, then soil subsoiling is closed according to Shanxi Agricultural portion
Lattice rate is
Present invention contrast prior art, improves subsoiling qualification rate computational accuracy.
Traditional soil subsoiling qualification rate computational methods are not fixed due to agricultural machinery operating rate, and sampling point distributions are uneven,
The problem of causing qualification rate computational accuracy relatively low, and there is " cheating " method.For sampling point distributions it is uneven the problem of, this hair
The bright thought using calculus carries out interpolation processing to the coordinate of upload, and operation soil can be so divided into several junior units
Lattice, and according to the ratio calculation subsoiling qualification rate of qualified cell number and total cell number, it is ensured that sampled point is uniformly distributed,
With higher computational accuracy.
Claims (4)
1. the soil subsoiling qualification rate acquisition methods with GPS track interpolation are scanned based on cell, it is characterised in that：Methods described
Detailed process is：
Step one：Data acquisition：By the GPS positioning device on agricultural equipment and subsoiling sensor every the set time
Gather agricultural machinery latitude and longitude coordinates and Subsoiler depth；
Step 2：Coordinate Conversion：The agricultural machinery latitude and longitude coordinates of collection are converted into plane rectangular coordinates according to gauss projection conversion
XY coordinates under system；
Step 3：Piecewise interpolation：Segment processing is carried out to each adjacent coordinates in coordinate points, k sections of coordinate point set C are obtainedi, i=
1 ..., k, k value be positive integer；Again by obtained k section coordinate point sets CiDistance is in the way of decile between middle adjacent coordinates point
The new coordinate points of insertion, make the distance between neighbor interpolation point take agricultural equipment cultivation width length less than or equal to set-point L, L；
Step 4：Build circumscribed rectangle：Judge the azimuthal point of each section of internal coordinate point, the horizontal seat of minimax is found according to azimuthal point
Four coordinate X of mark and minimax ordinatemin, Xmax, Ymin, Ymax, it is X to expand four coordinatesmin-L,Xmax+L,Ymin-L,
YmaxThis four coordinates are formed four boundary points by combining, are made the parallel lines of XY axles on the basis of this four boundary points by+L,
Build circumscribed rectangle Ri；
Step 5：Scanning：By the length of side of the cultivation width L of agricultural equipment in every section of coordinate point set CiCircumscribed rectangle RiIt is middle to build
Square shaped cells lattice, circumscribed rectangle R is paved with square shaped cells latticei, i=1 ..., k enter to each square shaped cells lattice successively
Row traversal, judges to whether there is the coordinate points after interpolation in this square shaped cells lattice, cumulative if carrying out counting in the presence of if, finally
To every section of coordinate point set CiTotal number A containing square shaped cells latticei, i=1 ..., k, to each section of AiAdded up
Overall operation square shaped cells lattice number in periodN values are positive integer；A depth is taken in each square shaped cells lattice
The maximum data point of loose depth of implements, judges whether the data point of depth maximum in square shaped cells lattice is qualified, if qualified
Square shaped cells lattice number is M, then obtains soil subsoiling qualification rate, and soil subsoiling qualification rate is M ÷ N × 100%.
2. the soil subsoiling qualification rate acquisition methods with GPS track interpolation are scanned based on cell according to claim 1, its
It is characterised by：Coordinate Conversion in the step 2：The latitude and longitude coordinates of coordinate acquisition are converted to according to gauss projection conversion flat
XY coordinates under the rectangular coordinate system of face；Detailed process is：
The latitude and longitude coordinates of agricultural machinery are uploaded every the set time by the GPS positioning device on agricultural equipment, and will be upper
The latitude and longitude coordinates of biography convert the XY coordinates be converted under plane right-angle coordinate according to gauss projection；Process is：
If the longitude and latitude that biography is engraved during GPS device is respectively longitude and latitude, unit is degree；
It is iPI=π/180.0 to calculate the radian corresponding to 1 °；
In formula, iPI is converted into the coefficient of radian for degree；
Using point method with projection：° band of ZoneWide=6, i.e., 6；
In formula, ZoneWide is through difference；
Coordinate system parameters were from 54 years Beijing coordinate system parameters：
Spheroid major semiaxis length：A=6378137.0
Spheroid semi-minor axis length：B=6356752.3142
Longitude and latitude is respectively converted into radian longitude1 and latitude1：
Longitude1=longitude*iPI
Latitude1=latitude*iPI
Search reel number：ProjNo=longitude/ZoneWide
Longitude0=ProjNo*ZoneWide+ZoneWide/2
In formula, ProjNo is longitude with being numbered through difference than obtained projection zone；Longitude0 is the longitude of central meridian；
Obtain central meridian：Longitude2=longitude0*iPI
In formula, longitude2 is the corresponding radian of longitude of central meridian；
The eccentricity of ellipsoid first：
The eccentricity of ellipsoid second：
Calculate radius of curvature in prime vertical：
Calculate Meridian arc length：
S=a* ((1-e2/4-3*e2 2/64-5*e2 3/256)*latitude1-(3*e2/8+3*e2 2/32+45*e2 3/1024)*sin
(2*latitude1)+(15*e2 2/256+45*e2 3/1024)*sin(4*latitude1)-(35*e2 3/3072)*sin(6*
latitude1))
In formula, * is multiplication sign；
With interior geodetic coordinates：
Xval=R* (A+ (1-T+C) * A3/6+(5-18*T+T2+72*C-58*e)*A5/120)
Yval=S+R*tan (latitude1) * (A2/2+(5-T+9*C+4*C2)*A4/24+(61-58*T+T2+600*C-330*
e)*A6/720)
In formula, xval is with interior the earth x-axis coordinate；Yval is with interior the earth y-axis coordinate；
T=tan2(latitude1)
C=e1*cos2(latitude1)
A=(longitude1-longitude0) * cos (latitude1)
In formula, T, C, A are intermediate variable；
Calculate the offset on abscissa direction：X0=1000000* (ProjNo+1)+500000；
Calculating obtains rectangular co-ordinate：X=xval+X0, Y=yval.
3. the soil subsoiling qualification rate acquisition methods with GPS track interpolation are scanned based on cell according to claim 2, its
It is characterised by：Piecewise interpolation in the step 3：Segment processing is carried out to each adjacent coordinates in coordinate points, k sections of coordinates are obtained
Point set Ci, i=1 ..., k, k value are positive integer；Again by obtained k section coordinate point sets CiBetween middle adjacent coordinates point distance according to
The mode of decile inserts new coordinate points, the distance between neighbor interpolation point is taken agricultural equipment farming less than or equal to set-point L, L
Breadth length；Detailed process is：
Piecewise interpolation：The distance in continuous coordinate points between each adjacent coordinates point, which is respectively less than, in temporal sequence is equal to M meters, then recognizes
It is a continuous active section for continuous coordinate points, is divided into one section；If the distance between adjacent coordinates point is more than M meters, then it is assumed that
It is not an active section, adjacent coordinates point is divided into two sections；K sections of coordinate point set C are obtained after segmented processingi, i=1 ..., k,
K values are positive integer；Again by obtained k section coordinate point sets CiDistance inserts new in the way of decile between middle adjacent coordinates point
Coordinate points, make the distance between neighbor interpolation point take agricultural equipment cultivation width length less than or equal to set-point L, L.
4. the soil subsoiling qualification rate acquisition methods with GPS track interpolation are scanned based on cell according to claim 3, its
It is characterised by：Judge whether the maximum data point of depth in square shaped cells lattice is qualified, and detailed process is in the step 5：
The depth of the maximum data point of depth and Ministry of Agriculture regulation Subsoiler depth are compared, when the data of depth maximum
When the depth of point is more than or equal to Ministry of Agriculture's regulation Subsoiler depth, the data point of the depth maximum in square shaped cells lattice is to close
Lattice, it is otherwise as unqualified.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710379117.9A CN106989664A (en) | 2017-05-25 | 2017-05-25 | Soil subsoiling qualification rate acquisition methods with GPS track interpolation are scanned based on cell |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710379117.9A CN106989664A (en) | 2017-05-25 | 2017-05-25 | Soil subsoiling qualification rate acquisition methods with GPS track interpolation are scanned based on cell |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106989664A true CN106989664A (en) | 2017-07-28 |
Family
ID=59420714
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710379117.9A Pending CN106989664A (en) | 2017-05-25 | 2017-05-25 | Soil subsoiling qualification rate acquisition methods with GPS track interpolation are scanned based on cell |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106989664A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109388773A (en) * | 2018-09-14 | 2019-02-26 | 北京金双狐油气技术有限公司 | A kind of figure coordinate operation method and device |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030112235A1 (en) * | 2001-12-19 | 2003-06-19 | Earth Science Associates, Inc. | Method and system for creating irregular three-dimensional polygonal volume models in a three-dimensional geographic information system |
| CN104977586A (en) * | 2014-04-10 | 2015-10-14 | 中国农业机械化科学研究院 | Suspension-type subsoiler and online tilling depth detection device and method |
| CN105652840A (en) * | 2016-02-03 | 2016-06-08 | 西北农林科技大学 | Mechanical deep scarification work quality monitoring and evaluation method and system |
| CN205827143U (en) * | 2016-06-20 | 2016-12-21 | 北京农业智能装备技术研究中心 | Subsoiling land preparation operation quality detecting system |
| CN106247926A (en) * | 2016-07-13 | 2016-12-21 | 哈尔滨工业大学 | Based on cell scanning and the farm machinery tillage Area computing method of GPS track interpolation |
-
2017
- 2017-05-25 CN CN201710379117.9A patent/CN106989664A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030112235A1 (en) * | 2001-12-19 | 2003-06-19 | Earth Science Associates, Inc. | Method and system for creating irregular three-dimensional polygonal volume models in a three-dimensional geographic information system |
| CN104977586A (en) * | 2014-04-10 | 2015-10-14 | 中国农业机械化科学研究院 | Suspension-type subsoiler and online tilling depth detection device and method |
| CN105652840A (en) * | 2016-02-03 | 2016-06-08 | 西北农林科技大学 | Mechanical deep scarification work quality monitoring and evaluation method and system |
| CN205827143U (en) * | 2016-06-20 | 2016-12-21 | 北京农业智能装备技术研究中心 | Subsoiling land preparation operation quality detecting system |
| CN106247926A (en) * | 2016-07-13 | 2016-12-21 | 哈尔滨工业大学 | Based on cell scanning and the farm machinery tillage Area computing method of GPS track interpolation |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109388773A (en) * | 2018-09-14 | 2019-02-26 | 北京金双狐油气技术有限公司 | A kind of figure coordinate operation method and device |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106247926B (en) | Farm machinery tillage Area computing method based on cell scanning and GPS track interpolation | |
| CN107063197B (en) | Reservoir characteristic curve extraction method based on spatial information technology | |
| CN105718751B (en) | Deep scarification operation area calculating method based on area grid point coverage | |
| CN105652840B (en) | Machinery subsoiling operation quality monitoring method and system | |
| CN107462208A (en) | A kind of agricultural machinery and agricultural machinery working area measuring device and measuring method | |
| CN106908415A (en) | A kind of big region crops time of infertility Soil Moisture Monitoring method based on amendment NDVI time serieses | |
| CN104424390A (en) | Irrigation area monitoring method and device | |
| Ferguson et al. | Soil sampling for site‐specific nitrogen management | |
| CN106990216A (en) | A kind of shallow lake wawter bloom risk analysis early warning system and its analysis and early warning method | |
| CN106485002B (en) | In the method for complicated landform climatic province estimation sugarcane potential production | |
| CN104376204A (en) | Method for inverting vegetation coverage by adopting improved pixel dichotomy | |
| CN103440648B (en) | The method for automatic measurement of corps canopy regularity and device | |
| CN104794336A (en) | Farmland straw resource spacial distribution estimation method | |
| CN103413172B (en) | Maize at Seedling Stage is short of seedling several measuring methods | |
| CN107941673A (en) | A kind of salt affected soil pore structure measures analysis method | |
| CN105181607A (en) | Soil quick-acting nutrient inversion method and device | |
| CN105371762A (en) | Image analysis-based fruit tree crown volume measurement method | |
| CN109410736B (en) | A kind of multi-source dem data seamless integration method and processing terminal | |
| JP2018066616A (en) | Soil diagnostic method and soil condition improvement method | |
| CN106960190B (en) | The irrigated area statistical method of medium resolution image is supervised based on high-definition picture | |
| CN114239756A (en) | Insect pest detection method and system | |
| CN106989664A (en) | Soil subsoiling qualification rate acquisition methods with GPS track interpolation are scanned based on cell | |
| CN107831516A (en) | Merge the method that GNSS and ground monitoring net obtain dam real-time high-precision displacement | |
| Farazmand et al. | Determining the factors affecting rangeland suitability for livestock and wildlife grazing. | |
| Catania et al. | Vineyard design supported by GPS application |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170728 |
|
| RJ01 | Rejection of invention patent application after publication |