Cultivated land planting state detection method based on multi-source remote sensing image

Abstract

The invention provides a cultivated land planting state detection method based on a multi-source remote sensing image, and the method comprises the following steps: S1, obtaining remote sensing satellite image video data shot by a remote sensing satellite, and downloading the remote sensing satellite image video data shot by the remote sensing satellite to a cultivated land planting state detection platform; s2, after the cultivated land planting state detection platform receives the remote sensing satellite image video data, making the remote sensing satellite image video data into remote sensing satellite image data; s3, after the remote sensing satellite image video data are made into remote sensing satellite image data, background fuzzification processing is carried out on the remote sensing satellite image data, and a remote sensing satellite image processing image is obtained; and S4, acquiring a cultivated land planting area in the remote sensing satellite image processing image, and judging whether crops are planted in the cultivated land planting area or not.

Description

Technical Field

[0001] The present invention relates to the field of image processing technology, and in particular to a method for detecting cultivated land planting status based on multi-source remote sensing images. Background Art

[0002] Traditional methods for monitoring the planting status of cultivated land rely primarily on manual field surveys. This method requires a significant investment of manpower, material resources, and time. Workers must delve into the fields to record and compile statistics on information such as the crops planted, the planting area, and the crop growth status. However, manual field surveys have many limitations. On the one hand, their monitoring range is limited, making it difficult to cover large areas of cultivated land. This is especially true in areas with complex terrain and inconvenient transportation, where manual surveys are more difficult and costly. On the other hand, manual surveys are less timely and cannot capture dynamic changes in the planting status of cultivated land in a timely manner, making it difficult to meet the needs of modern agriculture for real-time monitoring and rapid response. Summary of the Invention

[0003] The present invention aims to at least solve the technical problems existing in the prior art, and in particular innovatively proposes a method for detecting the planting status of cultivated land based on multi-source remote sensing images.

[0004] In order to achieve the above-mentioned object of the present invention, the present invention provides a method for detecting the planting status of cultivated land based on multi-source remote sensing images, comprising the following steps:

[0005] S1, obtaining remote sensing satellite image video data taken by a remote sensing satellite, and downloading the remote sensing satellite image video data taken by the remote sensing satellite to a cultivated land planting status detection platform;

[0006] S2, after receiving the remote sensing satellite image video data, the cultivated land planting status detection platform converts the remote sensing satellite image video data into remote sensing satellite image data;

[0007] S3, after converting the remote sensing satellite image video data into remote sensing satellite image data, performing background blurring processing on the remote sensing satellite image data to obtain a remote sensing satellite image processed image;

[0008] S4, obtaining the cultivated land planting area in the remote sensing satellite image processing image, and determining whether there is crop planting in the cultivated land planting area.

[0009] In a preferred embodiment of the present invention, step S2 includes the following steps:

[0010] S21, obtaining video parameters of remote sensing satellite image video data, the video parameters including resolution Res, frame rate FPS, duration Tim, video format Cod, color depth Col (pixel value bits) or any combination thereof;

[0011] The video format Cod includes a video encoding format and a video encapsulation format. The video encoding format includes one of H.264, H.265, VP9, and AV1. The video encapsulation format includes one of MP4, MKV, MOV, AVI, and FLV.

[0012] S22: Obtain the number of remote sensing satellite image data according to the acquired duration Tim video parameter of the remote sensing satellite image video data. The method for obtaining the number of remote sensing satellite image data according to the acquired duration Tim video parameter of the remote sensing satellite image video data is:

[0013] Qua=Tim*FPS,

[0014] Wherein, Qua is the number of remote sensing satellite image data;

[0015] Tim is the duration of remote sensing satellite image video data;

[0016] FPS is the frame rate FPS of remote sensing satellite image video data;

[0017] The resolution Res of each remote sensing satellite image data is consistent with the resolution Res video parameter of the remote sensing satellite image video data.

[0018] In a preferred embodiment of the present invention, step S1 includes the following steps:

[0019] S11, determining the size of the remote sensing satellite image video data to be downloaded to the cultivated land planting status detection platform, dividing the remote sensing satellite image video data to be downloaded into U0 parts, and calculating the number of remote sensing satellite image video data parts U0 as follows:

[0020]

[0021] U0 is the number of copies of remote sensing satellite image video data;

[0022] floor() represents the rounding function;

[0023] datasize() is the data size function;

[0024] spdata is the remote sensing satellite image video data to be downloaded;

[0025] Datasize0 is the size of each remote sensing satellite image video data;

[0026] {N *} is a set of integers;

[0027] S12, obtaining video parameters of each remote sensing satellite image video data to be downloaded, the video parameters including one or any combination of resolution Res, frame rate FPS, duration Tim′, video format Cod, and color depth Col;

[0028] The video format Cod includes a video encoding format and a video encapsulation format. The video encoding format includes one of H.264, H.265, VP9, and AV1. The video encapsulation format includes one of MP4, MKV, MOV, AVI, and FLV.

[0029] S13, obtaining the number of remote sensing satellite image data to be downloaded based on the obtained duration Tim′ video parameter of each remote sensing satellite image video data to be downloaded. The method for obtaining the number of remote sensing satellite image data to be downloaded based on the obtained duration Tim′ video parameter of each remote sensing satellite image video data to be downloaded is:

[0030] Qua′=Tim′*FPS,

[0031] Wherein, Qua′ is the number of each remote sensing satellite image data to be downloaded;

[0032] Tim′ is the duration Tim′ of each remote sensing satellite image video data to be downloaded;

[0033] FPS is the frame rate FPS of remote sensing satellite image video data;

[0034] The resolution Res of each remote sensing satellite image data to be downloaded is consistent with the resolution Res video parameter of each remote sensing satellite image video data to be downloaded; through the above method, the duration Tim′ of each remote sensing satellite image video data to be downloaded can be converted into Qua′ remote sensing satellite images to be downloaded;

[0035] S14, calculating the image 0-1 digital string of the remote sensing satellite image video data to be downloaded in step S11;

[0036] S15, cyclically calculating the image 0-1 digital string and the pixel values of each channel in each remote sensing satellite image to be downloaded, to obtain a new remote sensing satellite image to be downloaded:

[0037]

[0038] Where R(μ,ν) New is the new value of the νth pixel corresponding to the μth pixel in the red channel from left to right and from top to bottom;

[0039] R(μ,ν) is the ν-th bit value in the pixel value corresponding to the μ-th pixel in the red color channel from left to right and from top to bottom. Here, the pixel value of each pixel is a binary pixel value, and the number of bits of the pixel value is equal to the color depth Col.

[0040] Imstr 0-1 New (γ) is the value of the γth bit in the new image 0-1 digital string;

[0041] Will Connect the image 0-1 digital string front and back to get a new image 0-1 digital string;

[0042]

[0043] Among them, G(μ,ν) New is the new value of the nth pixel corresponding to the μth pixel in the green channel from left to right and from top to bottom;

[0044] G(μ,ν) is the ν-th pixel value corresponding to the μ-th pixel in the green color channel from left to right and from top to bottom;

