CN104360347A - Method and device for monitoring crop harvesting progress - Google Patents
Method and device for monitoring crop harvesting progress Download PDFInfo
- Publication number
- CN104360347A CN104360347A CN201410610674.3A CN201410610674A CN104360347A CN 104360347 A CN104360347 A CN 104360347A CN 201410610674 A CN201410610674 A CN 201410610674A CN 104360347 A CN104360347 A CN 104360347A
- Authority
- CN
- China
- Prior art keywords
- polarization
- scattering
- remote sensing
- plot
- sensing image
- 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
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9076—Polarimetric features in SAR
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9004—SAR image acquisition techniques
- G01S13/9005—SAR image acquisition techniques with optical processing of the SAR signals
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Image Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a method and a device for monitoring crop harvesting progress. The method includes acquiring complete polarimetric synthetic aperture radar (SAR) remote-sensing images in a monitoring area in the harvest season of crops, decomposing the remote-sensing images by polarization to obtain polarization parameters of each pixel; extracting the boundary of each plot in the remote-sensing images, and acquiring the average polarization parameter of all pixels in each plot according to the polarization parameter of each pixel; judging whether the crops on each plot are harvested or not according to the average polarization parameters. The method solves the problem that the acquired harvesting progress of the optical remote-sensing monitoring has an effect on the monitoring accuracy of the harvesting progress of the crops since the harvester crops are aired on the spots. The harvesting progress of the crops can be quickly and accurately monitored in large area, high-efficient harvesting is guaranteed to the greatest extent, risk in harvesting output due to weather and the like is reduced, and the method has great significance in realizing high yield, high quality and high efficiency of the crop harvesting.
Description
Technical field
The present invention relates to radar remote sensing applied technical field, be specifically related to a kind of method and device of monitoring crops harvesting progress.
Background technology
Crops harvesting is the key activities directly determining final effectively grain yield.Too early harvesting can reduce grain yield, and crossing harvesting in evening will increase the risk of adverse weather impact.But duration harvest time is very of short duration, considers the factors such as weather, mechanical manpower and succession crop sowing, many times need " rush-harvest ", thus very high to the ageing requirement of monitoring information.Therefore, monitoring the harvesting progress of crop harvest time on a large scale in time, will contributing to optimizing harvesting strategy, carrying out the more rationally Resource allocation and smoothing such as effective machinery, manpower, transport, following process, guarantee efficient reaping and ultimate capacity to greatest extent.
Traditional crop harvesting progress msg obtains and mainly relies on manpower on-site inspection, take time and effort, inefficiency, and monitoring range is limited, cannot obtains the macroscopic information in regional extent, and then cannot gather in arrangement to entirety and make overall planning and instruct.Macroscopical fast monitored of the crop harvesting progress that develops into of remote sensing technology provides a kind of effective means.Remote sensing has that coverage is large, detect cycle is short, currency strong, the feature such as expense cost is low, can large area repeatedly observe ground, can the crop dynamic change of monitored area.At present, utilize remote sensing technology monitoring crop harvesting aspect, researcher carries out some preliminary trials.Existing research utilizes optical remote sensing image more, according to the bare area after crop harvesting be the difference of standing crops in spectral signature, judges the harvesting situation in farmland.
But, adopt said method owing to continuing after crop harvesting in a lot of situation to stay original place airing, cause harvesting plot and do not gather in plot difference on remote optical sensing not significantly, having had a strong impact on remote optical sensing and having utilized spectral signature to monitor the validity of harvesting situation.
Summary of the invention
For defect of the prior art, the invention provides a kind of method and device of monitoring crops harvesting progress, solving remote optical sensing monitoring harvesting progress data obtains because continuing to stay original place airing for after harvesting, the problem of the accuracy of impact monitoring crops harvesting progress.
First aspect, the invention provides a kind of method of monitoring crops harvesting progress, comprising:
Obtain the polarimetric synthetic aperture radar remote sensing image of monitored area in the crop harvesting phase, polarization decomposing is carried out to described remote sensing image, obtains the polarization parameter of each pixel;
Extract the border in each plot in described remote sensing image, obtain the average polarization parameter of all pixels in described each plot according to the polarization parameter of each pixel described;
According to described average polarization parameter, judge whether the crops in described each plot gather in.
Optionally, the described polarimetric synthetic aperture radar remote sensing image obtaining monitored area in the crop harvesting phase, carries out polarization decomposing to described remote sensing image, obtains the polarization parameter of each pixel, comprising:
The polarimetric synthetic aperture radar remote sensing image of monitored area is obtained in the crop harvesting phase;
Pre-service is carried out to described remote sensing image;
Utilize Polarization target decomposition method to carry out polarization decomposing to described pretreated remote sensing image, obtain the polarization parameter of each pixel.
Optionally, described pre-service is carried out to described remote sensing image, comprising:
Radiation calibration is carried out to described SAR remote sensing image;
Carry out multiple look processing and spot to the image after radiation calibration to make an uproar removal;
By the polarization scattering matrix of POLARIZATION CHANNEL data transformations of four after denoising;
Geocoding, topographic correction and geometric accurate correction process are carried out to described polarization scattering matrix image.
Optionally, the described Polarization target decomposition method that utilizes carries out polarization decomposing to described pretreated remote sensing image, obtains the polarization parameter of each pixel, comprising:
Utilize Freeman-Durden three-component polarization decomposing method to carry out polarization decomposing to pretreated remote sensing image, obtain the power level P of surface scattering, even scattering and volume scattering three kinds of scattering components
s, P
dand P
v; Three-component divide solve an equation into:
T＝P
s*T
s+P
d*T
d+P
v*T
v
Wherein, T is the polarization scattering matrix represented by coherence matrix form, T
sthe coherence matrix of presentation surface scattering model, P
sthe power representing atural object surface scattering component; T
dthe coherence matrix representing even scattering model, P
dthe power representing atural object even scattering component; T
vthe coherence matrix representing volume scattering model, P
vthe power representing atural object volume scattering component;
Utilize Cloude decomposition method to carry out polarization decomposing to pretreated remote sensing image, obtain the eigenvalue of maximum λ of described polarization scattering matrix T.
Optionally, the border in each plot in the described remote sensing image of described extraction, comprising:
Utilize database or obtain monitored area map according to OO segmentation, sorting technique;
Vector quantization is carried out to described monitored area map, obtains the border in each plot in monitored area.
Optionally, described according to described average polarization parameter, judge whether the crops in described each plot gather in, comprising:
Obtain the scattering component of object fifty-fifty in described average polarization parameter
with described average eigenvalue of maximum
According to the described scattering component of object fifty-fifty
with described average eigenvalue of maximum
build by the described scattering component of object fifty-fifty
with described average eigenvalue of maximum
the two-dimensional feature space formed;
In described two-dimensional feature space, if the scattering component of object fifty-fifty of all pixels in the described plot obtained
with the average eigenvalue of maximum of pixels all in described plot
size be all less than pre-set threshold value, then represent that the crops in described plot are gathered in.
Second aspect, the invention provides a kind of device of monitoring crops harvesting progress, comprising:
First acquisition module, for obtaining the polarimetric synthetic aperture radar remote sensing image of monitored area in the crop harvesting phase, carries out polarization decomposing to described remote sensing image, obtains the polarization parameter of each pixel;
Second acquisition module, for extracting the border in each plot in described remote sensing image, obtains the average polarization parameter of all pixels in described each plot according to the polarization parameter of each pixel described;
Judge module, for according to described average polarization parameter, judges whether the crops in described each plot gather in.
Optionally, described first acquisition module, comprising:
Image capturing unit, for obtaining the polarimetric synthetic aperture radar remote sensing image of monitored area in the crop harvesting phase;
Pretreatment unit, for carrying out pre-service to described remote sensing image;
First parameter acquiring unit, for utilizing Polarization target decomposition method to carry out polarization decomposing to described pretreated remote sensing image, obtains the polarization parameter of each pixel.
Optionally, described first parameter acquiring unit, for:
Utilize Freeman-Durden three-component polarization decomposing method to carry out polarization decomposing to pretreated remote sensing image, obtain the power level P of surface scattering, even scattering and volume scattering three kinds of scattering components
s, P
dand P
v; Three-component divide solve an equation into:
T＝P
s*T
s+P
d*T
d+P
v*T
v
Wherein, T is the polarization scattering matrix represented by coherence matrix form, T
sthe coherence matrix of presentation surface scattering model, P
sthe power representing atural object surface scattering component; T
dthe coherence matrix representing even scattering model, P
dthe power representing atural object even scattering component; T
vthe coherence matrix representing volume scattering model, P
vthe power representing atural object volume scattering component;
Utilize Cloude decomposition method to carry out polarization decomposing to pretreated remote sensing image, obtain the eigenvalue of maximum λ of described polarization scattering matrix T.
Optionally, described judge module, for:
Obtain the scattering component of object fifty-fifty in described average polarization parameter
with described average eigenvalue of maximum
According to the described scattering component of object fifty-fifty
with described average eigenvalue of maximum
build by the described scattering component of object fifty-fifty
with described average eigenvalue of maximum
the two-dimensional feature space formed;
In described two-dimensional feature space, if the scattering component of object fifty-fifty of all pixels in the described plot obtained
with the average eigenvalue of maximum of pixels all in described plot
size be all less than pre-set threshold value, then represent that the crops in described plot are gathered in.
As shown from the above technical solution, a kind of method and device of monitoring crops harvesting progress provided by the invention, by the feature of Synthetic Aperture Radar satellite, solving remote optical sensing monitoring harvesting progress data obtains because continuing to stay original place airing for after harvesting, the problem of the accuracy of impact monitoring crops harvesting progress.The method achieve large area, quick and precisely monitor the harvesting progress of crops, ensure that efficient reaping to greatest extent.Reduce the risk of the reasons such as weather to crop, for realizing high crop yield, high-quality, efficiently significant.
In instructions of the present invention, describe a large amount of detail.But can understand, embodiments of the invention can be put into practice when not having these details.In some instances, be not shown specifically known method, structure and technology, so that not fuzzy understanding of this description.
Last it is noted that above each embodiment is only in order to illustrate technical scheme of the present invention, be not intended to limit; Although with reference to foregoing embodiments to invention has been detailed description, those of ordinary skill in the art is to be understood that: it still can be modified to the technical scheme described in foregoing embodiments, or carries out equivalent replacement to wherein some or all of technical characteristic; And these amendments or replacement, do not make the essence of appropriate technical solution depart from the scope of various embodiments of the present invention technical scheme, it all should be encompassed in the middle of the scope of claim of the present invention and instructions.
Accompanying drawing explanation
The schematic flow sheet of the method for the monitoring crops harvesting progress that Fig. 1 provides for one embodiment of the invention;
The monitoring result schematic diagram of the employing method monitoring of the present invention harvesting progress that Fig. 2 provides for the embodiment of the present invention;
The monitoring result schematic diagram of the method monitoring harvesting progress of the employing polarization decomposing that Fig. 3 provides for the embodiment of the present invention;
The structural representation of the device of the monitoring crops harvesting progress that Fig. 4 provides for one embodiment of the invention.
Embodiment
Below in conjunction with accompanying drawing, the embodiment of invention is further described.Following examples only for technical scheme of the present invention is clearly described, and can not limit the scope of the invention with this.
Of the present inventionly a kind ofly monitor the crops harvesting method of progress and device carries out the harvesting progress of monitoring crops based on SAR remote sensing image.
Synthetic Aperture Radar satellite has advantage all-time anf all-weather, has stronger data retrieval capabilities, can overcome optical satellite and often be subject to the shortcoming that the adverse weathers such as cloud, rain, mist cannot obtain data in time.Owing to continuing after crop harvesting in a lot of situation to stay original place airing, cause harvesting plot and not gather in plot difference on remote optical sensing not remarkable.Although remote optical sensing has superiority in detecting light spectrum change, but harvesting plot is be embodied in the change of vegetation structure with not gathering in the larger difference in plot, and the detection that radar remote sensing changes ground object structure is more responsive, therefore radar remote sensing has more advantage in crop monitoring harvesting situation.In addition, relative to optical satellite data, because radar remote sensing data are not by weather effect, can ensure that harvesting progress monitoring and measuring application is to the promptness required by data acquisition and real-time.
Fig. 1 shows a kind of schematic flow sheet of monitoring the method for crops harvesting progress that the embodiment of the present invention provides, and as shown in Figure 1, the method comprises the following steps:
101, obtain the polarimetric synthetic aperture radar remote sensing image of monitored area in the crop harvesting phase, polarization decomposing is carried out to described remote sensing image, obtains the polarization parameter of each pixel.
Concrete, above-mentioned steps 101 also comprises following sub-step:
1011, the polarimetric synthetic aperture radar remote sensing image of monitored area is obtained in the crop harvesting phase;
1012, pre-service is carried out to described remote sensing image;
For example, in specific implementation process, the present embodiment carries out pre-service to remote sensing image and can comprise the following steps:
1012A, radiation calibration is carried out to described SAR remote sensing image;
1012B, the image after radiation calibration carried out to multiple look processing and spot and to make an uproar removal;
1012C, by the polarization scattering matrix of POLARIZATION CHANNEL data transformations of four after denoising;
1012D, geocoding, topographic correction and geometric accurate correction process are carried out to described polarization scattering matrix image.
1013, utilize Polarization target decomposition method to carry out polarization decomposing to described pretreated remote sensing image, obtain the polarization parameter of each pixel.
Utilize Freeman-Durden three-component polarization decomposing method to carry out polarization decomposing to pretreated remote sensing image, obtain the power level P of surface scattering, even scattering and volume scattering three kinds of scattering components
s, P
dand P
v; Three-component divide solve an equation into:
T＝P
s*T
s+P
d*T
d+P
v*T
v
Wherein, T is the polarization scattering matrix represented by coherence matrix form, T
sthe coherence matrix of presentation surface scattering model, P
sthe power representing atural object surface scattering component; T
dthe coherence matrix representing even scattering model, P
dthe power representing atural object even scattering component; T
vthe coherence matrix representing volume scattering model, P
vthe power representing atural object volume scattering component;
Utilize Cloude decomposition method to carry out polarization decomposing to pretreated remote sensing image, obtain the eigenvalue of maximum λ of described polarization scattering matrix T.
102, extract the border in each plot in described remote sensing image, obtain the average polarization parameter of all pixels in described each plot according to the polarization parameter of each pixel described;
In above-mentioned steps 102, extract the border in each plot in described remote sensing image, comprise the following steps:
1021, utilize database or obtain monitored area map according to OO segmentation, sorting technique;
1022, vector quantization is carried out to described monitored area map, obtain the border in each plot in monitored area.
103, according to described average polarization parameter, judge whether the crops in described each plot gather in.
Concrete, above-mentioned steps 103 comprises following sub-step:
1031, the scattering component of object fifty-fifty in described average polarization parameter is obtained
with described average eigenvalue of maximum
1032, according to the described scattering component of object fifty-fifty
with described average eigenvalue of maximum
build by the described scattering component of object fifty-fifty
with described average eigenvalue of maximum
the two-dimensional feature space formed;
1033, in described two-dimensional feature space, if the scattering component of object fifty-fifty of all pixels in the described plot obtained
with the average eigenvalue of maximum of pixels all in described plot
size be all less than pre-set threshold value, then represent that the crops in described plot are gathered in.
Said method solves remote optical sensing monitoring harvesting progress data and obtains because continuing to stay original place airing for after harvesting, the problem of the accuracy of impact monitoring crops harvesting progress.The method achieve large area, quick and precisely monitor the harvesting progress of crops, ensure that efficient reaping to greatest extent.Reduce the risk of the reasons such as weather to crop, for realizing high crop yield, high-quality, efficiently significant.
In order to clearer explanation said method of the present invention, specific embodiment is adopted to be described in detail below.
The embodiment of the present invention proposes a kind of method utilizing SAR remote sensing imaging monitor rape to gather in progress, the harvesting situation in all rape plot on August 27th, 2013, Shang Kuli farm, Ergun City, Inner Mongolia Autonomous Region that utilized the method to monitor, comprises the following steps:
1011, the polarimetric synthetic aperture radar remote sensing image of monitored area is obtained in the crop harvesting phase;
For the harvesting progress in all rape plot between harvest time, Shang Kuli farm, Ergun City, Real-Time Monitoring Inner Mongolia Autonomous Region, obtain the complete polarization Radarsat-2 radar remote sensing image on August 27th, 2013.This image fabric width 25km × 25km, completely covers power farm, storehouse.This scape image obtains with Fine Quad pattern (four polarization fine patterns), distributes with haplopia complex data (SLC) product.During on August 27th, 2013, this farm rape is in harvest time, and part rape cuts (harvested rape slivering is dried in the air on the ground), and part rape is harvesting.
1012, pre-service is carried out to described remote sensing image;
Carry out data prediction to this scape Radarsat-2 image, removal that this preprocessing process comprises radiation calibration, spot is made an uproar, the generation of polarization matrix, geometry correction step, be specially:
1012A, the scaling parameter (Sigma nought) comprised in Radarsat-2 data product file is utilized to carry out radiation calibration to four POLARIZATION CHANNEL data;
1012B, on the basis of 1012A, carry out multiple look processing to the image after radiation calibration, re-use Boxcar wave filter and reduce the intrinsic speckle noise of radar image, filter window is 5 × 5;
1012C, on the basis of 1012B by the radar image of four POLARIZATION CHANNEL transform generate polarization coherence matrix T3 (each pixel is represented by a coherence matrix);
1012D, on the basis of 1012C, geocoding and topographic correction are carried out to coherence matrix T3 image: the dem data utilizing 30 meters of resolution in this region, in conjunction with the geographical location information that Radarsat-2 data product file carries, range-doppler algorithm is utilized to complete geocoding and topographic correction; The Ground Nuclear Magnetic Resonance dominating pair of vertices image of field acquisition is utilized to carry out further geometric accurate correction, by the geometric position precision controlling of every for image pixel within 1 pixel subsequently;
Said process completes under the support of the professional tools such as PolSARPro, ASF MapReady and ENVI;
1013, utilize Polarization target decomposition method to carry out polarization decomposing to described pretreated remote sensing image, obtain the polarization parameter of each pixel.
Utilize Freeman-Durden three-component polarization decomposing method to carry out polarization decomposing to the image after step 1012 process, obtain the power level P of volume scattering component
v; Three-component divide solve an equation into:
T＝P
s*T
s+P
d*T
d+P
v*T
v
Wherein, T is the polarization scattering matrix represented by coherence matrix form, T
sthe coherence matrix of presentation surface scattering model, P
sthe power representing atural object surface scattering component; T
dthe coherence matrix representing even scattering model, P
dthe power representing atural object even scattering component; T
vthe coherence matrix representing volume scattering model, P
vthe power representing atural object volume scattering component;
Recycling Cloude decomposition method carries out polarization decomposing to the image after step 1012 process, obtains the eigenvalue of maximum λ of coherence matrix;
Said process can complete under the support of PolSARPro professional software;
102, the border in each plot in described remote sensing image is extracted;
Obtain the crop-planting zoning figure on farm in this remote sensing image, in conjunction with remote sensing image data, vector quantization is carried out to crop-planting zoning figure, extract the border in each rape plot;
1031, the scattering component of object fifty-fifty in described average polarization parameter is obtained
with described average eigenvalue of maximum
1032, according to the described scattering component of object fifty-fifty
with described average eigenvalue of maximum
build by the described scattering component of object fifty-fifty
with described average eigenvalue of maximum
the two-dimensional feature space formed;
The average body scattering component of all pixels in each rape plot is calculated with plot
and the average eigenvalue of maximum of all pixels in plot
build
with
the two-dimensional feature space formed.
1033, in described two-dimensional feature space, if the scattering component of object fifty-fifty of all pixels in the described plot obtained
with the average eigenvalue of maximum of pixels all in described plot
size be all less than pre-set threshold value, then represent that the crops in described plot are gathered in.
?
with
in the two-dimensional feature space formed, adopt in the present embodiment
and
rape plot for gather in plot, other rape plot is not for gather in plot.
Above-mentioned
with
pre-set threshold value, the present embodiment does not specifically limit, and specifically on the basis of this numerical value, carries out accommodation according to region, described plot, but this value is judging that the degree of accuracy of harvesting progress result obtained produces a desired effect.
In order to verify the result of the present embodiment method, the harvesting situation in these 33 pieces of rape plot, farm field investigation on August 27 in 2013, also obtain high resolving power ZY-102C optical satellite data on August 27th, 2013 (resolution: 5 meters panchromatic+10 meters multispectral), field investigation accurately can identify the harvesting situation in 88 rape plot on the 27th in August in conjunction with high resolution image decipher, wherein 19 pieces of rape plot are not gathered in, and 69 pieces of rape plot are gathered in.Adopt the result of the present embodiment method monitoring as shown in Figure 2, horizontal ordinate represents average eigenvalue of maximum
ordinate represents object scattering component fifty-fifty
can see that 19 pieces are not gathered in the upper right portion that plot is distributed in two feature spaces,
value concentrates on 0.4 ~ 0.8,
value concentrates on 0.18 ~ 0.35; And remainingly gathered in the bottom left section that plot is distributed in two feature spaces,
value concentrates on 0.1 ~ 0.3,
value concentrates on 0.05 ~ 0.15.And predetermined threshold value
relatively, the plot of harvesting all in the present embodiment and do not gather in plot and can accurately be identified.Visible this method can monitor the rape harvesting situation in a certain concrete period preferably.
For verifying the present embodiment method further, also compare with the monitoring method of another kind based on polarization decomposing, the method monitoring parameter have employed the Polarization scattering entropy (Entropy) and two, average Polarization scattering angle (Alpha) parameter that are obtained by Cloude polarization decomposing, its to the monitoring result of the present embodiment identical data as Fig. 3, wherein, horizontal ordinate represents Polarization scattering entropy, ordinate represents average polarization scattering angle, can easily find, part has been gathered in plot and has not been gathered in plot and easily obscured, be difficult to distinguish, visible the method is far inferior to the inventive method for the monitoring result of harvesting situation.This contrast embodies the advantage of the inventive method, and the inventive method continues the situation of stand-down airing to the impact of monitoring effect after considering rape harvesting.
By adopting a kind of method utilizing Synthetic Aperture Radar images monitoring crop to gather in progress disclosed by the invention, achieve large area Real-Time Monitoring crop harvesting progress, increasing work efficiency, while alleviating working strength, effectively can improve the ageing and accuracy of crop harvesting progress monitoring, to contribute to optimizing harvesting strategy, on the whole to regional extent (as county, town, village etc.) harvesting arrangement carry out making overall planning and instructing, carry out more rationally effective machinery, manpower, transport, the Resource allocation and smoothings such as following process, ensure efficient reaping to greatest extent, reduce the reasons such as weather to the risk of crop, for realizing high crop yield, high-quality, efficiently significant.
Fig. 4 shows a kind of device of monitoring crops harvesting progress, and as shown in Figure 4, this device comprises: the first acquisition module 41, second acquisition module 42 and judge module 43.
First acquisition module 41, for obtaining the polarimetric synthetic aperture radar remote sensing image of monitored area in the crop harvesting phase, carries out polarization decomposing to described remote sensing image, obtains the polarization parameter of each pixel.
Concrete above-mentioned first acquisition module, comprising:
Image capturing unit 411, for obtaining the polarimetric synthetic aperture radar remote sensing image of monitored area in the crop harvesting phase;
Pretreatment unit 412, for carrying out pre-service to described remote sensing image;
First parameter acquiring unit 413, for utilizing Polarization target decomposition method to carry out polarization decomposing to described pretreated remote sensing image, obtains the polarization parameter of each pixel.
For example, utilize Freeman-Durden three-component polarization decomposing method to carry out polarization decomposing to pretreated remote sensing image, obtain the power level P of surface scattering, even scattering and volume scattering three kinds of scattering components
s, P
dand P
v; Three-component divide solve an equation into:
T＝P
s*T
s+P
d*T
d+P
v*T
v
Wherein, T is the polarization scattering matrix represented by coherence matrix form, T
sthe coherence matrix of presentation surface scattering model, P
sthe power representing atural object surface scattering component; T
dthe coherence matrix representing even scattering model, P
dthe power representing atural object even scattering component; T
vthe coherence matrix representing volume scattering model, P
vthe power representing atural object volume scattering component;
Utilize Cloude decomposition method to carry out polarization decomposing to pretreated remote sensing image, obtain the eigenvalue of maximum λ of described polarization scattering matrix T.
Second acquisition module 42, for extracting the border in each plot in described remote sensing image, obtains the average polarization parameter of all pixels in described each plot according to the polarization parameter of each pixel described.
Judge module 43, for according to described average polarization parameter, judges whether the crops in described each plot gather in.
For example, described judge module, for:
Obtain the scattering component of object fifty-fifty in described average polarization parameter
with described average eigenvalue of maximum
According to the described scattering component of object fifty-fifty
with described average eigenvalue of maximum
build by the described scattering component of object fifty-fifty
with described average eigenvalue of maximum
the two-dimensional feature space formed;
In described two-dimensional feature space, if the scattering component of object fifty-fifty of all pixels in the described plot obtained
with the average eigenvalue of maximum of pixels all in described plot
size be all less than pre-set threshold value, then represent that the crops in described plot are gathered in.
Apparatus and method of the present invention are one to one, and the computation process because of some parameters in the method is also applicable to the process calculated in this apparatus module, will no longer be described in detail in a device.
Claims (10)
1. monitor a method for crops harvesting progress, it is characterized in that, comprising:
Obtain the polarimetric synthetic aperture radar remote sensing image of monitored area in the crop harvesting phase, polarization decomposing is carried out to described remote sensing image, obtains the polarization parameter of each pixel;
Extract the border in each plot in described remote sensing image, obtain the average polarization parameter of all pixels in described each plot according to the polarization parameter of each pixel described;
According to described average polarization parameter, judge whether the crops in described each plot gather in.
2. method according to claim 1, is characterized in that, the described polarimetric synthetic aperture radar remote sensing image obtaining monitored area in the crop harvesting phase, carries out polarization decomposing, obtain the polarization parameter of each pixel, comprising described remote sensing image:
The polarimetric synthetic aperture radar remote sensing image of monitored area is obtained in the crop harvesting phase;
Pre-service is carried out to described remote sensing image;
Utilize Polarization target decomposition method to carry out polarization decomposing to described pretreated remote sensing image, obtain the polarization parameter of each pixel.
3. method according to claim 2, is characterized in that, describedly carries out pre-service to described remote sensing image, comprising:
Radiation calibration is carried out to described SAR remote sensing image;
Carry out multiple look processing and spot to the image after radiation calibration to make an uproar removal;
By the polarization scattering matrix of POLARIZATION CHANNEL data transformations of four after denoising;
Geocoding, topographic correction and geometric accurate correction process are carried out to described polarization scattering matrix image.
4. method according to claim 3, is characterized in that, the described Polarization target decomposition method that utilizes carries out polarization decomposing to described pretreated remote sensing image, obtains the polarization parameter of each pixel, comprising:
Utilize Freeman-Durden three-component polarization decomposing method to carry out polarization decomposing to pretreated remote sensing image, obtain the power level P of surface scattering, even scattering and volume scattering three kinds of scattering components
s, P
dand P
v; Three-component divide solve an equation into:
T＝P
s*T
s+P
d*T
d+P
v*T
v
Wherein, T is the polarization scattering matrix represented by coherence matrix form, T
sthe coherence matrix of presentation surface scattering model, P
sthe power representing atural object surface scattering component; T
dthe coherence matrix representing even scattering model, P
dthe power representing atural object even scattering component; T
vthe coherence matrix representing volume scattering model, P
vthe power representing atural object volume scattering component;
Utilize Cloude decomposition method to carry out polarization decomposing to pretreated remote sensing image, obtain the eigenvalue of maximum λ of described polarization scattering matrix T.
5. method according to claim 1, is characterized in that, the border in each plot in the described remote sensing image of described extraction, comprising:
Utilize database or obtain monitored area map according to OO segmentation, sorting technique;
Vector quantization is carried out to described monitored area map, obtains the border in each plot in monitored area.
6. method according to claim 5, is characterized in that, described according to described average polarization parameter, judges whether the crops in described each plot gather in, comprising:
Obtain the scattering component of object fifty-fifty in described average polarization parameter
with described average eigenvalue of maximum
;
According to the described scattering component Pv of object fifty-fifty and described average eigenvalue of maximum
build by the described scattering component Pv of object fifty-fifty and described average eigenvalue of maximum
the two-dimensional feature space formed;
In described two-dimensional feature space, if the scattering component of object fifty-fifty of all pixels in the described plot obtained
with the average eigenvalue of maximum of pixels all in described plot
size be all less than pre-set threshold value, then represent that the crops in described plot are gathered in.
7. monitor a device for crops harvesting progress, it is characterized in that, comprising:
First acquisition module, for obtaining the polarimetric synthetic aperture radar remote sensing image of monitored area in the crop harvesting phase, carries out polarization decomposing to described remote sensing image, obtains the polarization parameter of each pixel;
Second acquisition module, for extracting the border in each plot in described remote sensing image, obtains the average polarization parameter of all pixels in described each plot according to the polarization parameter of each pixel described;
Judge module, for according to described average polarization parameter, judges whether the crops in described each plot gather in.
8. device according to claim 7, is characterized in that, described first acquisition module, comprising:
Image capturing unit, for obtaining the polarimetric synthetic aperture radar remote sensing image of monitored area in the crop harvesting phase;
Pretreatment unit, for carrying out pre-service to described remote sensing image;
First parameter acquiring unit, for utilizing Polarization target decomposition method to carry out polarization decomposing to described pretreated remote sensing image, obtains the polarization parameter of each pixel.
9. device according to claim 8, is characterized in that, described first parameter acquiring unit, for:
Utilize Freeman-Durden three-component polarization decomposing method to carry out polarization decomposing to pretreated remote sensing image, obtain the power level P of surface scattering, even scattering and volume scattering three kinds of scattering components
s, P
dand P
v; Three-component divide solve an equation into:
T＝P
s*T
s+P
d*T
d+P
v*T
v
Wherein, T is the polarization scattering matrix represented by coherence matrix form, T
sthe coherence matrix of presentation surface scattering model, P
sthe power representing atural object surface scattering component; T
dthe coherence matrix representing even scattering model, P
dthe power representing atural object even scattering component; T
vthe coherence matrix representing volume scattering model, P
vthe power representing atural object volume scattering component;
Utilize Cloude decomposition method to carry out polarization decomposing to pretreated remote sensing image, obtain the eigenvalue of maximum λ of described polarization scattering matrix T.
10. device according to claim 7, is characterized in that, described judge module, for:
Obtain the scattering component of object fifty-fifty in described average polarization parameter
with described average eigenvalue of maximum
According to the described scattering component of object fifty-fifty
with described average eigenvalue of maximum
build by the described scattering component of object fifty-fifty
with described average eigenvalue of maximum
the two-dimensional feature space formed;
In described two-dimensional feature space, if the scattering component of object fifty-fifty of all pixels in the described plot obtained
with the average eigenvalue of maximum of pixels all in described plot
size be all less than pre-set threshold value, then represent that the crops in described plot are gathered in.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410610674.3A CN104360347B (en) | 2014-11-03 | 2014-11-03 | A kind of monitoring crop gathers in the method and device of progress |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410610674.3A CN104360347B (en) | 2014-11-03 | 2014-11-03 | A kind of monitoring crop gathers in the method and device of progress |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104360347A true CN104360347A (en) | 2015-02-18 |
| CN104360347B CN104360347B (en) | 2017-07-28 |
Family
ID=52527627
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410610674.3A Active CN104360347B (en) | 2014-11-03 | 2014-11-03 | A kind of monitoring crop gathers in the method and device of progress |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104360347B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108508416A (en) * | 2018-03-12 | 2018-09-07 | 西安电子工程研究所 | A kind of polarization reciprocity verification method rapidly and efficiently |
| CN109886142A (en) * | 2019-01-28 | 2019-06-14 | 中科光启空间信息技术有限公司 | A kind of crops decomposition method based on SAR technology |
| CN115372970A (en) * | 2022-08-19 | 2022-11-22 | 陕西省土地工程建设集团有限责任公司 | Remote sensing extraction method for crops SAR in mountainous and hilly areas |
| CN115578641A (en) * | 2022-11-08 | 2023-01-06 | 中化现代农业有限公司 | Crop shoveling progress monitoring method and device, electronic equipment and storage medium |
| WO2024019632A1 (en) * | 2022-07-22 | 2024-01-25 | Публичное Акционерное Общество "Сбербанк России" | Device and method for determining crop productivity |
| WO2024085780A1 (en) * | 2022-10-17 | 2024-04-25 | Публичное Акционерное Общество "Сбербанк России" | Device and method for identifying crop types |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103969632A (en) * | 2014-03-26 | 2014-08-06 | 北京农业信息技术研究中心 | Device and method of using radar remote sensing data for monitoring wheat lodging |
-
2014
- 2014-11-03 CN CN201410610674.3A patent/CN104360347B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103969632A (en) * | 2014-03-26 | 2014-08-06 | 北京农业信息技术研究中心 | Device and method of using radar remote sensing data for monitoring wheat lodging |
Non-Patent Citations (5)
| Title |
|---|
| H. MCNAIRN等: "The application of C-band polarimetric SAR for agriculture: a review", 《CANADIAN JOURNAL OF REMOTE SENSING》 * |
| 任俊英等: "基于中间层特性的全极化SAR监督地物分类", 《遥感技术与应用》 * |
| 化国强等: "基于全极化SAR数据散射机理的农作物分类", 《江苏农业学报》 * |
| 张淼等: "作物残茬覆盖度遥感监测研究进展", 《光谱学与光谱分析》 * |
| 郭华东等: "航空双波段全极化SAR信息分析", 《遥感学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108508416A (en) * | 2018-03-12 | 2018-09-07 | 西安电子工程研究所 | A kind of polarization reciprocity verification method rapidly and efficiently |
| CN109886142A (en) * | 2019-01-28 | 2019-06-14 | 中科光启空间信息技术有限公司 | A kind of crops decomposition method based on SAR technology |
| CN109886142B (en) * | 2019-01-28 | 2022-12-02 | 中科光启空间信息技术有限公司 | Crop interpretation method based on SAR technology |
| WO2024019632A1 (en) * | 2022-07-22 | 2024-01-25 | Публичное Акционерное Общество "Сбербанк России" | Device and method for determining crop productivity |
| CN115372970A (en) * | 2022-08-19 | 2022-11-22 | 陕西省土地工程建设集团有限责任公司 | Remote sensing extraction method for crops SAR in mountainous and hilly areas |
| WO2024085780A1 (en) * | 2022-10-17 | 2024-04-25 | Публичное Акционерное Общество "Сбербанк России" | Device and method for identifying crop types |
| CN115578641A (en) * | 2022-11-08 | 2023-01-06 | 中化现代农业有限公司 | Crop shoveling progress monitoring method and device, electronic equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104360347B (en) | 2017-07-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104360347A (en) | Method and device for monitoring crop harvesting progress | |
| Dong et al. | Estimating winter wheat biomass by assimilating leaf area index derived from fusion of Landsat-8 and MODIS data | |
| CN104656098B (en) | A kind of method of remote sensing forest biomass inverting | |
| CN111738165B (en) | Method for extracting individual plant canopy from high-resolution unmanned aerial vehicle visible light remote sensing image | |
| CN104463164B (en) | It is a kind of based on umbrella frame method and crown height than trees canopy structure information extracting method | |
| Baofeng et al. | Digital surface model applied to unmanned aerial vehicle based photogrammetry to assess potential biotic or abiotic effects on grapevine canopies | |
| González-Yebra et al. | Methodological proposal to assess plastic greenhouses land cover change from the combination of archival aerial orthoimages and Landsat data | |
| KR102366267B1 (en) | Analysis apparatus for present condition of vegetation, and control method thereof | |
| Zhao et al. | UAV-based individual shrub aboveground biomass estimation calibrated against terrestrial LiDAR in a shrub-encroached grassland | |
| CN110310246B (en) | Sugarcane planting area remote sensing information extraction method based on three-linear array image | |
| Millan et al. | Vineyard pruning weight assessment by machine vision: towards an on-the-go measurement system: This article is published in cooperation with the 21th GIESCO International Meeting, June 23-28 2019, Thessaloniki, Greece. Guests editors: Stefanos Koundouras and Laurent Torregrosa | |
| Panagiotidis et al. | Detection of fallen logs from high-resolution UAV images | |
| CN112434569A (en) | Thermal imaging system of unmanned aerial vehicle | |
| Wang et al. | A new object-class based gap-filling method for PlanetScope satellite image time series | |
| CN101788685B (en) | Remote sensing earthquake damage information extracting and digging method based on pixels | |
| Itakura et al. | Estimation of tree structure parameters from video frames with removal of blurred images using machine learning | |
| Xiao et al. | A Novel Image Fusion Method for Water Body Extraction Based on Optimal Band Combination. | |
| Demirkesen et al. | Hyperspectral images and lidar based DEM fusion: A multi-modal landuse classification strategy | |
| CN104330798A (en) | Synthetic aperture radar based crop seeding date monitoring method and device through remote sensing image | |
| Marcal et al. | Evaluation of satellite image segmentation using synthetic images | |
| Van Coillie et al. | Optimized feature fusion of LiDAR and hyperspectral data for tree species mapping in closed forest canopies | |
| Johansen et al. | Multi-temporal assessment of lychee tree crop structure using multi-spectral RPAS imagery | |
| Ichikawa et al. | Identification of paddy fields in Northern Japan using RapidEye images | |
| Angelini et al. | A Review and Test of Shoreline Extraction Techniques | |
| Ali-Sisto et al. | Comparison of 3D point clouds from aerial stereo images and lidar for forest change detection |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20210510 Address after: 100097 room 301-a315, 03 / F, 11 Shuguang Huayuan Middle Road, Haidian District, Beijing Patentee after: NONGXIN TECHNOLOGY (BEIJING) Co.,Ltd. Address before: 100097 room 206, block a, agricultural science building, No. 11, Shuguang Huayuan Middle Road, Haidian District, Beijing Patentee before: BEIJING RESEARCH CENTER FOR INFORMATION TECHNOLOGY IN AGRICULTURE |