This application claims the benefit of Korean Application No. 10-2015-0131899, filed on Sep. 17, 2015, the disclosures of which are incorporated by reference into the present application.
1. Field of the Invention
The present invention relates generally to a method and system for predicting a chlorophyll-a concentration and, more particularly, to a method and apparatus for predicting a chlorophyll-a concentration in a river using satellite image data and a nonlinear RANSAC method.
2. Description of the Related Art
Due to the effects of recent global warming, pollution levels in waters including rivers, reservoirs, lakes, etc., have been increasing every year. In particular, as water temperatures increase, eutrophication, caused primarily by an increase in chlorophyll-a, disturbs aquatic ecosystems as well as pollutes drinking water resources.
Current water quality, including measurements of chlorophyll-a concentrations in rivers, is managed in a manner whereby the Korea Environment Corporation and each local government install measurement devices at main water areas to remotely measure water quality and then deliver the measurement information to measurement stations for monitoring pollution levels. Such a method is advantageous in that accurate measurement values may be obtained from watersheds, but has limitations in that it is not plausible to install measurement devices at all positions of all rivers in the country. Further, frequent malfunctioning and maintenance of the measurement devices is costly. In order to make up for such limitations, recently, development has been ongoing into a method that estimates chlorophyll-a concentrations by extracting a recursive function based on satellite image data and inputting satellite image data values to the extracted recursive function. However, the recursive function extracted in an existing manner has a limitation of being largely influenced by outlier data. Regarding measurement of water quality of rivers, Korean Patent Laid-Open Application publication No. 10-2012-0098158, published on Sep. 5, 2012, discloses a system for measuring water pollution using a water quality measurement sensor and a gas concentration measurement sensor.
Accordingly, the present invention has been made keeping in mind the above problems occurring in the prior art, and an object of the present invention is to provide a method and apparatus for predicting a chlorophyll-a concentration in a river using satellite image data and a nonlinear random sample consensus (RANSAC method) method, which is capable of effectively predicting a wide range of chlorophyll-a concentrations in rivers by inputting satellite image data collected at a specific basin to a second-order function, which is extracted by applying, to the nonlinear RANSAC method, chlorophyll-a concentrations, which were actually measured at a measurement station at the specific basin, and image data, which was obtained from a satellite on the same date and at the same site as when and where the chlorophyll-a concentrations were measured at the corresponding measurement station.
Another object of the present invention is to provide a method and apparatus for predicting a chlorophyll-a concentration in a river using satellite image data and a nonlinear RANSAC method, which is capable of minimizing an effect of outlier data and predicting the chlorophyll-a concentration more accurately.
In order to accomplish the above object, the present invention provides a method for predicting a chlorophyll-a concentration in a river using satellite image data and a nonlinear RANSAC method. The method includes: (1) receiving chlorophyll-a concentration, which was actually measured at a measurement station at a prescribed basin and an image from a satellite on a same date and at a same site as when and where the chlorophyll-a concentration was measured at the measurement station; (2) correcting a distortion of the image received from the satellite; (3) applying the chlorophyll-a concentration actually measured at the measurement station and data of the corrected satellite image to a nonlinear RANSAC method to extract a second-order function; and (4) inputting corrected data of a satellite image collected at the prescribed basin to the extracted second-order function to predict a chlorophyll-a concentration at the prescribed basin.
In step (1), an operational land imager sensor image provided from LANDsat 8 satellite may be received.
Step (2) may include: (2-1) converting a value of a digital number in the image received from the satellite into a top of atmosphere (TOA) reflectance ratio by using the following Equation:
ρλ′=MρQcal+Ap
where ρλ′ denotes a TOA reflectance ratio to which a sun angle is not considered, Mρ denotes a multiplicative rescaling factor of a metadata file, Ap denotes an additive rescaling factor of a metadata file, and cal denotes a pixel of the image.
Step (2) may further include: (2-2) extracting a TOA reflectance ratio, for which a sun angle is considered, from the TOA reflectance ratio converted in step (2-1) by using the following Equation:
where ρλ denotes the TOA reflectance ratio for which the sun angle is considered, θSE denotes a sun angle at a measurement area, and θSZ denotes a solar zenith angle.
Step (3) may include (3-1) generating sample data in which the corrected data of the satellite image, which was received on the same date and at the same site as when and where the chlorophyll-a concentration was measured at the measurement station, is taken as x, and a value of the chlorophyll-a concentration, which was actually measured at the measurement station, is taken as y.
Step (3) may further include: (3-2) randomly extracting three sample data points of P1(x1,y1), P2(x2,y2), and P3(x3,y3) from among the sample data generated in step (3-1).
Step (3) may further include: (3-3) substituting the three points extracted in step (3-2) to the following Equation to obtain a second-order function crossing the three points:
y=ax
2
+bx+c.
Step (3) may further include (3-4) obtaining a number N of data of which distances to the second-order function obtained in step (3-3) are equal to or greater than a prescribed value.
Step (3) may further include: (3-5) storing N as a new maximum value N_max of N when N obtained in step (3-4) is greater than a stored maximum value N_max of N and storing a corresponding second-order function.
Step (3) may further include: (3-6) determining, as a final second-order function, a second-order function stored after repeating steps (3-2) to (3-5) a prescribed number of times.
In order to accomplish the above object, the present invention also provides a system for predicting a chlorophyll-a concentration in a river using satellite image data and a nonlinear RANSAC method. The system includes: an input module receiving a chlorophyll-a concentration at a predetermined basin, which was actually measured at a measurement station and an image from a satellite on a same date and at a same site as when and where the chlorophyll-a concentration was measured at the measurement station from a satellite; a correcting module correcting a distortion of an image received from the satellite; a function extracting module applying the chlorophyll-a concentration actually measured at the measurement station and data of the corrected satellite image to a nonlinear RANSAC method to extract a second-order function; and a chlorophyll-a concentration predicting module inputting corrected data of a satellite image collected at the prescribed basin to the extracted second-order function to predict a chlorophyll-a concentration at the prescribed basin.
The input module may receive an operational land imager sensor image provided from LANDsat 8 satellite.
The correcting module may include: a TOA converting unit converting a value of a digital number in the image received from the satellite into a TOA reflectance ratio by using the following Equation:
ρλ′=MρQcal+Ap
where ρλ′ denotes a TOA reflectance ratio for which a sun angle is not considered, Mρ denotes a multiplicative rescaling factor of a metadata file, Ap denotes an additive rescaling factor of a metadata file, and Qcal denotes a pixel of the image.
The TOA converting unit may extract a TOA reflectance ratio, for which a sun angle is considered, from the converted TOA reflectance ratio by using the following Equation:
where ρλ denotes the TOA reflectance ratio for which the sun angle is considered, θSE denotes a sun angle at a measurement area, and θSZ denotes a solar zenith angle.
The function extracting module may generate sample data in which the corrected data of the satellite image, which was received on the same date and at the same site as when and where the chlorophyll-a concentration was measured at the measurement station, is taken as x, and a value of the chlorophyll-a concentration, which was actually measured at the measurement station, is taken as y.
The function extracting module may randomly extract three sample data points of P1(x1,y1), P2(x2,y2), and P3(x3, y3) from among the generated data.
The function extracting module may substitute the extracted three points to the following Equation to obtain a second-order function crossing the three points:
y=ax
2
+bx+c.
The function extracting module may obtain a number N of data of which distances to the obtained second-order function are equal to or greater than a prescribed value from the sample data generated therein.
The function extracting module may store N as a new maximum value N_max of N when N is greater than a stored maximum value N_max of N and stores a corresponding second-order function.
The function extracting module may randomly extract three points from among the sample data to obtain a second-order function crossing the three points, obtains a number N of sample data of which distances to the second-order function are equal to or smaller than a prescribed value, and when N is greater than a maximum value N_max of N, stores N as a new maximum N_max of N, and stores, as a final second-order function, a second-order function stored after the process for storing a corresponding second-order function is repeated a prescribed number of times.
The above and other objects, features and advantages of the present invention will be more clearly understood from the following detailed description taken in conjunction with the accompanying drawings, in which:
Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings so that the present invention can be easily realized by those skilled in the art. Further, in the description of the present invention, when it is determined that the detailed description of the related functions and constructions would obscure the gist of the present invention, the description thereof will be omitted. In addition, like reference numerals refer to elements having like functions and operations throughout the drawings.
In addition, if certain parts are described as being connected to other parts, they are not only directly connected to the other parts, but also indirectly connected to the other parts with any other device intervened therebetween. In addition, when an element is referred to as “comprising” or “including” a component, it does not preclude another component but may further include the other component unless the context clearly indicates otherwise.
Hereinafter, each step of a method for predicting a chlorophyll-a concentration in a river using satellite image data and the nonlinear RANSAC method proposed in the present invention will be described in detail with reference to the accompanying drawings.
In step S210, the value of the digital number in the image, which is received from the satellite, may be converted into a TOA reflectance ratio by using the following Equation (1).
ρλ′=MρQcal+Ap (1)
where ρλ′ denotes a TOA reflectance ratio for which a sun angle is not considered, Mρ denotes a multiplicative rescaling factor of a metadata file, Ap denotes an additive rescaling factor of a metadata file, and Qcal denotes a pixel of the image.
In step S220, the TOA reflectance ratio, for which the sun angle is considered, may be extracted from the TOA reflectance ratio converted in step S210 by using the following the following Equation (2).
where ρλ denotes the TOA reflectance ratio for which the sun angle is considered, θSE denotes a sun angle at a measurement area, and θSZ denotes a solar zenith angle.
In step S310, the sample data may be generated by taking the data values, which are obtained through step S200 by correcting the satellite images received on the same date and at the same site as when and where the chlorophyll-a concentrations were measured at the measurement station, as x, and by taking values of chlorophyll-a concentrations, which were actually measured at the corresponding measurement station, as y. Here, according to an embodiment, the sample data may be generated on the basis of data measured for several days, not for one day.
In step S320, three data points P1(x1,y1) P2(x2,y2), and P3(x3,y3) may be randomly extracted from among the sample data generated in step S310.
In operation S330, the three points extracted in step S320 are substituted to the following Equation (3) to obtain a second-order function, which crosses the three points, as illustrated in
y=ax
2
+bx+c (3)
The nonlinear RANSAC method used herein may be used when data is nonlinearly configured. Typically, the RANSAC method is a parameter estimation method for selecting a model of which the number of supporting data is a greatest, and herein, the supporting data refers to data of which distances to an estimated model are equal to or smaller than the prescribed value. For example, in the RANSAC method, when the distance ri=|yi−f(xi)| between data (xi, yi) and a model f(x) is equal to or smaller than the prescribed value, a corresponding data may be considered as data for supporting the model. Here, when the prescribed value is too great, models are not discriminated from each other, and when too small, the RANSAC method may be unstable.
In step S340, from the sample data generated in step S310, the number N of data, of which distances to the second-order function obtained in step S330 are equal to or smaller than the prescribed value, may be obtained.
In step S350, when the N obtained in step S340 is greater than the stored maximum value N_max of N, the N may be stored as the new maximum value N_max of N and the corresponding second-order function may be stored. Here, the maximum value N_max of N may be initialized as 0 at first. A pre-stored maximum value N_max of N is compared with the N obtained in step S340, and when the N obtained in step S340 is greater, the N obtained in step S340 is stored as the maximum value N_max and the corresponding second-order function may be stored. On the other hand, when the N obtained in step S340 is smaller the pre-stored maximum value N_max of N, the pre-stored N and the second-order function corresponding thereto may be stored as they are.
As illustrated in
In the system 10 for predicting a chlorophyll-a concentration in a river using satellite image data and a nonlinear RANSAC method, the input module 100 receives chlorophyll-a concentrations, which were actually measured at a measurement station, and images from an OLE sensor provided at Landsat 8 satellite on the same date and at the same site as when and where the chlorophyll-a concentrations were measured at the corresponding measurement station. Next, the correcting module 200 corrects, through the TOA converting unit 200, distortions in the images received from the input module 100. In addition, the function extracting module 300 extracts a second-order function by applying, to the nonlinear RANSAC method, chlorophyll-a concentrations, which were actually measured at a measurement station, and data of the satellite image corrected through the correcting module 200. Finally, the chlorophyll-a concentration predicting module 400 inputs the corrected data of the satellite images, which are collected at the specific basin, to the second-order function extracted through the function extracting module 300 to predict the chlorophyll-a concentrations at the corresponding basin.
According to a method and apparatus for predicting a chlorophyll-a concentration in a river using satellite image data and a nonlinear RANSAC method proposed in the present invention, a wide range of chlorophyll-a concentrations may be effectively predicted in rivers by inputting satellite image data collected at a specific basin to a second-order function, which is extracted by applying, to a RANSAC method, chlorophyll-a concentrations at the specific basin, which were actually measured at a measurement station, and data of images that were obtained from the satellite on the same date and at the same site as when and where the chlorophyll-a concentrations were measured at the corresponding measurement station.
In addition, the present invention enables an effect of outlier data to be minimized and chlorophyll-a concentrations to be more accurately predicted by using a nonlinear RANSAC method.
Although the preferred embodiments of the present invention have been disclosed for illustrative purposes, those skilled in the art will appreciate that various modifications, additions and substitutions are possible, without departing from the scope and spirit of the invention as disclosed in the accompanying claims.
Number | Date | Country | Kind |
---|---|---|---|
10-2015-0131899 | Sep 2015 | KR | national |