[0045] Imstr 0-1 New (γ) is the value of the γth bit in the new image 0-1 digital string;

[0046] The relationship between γ, μ and ν is:

[0047]

[0048] Among them, γ is the position of the new image 0-1 digital string from left to right;

[0049] μ is the pixel position from left to right and from top to bottom in the red / green / blue color channel; μ = 1, 2, 3, ..., O K *O G ;

[0050] is the color depth Col;

[0051] ν is the position of the pixel value from left to right,

[0052]

[0053] Among them, B(μ,ν) New is the new value of the νth pixel corresponding to the μth pixel in the blue color channel from left to right and from top to bottom;

[0054] B(μ,ν) is the ν-th pixel value corresponding to the μ-th pixel in the blue color channel from left to right and from top to bottom;

[0055] Imstr 0-1 New (γ) is the value of the γth bit in the new image 0-1 digital string;

[0056] Through the operation of step S15, a new remote sensing satellite image to be downloaded can be obtained at this time. All Qua' remote sensing satellite images to be downloaded are subjected to the operation of step S15 to obtain all new Qua' remote sensing satellite images to be downloaded;

[0057] S16, all the new remote sensing satellite images to be downloaded are combined into a new remote sensing satellite image video data to be downloaded in the order in which they were played. By operating in step S16, a new set of remote sensing satellite image video data to be downloaded can be obtained; U0 new sets of remote sensing satellite image video data to be downloaded can be obtained by simply operating in steps S15 and S16;

[0058] S17, transmit the image 0-1 digital string and U0 copies of new remote sensing satellite image video data to be downloaded to the cultivated land planting status detection platform.

[0059] In a preferred embodiment of the present invention, step S3 includes the following steps:

[0060] S31, dividing the remote sensing satellite image data into K0 pixel groups of uniform size. The number of groups is calculated as follows:

[0061]

[0062] K0 represents the number of pixel groups into which the remote sensing satellite image data is divided;

[0063] O K is the horizontal pixel number of remote sensing satellite image data;

[0064] O G is the number of vertical pixels of remote sensing satellite image data;

[0065] O K A factor of , and That is O K Can be divisibility;

[0066] O G A factor of , and That is O G Can be divisibility;

[0067] S32, performing pixel group processing on each pixel group using an image pixel group processing algorithm;

[0068] S33, determining the background judgment range of the remote sensing satellite processed image;

[0069] S34, obtaining the grayscale pixel value of each pixel in the remote sensing satellite processed image, and determining the pixel area where the pixels belong to the background judgment range are located as the background area of the remote sensing satellite processed image;

[0070] S35, performing fuzzy processing on the background area of the remote sensing satellite processed image to obtain a remote sensing satellite image processed image.

[0071] In a preferred embodiment of the present invention, the method for obtaining the cultivated land planting area in step S4 includes the following steps:

[0072] S41, acquiring historical image data containing rice crops, calculating the grayscale pixel value variance and the grayscale pixel value average of the historical image data containing rice crops, and obtaining a moving frame determination threshold based on the grayscale pixel value variance and the grayscale pixel value average;

[0073] S42, setting a moving frame, performing moving scanning on the remote sensing satellite image processing image, and calculating the moving determination value of each pixel in the remote sensing satellite image processing image according to the moving frame;

[0074] S43, determining the cultivated land and non-cultivated land areas in the remote sensing satellite image processing image according to the relationship between the moving frame judgment threshold and the moving determination value:

[0075] If f0≥d0, f0 is the moving frame judgment threshold, d0 is the moving determination value; then the cultivated land planting area in the remote sensing satellite image processing image;

[0076] If f0 ≥ d0, the non-cultivated land in the remote sensing satellite image processing image is the planting area.

[0077] In a preferred embodiment of the present invention, the method for determining whether there is crop planting in the cultivated land planting area in step S4 is:

[0078] S4-1, obtaining the date on which the remote sensing satellite image video data was taken, and after obtaining the date on which the remote sensing satellite image video data was taken, obtaining the grayscale pixel value range of the rice crop on that date;

[0079] S4-2, arrange the grayscale pixel values in the cultivated area in ascending order, remove the first and last W pixel values, and average the remaining grayscale pixel values:

[0080]

[0081] Among them, p final is the grayscale pixel value of rice crops in the cultivated land area;

[0082] O K is the horizontal pixel number of remote sensing satellite image data;

[0083] O G is the number of vertical pixels of remote sensing satellite image data;

[0084] W is the number of pixels removed at the first position;

[0085] p i ′ is the grayscale pixel value of each pixel in the cultivated land planting area;

[0086] p i ″ is the grayscale pixel value of each pixel removed from the front of the cultivated land planting area;

[0087] p i ″′ is the grayscale pixel value of each pixel that is finally removed in the cultivated land planting area;

[0088] S4-3, the relationship between the grayscale pixel value of the rice crop in the cultivated land area and the grayscale pixel value range of the rice crop on that date is:

[0089] If the grayscale pixel value of the rice crop in the cultivated land area is within the grayscale pixel value range of the rice crop on that date, then the rice crop exists on that date;

[0090] Otherwise, there is no rice crop planted on that date.

[0091] The present invention also discloses a computer system, comprising:

[0092] processor;

[0093] a memory for storing processor-executable instructions;

[0094] Wherein, the processor is configured to implement the method for detecting cultivated land planting status based on multi-source remote sensing images when executing the executable instructions.

[0095] The present invention also discloses a computer-readable storage medium, comprising:

[0096] a memory having a computer program stored thereon;

[0097] A processor is used to execute the program in the memory to implement the method for detecting the planting status of cultivated land based on multi-source remote sensing images.

[0098] In summary, due to the adoption of the above technical solution, the present invention can efficiently detect the planting status of cultivated land and can control the planting status of cultivated land at any time throughout the day.

[0099] Additional aspects and advantages of the present invention will be set forth in part in the description which follows and, in part, will be obvious from the description which follows, or may be learned by practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0100] The above and / or additional aspects and advantages of the present invention will become apparent and readily understood from the description of the embodiments with reference to the following drawings, in which:

[0101] Figure 1 It is a schematic block diagram of the process of the present invention.

[0102] Figure 2 It is a schematic diagram of the image of the present invention.

[0103] Figure 3 It is a schematic diagram of the image of the present invention.

[0104] Figure 4 It is a schematic diagram of the image of the present invention.

[0105] Figure 5 It is a schematic diagram of the image of the present invention.

[0106] Figure 6 It is a schematic diagram of the image of the present invention.

[0107] Figure 7 It is a schematic diagram of the image of the present invention.

[0108] Figure 8 It is a schematic diagram of the image of the present invention.

[0109] Figure 9 It is a schematic diagram of the image of the present invention.

[0110] Figure 10 It is a schematic diagram of the image of the present invention.

[0111] Figure 11 It is a schematic diagram of the image of the present invention. DETAILED DESCRIPTION

[0112] The following describes embodiments of the present invention in detail. Examples of the embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are intended only to explain the present invention and are not to be construed as limiting the present invention.

