Patent Application
20190114744
The present invention relates to satellite images, and more specifically, to enhancing observation resolution using continuous learning.
Images obtained with orbiting satellites are used in numerous applications such as mapping, border control, infrastructure and regional planning, agricultural and environmental monitoring, intelligence gathering, and biodiversity conservation, for example. In satellite imagery, spatial resolution refers to the pixel size of an image representing the surface area being measured, and temporal resolution refers to the amount of time between imagery collection periods for a given surface location. Some satellites can obtain images at the same viewing angle every 5 or 15 minutes while others can only revisit the same viewing angle every 5 or more days, for example. A satellite with relatively high temporal resolution (i.e., one that revisits the same surface location relatively more often) can have relatively low spatial resolution compared with a satellite with relatively lower temporal resolution. Enhancing the observation resolution of the satellite with relatively higher temporal resolution but relatively lower spatial resolution using continuous learning would provide higher spatial resolution at higher temporal resolution.
Embodiments of the present invention are directed to a method, system, and computer program product to enhance observation resolution using continuous learning. The method includes obtaining a first spatial resolution image of a surface area from a first satellite, and obtaining a second spatial resolution image of the surface area from a second satellite. The first spatial resolution image and the second spatial resolution image are obtained within a specified time period of each other, the first spatial resolution image has a lower spatial resolution than the second spatial resolution image, and temporal resolution of the first spatial resolution images obtained by the first satellite is higher than temporal resolution of the second spatial resolution images obtained by the second satellite. The method also includes determining a convolution matrix A or training a neural network using the first spatial resolution image and the second spatial resolution image, obtaining additional one or more of the first spatial resolution images prior to obtaining an additional one of the second spatial resolution images, and generating a new image from each of the one or more of the first spatial resolution images using the convolution matrix A or the neural network. The new image corresponding with each of the one or more of the first spatial resolution images has a higher spatial resolution than the one or more of the first spatial resolution images.
As previously noted, satellites that obtain surface images can have different temporal and spatial resolutions. For example, a satellite that obtains images of the same surface area more frequently (i.e., a satellite with higher temporal resolution) can have lower spatial resolution than a satellite that obtains images less frequently. Yet, higher spatial resolution images at a higher temporal resolution can be beneficial in each of the technical areas in which satellite images are used. Prior approaches to obtaining higher resolution images involve interpolation and the application of an optical flow method. These approaches provide solutions that are fixed rather than continuously improved and that are not necessarily applicable to every type of surface area.
Embodiments of the systems and methods detailed herein relate to enhancing observation resolution using continuous learning. Unlike prior solutions, the embodiments facilitate combining the higher temporal resolution obtained with one satellite coupled and the higher spatial resolution obtained with another satellite. Specifically, continuous learning is performed using lower and higher spatial resolution images of the same surface location. The learning is applied to increase the resolution of the lower spatial resolution images that are obtained more frequently than the higher spatial resolution images. Additional data (e.g., information regarding precipitation accumulation in the hour prior to the satellite image being obtained) can be used in the process of obtaining the higher spatial resolution image. This additional data can be obtained by a different source that either of the satellites.
By increasing the spatial resolution of images that are obtained more frequently (i.e., higher temporal resolution images), the analysis that is facilitated by the images is improved. Thus, the one or more relevant technical areas that analyze the images—e.g., mapping, border control, infrastructure and regional planning, agricultural and environmental monitoring, intelligence gathering, and biodiversity conservation—are improved based on the increased resolution.
The system 100 can be implemented as a computer system as shown in the exemplary embodiment of
In machine learning and, in particular, in convolutional neural networks (CNNs), an input image or, more generally, a multidimensional data array or matrix, is convolved with a kernel (i.e., convolutional kernel or convolutional matrix) to obtain an output, which can be an output image. A convolution matrix is generated from the kernel as detailed below. The process of convolution can also be referred to as filtering. According to the present embodiment, as
A·F=G [EQ. 1]
In EQ. 1, the high spatial resolution image 140b-1 is indicated as F, and the low spatial resolution image 140a-1 is indicated as G. Thus, EQ. 1 indicates the convolution matrix A that filters the high spatial resolution image 140b-1 to provide the low spatial resolution image 140a-1. When the next low spatial resolution image 140a-2 is obtained before the next high spatial resolution image 140b-2 is obtained, the convolution matrix A can be used to enhance the spatial resolution of the low spatial resolution image 140a-2. Specifically, a deconvolution of the low spatial resolution image 140a-2 with the convolution matrix A will provide the enhanced spatial resolution image 140c-1, and a deconvolution of the low spatial resolution image 140a-3 with the convolution matrix A will provide the enhanced spatial resolution image 140c-2.
When the next high spatial resolution image 140b-2 is obtained, the convolution matrix A is refined, as indicated in
The determination of the convolution matrix A is explained with reference to the image 140a-1 being a 4×4 pixel image, the image 140b-1 being a 12×12 pixel image. The image 140a-1 is written as a vector of 16 elements (4*4), and the image 140b-1 is written as a vector of 144 elements (12*12). Thus, the convolution matrix A according to EQ. 1 is 144×16 such that each of the 16 rows of 144 elements of the convolution matrix A are convolved (i.e., linear combination is performed) with the 144 elements of the image 140b-1 to provide one of the 16 elements of the image 140a-1. Of the 144 elements of each of the 16 rows of the convolution matrix A, most of the 144 elements are 0. Only a subset of the 144 elements, corresponding with a kernel H of k×k elements, is non-zero.
The location of the kernel H among the other elements differs for each of the 16 rows (i.e., the kernel H slides along the image 140b-1). In the exemplary case, the kernel H is a 3×3 matrix (i.e., k=3). In this case, 16 equations corresponding with the 16 rows of the convolution matrix A are used to determine 9 elements (3*3) of the kernel H. This greater number of equations (16) than unknowns (9) means that EQ. 1 is an overdetermined equation. Solving for kernel H and, thus, convolution matrix A using images 140a-1 and 140b-1 can be accomplished with known techniques such as a least square fit or QR factorization.
Once the convolution matrix A is obtained or refined, a low resolution image such as image 140a-2 can be deconvolved with the convolution matrix A to obtain a higher resolution image such as image 140c-1. Deconvolution according to a mathematical inversion of EQ. 1 leads to an under-determined set of linear equations. That is, the number of unknowns (144 in the exemplary case discussed herein) is greater than the number of equations (16 in the exemplary case). In addition to known mathematical regularization techniques to solve such equations, the additional data (e.g., precipitation information for the surface area corresponding with the spatial area sensed by the satellites 105a, 105b) can impose additional constraints on the solution of the equations. Accordingly, a physically meaningful objective function can be constructed that constrains EQ. 1 in order to obtain a higher spatial resolution image 140c from a lower spatial resolution image 140a. For example, precipitation information can be used to restrict the variance, absolute value, or both of the values of pixels of the higher resolution image 140c.
In an additional or alternate embodiment, multiple convolution matrices A can be determined and bucketized to correspond with the additional data. For example, the variance of pixel values can be lower during rain than in sun. Thus, a convolution matrix A can be determined, refined, and maintained for rainy weather while other convolution matrices A are maintained for other weather conditions.
In additional or alternate embodiments, additional data (e.g., weather-related data) can be used in the process by inputting a three-dimensional matrix to the decoder rather than a two-dimensional image alone. The first two dimensions of the input matrix are unchanged and refer to the geo-spatial dimensions such as, for example, latitude and longitude coordinates defined by the center of the pixels of the lower spatial resolution image 140a G, for example. The third dimension can index the geo-referenced additional data per pixel. For example, when the satellite image (e.g., lower spatial resolution image 140a) has multiple bands of the electromagnetic spectrum (e.g., an image with red, green, and blue information), these bands define the third dimension having three distinct values red, green, and blue. When other data, such as precipitation P or temperature T is available, this data (P and T) can define the third dimension of the input matrix either instead of or in addition to the index values red, green, and blue.
Unlike the process discussed with reference to
The present invention may be a system, a method, and/or a computer program product at any possible technical detail level of integration. The computer program product may include a computer readable storage medium (or media) having computer readable program instructions thereon for causing a processor to carry out aspects of the present invention.
The computer readable storage medium can be a tangible device that can retain and store instructions for use by an instruction execution device. The computer readable storage medium may be, for example, but is not limited to, an electronic storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the foregoing. A non-exhaustive list of more specific examples of the computer readable storage medium includes the following: a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), a static random access memory (SRAM), a portable compact disc read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, a mechanically encoded device such as punch-cards or raised structures in a groove having instructions recorded thereon, and any suitable combination of the foregoing. A computer readable storage medium, as used herein, is not to be construed as being transitory signals per se, such as radio waves or other freely propagating electromagnetic waves, electromagnetic waves propagating through a waveguide or other transmission media (e.g., light pulses passing through a fiber-optic cable), or electrical signals transmitted through a wire.
Computer readable program instructions described herein can be downloaded to respective computing/processing devices from a computer readable storage medium or to an external computer or external storage device via a network, for example, the Internet, a local area network, a wide area network and/or a wireless network. The network may comprise copper transmission cables, optical transmission fibers, wireless transmission, routers, firewalls, switches, gateway computers and/or edge servers. A network adapter card or network interface in each computing/processing device receives computer readable program instructions from the network and forwards the computer readable program instructions for storage in a computer readable storage medium within the respective computing/processing device.
Computer readable program instructions for carrying out operations of the present invention may be assembler instructions, instruction-set-architecture (ISA) instructions, machine instructions, machine dependent instructions, microcode, firmware instructions, state-setting data, configuration data for integrated circuitry, or either source code or object code written in any combination of one or more programming languages, including an object oriented programming language such as Smalltalk, C++, or the like, and procedural programming languages, such as the “C” programming language or similar programming languages. The computer readable program instructions may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the latter scenario, the remote computer may be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet Service Provider). In some embodiments, electronic circuitry including, for example, programmable logic circuitry, field-programmable gate arrays (FPGA), or programmable logic arrays (PLA) may execute the computer readable program instruction by utilizing state information of the computer readable program instructions to personalize the electronic circuitry, in order to perform aspects of the present invention.
Aspects of the present invention are described herein with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and/or block diagrams, and combinations of blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer readable program instructions.
These computer readable program instructions may be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks. These computer readable program instructions may also be stored in a computer readable storage medium that can direct a computer, a programmable data processing apparatus, and/or other devices to function in a particular manner, such that the computer readable storage medium having instructions stored therein comprises an article of manufacture including instructions which implement aspects of the function/act specified in the flowchart and/or block diagram block or blocks.
The computer readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable apparatus or other device to produce a computer implemented process, such that the instructions which execute on the computer, other programmable apparatus, or other device implement the functions/acts specified in the flowchart and/or block diagram block or blocks.
The flowchart and block diagrams in the Figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of instructions, which comprises one or more executable instructions for implementing the specified logical function(s). In some alternative implementations, the functions noted in the blocks may occur out of the order noted in the Figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems that perform the specified functions or acts or carry out combinations of special purpose hardware and computer instructions.
The descriptions of the various embodiments of the present invention have been presented for purposes of illustration, but are not intended to be exhaustive or limited to the embodiments described. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments. The terminology used herein was chosen to best explain the principles of the embodiments, the practical application or technical improvement over technologies found in the marketplace, or to enable others of ordinary skill in the art to understand the embodiments described herein.
