CN107369175A - A kind of oyster for aquaculture arranges Area computing method and system - Google Patents
A kind of oyster for aquaculture arranges Area computing method and system Download PDFInfo
- Publication number
- CN107369175A CN107369175A CN201710458345.5A CN201710458345A CN107369175A CN 107369175 A CN107369175 A CN 107369175A CN 201710458345 A CN201710458345 A CN 201710458345A CN 107369175 A CN107369175 A CN 107369175A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- image
- mtd
- remote sensing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 241000237502 Ostreidae Species 0.000 title claims abstract description 67
- 235000020636 oyster Nutrition 0.000 title claims abstract description 67
- 238000009360 aquaculture Methods 0.000 title claims abstract description 21
- 244000144974 aquaculture Species 0.000 title claims abstract description 21
- 238000004364 calculation method Methods 0.000 title abstract description 21
- 238000000034 method Methods 0.000 claims abstract description 34
- 238000012545 processing Methods 0.000 claims abstract description 31
- 238000012937 correction Methods 0.000 claims abstract description 20
- 238000009499 grossing Methods 0.000 claims description 20
- 238000012952 Resampling Methods 0.000 claims description 10
- 230000009467 reduction Effects 0.000 claims description 10
- 230000014509 gene expression Effects 0.000 claims description 7
- 238000003384 imaging method Methods 0.000 claims description 5
- 238000010835 comparative analysis Methods 0.000 claims description 4
- 239000003550 marker Substances 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 238000013139 quantization Methods 0.000 abstract 1
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 9
- 238000012364 cultivation method Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 238000005259 measurement Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 230000000007 visual effect Effects 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000009395 breeding Methods 0.000 description 3
- 230000001488 breeding effect Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000012136 culture method Methods 0.000 description 2
- 238000000556 factor analysis Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 241000894007 species Species 0.000 description 2
- 235000002566 Capsicum Nutrition 0.000 description 1
- 241000548230 Crassostrea angulata Species 0.000 description 1
- 241001522252 Crassostrea rivularis Species 0.000 description 1
- 241000237852 Mollusca Species 0.000 description 1
- 239000006002 Pepper Substances 0.000 description 1
- 235000016761 Piper aduncum Nutrition 0.000 description 1
- 235000017804 Piper guineense Nutrition 0.000 description 1
- 244000203593 Piper nigrum Species 0.000 description 1
- 235000008184 Piper nigrum Nutrition 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 235000015170 shellfish Nutrition 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000004304 visual acuity Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10036—Multispectral image; Hyperspectral image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20182—Noise reduction or smoothing in the temporal domain; Spatio-temporal filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30188—Vegetation; Agriculture
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Processing (AREA)
Abstract
The application is related to a kind of oyster for aquaculture and arranges Area computing method and system, the measuring method obtains view data corresponding to oyster row's area by way of multispectral remote sensing, described image data are carried out with geometric correction, correction, then image enhancement processing is carried out, with the prominent useful information related from oyster row's area and expand the different direct characteristic differences of atural object, the image that grid is obtained to processing carries out vector quantization, chooses the related feature of oyster row's area using the spatial database being connected to calculate oyster row's area.By such a mode, the application is handled the view data of remote sensing images, it is possible to achieve automatic interpretation and calculating, is reduced cost of labor, and can be improved computational accuracy by geometric correction and correction.
Description
Technical Field
The application relates to the technical field of measurement and calculation, in particular to an oyster arrangement area measurement and calculation method and system for aquaculture.
Background
Oyster, named oyster by its scholar name, belongs to the phylum mollusca, the family of oysters. Is the most common bivalve shellfish along the sea and has various varieties. More than 200 species of oysters exist in the world, and China has species of crassostrea gigas, ostrea rivularis and the like. The offshore oyster cultivation method mainly comprises a stone shell attachment cultivation method, a pile cultivation method, a bridge cultivation method, a felling cultivation method, a pile-connected cultivation method and the like. Along with the development of oyster cultivation cardinality and the improvement of social environment protection consciousness, the fell cultivation method is more and more popular. The raft culture method is one of oyster three-dimensional culture methods, is of a multi-purpose longline type, and the cultured oysters grow fast and have high yield per unit area. The most of the cultivated oysters are the seeds collected by a shell seedling collector, the shells are connected in series by ropes, the length of the shells is determined according to the depth of water, and the shells are hung under a raft.
The oyster cultivation area in China is increased from 112818 hectares in 2009 to 141498 hectares in 2015, and the oyster industry yield is increased from 350.38 million tons in 2009 to 457.34 million tons in 2015. Because the oyster cultivation scale is continuously enlarged, the offshore excessive cultivation and the cultivation environment are deteriorated, the coastal oyster cultivation area industry needs to be intensively monitored and managed, and how to detect that the coastal oyster cultivation area in each area is more important. The traditional monitoring means can not meet the real requirement, and has great limitation and risk, so that a new monitoring means and a new monitoring method are required.
Remote sensing has the advantages of wide detection range, high data acquisition speed, short period, strong timeliness, low cost, great economic benefit and the like. The remote sensing image can be used for quickly extracting the area of the needed aquaculture, helping the farm to select sites, determining breeding varieties and monitoring the breeding density and the pollution (red tide, water quality and the like) of the breeding water body; by combining with GIS (Geographic Information system) technology, the method can also plan and manage the culture area and evaluate the influence of the aquaculture area on the environment.
There is no method for identifying the most common and optimal aquaculture area, limited by objective factors such as study time, study area and data source. The current commonly used aquaculture area identification methods mainly include visual interpretation, information extraction based on ratio index analysis, information extraction based on neighborhood analysis, information extraction based on corresponding analysis, and the like.
The visual interpretation is to combine with various non-remote sensing information data according to the visual interpretation mark and interpretation experience of the remote sensing image, and to use the relevant knowledge to carry out comprehensive analysis and logical reasoning, so as to obtain the required thematic information from the remote sensing image.
The ratio index is in the multispectral wave band of the same image, the ratio of the brightness value of each pixel in different wave bands is obtained, a new image is formed, the influence of certain factors or backgrounds causing illumination difference is reduced, and the radiation characteristic of a target ground object is highlighted.
The corresponding analysis is an analysis method developed on the basis of factor analysis, also called R-Q type factor analysis.
Neighborhood analysis is a method for spatially analyzing each pixel of a wave band according to the neighboring pixels around, and the number and position of the analyzed and operated pixels are determined by a scanning window.
However, oyster cultivation land has very similar spectral characteristics with other water body types, and conventional multispectral remote sensing can only provide discontinuous spectral band information with spectral resolution of more than 100nm (nanometers). In the oyster cultivation area identification method, the problems of foreign matter co-spectrum cannot be well solved by using visual interpretation, ratio index analysis and corresponding analysis methods, and pepper salt noise can be generated in classification results.
Therefore, there is a need in the art for a new technique to solve the above-mentioned problems.
Disclosure of Invention
Therefore, it is necessary to provide an oyster rank area calculation method and system for aquaculture, which can reduce labor cost and improve calculation accuracy.
An oyster row area measuring and calculating method for aquaculture comprises the following steps:
acquiring image data corresponding to the oyster arrangement area in a multispectral remote sensing mode;
performing geometric correction on the image data, and correcting the deformed image data;
carrying out image enhancement processing on the corrected and rectified image data so as to highlight useful information related to oyster arrangement areas and enlarge direct characteristic differences of different objects;
and processing the image data subjected to image enhancement processing to obtain a grid-form image, carrying out vectorization on the image, and selecting features related to the oyster arrangement area by using a connected spatial database to calculate the oyster arrangement area.
An oyster placement area measurement and calculation system for aquaculture, the measurement and calculation system comprising a processor, the processor when executing program instructions:
the processor is used for acquiring image data corresponding to the oyster arrangement area in a multispectral remote sensing mode;
the processor is used for carrying out geometric correction on the image data and rectifying the deviation of the deformed image data;
the processor is used for carrying out image enhancement processing on the corrected and rectified image data so as to highlight useful information related to oyster arrangement areas and enlarge direct characteristic differences of different objects;
the processor is used for processing the image data subjected to image enhancement processing to obtain a grid-form image, vectorizing the image, and selecting features related to the oyster arrangement area by using a connected spatial database to calculate the oyster arrangement area.
According to the oyster rank area measuring and calculating method and system for aquaculture, image data corresponding to the oyster rank area is obtained through a multispectral remote sensing mode, geometric correction and deviation rectification are carried out on the image data, image enhancement processing is carried out, useful information related to the oyster rank area is highlighted, direct characteristic differences of different objects are enlarged, vectorization is carried out on the processed image in a grid mode, and features related to the oyster rank area are selected through a connected spatial database to measure and calculate the oyster rank area. By the method, the image data of the remote sensing image is processed, automatic interpretation and calculation can be achieved, labor cost is reduced, and calculation accuracy can be improved through geometric correction and deviation rectification.
Drawings
FIG. 1 is a flowchart illustrating an oyster rank area estimation method for aquaculture according to an embodiment;
FIG. 2 is a diagram illustrating a gray scale value weighted interpolation of pixels during resampling according to an embodiment;
FIG. 3 is a schematic diagram of a coordinate system used in the linear stretch calculation of the present application;
fig. 4 is a block diagram illustrating an oyster rank area estimation system for aquaculture according to an embodiment of the present disclosure;
FIG. 5 is a general flow chart of the oyster rank area calculation system for aquaculture according to the present application
Detailed Description
In an embodiment, referring to fig. 1, the oyster rank area calculation method for aquaculture of this embodiment includes, but is not limited to, the following steps.
And S101, acquiring image data corresponding to the oyster arrangement area in a multispectral remote sensing mode.
In this embodiment S101, multispectral resolution remote sensing may be adopted, and two or more sensors of spectrum channels are used to perform synchronous imaging on the surface feature including the oyster row, so as to divide electromagnetic wave information of radiation reflected by the surface feature into a plurality of spectrum segments for receiving and recording.
S102, performing geometric correction on the image data, and correcting the deformed image data.
In this embodiment S102, the performing geometric correction on the image data may include: and performing geometric correction by selecting ground control points and image resampling.
It should be noted that the geometric correction by selecting the ground control point and resampling the image may include: taking the vectorized topographic map as a base map, carrying out comparative analysis on a remote sensing image of the image data acquired by multispectral remote sensing, and selecting an obvious and clear positioning marker on the base map and the remote sensing image as a ground control point; and resampling the remote sensing image according to the position of each pixel on the output image in the input image, and establishing a new remote sensing image.
It should be noted that, the remote sensing image is resampled according to the position of each pixel on the output image in the input image to create a new remote sensing image, and specifically, a cubic interpolation method may be adopted, and a cubic polynomial is used to solve a manner of approximating a theoretically optimal interpolation function.
Specifically, with a cubic polynomial s (x), an optimal interpolation function sin (x)/x, and the cubic interpolation method of the present embodiment, the expressions include, but are not limited to, the following:
it can be seen that the gray value of the pixel (x, y) to be determined is obtained by weighted interpolation of 16 gray values around the pixel, as shown in fig. 2.
Then the gray scale of the pixel to be solved is calculated as follows:
f(x,y)=f(i+u,j+v)=ABC
wherein:
s103, carrying out image enhancement processing on the corrected and rectified image data so as to highlight useful information related to oyster arrangement areas and enlarge direct characteristic differences of different objects.
In S103, it should be noted that there are many noises in the image data of the remote sensing image, and these noises not only limit the radiation resolution of the remote sensing image and affect the resolving power for different signal intensities, but also reduce the recognition capability for ground targets and structures, so in this embodiment, the image data after being corrected and rectified is subjected to image enhancement, and specifically, a Frost filter may be used to perform noise reduction on the image data.
It should be noted that the impulse response of the Frost filter employed in this embodiment is a bilateral exponential function, which is approximated to a low-pass filter, the filter parameters of which are determined by the image local variance coefficient, and the attenuation speed of the impulse response depends on the size of the local variance coefficient and is in direct proportion thereto. The front filter is a circular symmetric filter with a weight M value as an adaptive adjustment parameter, and the mathematical expression of the front filter in this embodiment specifically includes the following:
wherein,
in the above mathematical expression of this embodiment, g'(i,j)The pixel gray value is the pixel gray value after smoothing treatment; g(i,j)The original gray value of each pixel in the smooth window;the average value of the pixel gray levels in the window is obtained; m(i,j)The weight index of each corresponding pixel in the smoothing window; t is(i,j)The absolute distance from the central pixel to the adjacent pixels in the smoothing window is calculated; sigma(i,j)Is the variance of the pixel values in the smoothing window;is the size of the smoothing window; l is the imaging view.
Furthermore, after performing noise reduction processing on image data by using a Frost filter, the embodiment may further include: and the dynamic range of the gray value of the remote sensing image is stretched to a specified range according to a linear relation formula in a linear stretching mode, so that the contrast of the remote sensing image is increased. In the embodiment, original pixels with different brightness values in the image are converted into the same brightness value through linear stretching, and the original similar brightness values become dissimilar, so that the contrast of the whole image is enhanced, and the brightness of the image is improved.
It should be noted that, as shown in fig. 3, in the schematic diagram of the coordinate system shown in fig. 3, each grid represents a coordinate point, such as an intermediate pixel unit, whose coordinate is (i, j), the transformation formula of the linear stretching includes:
wherein, G'(i,j)Is gray value with pixel coordinate of (i, j) after linear stretching treatment, G(i,j)Is the gray value with the coordinate of (i, j) of the pixel after the noise reduction processing, C0Is a constant value of C0May be determined from the overall actual image.
And S104, processing the image data subjected to image enhancement to obtain a grid-form image, vectorizing the image, and selecting features related to the oyster arrangement area by using a connected spatial database to calculate the oyster arrangement area.
In this embodiment S104, a professional software may be specifically used to perform vectorization on the processed image in the grid form, and the scope understood by those skilled in the art is not limited thereto.
According to the embodiment of the oyster row area measuring and calculating method for aquaculture, automatic interpretation and calculation can be achieved by processing the image data of the remote sensing image, labor cost is reduced, and calculation accuracy can be improved through geometric correction and deviation correction.
In an embodiment, referring to fig. 3, the oyster rank area measuring and calculating system for aquaculture of the present embodiment includes a processor 30, and when the processor 30 executes a program instruction, the following processes may be implemented, but not limited to. The measuring and calculating system of the present embodiment may further include a memory, which may be an external memory connected to a network, or a memory installed in a wired connection, and is not limited herein.
Specifically, the processor 30 is configured to obtain image data corresponding to the oyster arrangement area in a multispectral remote sensing manner;
the processor 30 is configured to perform geometric correction on the image data and perform deviation rectification on the deformed image data;
the processor 30 is used for performing image enhancement processing on the corrected and rectified image data so as to highlight useful information related to oyster arrangement areas and enlarge direct characteristic differences of different objects;
the processor 30 is configured to process the image data subjected to the image enhancement processing to obtain an image in a grid form, perform vectorization on the image, and select features related to the oyster allocation area by using a connected spatial database to calculate the oyster allocation area.
In a specific implementation process, the processor 30 is specifically configured to perform comparative analysis on a remote sensing image of the image data obtained by multispectral remote sensing with the vectorized topographic map as a base map, and select an obvious and clear positioning marker as a ground control point on the base map and the remote sensing image; the processor 30 is specifically configured to resample the remote sensing image by using a cubic interpolation method according to the position of each pixel on the output image in the input image, and establish a new remote sensing image.
For a specific process of interpolating three times in this embodiment, please refer to the related description of the previous method embodiment, which is not described herein again.
The processor 30 is specifically configured to perform noise reduction processing on the image data by using a Frost filter; the processor 30 is specifically configured to stretch the dynamic range of the gray value of the remote sensing image to a specified range according to a linear relation formula in a linear stretching manner, so as to increase the contrast of the remote sensing image.
The mathematical expression of the Frost filter in this embodiment specifically includes the following:
wherein,
in the above mathematical expression of this embodiment, g'(i,j)The pixel gray value is the pixel gray value after smoothing treatment; g(i,j)The original gray value of each pixel in the smooth window;the average value of the pixel gray levels in the window is obtained; m(i,j)The weight index of each corresponding pixel in the smoothing window; t is(i,j)The absolute distance from the central pixel to the adjacent pixels in the smoothing window is calculated; sigma(i,j)Is the variance of the pixel values in the smoothing window;is the size of the smoothing window; l is the imaging view.
As mentioned above, the transformation formula of the linear stretching includes:
wherein, G'(i,j)Is gray value with pixel coordinate of (i, j) after linear stretching treatment, G(i,j)For pixel seat after noise reduction treatmentGray value denoted (i, j), C0Is a constant value of C0May be determined from the overall actual image.
Please refer to fig. 5, fig. 5 is a general flowchart of the oyster rank area measurement and calculation system for aquaculture according to the present invention, where the oyster rank area measurement and calculation system of the present embodiment includes the processes of ground point selection, resampling, image denoising, linear stretching, rasterization, area calculation, etc. as described above, the calculation processing method related to any of the above embodiments may be specifically adopted, and is not described herein again.
According to the embodiment, automatic interpretation and calculation can be realized by processing the image data of the remote sensing image, the labor cost is reduced, and the calculation precision can be improved by geometric correction and deviation rectification.
The above description is only an example of the present application and is not intended to limit the scope of the present application, and all modifications of equivalent structures and equivalent processes, which are made by the contents of the specification and the drawings, or which are directly or indirectly applied to other related technical fields, are intended to be included within the scope of the present application.
Claims (10)
1. An oyster row area measuring and calculating method for aquaculture is characterized by comprising the following steps:
acquiring image data corresponding to the oyster arrangement area in a multispectral remote sensing mode;
performing geometric correction on the image data, and correcting the deformed image data;
carrying out image enhancement processing on the corrected and rectified image data so as to highlight useful information related to oyster arrangement areas and enlarge direct characteristic differences of different objects;
and processing the image data subjected to image enhancement processing to obtain a grid-form image, carrying out vectorization on the image, and selecting features related to the oyster arrangement area by using a connected spatial database to calculate the oyster arrangement area.
2. The method of claim 1, wherein the geometrically correcting the image data comprises: and performing geometric correction by selecting ground control points and image resampling.
3. The method for estimation according to claim 2, wherein the geometric correction by means of selection of ground control points and image resampling comprises:
taking the vectorized topographic map as a base map, carrying out comparative analysis on a remote sensing image of the image data acquired by multispectral remote sensing, and selecting an obvious and clear positioning marker on the base map and the remote sensing image as a ground control point;
and resampling the remote sensing image according to the position of each pixel on the output image in the input image, and establishing a new remote sensing image.
4. The method for measuring and calculating according to claim 3, wherein the resampling the remote sensing image according to the positions of the image elements on the output image in the input image to establish a new remote sensing image comprises: and processing by adopting a cubic interpolation method and utilizing a cubic polynomial to solve a mode of approximating to a theoretically optimal interpolation function.
5. The method according to claim 3 or 4, wherein the image enhancement processing on the corrected and rectified image data specifically comprises:
performing noise reduction processing on image data by using a Frost filter, wherein a mathematical expression of the Frost filter comprises the following steps:
<mrow> <msubsup> <mi>g</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>&prime;</mo> </msubsup> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mover> <msub> <mi>g</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>&OverBar;</mo> </mover> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msub> <mi>C</mi> <mi>I</mi> </msub> <mo><</mo> <msub> <mi>C</mi> <mi>u</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>g</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>&times;</mo> <msub> <mi>M</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> </mrow> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>M</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> </mrow> </mfrac> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msub> <mi>C</mi> <mi>u</mi> </msub> <mo>&le;</mo> <msub> <mi>C</mi> <mi>I</mi> </msub> <mo>&le;</mo> <msub> <mi>C</mi> <mi>max</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>g</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msub> <mi>C</mi> <mi>I</mi> </msub> <mo>></mo> <msub> <mi>C</mi> <mi>max</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
wherein,
<mrow> <msub> <mi>A</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mfrac> <msub> <mi>&sigma;</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <msup> <mover> <msub> <mi>g</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>&OverBar;</mo> </mover> <mn>2</mn> </msup> </mfrac> </mrow>
g′(i,j)the pixel gray value is the pixel gray value after smoothing treatment; g(i,j)The original gray value of each pixel in the smooth window;the average value of the pixel gray levels in the window is obtained; m(i,j)The weight index of each corresponding pixel in the smoothing window; t is(i,j)The absolute distance from the central pixel to the adjacent pixels in the smoothing window is calculated; sigma(i,j)For smoothing windowsVariance of pixel values in the mouth;is the size of the smoothing window;l is the imaging view.
6. The method of claim 5, wherein after performing noise reduction processing on the image data by using a Frost filter, the method further comprises:
and the dynamic range of the gray value of the remote sensing image is stretched to a specified range according to a linear relation formula in a linear stretching mode, so that the contrast of the remote sensing image is increased.
7. The method of claim 6, wherein the transformation formula of the linear stretching comprises:
<mrow> <msub> <msup> <mi>G</mi> <mo>,</mo> </msup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mfrac> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mn>2</mn> </mfrac> <mo>+</mo> <mfrac> <mrow> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>+</mo> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>+</mo> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>+</mo> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>+</mo> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>+</mo> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>+</mo> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>+</mo> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </mrow> <mn>16</mn> </mfrac> <mo>+</mo> <msub> <mi>C</mi> <mn>0</mn> </msub> </mrow>
wherein, G'(i,j)Is gray value with pixel coordinate of (i, j) after linear stretching treatment, G(i,j)Is the gray value with the coordinate of (i, j) of the pixel after the noise reduction processing, C0Is a constant.
8. The oyster arrangement area measuring and calculating system for aquaculture is characterized by comprising a processor, wherein when the processor executes program instructions:
the processor is used for acquiring image data corresponding to the oyster arrangement area in a multispectral remote sensing mode;
the processor is used for carrying out geometric correction on the image data and rectifying the deviation of the deformed image data;
the processor is used for carrying out image enhancement processing on the corrected and rectified image data so as to highlight useful information related to oyster arrangement areas and enlarge direct characteristic differences of different objects;
the processor is used for processing the image data subjected to image enhancement processing to obtain a grid-form image, vectorizing the image, and selecting features related to the oyster arrangement area by using a connected spatial database to calculate the oyster arrangement area.
9. The system of claim 8, wherein:
the processor is specifically used for carrying out comparative analysis on the remote sensing image of the image data acquired by multispectral remote sensing by taking the vectorized topographic map as a base map, and selecting an obvious and clear positioning marker as a ground control point on the base map and the remote sensing image; the processor is specifically used for resampling the remote sensing image by adopting a cubic interpolation method according to the position of each pixel on the output image in the input image, and establishing a new remote sensing image;
the processor is specifically configured to perform noise reduction processing on the image data by using a Frost filter; the processor is specifically used for stretching the dynamic range of the gray value of the remote sensing image to a specified range according to a linear relation formula in a linear stretching mode, and increasing the contrast of the remote sensing image;
wherein, the mathematical expression of the Frost filter comprises the following:
<mrow> <msubsup> <mi>g</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>&prime;</mo> </msubsup> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mover> <msub> <mi>g</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>&OverBar;</mo> </mover> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msub> <mi>C</mi> <mi>I</mi> </msub> <mo><</mo> <msub> <mi>C</mi> <mi>u</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>g</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>&times;</mo> <msub> <mi>M</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> </mrow> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>M</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> </mrow> </mfrac> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msub> <mi>C</mi> <mi>u</mi> </msub> <mo>&le;</mo> <msub> <mi>C</mi> <mi>I</mi> </msub> <mo>&le;</mo> <msub> <mi>C</mi> <mi>max</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>g</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msub> <mi>C</mi> <mi>I</mi> </msub> <mo>></mo> <msub> <mi>C</mi> <mi>max</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
wherein,
<mrow> <msub> <mi>A</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mfrac> <msub> <mi>&sigma;</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <msup> <mover> <msub> <mi>g</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>&OverBar;</mo> </mover> <mn>2</mn> </msup> </mfrac> </mrow>
g′(i,j)the pixel gray value is the pixel gray value after smoothing treatment; g(i,j)The original gray value of each pixel in the smooth window;the average value of the pixel gray levels in the window is obtained; m(i,j)The weight index of each corresponding pixel in the smoothing window; t is(i,j)The absolute distance from the central pixel to the adjacent pixels in the smoothing window is calculated; sigma(i,j)Is the variance of the pixel values in the smoothing window;is the size of the smoothing window;l is the imaging view.
10. The system of claim 9, wherein the linear stretch transformation equation comprises:
<mrow> <msub> <msup> <mi>G</mi> <mo>,</mo> </msup> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mfrac> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mn>2</mn> </mfrac> <mo>+</mo> <mfrac> <mrow> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>+</mo> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>+</mo> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>+</mo> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>+</mo> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </msub> <mo>+</mo> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>+</mo> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>+</mo> <msub> <mi>G</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </mrow> <mn>16</mn> </mfrac> <mo>+</mo> <msub> <mi>C</mi> <mn>0</mn> </msub> </mrow>
wherein, G'(i,j)Is gray value with pixel coordinate of (i, j) after linear stretching treatment, G(i,j)Is the gray value with the coordinate of (i, j) of the pixel after the noise reduction processing, C0Is a constant.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710458345.5A CN107369175A (en) | 2017-06-16 | 2017-06-16 | A kind of oyster for aquaculture arranges Area computing method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710458345.5A CN107369175A (en) | 2017-06-16 | 2017-06-16 | A kind of oyster for aquaculture arranges Area computing method and system |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107369175A true CN107369175A (en) | 2017-11-21 |
Family
ID=60305364
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710458345.5A Pending CN107369175A (en) | 2017-06-16 | 2017-06-16 | A kind of oyster for aquaculture arranges Area computing method and system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107369175A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109063605A (en) * | 2018-07-16 | 2018-12-21 | 中国水产科学研究院南海水产研究所 | A kind of oyster row's oyster frame cultivation measuring method, system and device for aquaculture |
CN113160262A (en) * | 2021-03-23 | 2021-07-23 | 珠海欧比特宇航科技股份有限公司 | Oyster row extraction method, system and medium based on high-resolution satellite remote sensing image |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103363962A (en) * | 2013-07-23 | 2013-10-23 | 中国科学院遥感与数字地球研究所 | Remote sensing evaluation method of lake water reserves based on multispectral images |
CN106791320A (en) * | 2016-12-30 | 2017-05-31 | 北京农业信息技术研究中心 | A kind of cultured fishes image intensifier device and method |
CN106815559A (en) * | 2016-12-21 | 2017-06-09 | 中国科学院深圳先进技术研究院 | A kind of utilization SAR data monitoring oyster arranges method and device, the user equipment in region |
-
2017
- 2017-06-16 CN CN201710458345.5A patent/CN107369175A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103363962A (en) * | 2013-07-23 | 2013-10-23 | 中国科学院遥感与数字地球研究所 | Remote sensing evaluation method of lake water reserves based on multispectral images |
CN106815559A (en) * | 2016-12-21 | 2017-06-09 | 中国科学院深圳先进技术研究院 | A kind of utilization SAR data monitoring oyster arranges method and device, the user equipment in region |
CN106791320A (en) * | 2016-12-30 | 2017-05-31 | 北京农业信息技术研究中心 | A kind of cultured fishes image intensifier device and method |
Non-Patent Citations (2)
Title |
---|
刘德长: "《后遥感应用技术研究》", 31 May 2007 * |
李天华 等: "《"3S"技术与高原机场建设工程》", 31 May 2011, 电子科技大学出版社 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109063605A (en) * | 2018-07-16 | 2018-12-21 | 中国水产科学研究院南海水产研究所 | A kind of oyster row's oyster frame cultivation measuring method, system and device for aquaculture |
CN113160262A (en) * | 2021-03-23 | 2021-07-23 | 珠海欧比特宇航科技股份有限公司 | Oyster row extraction method, system and medium based on high-resolution satellite remote sensing image |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Malambo et al. | Automated detection and measurement of individual sorghum panicles using density-based clustering of terrestrial lidar data | |
Roelfsema et al. | Multi-temporal mapping of seagrass cover, species and biomass: A semi-automated object based image analysis approach | |
CN103914678B (en) | Abandoned land remote sensing recognition method based on texture and vegetation indexes | |
CN102800074B (en) | Synthetic aperture radar (SAR) image change detection difference chart generation method based on contourlet transform | |
CN107423537B (en) | Surface temperature downscaling method based on self-adaptive threshold | |
CN109146889A (en) | A kind of field boundary extracting method based on high-resolution remote sensing image | |
CN107247927B (en) | Method and system for extracting coastline information of remote sensing image based on tassel cap transformation | |
CN103226832B (en) | Based on the multi-spectrum remote sensing image change detecting method of spectral reflectivity mutation analysis | |
CN108918820B (en) | Method and device for acquiring graded distribution of salinization degree of farmland soil | |
Delenne et al. | From pixel to vine parcel: A complete methodology for vineyard delineation and characterization using remote-sensing data | |
CN114740180A (en) | Soil organic carbon estimation method and device based on multi-source remote sensing data | |
CN116188465B (en) | Crop growth state detection method based on image processing technology | |
CN115861409B (en) | Soybean leaf area measuring and calculating method, system, computer equipment and storage medium | |
CN104952070A (en) | Near-rectangle guide based remote-sensing cornfield image segmentation method | |
CN106156758A (en) | A kind of tidal saltmarsh method in SAR coast chart picture | |
Kamal et al. | Comparison of Google Earth Engine (GEE)-based machine learning classifiers for mangrove mapping | |
CN117975282B (en) | Remote sensing extraction method and system for laver cultivation area based on multi-element optical feature fusion | |
CN107369175A (en) | A kind of oyster for aquaculture arranges Area computing method and system | |
CN117058522A (en) | Remote sensing earth surface ecological change detection method integrating spectral neural network and leaf area index | |
Wang et al. | Automatic estimation of trunk cross sectional area using deep learning | |
Prerna et al. | Evaluation of LiDAR and image segmentation based classification techniques for automatic building footprint extraction for a segment of Atlantic County, New Jersey | |
CN109063577A (en) | Method is determined based on the satellite image segmentation best segmental scale of information gain-ratio | |
Xu et al. | Spatiotemporal distribution of cage and raft aquaculture in China's offshore waters using object-oriented random forest classifier | |
CN113129320B (en) | Method and system for extracting culture pond boundary from aerial orthographic image | |
Kibele | Submerged habitats from space: Increasing map production capacity with new methods and software |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171121 |