[0113] The present invention discloses a method for detecting the planting status of cultivated land based on multi-source remote sensing images. Figures 1 to 11As shown, the following steps are included:

[0114] S1, obtaining remote sensing satellite image video data taken by a remote sensing satellite, and downloading the remote sensing satellite image video data taken by the remote sensing satellite to a cultivated land planting status detection platform;

[0115] S2, after receiving the remote sensing satellite image video data, the cultivated land planting status detection platform converts the remote sensing satellite image video data into remote sensing satellite image data;

[0116] S3, after converting the remote sensing satellite image video data into remote sensing satellite image data, performing background blurring processing on the remote sensing satellite image data to obtain a remote sensing satellite image processed image;

[0117] S4, obtaining the cultivated land planting area in the remote sensing satellite image processing image, and determining whether there is crop planting in the cultivated land planting area.

[0118] In a preferred embodiment of the present invention, step S2 includes the following steps:

[0119] S21, obtaining video parameters of remote sensing satellite image video data, the video parameters including resolution Res, frame rate FPS, duration Tim, video format Cod, color depth Col (pixel value bits) or any combination thereof;

[0120] The video format Cod includes a video encoding format and a video encapsulation format. The video encoding format includes one of H.264, H.265, VP9, and AV1. The video encapsulation format includes one of MP4, MKV, MOV, AVI, and FLV.

[0121] S22: Obtain the number of remote sensing satellite image data according to the acquired duration Tim video parameter of the remote sensing satellite image video data. The method for obtaining the number of remote sensing satellite image data according to the acquired duration Tim video parameter of the remote sensing satellite image video data is:

[0122] Qua=Tim*FPS,

[0123] Wherein, Qua is the number of remote sensing satellite image data;

[0124] Tim is the duration of remote sensing satellite image video data;

[0125] FPS is the frame rate FPS of remote sensing satellite image video data;

[0126] The resolution Res of each remote sensing satellite image data is consistent with the resolution Res video parameter of the remote sensing satellite image video data.

[0127] In a preferred embodiment of the present invention, step S1 includes the following steps:

[0128] S11, determining the size of the remote sensing satellite image video data to be downloaded to the cultivated land planting status detection platform, and dividing the remote sensing satellite image video data to be downloaded into U0 parts, which are respectively the first remote sensing satellite image video data to be downloaded, the second remote sensing satellite image video data to be downloaded, the third remote sensing satellite image video data to be downloaded, ..., the U0th remote sensing satellite image video data to be downloaded in order of playback; wherein the calculation formula of the number of remote sensing satellite image video data parts U0 is:

[0129]

[0130] U0 is the number of copies of remote sensing satellite image video data;

[0131] floor() represents the rounding function;

[0132] datasize() is the data size function;

[0133] spdata is the remote sensing satellite image video data to be downloaded;

[0134] Datasize0 is the size of each remote sensing satellite image video data;

[0135] {N *} is a set of integers;

[0136] S12, obtaining video parameters of each remote sensing satellite image video data to be downloaded, the video parameters including one or any combination of resolution Res, frame rate FPS, duration Tim′, video format Cod, and color depth Col (number of pixel value bits);

[0137] The video format Cod includes a video encoding format and a video encapsulation format. The video encoding format includes one of H.264, H.265, VP9, and AV1. The video encapsulation format includes one of MP4, MKV, MOV, AVI, and FLV.

[0138] S13, obtaining the number of remote sensing satellite image data to be downloaded based on the obtained duration Tim′ video parameter of each remote sensing satellite image video data to be downloaded. The method for obtaining the number of remote sensing satellite image data to be downloaded based on the obtained duration Tim′ video parameter of each remote sensing satellite image video data to be downloaded is:

[0139] Qua′=Tim′*FPS,

[0140] Wherein, Qua′ is the number of each remote sensing satellite image data to be downloaded;

[0141] Tim′ is the duration Tim′ of each remote sensing satellite image video data to be downloaded;

[0142] FPS is the frame rate FPS of remote sensing satellite image video data;

[0143] The resolution Res of each remote sensing satellite image data to be downloaded is consistent with the resolution Res video parameter of each remote sensing satellite image video data to be downloaded; through the above method, the duration Tim′ of each remote sensing satellite image video data to be downloaded can be converted into Qua′ remote sensing satellite images to be downloaded;

[0144] S14, calculating the image 0-1 digital string of the remote sensing satellite image video data to be downloaded in step S11, the calculation method of the image 0-1 digital string of the remote sensing satellite image video data to be downloaded is:

[0145] Imstr 0-1 =MD5Imstr 0-1 (spingdata),

[0146] Imstr 0-1 It is a 0-1 digital string of images, and its length is 128 bits;

[0147] MD5Imstr 0-1 () is the calculation function of the image 0-1 digital string, preferably using the MD5 digest function;

[0148] spingdata is the remote sensing satellite image video data to be downloaded;

[0149] S15, cyclically calculating the image 0-1 digital string and the pixel values of each channel in each remote sensing satellite image to be downloaded, to obtain a new remote sensing satellite image to be downloaded:

[0150]

[0151] Among them, R(μ,ν) New is the new value of the νth pixel corresponding to the μth pixel in the red channel from left to right and from top to bottom;

[0152] R(μ,ν) is the ν-th bit value in the pixel value corresponding to the μ-th pixel in the red color channel from left to right and from top to bottom. Here, the pixel value of each pixel is a binary pixel value, and the number of bits of the pixel value is equal to the color depth Col.

[0153] Imstr 0-1 New (γ) is the value of the γth bit in the new image 0-1 digital string;

[0154] Will The image 0-1 digital string is connected front and back to get a new image 0-1 digital string:

[0155]

[0156] Imstr 0-1 New New image 0-1 digital string;

[0157]

[0158] in, Indicates the number of 0-1 digital strings in the image;

[0159] is the color depth Col;

[0160] floor() represents the rounding function;

[0161] O K is the horizontal pixel number of remote sensing satellite image data;

[0162] O G is the number of vertical pixels of remote sensing satellite image data;

[0163]

[0164] Among them, G(μ,ν) New is the new value of the nth pixel corresponding to the μth pixel in the green channel from left to right and from top to bottom;

[0165] G(μ,ν) is the ν-th pixel value corresponding to the μ-th pixel in the green color channel from left to right and from top to bottom;

[0166] Imstr 0-1 New (γ) is the value of the γth bit in the new image 0-1 digital string;

[0167] The relationship between γ, μ and ν is:

[0168]

[0169] Among them, γ is the position of the new image 0-1 digital string from left to right;

[0170] μ is the pixel position from left to right and from top to bottom in the red / green / blue color channel; μ = 1, 2, 3, ..., O K *O G ;

[0171] is the color depth Col;

[0172] ν is the position of the pixel value from left to right,

[0173]

[0174] Among them, B(μ,ν) New is the new value of the νth pixel corresponding to the μth pixel in the blue color channel from left to right and from top to bottom;

[0175] B(μ,ν) is the ν-th pixel value corresponding to the μ-th pixel in the blue color channel from left to right and from top to bottom;

[0176] Imstr 0-1 New (γ) is the value of the γth bit in the new image 0-1 digital string;

[0177] Through the operation of step S15, a new remote sensing satellite image to be downloaded can be obtained at this time. All Qua' remote sensing satellite images to be downloaded are subjected to the operation of step S15 to obtain all new Qua' remote sensing satellite images to be downloaded;

[0178] S16, all the new remote sensing satellite images to be downloaded are combined into a new remote sensing satellite image video data to be downloaded in the order in which they were played. By operating in step S16, a new set of remote sensing satellite image video data to be downloaded can be obtained; U0 new sets of remote sensing satellite image video data to be downloaded can be obtained by simply operating in steps S15 and S16;

[0179] S17, transmit the image 0-1 digital string and U0 copies of new remote sensing satellite image video data to be downloaded to the cultivated land planting status detection platform.

[0180] In a preferred embodiment of the present invention, step S2 includes the following steps:

[0181] S2-1. After the cultivated land planting status detection platform receives the image 0-1 digital string and U0 downloaded remote sensing satellite image video data, the cultivated land planting status detection platform obtains the number of downloaded remote sensing satellite image data according to the duration Tim′ video parameter of each received downloaded remote sensing satellite image video data. The method for obtaining the number of downloaded remote sensing satellite image data according to the duration Tim″ video parameter of each received remote sensing satellite image video data to be downloaded is as follows:

[0182] Qua″=Tim″*FPS,

[0183] Wherein, Qua″ is the number of remote sensing satellite image data to be downloaded;

[0184] Tim″ is the duration Tim″ of each downloaded remote sensing satellite image video data;

[0185] FPS is the frame rate FPS of remote sensing satellite image video data;

[0186] The resolution Res of each remote sensing satellite image data is consistent with the resolution Res video parameter of each downloaded remote sensing satellite image video data; through the above method, the duration Tim′ of each downloaded remote sensing satellite image video data can be converted into Qua′ downloaded remote sensing satellite images;

[0187] S2-2, the received image 0-1 digital string is cyclically calculated with the pixel value of each channel in each downloaded remote sensing satellite image to obtain the cultivated land platform remote sensing satellite image:

[0188]

[0189] Among them, R′(μ′,ν′) New is the new value of the ν′th platform of the pixel value corresponding to the μ′th pixel from left to right and from top to bottom in the red color channel;

[0190] R′(μ′,ν′) is the ν′th platform value in the pixel value corresponding to the μ′th pixel from left to right and from top to bottom in the red color channel. Here, the pixel value of each pixel is a binary pixel value, and the number of bits of the pixel value is equal to the color depth Col.

[0191] Im′str 0-1 New (γ′) is the value of the γ′th position in the new platform image 0-1 digital string;

[0192] Will The received image 0-1 digital strings are connected front to back to get the new platform image 0-1 digital string:

[0193]

[0194] Im′str 0-1 New New platform image 0-1 digital string;

[0195]

[0196] in, Indicates the number of connected 0-1 digital strings of the received image;

[0197] is the color depth Col;

[0198] floor() represents the rounding function;

[0199] O K is the horizontal pixel number of remote sensing satellite image data;

[0200] O G is the number of vertical pixels of remote sensing satellite image data;

[0201]

[0202] Among them, G′(μ′,ν′) New is the new value of the ν′th platform of the pixel value corresponding to the μ′th pixel from left to right and from top to bottom in the green color channel;

[0203] G′(μ′,ν′) is the ν′th platform value of the pixel value corresponding to the μ′th pixel from left to right and from top to bottom in the green color channel;

[0204] Im′str 0-1 New (γ′) is the value of the γ′th position in the new platform image 0-1 digital string;

[0205] The relationship between γ′, μ′ and ν′ is:

[0206]

[0207] Where γ′ is the position of the new platform image 0-1 digital string from left to right;

[0208] μ′ is the platform pixel position from left to right and from top to bottom in the red / green / blue color channel; μ′=1, 2, 3, ..., O K *O G ;

[0209] is the color depth Col;

[0210] ν′ is the position of the platform pixel value from left to right,

[0211]

[0212] Among them, B′(μ′,ν′) New is the new value of the ν′th platform of the pixel value corresponding to the μ′th pixel from left to right and from top to bottom in the blue color channel;

[0213] B′(μ′,ν′) is the ν′th platform value of the pixel value corresponding to the μ′th pixel from left to right and from top to bottom in the blue color channel;

[0214] Im′str 0-1 New (γ′) is the value of the γ′th position in the new platform image 0-1 digital string;

[0215] Through the operation of step S2-2, a new downloaded remote sensing satellite image can be obtained at this time. All Quad″ downloaded remote sensing satellite images can be obtained by performing the operation of step S2-2;

[0216] S2-3, all newly downloaded remote sensing satellite image images are combined into new remote sensing satellite image video data to be downloaded in the order of playback. By performing step S2-3, a new set of downloaded remote sensing satellite image video data can be obtained. U0 new sets of downloaded remote sensing satellite image video data can be obtained by simply performing steps S2-2 and S2-3.

[0217] S2-4, composing the U0 copies of newly downloaded remote sensing satellite image video data into a new entire downloaded remote sensing satellite image video data in the order of playback.

[0218] In a preferred embodiment of the present invention, step S2-5 is included to calculate the new image 0-1 digital string of the entire downloaded remote sensing satellite image video data in step S2-4. The calculation method of the new image 0-1 digital string of the entire downloaded remote sensing satellite image video data is:

[0219] Im″str 0-1 =MD5Imstr 0-1 (spdata′),

[0220] Im″str 0-1 It is the image 0-1 digital string of the newly downloaded remote sensing satellite image video data, and its length is 128 bits;

[0221] MD5Imstr 0-1 () is the calculation function of the image 0-1 digital string, preferably using the MD5 digest function;

[0222] spdata′ is the newly downloaded remote sensing satellite image video data;

[0223] If the image 0-1 digital string of the new entire downloaded remote sensing satellite image video data is consistent with the image 0-1 digital string received in step S2-1, then the new entire downloaded remote sensing satellite image video data is the remote sensing satellite image video data obtained in step S1; otherwise, the new entire downloaded remote sensing satellite image video data is not the remote sensing satellite image video data obtained in step S1.

[0224] In a preferred embodiment of the present invention, step S3 includes the following steps:

[0225] S31, dividing the remote sensing satellite image data into K0 pixel groups of uniform size. The number of groups is calculated as follows:

[0226]

[0227] K0 represents the number of pixel groups into which the remote sensing satellite image data is divided;

[0228] O K is the horizontal pixel number of remote sensing satellite image data;

[0229] O G is the number of vertical pixels of remote sensing satellite image data;

[0230] O K A factor of , and That is O K Can be divisibility;

[0231] O G A factor of , and That is O G Can be divisibility;

[0232] S32, performing pixel group processing on each pixel group using an image pixel group processing algorithm to obtain a remote sensing satellite processed image; the calculation formula of the image pixel group processing algorithm is:

[0233]

[0234] Among them, K υ Represents the grayscale pixel value of the υth pixel group in the remote sensing satellite processed image;

[0235] Represents the grayscale pixel value of the υth pixel group in the remote sensing satellite image;

[0236] m υ Represents the mean grayscale pixel value of the υth pixel group in the remote sensing satellite image; p υ,i is the grayscale pixel value of the i-th pixel in the υ-th pixel group, p υ,1 is the grayscale pixel value of the first pixel in the υth pixel group, p υ,2 is the grayscale pixel value of the second pixel in the υth pixel group, p υ,3 is the grayscale pixel value of the 3rd pixel in the υth pixel group, is the first pixel in the υth pixel group Grayscale pixel value of the pixel;

[0237] K ε Indicates the group adjustment threshold coefficient, which can be 2 here;

[0238] K0 represents the number of pixel groups into which the remote sensing satellite image data is divided;

[0239] S33, determine the background judgment range of the remote sensing satellite processed image, the background judgment range of the remote sensing satellite processed image (p min ,p max ) is determined as follows:

[0240]

[0241] Among them, p min Indicates the minimum value on the left side of the background judgment range;

[0242] p max Indicates the maximum value on the right side of the background judgment range;

[0243] O represents the number of pixels in the image processed by the remote sensing satellite; O=O K *O G , O K is the number of horizontal pixels of remote sensing satellite image data; O G is the number of vertical pixels of remote sensing satellite image data;

[0244] p j Represents the pixel value in the grayscale pixel value set of the jth row in the remote sensing satellite processed image; j = 1, 2, 3, ..., O G ;

[0245] Represents the average value of all grayscale pixel values in the remote sensing satellite processed image; p ji is the grayscale pixel value of the j-th row and i-th column in the image;

[0246] p0 is the pixel value adjustment threshold coefficient, which is 2 here;

[0247] min{} means taking the minimum value function;

[0248] max{} represents the minimum value function;

[0249] S34, obtaining the grayscale pixel value of each pixel in the remote sensing satellite processed image, and determining the pixel area where the pixels belong to the background judgment range are located as the background area of the remote sensing satellite processed image;

[0250] S35, performing fuzzy processing on the background area of the remote sensing satellite processed image to obtain a remote sensing satellite image processed image.

[0251] In a preferred embodiment of the present invention, the method for obtaining the cultivated land planting area in step S4 includes the following steps:

[0252] S41, obtaining historical image data containing rice crops, calculating the grayscale pixel value variance and the grayscale pixel value average of the historical image data containing rice crops, and obtaining a moving frame judgment threshold according to the grayscale pixel value variance and the grayscale pixel value average:

[0253]

[0254] Among them, f0 is the moving frame judgment threshold;

[0255] is the average grayscale pixel value of historical image data;

[0256] δ0 is the variance of the grayscale pixel value of the historical image data;

[0257] S42, setting a moving frame, performing moving scanning on the remote sensing satellite image processing image, and calculating the moving determination value of each pixel in the remote sensing satellite image processing image according to the moving frame;

[0258]

[0259] Among them, d0 is the movement determination value;

[0260] l0 is the moving step of the moving frame;

[0261] K0 is the number of horizontal pixels of the set moving frame;

[0262] G0 is the number of vertical pixels of the set moving frame;

[0263] O K is the horizontal pixel number of remote sensing satellite image data;

[0264] O G is the number of vertical pixels of remote sensing satellite image data;

[0265] p τ is the pixel value of the τth pixel in the remote sensing satellite image processing image;

[0266] p(0,0) is the pixel value of the center point of the remote sensing satellite image processing image;

[0267] S43, determining the cultivated land and non-cultivated land areas in the remote sensing satellite image processing image according to the relationship between the moving frame judgment threshold and the moving determination value:

[0268] If f0 ≥ d0, f0 is the moving frame judgment threshold, d0 is the moving determination value; then the cultivated land planting area in the remote sensing satellite image processing image;

[0269] If f0 ≥ d0, the non-cultivated land in the remote sensing satellite image processing image is the planting area.

[0270] In a preferred embodiment of the present invention, the method for determining whether there is crop planting in the cultivated land planting area in step S4 is:

[0271] S4-1, obtaining the date on which the remote sensing satellite image video data was taken, and after obtaining the date on which the remote sensing satellite image video data was taken, obtaining the grayscale pixel value range of the rice crop on that date;

[0272] S4-2, arrange the grayscale pixel values in the cultivated area in ascending order, remove the first and last W pixel values, and average the remaining grayscale pixel values:

[0273]

[0274] Among them, p final is the grayscale pixel value of rice crops in the cultivated land area;

[0275] O K is the horizontal pixel number of remote sensing satellite image data;

[0276] O G is the number of vertical pixels of remote sensing satellite image data;

[0277] W is the number of pixels removed at the first position;

[0278] p i ′ is the grayscale pixel value of each pixel in the cultivated land planting area;

[0279] p i ″ is the grayscale pixel value of each pixel removed from the front of the cultivated land planting area;

[0280] p i ″′ is the grayscale pixel value of each pixel that is finally removed in the cultivated land planting area;

[0281] S4-3, the relationship between the grayscale pixel value of the rice crop in the cultivated land area and the grayscale pixel value range of the rice crop on that date is:

[0282] If the grayscale pixel value of the rice crop in the cultivated land area is within the grayscale pixel value range of the rice crop on that date, then the rice crop exists on that date;

[0283] Otherwise, there is no rice crop planted on that date.

[0284] The present invention also discloses a computer system, comprising:

[0285] processor;

[0286] a memory for storing processor-executable instructions;

[0287] Wherein, the processor is configured to implement the method for detecting cultivated land planting status based on multi-source remote sensing images when executing the executable instructions.

[0288] In a preferred embodiment of the present invention, step S0 is further included before step S1, using a mobile phone or an office computer to enter the cultivated land planting status detection platform; after using a mobile phone or an office computer to enter the cultivated land planting status detection platform, the abnormal land use in the cultivated land planting area is checked.

[0289] In a preferred embodiment of the present invention, in step S0, the method of using a mobile phone or an office computer to access the cultivated land planting status detection platform includes the following steps:

[0290] S01, enter the account number and password on the login interface of the cultivated land planting status detection platform; the account number and password must be no less than 6 digits;

[0291] S02, after entering the account number and password on the farmland planting status detection platform login interface, the input account number and password are securely encrypted;

[0292] S03, after securely encrypting the input account number and password, the securely encrypted account number and password are sent to the cultivated land planting status detection platform verification platform;

[0293] S04. After the cultivated land planting status detection platform verification platform receives the securely encrypted account and password, it decrypts the received securely encrypted account and password; it can also use the public key held by the login end to securely encrypt the input account and password; the corresponding cultivated land planting status detection platform verification platform uses the private key it holds to decrypt the received securely encrypted account and password, where the private key is held by the platform and the public key is held by the login end, and the private key and public key are a pair of asymmetric encryption and decryption keys.

[0294] S05, if the decrypted account and password are both present in the farmland planting status detection platform verification platform, then enter the farmland planting status detection platform;

[0295] Otherwise, the account and password entered in the login interface of the cultivated land planting status detection platform are incorrect, and you need to re-enter the correct account and password.

[0296] In a preferred embodiment of the present invention, the method for securely encrypting the input account number and password in step S02 includes the following steps:

[0297] S021, obtain the current time and convert it into a standard ten-digit format. The current time includes year, month, day, and hour. For example, if the obtained time is 17:00 on March 26, 2025, the standard ten-digit format is 2025032617. If the obtained time is 21:00 on April 16, 2024, the standard ten-digit format is 2024041621.

[0298] S022, determine the relationship between the number of digits of the input account number and password and 10:

[0299] S0221, if the number of digits in the account number is less than 10, remove the leading e1 digits from the tens-digit standard time, until the digits after removal equal the number of digits in the account number; then proceed to the next step; e1 = 10 - e0, where e0 is the number of digits in the account number and e1 is the number of digits removed from the time;

[0300] If the number of digits in the account number is greater than 10, then add digits e2 and e3 to the front of the tens-digit standard time, and the added digits are equal to the number of digits in the account number; then proceed to the next step; e2 = e0 - 10, where e0 is the number of digits in the account number, e2 is the number of digits in the added time, and e3 is the last digit in the tens-digit standard format;

[0301] If the number of digits in the account number is equal to 10, proceed to the next step;

[0302] S0222, if the number of digits in the password is less than 10, remove the f1 digits in front of the ten-digit standard time, and the number after removal is equal to the number of digits in the password; then proceed to the next step; f1 = 10 - f0, where f0 is the number of digits in the password and f1 is the number of digits in the time removed;

[0303] If the number of digits in the password is greater than 10, then add f2 and f3 to the front of the ten-digit standard format time, and the added value is equal to the number of digits in the password; then proceed to the next step; f2 = f0-10, where f0 is the number of digits in the password, f2 is the number of digits in the added time, and f3 is the last digit in the ten-digit standard format;

[0304] If the number of digits in the password is 10, proceed to the next step. For example, if the account number and password are Zhanghao123 and miMa3211, and the time in the standard format of tens digits is 2025050515, then the times after processing are 52025050515 and 25050515 respectively. If the account number and password are Zhanghao12332 and miMa3211222, and the time in the standard format of tens digits is 2025050515, then the times after processing are 5552025050515 and 552025050515 respectively.

[0305] The account and password remain unchanged, but the time corresponding to the account and password may change;

[0306] S023, shift the first digit of the account number and password backward by f4 digits, f4 = floor(f3 / 2), where f4 is the number of digits to be shifted backward. For example, if the account number and password are Zhang1hao3 and mi2Ma2, respectively, and the time in the tens-digit standard format is 2025060715, then the processed account number and password are haZng1hao3 and i2mMa2, respectively. If the account number and password are Zhang1hao3 and mi2Ma2, respectively, and the time in the tens-digit standard format is 2025060716, then the processed account number and password are hanZg1hao3 and i2Mma2, respectively.

[0307] S024, convert each character in the account number, password, and time into 7-bit binary ASCII form, calculate the account number in ASCII form with the time in ASCII form, and calculate the password in ASCII form with the time in ASCII form to obtain the encrypted account number and encrypted password respectively.

[0308]

[0309] g′2(φ) represents the value of the φth digit in the ASCII encrypted account number;

[0310] g2(φ) represents the value of the φth digit in the ASCII account number;

[0311] t2(φ) represents the value of the φth digit in the ASCII format time;

[0312] φ=1, 2, 3, ..., φ1; φ1 is the number of digits of the account number in ASCII format;

[0313] φ1=7φ2, φ2 is the number of digits in the account number;

[0314] ◇ represents the exchange value operator;

[0315]

[0316] v′2(θ) represents the value of the θth bit in the ASCII encrypted password;

[0317] v2(θ) represents the value of the θth bit in the ASCII password;

[0318] t′2(θ) represents the value of the θth digit in the ASCII form of time;

[0319] θ=1, 2, 3, ..., θ1; θ1 is the number of digits of the ASCII password;

[0320] θ1=7θ2, θ2 is the number of digits in the password;

[0321] Convert the ASCII encrypted account and ASCII encrypted password into character encrypted account and encrypted password.

[0322] For example, if the account and password entered are ZHangHao12234 and MIma912 respectively, the obtained time is 23:00 on January 22, 2025;

[0323] In step S021, the time in the tens-digit standard format is 2025012223;

[0324] In step S022, the account number is ZHangHao12234, and the time corresponding to the account number ZHangHao12234 is 3332025012223;

[0325] The password is MIma912, and the time corresponding to the password MIma912 is 5012223;

[0326] In step S023, the account number is HanZgHao12234, and the time corresponding to the account number HanZgHao12234 is 3332025012223;

[0327] The password is ImaM912, and the time corresponding to the password ImaM912 is 5012223;

[0328] In step S024, the ASCII format of the account number HanZgHao12234 is 1001000110000111011101011010110011110010001100001110111101100010110010011001100110110100, and the ASCII format of the time corresponding to the account number HanZgHao12234, 3332025012223, is 01100110110011011001 1011001001100000110010011010101100000110001011001001100100110010011; the encrypted account number in ASCII form is 111101110100101011101110100010101011111101010101010101011111000000000000000000000000000000010000111, and the corresponding encrypted account number is {R]hWzT_NUT NUT NUTSOH BEL, see Table 1 for details.

[0329] Table 1 Account and time and the ASCII code corresponding to each character in the encrypted account

[0330] character ASCII character ASCII character ASCII H 1001000 3 0110011 { 1111011 a 1100001 3 0110011 R 1010010 n 1101110 3 0110011 ] 1011101 Z 1011010 2 0110010 h 1101000 g 1100111 0 0110000 W 1010111 H 1001000 2 0110010 z 1111010 a 1100001 5 0110101 T 1010100 o 1101111 0 0110000 _ 1011111 1 0110001 1 0110001 NUT 0000000 2 0110010 2 0110010 NUT 0000000 2 0110010 2 0110010 NUT 0000000 3 0110011 2 0110010 SOH 0000001 4 0110100 3 0110011 BEL 0000111

[0331] The account number in ASCII form of the password ImaM912 is 1001001110110111000011001101011100101100010110010, and the time 5012223 corresponding to the password ImaM912 is 0110101011000001100010110010011001001100110011; the encrypted password in ASCII form is 1111100101110110100001111111000101100000110000001, and the corresponding encrypted password is |]P DEL VT ETX SOH, see Table 2 for details.

[0332] Table 2 Password and time and ASCII code corresponding to each character in the encrypted password

[0333] character ASCII character ASCII character ASCII I 1001001 5 0110101 | 1111100 m 1101101 0 0110000 ] 1011101 a 1100001 1 0110001 P 1010000 M 1001101 2 0110010 DEL 1111111 9 0111001 2 0110010 VT 0001011 1 0110001 2 0110010 ETX 0000011 2 0110010 3 0110011 SOH 0000001

[0334] In a preferred embodiment of the present invention, the method for decrypting the received securely encrypted account and password in step S04 includes the following steps:

[0335] S041, obtain the current time and convert it into a standard decimal format;

[0336] S042, determine the relationship between the number of digits of the received encrypted account number and the encrypted password and 10:

[0337] S0221, if the number of digits in the encrypted account number is less than 10, remove the leading e1′ digits from the ten-digit standard time format until the digits after removal equal the number of digits in the encrypted account number; then proceed to the next step; e1′ = 10-e′0, where e′0 is the number of digits in the encrypted account number and e1′ is the number of digits removed from the time;

[0338] If the number of digits in the encrypted account number is greater than 10, then add e'2 digits e3' to the front of the time in the ten-digit standard format, and the added digits equal to the number of digits in the encrypted account number; then proceed to the next step; e'2 = e'0 - 10, where e'0 is the number of digits in the encrypted account number, e'2 is the number of digits in the added time, and e'3 is the last digit in the ten-digit standard format;

[0339] If the number of digits in the encrypted account is equal to 10, proceed to the next step;

[0340] S0222, if the number of digits in the encrypted password is less than 10, remove the leading f1′ digits from the ten-digit standard format time until the number after removal equals the number of digits in the encrypted password; then proceed to the next step; f1′ = 10-f0′, where f0′ is the number of digits in the encrypted password and f1′ is the number of digits in the removed time;

[0341] If the number of digits in the encrypted password is greater than 10, then add f2′ and f3′ to the front of the time in the ten-digit standard format, and the added value is equal to the number of digits in the encrypted password; then proceed to the next step; f2′=f0′-10, where f0′ is the number of digits in the encrypted password, f2′ is the number of digits in the added time, and f3′ is the last digit in the ten-digit standard format;

[0342] If the number of bits in the encrypted password is equal to 10, proceed to the next step;

[0343] The encryption account and encryption password remain unchanged, but the time corresponding to the encryption account and encryption password may change;

[0344] S043, shift the first f′4 digits of the encrypted account number and the encrypted password backward by one digit, f′4=floor(f3′ / 2), where f′4 is the total number of digits;

[0345] S044, converting each character in the encrypted account, encrypted password, and time into 7-bit binary ASCII form, calculating the encrypted account in ASCII form with the time in ASCII form, and calculating the encrypted password in ASCII form with the time in ASCII form, to obtain the decrypted account and decrypted password respectively.

[0346]

[0347] g′2″(φ′) represents the value of the φ′th digit in the decrypted account number in ASCII form;

[0348] g′2′(φ′) represents the value of the φ′th digit in the ASCII encrypted account number;

[0349] t2″(φ′) represents the value of the φ′th digit in the ASCII format time;

[0350] φ′=1, 2, 3, ..., φ1′; φ1′ is the number of digits of the encrypted account number in ASCII form;

[0351] φ1′=7φ2′, φ2′ is the number of digits of the encrypted account number;

[0352] ◇ represents the exchange value operator;

[0353]

[0354] v″′2(θ′) represents the value of the θ′th digit in the ASCII decrypted password;

[0355] v″2(θ′) represents the value of the θ′th digit in the ASCII encrypted password;

[0356] t″2(θ′) represents the value of the θ′th digit in the ASCII form of the time;

[0357] θ′=1, 2, 3, ..., θ′1; θ′1 is the number of bits of the ASCII encrypted password;

[0358] θ′1=7θ′2, θ′2 is the number of bits of the encryption code;

[0359] Convert the decrypted account number and decrypted password in ASCII format into the decrypted account number and decrypted password in character format.

[0360] The present invention also discloses a computer-readable storage medium, comprising:

[0361] a memory having a computer program stored thereon;

[0362] A processor is used to execute the program in the memory to implement the method for detecting the planting status of cultivated land based on multi-source remote sensing images.

[0363] While embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions, and variations may be made to the embodiments without departing from the principles and spirit of the invention, and that the scope of the invention is defined by the claims and their equivalents.

Claims

1. A method for detecting cultivated land planting status based on multi-source remote sensing images, characterized in that: The following steps are involved: S1, obtaining remote sensing satellite image video data taken by a remote sensing satellite, and downloading the remote sensing satellite image video data taken by the remote sensing satellite to a cultivated land planting status detection platform; S2, after receiving the remote sensing satellite image video data, the cultivated land planting status detection platform converts the remote sensing satellite image video data into remote sensing satellite image data; S3, after converting the remote sensing satellite image video data into remote sensing satellite image data, performing background blurring processing on the remote sensing satellite image data to obtain a remote sensing satellite image processed image; S4, obtaining the cultivated land planting area in the remote sensing satellite image processing image, and determining whether there is crop planting in the cultivated land planting area.

2. The method for detecting cultivated land planting status based on multi-source remote sensing images according to claim 1, characterized in that: Step S2 includes the following steps: S21, obtaining video parameters of remote sensing satellite image video data, where the video parameters include one or any combination of resolution Res, frame rate FPS, duration Tim, video format Cod, and color depth Col; The video format Cod includes a video encoding format and a video encapsulation format. The video encoding format includes one of H.264, H.265, VP9, and AV1. The video encapsulation format includes one of MP4, MKV, MOV, AVI, and FLV. S22: Obtain the number of remote sensing satellite image data according to the acquired duration Tim video parameter of the remote sensing satellite image video data. The method for obtaining the number of remote sensing satellite image data according to the acquired duration Tim video parameter of the remote sensing satellite image video data is: Qua=Tim*FPS, Wherein, Qua is the number of remote sensing satellite image data; Tim is the duration of remote sensing satellite image video data; FPS is the frame rate FPS of remote sensing satellite image video data; The resolution Res of each remote sensing satellite image data is consistent with the resolution Res video parameter of the remote sensing satellite image video data.

3. The method for detecting cultivated land planting status based on multi-source remote sensing images according to claim 1, characterized in that: Step S1 includes the following steps: S11, determining the size of the remote sensing satellite image video data to be downloaded to the cultivated land planting status detection platform, dividing the remote sensing satellite image video data to be downloaded into U0 parts, and calculating the number of remote sensing satellite image video data parts U0 as follows: U0 is the number of copies of remote sensing satellite image video data; floor() represents the rounding function; datasize() is the data size function; spdata is the remote sensing satellite image video data to be downloaded; Datasize0 is the size of each remote sensing satellite image video data; {N * } is a set of integers; S12, obtaining video parameters of each remote sensing satellite image video data to be downloaded, the video parameters including one or any combination of resolution Res, frame rate FPS, duration Tim′, video format Cod, and color depth Col; The video format Cod includes a video encoding format and a video encapsulation format. The video encoding format includes one of H.264, H.265, VP9, and AV1. The video encapsulation format includes one of MP4, MKV, MOV, AVI, and FLV. S13, obtaining the number of remote sensing satellite image data to be downloaded based on the obtained duration Tim′ video parameter of each remote sensing satellite image video data to be downloaded. The method for obtaining the number of remote sensing satellite image data to be downloaded based on the obtained duration Tim′ video parameter of each remote sensing satellite image video data to be downloaded is: Qua′=Tim′*FPS, Wherein, Qua′ is the number of each remote sensing satellite image data to be downloaded; Tim′ is the duration Tim′ of each remote sensing satellite image video data to be downloaded; FPS is the frame rate FPS of remote sensing satellite image video data; The resolution Res of each remote sensing satellite image data to be downloaded is consistent with the resolution Res video parameter of each remote sensing satellite image video data to be downloaded; through the above method, the duration Tim′ of each remote sensing satellite image video data to be downloaded can be converted into Qua′ remote sensing satellite images to be downloaded; S14, calculating the image 0-1 digital string of the remote sensing satellite image video data to be downloaded in step S11; S15, cyclically calculating the image 0-1 digital string and the pixel values of each channel in each remote sensing satellite image to be downloaded, to obtain a new remote sensing satellite image to be downloaded: Among them, R(μ,ν) New is the new value of the νth pixel corresponding to the μth pixel in the red channel from left to right and from top to bottom; R(μ,ν) is the ν-th bit value in the pixel value corresponding to the μ-th pixel in the red color channel from left to right and from top to bottom. Here, the pixel value of each pixel is a binary pixel value, and the number of bits of the pixel value is equal to the color depth Col. Imstr 0-1 New (γ) is the value of the γth bit in the new image 0-1 digital string; Will Connect the image 0-1 digital string front and back to get a new image 0-1 digital string; Among them, G(μ,ν) New is the new value of the nth pixel corresponding to the μth pixel from left to right and from top to bottom in the green color channel; G(μ,ν) is the ν-th pixel value corresponding to the μ-th pixel in the green color channel from left to right and from top to bottom; Imstr 0-1 New (γ) is the value of the γth bit in the new image 0-1 digital string; The relationship between γ, μ and ν is: Among them, γ is the position of the new image 0-1 digital string from left to right; μ is the pixel position from left to right and from top to bottom in the red / green / blue color channel; μ = 1, 2, 3, ..., O K *O G ; is the color depth Col; ν is the position of the pixel value from left to right, Among them, B(μ,ν) New is the new value of the νth pixel corresponding to the μth pixel in the blue color channel from left to right and from top to bottom; B(μ,ν) is the ν-th pixel value corresponding to the μ-th pixel in the blue color channel from left to right and from top to bottom; Imstr 0-1 New (γ) is the value of the γth bit in the new image 0-1 digital string; Through the operation of step S15, a new remote sensing satellite image to be downloaded can be obtained at this time. All Qua' remote sensing satellite images to be downloaded are subjected to the operation of step S15 to obtain all new Qua' remote sensing satellite images to be downloaded; S16, all the new remote sensing satellite images to be downloaded are combined into a new remote sensing satellite image video data to be downloaded in the order in which they were played. By operating in step S16, a new set of remote sensing satellite image video data to be downloaded can be obtained; U0 new sets of remote sensing satellite image video data to be downloaded can be obtained by simply operating in steps S15 and S16; S17, transmit the image 0-1 digital string and U0 copies of new remote sensing satellite image video data to be downloaded to the cultivated land planting status detection platform.

4. The method for detecting cultivated land planting status based on multi-source remote sensing images according to claim 1, characterized in that: Step S3 includes the following steps: S31, dividing the remote sensing satellite image data into K0 pixel groups of uniform size. The number of groups is calculated as follows: K0 represents the number of pixel groups into which the remote sensing satellite image data is divided; O K is the horizontal pixel number of remote sensing satellite image data; O G is the number of vertical pixels of remote sensing satellite image data; O K A factor of , and That is O K Can be divisibility; O G A factor of , and That is O G Can be divisibility; S32, performing pixel group processing on each pixel group using an image pixel group processing algorithm; S33, determining the background judgment range of the remote sensing satellite processed image; S34, obtaining the grayscale pixel value of each pixel in the remote sensing satellite processed image, and determining the pixel area where the pixels belong to the background judgment range are located as the background area of the remote sensing satellite processed image; S35, performing fuzzy processing on the background area of the remote sensing satellite processed image to obtain a remote sensing satellite image processed image.

5. The method for detecting cultivated land planting status based on multi-source remote sensing images according to claim 1, characterized in that: The method for obtaining the cultivated land planting area in step S4 includes the following steps: S41, acquiring historical image data containing rice crops, calculating the grayscale pixel value variance and the grayscale pixel value average of the historical image data containing rice crops, and obtaining a moving frame determination threshold based on the grayscale pixel value variance and the grayscale pixel value average; S42, setting a moving frame, performing moving scanning on the remote sensing satellite image processing image, and calculating the moving determination value of each pixel in the remote sensing satellite image processing image according to the moving frame; S43, determining the cultivated land and non-cultivated land areas in the remote sensing satellite image processing image according to the relationship between the moving frame judgment threshold and the moving determination value: If f0 ≥ d0, f0 is the moving frame judgment threshold, d0 is the moving determination value; then the cultivated land planting area in the remote sensing satellite image processing image; If f0 ≥ d0, the non-cultivated land in the remote sensing satellite image processing image is the planting area.

6. The method for detecting cultivated land planting status based on multi-source remote sensing images according to claim 1, characterized in that: In step S4, the method for determining whether there is crop planting in the cultivated land planting area is: S4-1, obtaining the date on which the remote sensing satellite image video data was taken, and after obtaining the date on which the remote sensing satellite image video data was taken, obtaining the grayscale pixel value range of the rice crop on that date; S4-2, arrange the grayscale pixel values in the cultivated area in ascending order, remove the first and last W pixel values, and average the remaining grayscale pixel values: Among them, p final is the grayscale pixel value of rice crops in the cultivated land area; O K is the horizontal pixel number of remote sensing satellite image data; O G is the number of vertical pixels of remote sensing satellite image data; W is the number of pixels removed at the first position; p i ′ is the grayscale pixel value of each pixel in the cultivated land planting area; p i ″ is the grayscale pixel value of each pixel removed from the front of the cultivated land planting area; p i ″′ is the grayscale pixel value of each pixel that is finally removed in the cultivated land planting area; S4-3, the relationship between the grayscale pixel value of the rice crop in the cultivated land area and the grayscale pixel value range of the rice crop on that date is: If the grayscale pixel value of the rice crop in the cultivated land area is within the grayscale pixel value range of the rice crop on that date, then the rice crop exists on that date; Otherwise, there is no rice crop planted on that date.

7. A computer system, characterized in that: include: processor; a memory for storing processor-executable instructions; Wherein, the processor is configured to implement the method for detecting cultivated land planting status based on multi-source remote sensing images as described in one of claims 1 to 6 when executing the executable instructions.

8. A computer-readable storage medium, characterized in that include: a memory having a computer program stored thereon; A processor is used to execute the program in the memory to implement the method for detecting cultivated land planting status based on multi-source remote sensing images as described in any one of claims 1 to 6.

Images