# US 2014/0153785 A1 · Methods and systems for detection and characterization of channel systems > Agents: this data is available as JSON from the patentgrep API and MCP server — see https://patentgrep.com/llms.txt - Number: US 2014/0153785 A1 - Document ID: US20140153785A1 - Status: Published application - Kind: A1 - Application number: 13/706,148 - Filed: 2012-12-05 - Published: 2014-06-05 - Earliest claimed: 2012-12-05 - Assignees: Chevron U.S.A. Inc. - Inventors: Eitan Shelef; Jesse Lomask; Andrea Fildani - Web page: https://patentgrep.com/patent/US20140153785A1 ## CPC - G06V 10/443: {by matching or filtering} - G06T 7/12: Edge-based segmentation - G06T 7/168: involving transform domain methods - G06T 2207/10032: Satellite or aerial image; Remote sensing - G06T 2207/20064: Wavelet transform [DWT] - G06T 2207/30181: Earth observation - G06V: IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING - G06T: IMAGE DATA PROCESSING OR GENERATION, IN GENERAL ## Abstract A computer implemented method for detecting a channel system comprises importing channel data, wherein the channel data includes intensity measurements associated with locations in the channel system. The method further comprises calculating, by a processor, directional first derivative data of the intensity measurements; selecting a plurality of localized test wavelets; calculating, by the processor, a plurality of fit-measures, wherein the plurality of fit-measures indicate correlations between the directional first derivatives and the plurality of localized test wavelets; and determining a plurality of selected wavelets from the plurality of localized test wavelets based on the plurality of fit-measures, wherein the plurality of selected wavelets model the channel system. ## Claims ### Claim 1 A computer implemented method for detecting a channel system, the method comprising: importing channel data, wherein the channel data includes intensity measurements associated with locations in the channel system; calculating, by a processor, directional first derivative data of the intensity measurements; selecting a plurality of localized test wavelets; calculating, by the processor, a plurality of fit-measures, wherein the plurality of fit-measures indicate correlations between the directional first derivatives and the plurality of localized test wavelets; and determining a plurality of selected wavelets from the plurality of localized test wavelets based on the plurality of fit-measures, wherein the plurality of selected wavelets model the channel system. ### Claim 2 (depends on 1) The method of claim 1, wherein each of the plurality of selected wavelets is centered at a corresponding point of the channel system and has an orientation that is perpendicular to a direction of a corresponding directional first derivative at the corresponding point. ### Claim 3 (depends on 1) The method of claim 1, wherein the plurality of localized test wavelets include a plurality of wavelets having different orientations. ### Claim 4 (depends on 1) The method of claim 1, wherein the plurality of localized test wavelets include a plurality of wavelets having different widths. ### Claim 5 (depends on 1) The method of claim 1, wherein the channel data correspond to a plurality of grid points and the plurality of selected wavelets are assigned to the plurality of grid points, the method further comprising calculating a plurality of signal to noise values for each of the plurality of grid points, wherein a signal to noise value for one grid point is an increasing function of the calculated fit-measure of the selected wavelet assigned to the one grid point. ### Claim 6 (depends on 1) The method of claim 1, wherein calculating the plurality of fit-measures include calculating a plurality of convolutions. ### Claim 7 (depends on 6) The method of claim 6, wherein calculating the plurality of convolutions is performed in a Fourier space. ### Claim 8 (depends on 1) The method of claim 1, wherein calculating a plurality of fit-measures includes: calculating a first fit-measure between the directional first derivatives and a first test wavelet localized at one location in the channel system; and calculating a second fit-measure between the directional first derivatives and a second test wavelet that is localized at the one location and is different from the first wavelet; and wherein determining the plurality of selected wavelets includes selecting, as a selected wavelet, one of the first test wavelet and the second test wavelet for which the corresponding calculated fit-measure is larger. ### Claim 9 A computer implemented method for analyzing data for a channel system, wherein the data correspond to a plurality of grid points that represent a surface containing the channel system, the method comprising: storing a plurality of fit-measures assigned to the plurality of grid points, wherein the plurality of fit-measures indicate correlations between a plurality of wavelets and directional first derivates of the data for the channel system; determining a plurality of links connecting a subset of the plurality of grid points; assigning to a subset of the plurality of links a plurality of link costs, wherein the plurality of link costs are based on the plurality of fit-measures; selecting a focal grid point; and determining, by a processor, a plurality of channel paths each connecting the focal grid point to one of a subset of the plurality of grid points, wherein the plurality of channel paths is determined based on the plurality of link costs. ### Claim 10 (depends on 9) The method of claim 9, wherein determining the plurality of channel paths from the focal grid point to one grid point of the plurality of grid points includes: determining a plurality of paths from the focal grid point to the one grid point; assigning to each of the plurality of paths a path cost, wherein the path cost for one path is an increasing function of a subset of the link costs for links covered by the one path; based on the path costs, selecting a path among the plurality of paths as the determined channel path from the focal grid point to the one grid point. ### Claim 11 (depends on 10) The method of claim 10, wherein the selected path has a minimum path cost among the plurality of paths. ### Claim 12 (depends on 9) The method of claim 9 further comprising: assigning a plurality of wavelet parameters to the plurality of grid points; and increasing a link cost for one link of the plurality of links based on a difference between values of wavelet parameters for the grid points connected by the one link. ### Claim 13 A system for detecting a channel system, the system comprising: a detector module configured to: import channel data, wherein the channel data includes intensity measurements associated with locations in the channel system; calculate directional first derivative data of the intensity measurements; selecting a plurality of localized test wavelets; calculate a plurality of fit-measures, wherein the plurality of fit-measures indicate correlations between the directional first derivatives and the plurality of localized test wavelets; and determine a plurality of selected wavelets from the plurality of localized test wavelets based on the plurality of fit-measures, wherein the plurality of selected wavelets model the channel system. ### Claim 14 (depends on 13) The system of claim 13, further comprising a connector module for analyzing data for the channel system, wherein the data correspond to a plurality of grid points that represent a surface containing the channel system, wherein the connector module is configured to: store a plurality of fit-measures assigned to the plurality of grid points, wherein the plurality of fit-measures indicate correlations between a plurality of wavelets and directional first derivates of the data for the channel system; determine a plurality of links connecting a subset of the plurality of grid points; assign to a subset of the plurality of links a plurality of link costs, wherein the plurality of link costs are based on the plurality of fit-measures; select a focal grid point; and determine a plurality of channel paths each connecting the focal grid point to one of a subset of the plurality of grid points, wherein the plurality of channel paths is determined based on the plurality of link costs. ### Claim 15 (depends on 13) The method of claim 13, wherein one wavelet is assigned to one grid point of the plurality of grid points and wherein one wavelet parameter assigned to the one grid point includes a width of the one wavelet, a standard deviation of the one wavelet, or an orientation of the one wavelet. ### Claim 16 (depends on 13) The method of claim 13 further comprising setting a constraint on a subset of the plurality of grid points or on a subset of the plurality of links based on the plurality of wavelet parameters, wherein one or more of the plurality of channel paths satisfies the constraint. ### Claim 17 (depends on 16) The method of claim 16, wherein the constraint includes blocking the subset from being included in a subset of the plurality of channel paths. ### Claim 18 (depends on 17) The method of claim 17, wherein a link connecting a first grid point and a second grid point neighboring the first grid point is blocked from being included in a channel path if a difference between values of one wavelet parameter for the first and second grid points exceeds a threshold difference value. ### Claim 19 (depends on 9) The method of claim 9, wherein determining the plurality of channel paths includes removing a subset of the plurality of links from being included in the plurality of channel paths, wherein for each link in the subset of the plurality of links the link cost exceeds a link cost maximum. ### Claim 20 (depends on 9) The method of claim 9, wherein determining one of the plurality of channel paths from the focal point to one grid point of the plurality of grid points includes determining a lowest cost path from the focal point to the one grid point using a Dijkstra algorithm. ### Claim 21 (depends on 9) The method of claim 9, further comprising: assigning to a subset of the plurality of links a plurality of focal link costs, wherein the plurality of focal link costs are below a threshold cost value; and determining, by a processor, the plurality of channel paths such that a subset of the plurality of channel paths overlap one or more of the low-cost subset of the plurality of links. ### Claim 22 (depends on 9) The method of claim 9, further comprising assigning statistical measures to the channel system. ### Claim 23 (depends on 22) The method of claim 22, wherein the statistical measures comprise channel orders, channel widths, channel lengths, or number of first order basins. ### Claim 24 (depends on 22) The method of claim 22, wherein assigning statistical measures comprises assigning stream directions to a subset of the plurality of channel paths and drainage areas to a subset of the plurality of grid points, wherein a drainage area assigned to one grid point is proportional to a number of grid points that are located on the same channel path as the one grid point and are upstream on that channel path compared to the one grid point. ### Claim 25 (depends on 22) The method of claim 22, wherein a test wavelet is based on a diffusion function and wherein the statistical measures include a morphologic age of a wavelet included in the channel system. ### Claim 26 (depends on 9) The method of claim 9, wherein the plurality of fit-measures are a plurality of signal to noise ratios. ### Claim 27 A system for analyzing data for a channel system, wherein the data correspond to a plurality of grid points that represent a surface containing the channel system, the system comprising a connector module, wherein the connector module is configured to: store a plurality of fit-measures assigned to the plurality of grid points, wherein the plurality of fit-measures indicate correlations between a plurality of wavelets and directional first derivates of the data for the channel system; determine a plurality of links connecting a subset of the plurality of grid points; assign to a subset of the plurality of links a plurality of link costs, wherein the plurality of link costs are based on the plurality of fit-measures; select a focal grid point; and determine a plurality of channel paths each connecting the focal grid point to one of a subset of the plurality of grid points, wherein the plurality of channel paths is determined based on the plurality of link costs. ### Claim 28 A non-transitory computer-readable medium storing instructions that, when executed by one or more computers, causes the one or more computers to perform a method for detecting a channel system, the method comprising: importing channel data, wherein the channel data includes intensity measurements associated with locations in the channel system; calculating, by a processor, directional first derivative data of the intensity measurements; selecting a plurality of localized test wavelets; calculating, by the processor, a plurality of fit-measures, wherein the plurality of fit-measures indicate correlations between the directional first derivatives and the plurality of localized test wavelets; and determining a plurality of selected wavelets from the plurality of localized test wavelets based on the plurality of fit-measures, wherein the plurality of selected wavelets model the channel system. ### Claim 29 A non-transitory computer-readable medium storing instructions that, when executed by one or more computers, causes the one or more computers to perform a method for analyzing data for a channel system, wherein the data correspond to a plurality of grid points that represent a surface containing the channel system, the method comprising: storing a plurality of fit-measures assigned to the plurality of grid points, wherein the plurality of fit-measures indicate correlations between a plurality of wavelets and directional first derivates of the data for the channel system; determining a plurality of links connecting a subset of the plurality of grid points; assigning to a subset of the plurality of links a plurality of link costs, wherein the plurality of link costs are based on the plurality of fit-measures; selecting a focal grid point; and determining, by a processor, a plurality of channel paths each connecting the focal grid point to one of a subset of the plurality of grid points, wherein the plurality of channel paths is determined based on the plurality of link costs.