CN106323670A - Method and device for determining soil sampling points - Google Patents
Method and device for determining soil sampling points Download PDFInfo
- Publication number
- CN106323670A CN106323670A CN201610829349.5A CN201610829349A CN106323670A CN 106323670 A CN106323670 A CN 106323670A CN 201610829349 A CN201610829349 A CN 201610829349A CN 106323670 A CN106323670 A CN 106323670A
- Authority
- CN
- China
- Prior art keywords
- grid
- soil sampling
- sampled point
- district
- sampling district
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/02—Devices for withdrawing samples
- G01N1/04—Devices for withdrawing samples in the solid state, e.g. by cutting
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
The invention discloses a method and device for determining soil sampling points. The method includes: determining a soil sampling area; determining the interval among initial soil sampling points of the soil sampling area according to the area of the soil sampling area and appointed sampling point quantity; converting the soil sampling area into grid data according to the interval among the initial soil sampling points, and determining effective grid quantity of grids of the current soil sampling area; regulating the size of the grids of the current soil sampling area according to the effective grid quantity to allow the regulated effective grid quantity to be equal to the sampling point quantity or the difference between the effective grid quantity and the sampling point quantity is smaller than a preset difference threshold; determining the soil sampling points of the soil sampling area according to the regulated grids of the soil sampling area. By the method, the soil sampling points of appointed quantity can be automatically determined in an appointed range, the interval among the soil sampling points is strictly identical, and the distribution efficiency of the soil sampling points is increased.
Description
Technical field
The present invention relates to soil science technology, determine method and device in particular to a kind of soil sampling point.
Background technology
Determine that soil sampling point is to obtain soil physics and the important technical of chemical characteristic.At present, net is generally used
The method of formatting determines soil sampling point, and it is more common in space interpolation analysis and other statistical methods.
In soil sampling is put into practice, it is often necessary to the soil designing specified quantity according to sampling funds within the specified range is adopted
Sampling point, but due to the irregular characteristic of sample area shape, cause the spacing being difficult to accurately determine sampling network lattice point.At present, logical
Sampling network lattice point is tentatively generated, then according to sample region shape and fall into sampling usually through rough estimate mesh point spacing mode
The quantity of sampled point grid in district, carries out manual setting fractional-sample dot grid position, supplements or deletes fractional-sample dot grid,
Thus generate specified quantity or the sampled point grid of approximation specified quantity.But the method efficiency is low, need a large amount of manual operations, raw
The sampling grid dot spacing become can not be essentially equal, it is impossible to the automatization realizing sampled point grid generates.
Be only used for strengthening the understanding to the background of the present invention in information above-mentioned disclosed in described background section, therefore it
Can include not constituting the information to prior art known to persons of ordinary skill in the art.
Summary of the invention
In view of this, the present invention provides a kind of soil sampling point to determine method and device, it is possible in the range of realizing specifying certainly
The dynamic soil sampling point determining specified quantity, and ensure that the spacing between each soil sampling dot grid is the most equal, improve soil
The efficiency that sampling is layouted.
Other characteristics of the present invention and advantage will be apparent from by detailed description below, or partially by the present invention
Practice and acquistion.
According to an aspect of the present invention, it is provided that a kind of soil sampling point determines method, comprises determining that soil sampling district；
Area according to described soil sampling district and the sampled point quantity specified, determine the initial soil sampled point in described soil sampling district
Spacing；According to described initial soil sampled point spacing, described soil sampling district is converted to raster data, and determines presently described
Effective grid quantity of soil sampling district grid；According to effective grid quantity of presently described soil sampling district grid, adjust and work as
The size of front described soil sampling district grid, so that effective grid quantity of the described soil sampling district grid after Tiao Zhenging is with described
The difference threshold that equal or described effective grid quantity and described sampled point quantity the difference of sampled point quantity is preset less than one；
And according to the described soil sampling district grid after described adjustment, determine each soil sampling point in described soil sampling district.
According to an embodiment of the present invention, described soil sampling district includes at least one piece of irregular sampling district.
According to an embodiment of the present invention, determine that soil sampling district includes: obtain described according to a GIS-Geographic Information System
The vectorial geographical information data in soil sampling district.
According to an embodiment of the present invention, described initial soil sampled point spacing calculates according to following formula:
Wherein, the described soil sampling that according to S_Area, the vectorial geographical information data in described soil sampling district calculates
The area in district, P_Count is described sampled point quantity, and D_CellSize is described initial soil sampled point spacing.
According to an embodiment of the present invention, according to effective grid quantity of presently described soil sampling district grid, adjust
The size of presently described soil sampling district grid, so that effective grid quantity of the described soil sampling district grid after Tiao Zhenging and institute
State the difference threshold that equal or described effective grid quantity and described sampled point quantity the difference of sampled point quantity is preset less than one
Value includes: step a, compares effective grid quantity of presently described soil sampling district grid and described sampled point quantity, if worked as
Before described effective grid quantity is equal with described sampled point quantity or presently described effective grid quantity and described sampled point quantity
Difference less than described difference threshold, then perform step b；Otherwise, step c is performed；Step b, determines presently described soil sampling
District's grid is the described soil sampling district grid after adjusting；Step c, if presently described effective grid quantity is more than described sampling
Put quantity, or presently described effective grid quantity is more than described sampled point quantity and presently described effective grid quantity with described
The difference of sampled point quantity is more than described difference threshold, then increase the size of presently described soil sampling district grid；If it is current
Described effective grid quantity is less than described sampled point quantity, or presently described effective grid quantity is less than described sampled point quantity
And the difference of presently described effective grid quantity and described sampled point quantity is more than described difference threshold, then reduce presently described soil
The size of earth sample region grid；Described increase or reduce presently described soil sampling district grid size adjustment size according under
State formula to calculate:
CellSize_change=sign (T_Count-P_Count) * Step_Change；
Wherein, CellSize_change is described adjustment size, and P_Count is described sampled point quantity, and CellSize is
The size of presently described soil sampling district grid, T_Count is presently described effective grid quantity, and Step_Change is one pre-
If adjustment step-length；And
Step d, according to the size of the described soil sampling district grid after adjusting, determines new described soil sampling district grid
Lattice, and to arrange described new described soil sampling district grid be presently described soil sampling district grid, again performs step b；
According to an embodiment of the present invention, according to the described soil sampling district grid after described adjustment, described soil is determined
Each soil sampling point of earth sample region includes: the described soil sampling district grid after described adjustment is converted to vectorial geographical information
Data, the center position coordinates of the most each grid is that each soil sampling point in described soil sampling district is believed at described vectorial geographical
Coordinate in breath data.
According to a further aspect in the invention, it is provided that a kind of soil sampling point determines device, including: sample region determines module,
For determining soil sampling district；Spacing determines module, is used for the area according to described soil sampling district and the sampling number specified
Amount, determines the initial soil sampled point spacing in described soil sampling district；Initial determines module, for according to described initial soil
Earth sampled point spacing, is converted to raster data by described soil sampling district, and determines having of presently described soil sampling district grid
Effect grid quantity；Grid adjusting module, for the effective grid quantity according to presently described soil sampling district grid, adjusts current
The size of described soil sampling district grid, so that effective grid quantity of the described soil sampling district grid after Tiao Zhenging is adopted with described
The difference threshold that equal or described effective grid quantity and described sampled point quantity the difference of sampling point quantity is preset less than one；With
And sampled point determines module, for according to the described soil sampling district grid after described adjustment, determine described soil sampling district
Each soil sampling point.
According to an embodiment of the present invention, described soil sampling district includes at least one piece of irregular sampling district.
According to an embodiment of the present invention, described grid adjusting module includes: quantity comparison sub-module, is used for comparing and works as
Effective grid quantity of front described soil sampling district grid and described sampled point quantity；Grid determines submodule, for when described
Quantity comparison sub-module more the currently active grid quantity is equal with described sampled point quantity or presently described effective grid quantity
When being less than described difference threshold with the difference of described sampled point quantity, after determining that presently described soil sampling district grid is adjustment
Described soil sampling district grid；Spacing adjusts submodule, for when the more presently described effective grid of described quantity comparison sub-module
Lattice quantity is more than described sampled point quantity, or presently described effective grid quantity is more than described sampled point quantity and presently described
When effectively the difference of grid quantity and described sampled point quantity is more than described difference threshold, increase presently described soil sampling district grid
The size of lattice；Or when the more presently described effective grid quantity of described quantity comparison sub-module is less than described sampled point quantity,
Or presently described effective grid quantity is less than described sampled point quantity and presently described effective grid quantity and described sampled point
When the difference of quantity is more than described difference threshold, reduce the size of presently described soil sampling district grid；Described increase or reduction
The adjustment size of the size of presently described soil sampling district grid calculates according to following formula:
CellSize_change=sign (T_Count-P_Count) * Step_Change；
Wherein, CellSize_change is described adjustment size, and P_Count is described sampled point quantity, and CellSize is
The size of presently described soil sampling district grid, T_Count is presently described effective grid quantity, and Step_Change is one pre-
If adjustment step-length；And
New grid determines submodule, for the size according to the described soil sampling district grid after adjusting, determines new institute
State soil sampling district grid, and to arrange described new described soil sampling district grid be presently described soil sampling district grid.
According to an embodiment of the present invention, described sampled point determines that module includes: raster data transform subblock, is used for
Described soil sampling district grid after described adjustment is converted to vectorial geographical information data, the center of the most each grid
Coordinate is each soil sampling point in described soil sampling district coordinate in described vectorial geographical information data
Soil sampling point according to the present invention determines method, it is possible to achieve automatically generate specified quantity within the specified range
Spacing between soil sampling dot grid, and each grid is the most equal, substantially increases the efficiency that soil sampling is layouted.Additionally,
The soil sampling point that the present invention provides determines that method applies also for irregular sampling district, solves in prior art for irregularly
Sample region is arranged specifies sampled point and the needs that produce are the most hand-manipulated and the unequal problem of spacing between each sampled point.
It should be appreciated that it is only exemplary that above general description and details hereinafter describe, can not be limited this
Invention.
Accompanying drawing explanation
Describing its example embodiment in detail by referring to accompanying drawing, above and other target of the present invention, feature and advantage will
Become more fully apparent.
Fig. 1 is the flow chart determining method according to a kind of soil sampling point shown in an illustrative embodiments.
Fig. 2 is the flow chart determining method according to the another kind of soil sampling point shown in an illustrative embodiments.
Fig. 3 A is the schematic diagram according to the soil sampling district shown in an example.
Fig. 3 B is the distribution schematic diagram according to the soil sampling district initial samples point grid shown in an example.
Fig. 4 is the block diagram determining device according to a kind of soil sampling point shown in an illustrative embodiments.
Fig. 5 is the block diagram determining device according to the another kind of soil sampling point shown in an illustrative embodiments.
Detailed description of the invention
It is described more fully with example embodiment referring now to accompanying drawing.But, example embodiment can be with multiple shape
Formula is implemented, and is not understood as limited to example set forth herein；On the contrary, it is provided that these embodiments make the present invention will more
Fully and completely, and by the design of example embodiment those skilled in the art is conveyed to all sidedly.Accompanying drawing is only the present invention
Schematic illustrations, be not necessarily drawn to scale.Reference identical in figure represents same or similar part, thus
Repetition thereof will be omitted.
Additionally, described feature, structure or characteristic can be combined in one or more enforcement in any suitable manner
In mode.In the following description, it is provided that many details thus be given and embodiments of the present invention fully understood.So
And, it will be appreciated by persons skilled in the art that and can put into practice technical scheme and omit in described specific detail one
Or more, or other method, constituent element, device, step etc. can be used.In other cases, it is not shown in detail or describes
Known features, method, device, realize or operate avoiding that a presumptuous guest usurps the role of the host and each aspect of the present invention is thickened.
Fig. 1 is the flow chart determining method according to a kind of soil sampling point shown in an illustrative embodiments.Such as Fig. 1 institute
Showing, soil sampling point determines that method 10 includes:
In step s 102, soil sampling district is determined.
This soil sampling district can include at least one piece of irregular sampling district.
In step S104, according to the area in soil sampling district and the sampled point quantity specified, determine that initial soil is sampled
Dot spacing.
In step s 106, to specify lattice dimensions size, soil sampling district is converted to raster data.
During first operating procedure S106, the initial soil sampled point spacing determined using step S104 is as indication in this step
" specify lattice dimensions size ".
In step S108, determine effective grid quantity of current soil sample region grid.
In step s 110, according to effective grid quantity of current soil sample region grid, adjust current grid size big
Little so that after Tiao Zhenging lattice dimensions generate sample region raster data effective grid quantity with the sampled point quantity specified more
Close or equal.
Wherein close to such as referring to that effective grid quantity difference with sampled point quantity is less than a difference threshold specified.
In step S112, according to adjusting post-sampling point grid, determine the soil sampling point in this soil sampling district.
Step S106, step S108 and step S110 may need repeatedly, until the sampling number of actual arrangement with
The sampled point quantity specified is consistent, or the difference between the two is less than the above-mentioned difference threshold specified.
The soil sampling point that the present invention provides determines method, it is possible to achieve automatically generate specified quantity within the specified range
Spacing between soil sampling dot grid, and each grid is the most equal, substantially increases the efficiency that soil sampling is layouted.Additionally,
The soil sampling point that the present invention provides determines that method applies also for irregular sampling district, solves in prior art for irregularly
Sample region is arranged specifies sampled point and the needs that produce are the most hand-manipulated and the unequal problem of spacing between each sampled point.
It will be clearly understood that present disclosure describe how to be formed and use particular example, but the principle of the present invention is not limited to
Any details of these examples.On the contrary, teaching based on present disclosure, these principles can be applied to many other
Embodiment.
Fig. 2 is the flow chart determining method according to the another kind of soil sampling point shown in an illustrative embodiments.Such as Fig. 2
Shown in, soil sampling point determines that method 20 includes:
In step S202, the geography information obtained according to GIS-Geographic Information System, determine soil sampling district.
In the present embodiment, can be according to the geography information in the data acquisition soil sampling district in GIS-Geographic Information System, its
In the geography information in soil sampling district that gets be vector data.
Fig. 3 A is the schematic diagram according to the soil sampling district shown in an example.As shown in Figure 3A, soil sampling district A is (in figure
Dash area) include at least one piece of irregular sampling district A1-A3.
In step S204, determine the area in soil sampling district.
According to the vectorial geographical information in the soil sampling district obtained from GIS-Geographic Information System, calculate this soil sampling district
Area.
In step S206, according to the area in the soil sampling district determined and the sampled point quantity specified, determine initial soil
The spacing of earth sampled point.
For example, it is assumed that the area in the soil sampling district determined is S_Area, it is intended that sampled point quantity be P_Count, then
May determine that the spacing of initial soil sampled point:
In step S208, according to the lattice dimensions value (being designated as CellSize) specified, by the vector number in soil sampling district
The raster data that the grid length of side is CellSize is converted to according to geography information, so that it is determined that initial samples in the range of soil sampling district
The distribution of some grid.
It should be noted that when vector data is converted to raster data, it should be ensured that grid central point falls in vector data
Portion.For the first time during operating procedure S208 CellSize value determine in being step S206 initial soil sampled point space D _
CellSize。
Fig. 3 B is the distribution schematic diagram according to the soil sampling district initial samples point grid shown in an example.Such as Fig. 3 B institute
Showing, the grid wherein comprising point " " is the raster data that the vector data conversion of soil sampling district generates, and wherein " " is sampling
The central point of district's grid, other are then non-sampled district grid.Non-sampled district grid is typically expressed as " no data grid " or invalid
Grid.Therefore press effective grid quantity that sampled point spacing CellSize is generated by the conversion of sample region vector data, be spacing
The quantity layouted for the grid sampling of CellSize, the position of " " is the position of grid sampled point.
In step S210, the quantity of the effective grid generated in statistic procedure S208, it is designated as T_Count, this quantity is i.e.
For initial samples point grid point quantity.
In step S212, according to T_Count calculated in step S210, adjust current soil sample region grid
Size, so that before after Tiao Zhenging, the quantity of effective grid is compared to adjusting, being more nearly with the sampled point quantity specified, or
Person make to adjust after the quantity of effective grid equal to the sampled point quantity specified.
Such as, compare the quantity of the currently active grid and the sampled point quantity specified, if the quantity of the currently active grid
More than the sampled point quantity specified, then increase the size of current grid, if the quantity of the currently active grid is less than adopting of specifying
Sampling point quantity, then reduce the size of current grid.Adjustment size CellSize_change of lattice dimensions can be such as:
CellSize_change=sign (T_Count-P_Count) * Step_Change
The grid sampling number generated when the lattice dimensions obtained during wherein T_Count is step S210 is CellSize
Amount；Sign (T_Count-P_Count) is a discriminant function, and when T_Count-P_Count is more than zero, this function value is 1,
When T_Count-P_Count is less than zero, this function value is-1, when T_Count-P_Count is equal to zero, and this function value
It is 0；Step_Change is an adjustment step-length preset, and its length can be configured according to the actual requirements, and the present invention is not with this
It is limited.
Lattice dimensions CellSize after adjustment adds CellSize_change equal to current grid size.
Repeat S208 to step S212, until the quantity of the currently active grid is equal to the sampled point quantity specified, or
The quantity of the currently active grid of person is less than, with the difference specifying sampled point quantity, the difference threshold preset, and stops above-mentioned adjustment operation,
And generate the sample region raster data meeting quantitative requirement.
In step S214, according to meeting the sample region raster data of quantitative requirement, determine the soil in this soil sampling district
Sampled point.
Stopping repeating S208 to after step S212, the soil sampling district raster data of generation, this raster data effective
Grid quantity is equal with the sampled point quantity specified, or the two difference is less than an above-mentioned difference threshold specified.
Further according to this raster data, determine each soil sampling point in final soil sampling district.
Specifically, above-mentioned raster data is converted to again vector point data, the center position coordinates of the most each grid
Corresponding point is each soil sampling point in soil sampling district.
It will be appreciated by those skilled in the art that all or part of step realizing above-mentioned embodiment is implemented as being held by CPU
The computer program of row.When this computer program is performed by CPU, it is above-mentioned that the said method that performing the present invention provides is limited
Function.Described program can be stored in a kind of computer-readable recording medium, and this storage medium can be read only memory,
Disk or CD etc..
Further, it should be noted that above-mentioned accompanying drawing is only according to included by the method for exemplary embodiment of the invention
Process schematically illustrates rather than limits purpose.It can be readily appreciated that above-mentioned process shown in the drawings is not intended that or limits these
The time sequencing processed.It addition, be also easy to understand, these process can be such as either synchronously or asynchronously to perform in multiple modules
's.
Following for apparatus of the present invention embodiment, may be used for performing the inventive method embodiment.Real for apparatus of the present invention
Execute the details not disclosed in example, refer to the inventive method embodiment.
Fig. 4 is the block diagram determining device according to a kind of soil sampling point shown in an illustrative embodiments.Such as Fig. 4 institute
Showing, soil sampling point determines that device 30 includes: sample region determines that module 302, spacing determine that module 304, initial determine mould
Block 306, grid adjusting module 308 and sampled point determine module 310.
Wherein, sample region determines that module 302 is for determining soil sampling district.
This soil sampling district such as can include at least one piece of irregular sampling district as shown in Figure 3A.
Spacing determines that module 304, for the area according to described soil sampling district and the sampled point quantity specified, determines institute
State the initial soil sampled point spacing in soil sampling district.
Described initial soil sampled point spacing can calculate according to following formula:
Wherein, the described soil sampling that according to S_Area, the vectorial geographical information data in described soil sampling district calculates
The area in district, P_Count is described sampled point quantity, and D_CellSize is described initial soil sampled point spacing.
Initial determines that module 306 is for according to described initial soil sampled point spacing, turning described soil sampling district
It is changed to raster data, and determines effective grid quantity of presently described soil sampling district grid.
Grid adjusting module 308, for the effective grid quantity according to presently described soil sampling district grid, adjusts current
The size of described soil sampling district grid, so that effective grid quantity of the described soil sampling district grid after Tiao Zhenging is adopted with described
The difference threshold that equal or described effective grid quantity and described sampled point quantity the difference of sampling point quantity is preset less than one.
Sampled point determines that module 310 is for according to the described soil sampling district grid after described adjustment, determining described soil
Each soil sampling point of sample region.
The soil sampling point that the present invention provides determines device, it is possible to achieve automatically generate specified quantity within the specified range
Soil sampling point grid, and the spacing between each grid is the most equal, substantially increases the efficiency that soil sampling is layouted.Additionally,
The soil sampling point that the present invention provides determines that method applies also for irregular sampling district, solves in prior art for irregularly
Sample region is arranged specifies sampled point and the needs that produce are the most hand-manipulated and the unequal problem of spacing between each sampled point.
Fig. 5 is the block diagram determining device according to the another kind of soil sampling point shown in an illustrative embodiments.Such as Fig. 5 institute
Showing, soil sampling point determines with the soil sampling point shown in Fig. 4, device 40 determines that device 30 difference is:
Sample region determines that module 302 may include that vector data obtains submodule 3022 for according to a geography information system
System obtains the vectorial geographical information data in described soil sampling district.
Grid adjusting module 308 may include that quantity comparison sub-module 3084, grid determine that submodule 3086, spacing are adjusted
Whole submodule 3088 and new grid determine submodule 3090.
Wherein, quantity comparison sub-module 3084 for relatively presently described soil sampling district grid effective grid quantity and
Described sampled point quantity；
Grid determines that submodule 3086 is for when described quantity comparison sub-module more the currently active grid quantity is with described
When equal or presently described effective grid quantity and described sampled point quantity the difference of sampled point quantity is less than described difference threshold,
Determine that presently described soil sampling district grid is the described soil sampling district grid after adjusting.
Spacing adjusts submodule 3088 and is used for when quantity comparison sub-module 3084 relatively presently described effective grid quantity is big
In described sampled point quantity, or presently described effective grid quantity is more than described sampled point quantity and presently described effective grid
When the difference of quantity and described sampled point quantity is more than described difference threshold, increase the chi of presently described soil sampling district grid
Very little；Or when the more presently described effective grid quantity of described quantity comparison sub-module is less than described sampled point quantity, or work as
Front described effective grid quantity is less than described sampled point quantity and presently described effective grid quantity and described sampled point quantity
When difference is more than described difference threshold, reduce the size of presently described soil sampling district grid；Described increase or reduce current institute
The adjustment size of the size stating soil sampling district grid calculates according to following formula:
CellSize_change=sign (T_Count-P_Count) * Step_Change
Wherein CellSize_change is described adjustment size, and P_Count is described sampled point quantity, and CellSize is
The size of presently described soil sampling district grid, T_Count is presently described effective grid quantity, and Step_Change is one pre-
If adjustment step-length.New grid determine submodule 3090 for the size according to the described soil sampling district grid after adjusting, really
The newest fixed described soil sampling district grid, and to arrange described new described soil sampling district grid be presently described soil sampling district
Grid.
Sampled point determines that module 310 may include that raster data transform subblock 3102, for by the institute after described adjustment
Stating soil sampling district grid and be converted to vectorial geographical information data, the center position coordinates of the most each grid is that described soil is adopted
The each soil sampling point in sample district coordinate in described vectorial geographical information data.
It should be noted that the block diagram shown in above-mentioned accompanying drawing is functional entity, it is not necessary to must with physically or logically
Independent entity is corresponding.Can use software form to realize these functional entitys, or in one or more hardware modules or
Integrated circuit realizes these functional entitys, or realizes in heterogeneous networks and/or processor device and/or microcontroller device
These functional entitys.
Through the above description of the embodiments, those skilled in the art is it can be readily appreciated that example described herein is implemented
Mode can be realized by software, it is also possible to realizes by the way of software combines necessary hardware.Therefore, according to the present invention
The technical scheme of embodiment can embody with the form of software product, this software product can be stored in one non-volatile
In property storage medium (can be CD-ROM, USB flash disk, portable hard drive etc.) or on network, including some instructions so that a calculating
Equipment (can be personal computer, server, mobile terminal or the network equipment etc.) performs according to embodiment of the present invention
Method.
More than it is particularly shown and described the illustrative embodiments of the present invention.It should be appreciated that the present invention does not limits
In detailed construction described herein, set-up mode or implementation method；It is included in claims on the contrary, it is intended to contain
Spirit and scope in various amendments and equivalence arrange.
Claims (10)
1. a soil sampling point determines method, it is characterised in that including:
Determine soil sampling district；
Area according to described soil sampling district and the sampled point quantity specified, determine that the initial soil in described soil sampling district is adopted
Sampling point spacing；
According to described initial soil sampled point spacing, described soil sampling district is converted to raster data, and determines presently described
Effective grid quantity of soil sampling district grid；
According to effective grid quantity of presently described soil sampling district grid, adjust the size of presently described grid, so that adjusting
After effective grid quantity of described soil sampling district grid is equal with described sampled point quantity or described effective grid quantity
The difference threshold preset less than one with the difference of described sampled point quantity；And
According to the described soil sampling district grid after described adjustment, determine each soil sampling point in described soil sampling district.
Method the most according to claim 1, it is characterised in that described soil sampling district includes at least one piece of irregular sampling
District.
Method the most according to claim 1 and 2, it is characterised in that determine that soil sampling district includes: according to a geography information
System obtains the vectorial geographical information data in described soil sampling district.
Method the most according to claim 3, it is characterised in that described initial soil sampled point spacing is according to following formula meter
Calculate:
Wherein, the described soil sampling district that according to S_Area, the vectorial geographical information data in described soil sampling district calculates
Area, P_Count is described sampled point quantity, and D_CellSize is described initial soil sampled point spacing.
Method the most according to claim 3, it is characterised in that according to effective grid of presently described soil sampling district grid
Quantity, adjusts the size of presently described soil sampling district grid, so that effective grid of the described soil sampling district grid after Tiao Zhenging
Lattice quantity is equal with described sampled point quantity or described effective grid quantity is pre-less than one with the difference of described sampled point quantity
If difference threshold include:
Effective grid quantity of step a, relatively presently described soil sampling district grid and described sampled point quantity, if current institute
State that effective grid quantity is equal with described sampled point quantity or the difference of presently described effective grid quantity and described sampled point quantity
Value less than described difference threshold, then performs step b；Otherwise, step c is performed；
Step b, determines that presently described soil sampling district grid is the described soil sampling district grid after adjusting；
Step c, if presently described effective grid quantity is more than described sampled point quantity, or presently described effective grid quantity
It is more than described difference threshold more than the difference of described sampled point quantity and presently described effective grid quantity and described sampled point quantity
Value, then increase the size of presently described soil sampling district grid；If presently described effective grid quantity is less than described sampled point
Quantity, or presently described effective grid quantity is less than described sampled point quantity and presently described effective grid quantity is adopted with described
The difference of sampling point quantity is more than described difference threshold, then reduce the size of presently described soil sampling district grid；Described increase or
The adjustment size of the size reducing presently described soil sampling district grid calculates according to following formula:
CellSize_change=sign (T_Count-P_Count) * Step_Change；
Wherein, CellSize_change is described adjustment size, and P_Count is described sampled point quantity, and CellSize is current
The size of described soil sampling district grid, T_Count is presently described effective grid quantity, Step_Change be one preset
Adjust step-length；And
Step d, according to the size of the described soil sampling district grid after adjusting, determines new described soil sampling district grid, and
Arranging described new described soil sampling district grid is presently described soil sampling district grid, again performs step b.
Method the most according to claim 3, it is characterised in that according to the described soil sampling district grid after described adjustment,
Determine that each soil sampling point in described soil sampling district includes: be converted to the described soil sampling district grid after described adjustment vow
Amount geographic information data, the center position coordinates of the most each grid is that each soil sampling point in described soil sampling district is described
Coordinate in vectorial geographical information data.
7. a soil sampling point determines device, it is characterised in that including:
Sample region determines module, is used for determining soil sampling district；
Spacing determines module, is used for the area according to described soil sampling district and the sampled point quantity specified, determines described soil
The initial soil sampled point spacing of sample region；
Initial determines module, for according to described initial soil sampled point spacing, described soil sampling district is converted to grid
Lattice data, and determine effective grid quantity of presently described soil sampling district grid；
Grid adjusting module, for the effective grid quantity according to presently described soil sampling district grid, adjusts presently described soil
The size of earth sample region grid, so that effective grid quantity of the described soil sampling district grid after Tiao Zhenging and described sampling number
Measure the difference threshold that equal or described effective grid quantity is default less than with the difference of described sampled point quantity；And
Sampled point determines module, for according to the described soil sampling district grid after described adjustment, determines described soil sampling district
Each soil sampling point.
Device the most according to claim 7, it is characterised in that described soil sampling district includes at least one piece of irregular sampling
District.
9. according to the device described in claim 7 or 8, it is characterised in that described grid adjusting module includes:
Quantity comparison sub-module, for effective grid quantity and the described sampling number of relatively presently described soil sampling district grid
Amount；
Grid determines submodule, for when described quantity comparison sub-module more the currently active grid quantity and described sampling number
When measuring the equal or presently described effective grid quantity difference with described sampled point quantity less than described difference threshold, determine current
Described soil sampling district grid is the described soil sampling district grid after adjusting；
Spacing adjusts submodule, adopts described in being more than when the more presently described effective grid quantity of described quantity comparison sub-module
Sampling point quantity, or presently described effective grid quantity is more than described sampled point quantity and presently described effective grid quantity and institute
When stating the difference of sampled point quantity more than described difference threshold, increase the size of presently described soil sampling district grid；Or work as
The more presently described effective grid quantity of described quantity comparison sub-module be less than described sampled point quantity, or presently described effectively
Grid quantity is more than institute less than the difference of described sampled point quantity and presently described effective grid quantity and described sampled point quantity
When stating difference threshold, reduce the size of presently described soil sampling district grid；Described increase or reduce presently described soil sampling
The adjustment size of the size of district's grid calculates according to following formula:
CellSize_change=sign (T_Count-P_Count) * Step_Change；
Wherein, CellSize_change is described adjustment size, and P_Count is described sampled point quantity, and CellSize is current
The size of described soil sampling district grid, T_Count is presently described effective grid quantity, Step_Change be one preset
Adjust step-length；And
New grid determines submodule, for the size according to the described soil sampling district grid after adjusting, determines new described soil
Earth sample region grid, and to arrange described new described soil sampling district grid be presently described soil sampling district grid.
10. according to the device described in claim 7 or 8, it is characterised in that described sampled point determines that module includes: raster data
Transform subblock, for being converted to vectorial geographical information data, the most often by the described soil sampling district grid after described adjustment
The center position coordinates of individual grid is that each soil sampling point in described soil sampling district is in described vectorial geographical information data
Coordinate.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610829349.5A CN106323670B (en) | 2016-09-18 | 2016-09-18 | Soil sampling point determines method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610829349.5A CN106323670B (en) | 2016-09-18 | 2016-09-18 | Soil sampling point determines method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106323670A true CN106323670A (en) | 2017-01-11 |
| CN106323670B CN106323670B (en) | 2019-08-09 |
Family
ID=57787282
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610829349.5A Expired - Fee Related CN106323670B (en) | 2016-09-18 | 2016-09-18 | Soil sampling point determines method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106323670B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110411778A (en) * | 2019-06-17 | 2019-11-05 | 淮北矿业股份有限公司 | A kind of coal slime method of sampling |
| CN110457422A (en) * | 2019-08-20 | 2019-11-15 | 重庆壤科农业数据服务有限公司 | Collecting soil sample auto-distribution dot system and method |
| CN110623555A (en) * | 2019-09-19 | 2019-12-31 | 杭州九阳小家电有限公司 | Food material size calculation method in oven, size identification device and oven |
| CN112200361A (en) * | 2020-09-30 | 2021-01-08 | 杭州鲁尔物联科技有限公司 | Landslide prediction method, device, equipment and storage medium |
| CN112200356A (en) * | 2020-09-30 | 2021-01-08 | 杭州鲁尔物联科技有限公司 | Landslide prediction method, device, equipment and storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102645350A (en) * | 2012-03-08 | 2012-08-22 | 浙江大学 | Soil-sampling method based on high-resolution satellite remote sensing data |
| CN103196698A (en) * | 2013-03-20 | 2013-07-10 | 浙江大学 | Soil sampling method based on near-earth sensor technology |
-
2016
- 2016-09-18 CN CN201610829349.5A patent/CN106323670B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102645350A (en) * | 2012-03-08 | 2012-08-22 | 浙江大学 | Soil-sampling method based on high-resolution satellite remote sensing data |
| CN103196698A (en) * | 2013-03-20 | 2013-07-10 | 浙江大学 | Soil sampling method based on near-earth sensor technology |
Non-Patent Citations (3)
| Title |
|---|
| WANG HONGBIN等: "Determining optimal density of grid soil-sampl ing", 《农业工程学报》 * |
| 国家环境保护总局: "《中华人民共和国环境保护行业标准 HJ/T166-2004》", 31 December 2004 * |
| 王秀等: "土壤采样密度及栅格尺度对处方图结果的影响", 《中墨农业工程学会2005年学术年会论文集》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110411778A (en) * | 2019-06-17 | 2019-11-05 | 淮北矿业股份有限公司 | A kind of coal slime method of sampling |
| CN110457422A (en) * | 2019-08-20 | 2019-11-15 | 重庆壤科农业数据服务有限公司 | Collecting soil sample auto-distribution dot system and method |
| CN110623555A (en) * | 2019-09-19 | 2019-12-31 | 杭州九阳小家电有限公司 | Food material size calculation method in oven, size identification device and oven |
| CN112200361A (en) * | 2020-09-30 | 2021-01-08 | 杭州鲁尔物联科技有限公司 | Landslide prediction method, device, equipment and storage medium |
| CN112200356A (en) * | 2020-09-30 | 2021-01-08 | 杭州鲁尔物联科技有限公司 | Landslide prediction method, device, equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106323670B (en) | 2019-08-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106323670A (en) | Method and device for determining soil sampling points | |
| Christen et al. | A general purpose sampling algorithm for continuous distributions (the t-walk) | |
| Bhend et al. | An ensemble-based approach to climate reconstructions | |
| Tabak et al. | A family of nonparametric density estimation algorithms | |
| Castro et al. | Minimax bounds for active learning | |
| CN109886388A (en) | A kind of training sample data extending method and device based on variation self-encoding encoder | |
| CN106407682B (en) | Soil sampling point determines method and device | |
| DE102016221283A1 (en) | REQUESTING EMISSION SOURCES USING CFD AND SATELLITE DATA | |
| US11954115B2 (en) | Methods and systems for wavelet based representation | |
| US20090037147A1 (en) | Fast intrinsic mode decomposition of time series data with sawtooth transform | |
| Chen et al. | Quantile and tolerance-interval estimation in simulation | |
| Butucea et al. | Spectral thresholding quantum tomography for low rank states | |
| US20130185238A1 (en) | Splitting device, splitting method, and recording medium | |
| Obregón et al. | A deterministic geometric representation of temporal rainfall: sensitivity analysis for a storm in Boston | |
| US8457935B2 (en) | Data processing method for sampling data from sets of original data | |
| Lazić et al. | On sub-geometric ergodicity of diffusion processes | |
| KR20190125575A (en) | Leakage distribution estimating system of semiconductor and analysis method of the same | |
| Kent et al. | Simulation for the complex Bingham distribution | |
| CN114492809A (en) | Uncertainty analysis method for marine environment design parameter calculation model | |
| Opitz et al. | The spectrogram: A threshold-based inferential tool for extremes of stochastic processes | |
| Ali et al. | Statistical downscaling of precipitation and temperature using gene expression programming | |
| Ma et al. | Calibrated Uncertainty Quantification for Operator Learning via Conformal Prediction | |
| Weir et al. | Spatial modelling for binary data using␣ a␣ hidden conditional autoregressive Gaussian process: a multivariate extension of the probit model | |
| Beranger et al. | Exploratory data analysis for moderate extreme values using non-parametric Kernel methods | |
| Xiong et al. | Ensemble particle filter with posterior Gaussian resampling |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190809 Termination date: 20200918 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |