# Patent

## Title
Electromagnetic (EM) power density and field characterization technique

## Abstract

An apparatus and method for characterization of a directed beam of electromagnetic radiation is provided. An exemplary embodiment of the invention can include an apparatus and measuring technique method which uses a model for blackbody radiation that includes consideration all the degrees of freedom due to translation, vibration, and rotation of molecules or atoms that make up the absorber and a heat transfer term which averages the behavior of all the atoms of the material as a function of temperature. This apparatus and method provides an advantage of increased accuracy, substantial reductions of time required for processing, simplification of measuring processes, and reduction required equipment.

## Background
CROSS-REFERENCE TO RELATED APPLICATION
This application is a continuation of and claims priority to U.S. patent application Ser. No. 13/967,823, filed Aug. 15, 2013, entitled “ELECTROMAGNETIC (EM) POWER DENSITY AND FIELD CHARACTERIZATION TECHNIQUE”, which claims priority to U.S. patent application Ser. No. 12/983,230, filed Dec. 31, 2010, entitled “RF POWER DENSITY AND FIELD CHARACTERIZATION TECHNIQUE”, now U.S. Pat. No. 8,543,357, which claims priority to U.S. Provisional Patent Application Ser. No. 61/362,823, entitled “ELECTROMAGNETIC (EM) POWER DENSITY AND FIELD CHARACTERIZATION TECHNIQUE”, having a filing date of Jul. 30, 2010, the disclosures of which are expressly incorporated by reference herein.


STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT
The invention described herein includes contributions by one or more employees of the Department of the Navy made in performance of official duties and may be manufactured, used and licensed by or for the United States Government for any governmental purpose without payment of any royalties thereon. This invention (Navy Case 200,347) is assigned to the United States Government and is available for licensing for commercial purposes. Licensing and technical inquiries may be directed to the Technology Transfer Office, Naval Surface Warfare Center Crane, email: Cran_CTO@navy.mil.


BACKGROUND
The present invention relates to an apparatus, system and method for characterizing and measuring a directed electromagnetic (EM) field characterization, e.g., a radio frequency (RF) beam. Previous approaches to RF beam measurement employed a model of heat transfer from one macroscopic body to another and blackbody radiation of a particle with only one degree of freedom.
Other techniques can use a horn antenna, attenuators, and power meter to collect data at different points along a plane perpendicular to the beam. This technique is much more time consuming the apparatus is more complicated and in many cases must be automated with translational stages in order to make an accurate measurement.
Additional features and advantages of the present invention will become apparent to those skilled in the art upon consideration of the following detailed description of the illustrative embodiment exemplifying the best mode of carrying out the invention as presently perceived.
SUMMARY OF THE INVENTION
According to an illustrative embodiment of the present disclosure, an embodiment of the invention can be used to characterize a directed beam of EM radiation. In particular, an embodiment of the invention can include an apparatus and measuring technique method which uses a model for blackbody radiation which includes consideration of all the degrees of freedom due to translation, vibration, and rotation of molecules or atoms that make up the absorber and a heat transfer term which averages the behavior of all the atoms of the material as a function of temperature. This apparatus and method provides an advantage of increased accuracy, substantial reductions of time required for processing, simplification of measuring processes, and reduction required equipment.



BRIEF DESCRIPTION OF THE DRAWINGS
The detailed description of the drawings particularly refers to the accompanying figures in which:
FIG. 1 shows a description of an exemplary system in accordance with one embodiment of the invention.
FIG. 2 shows an exemplary functional description of software modules describing one way of organizing processing sequences or software in accordance with one embodiment of the invention.
FIG. 3 shows an exemplary graphical user interface for controlling processing, inputs, and producing outputs in accordance with one embodiment of the invention.
FIG. 4 shows an exemplary description of data arrays used with software used for processing in accordance with one embodiment of the invention.
FIG. 5A shows an exemplary output associated with one type of EM source produced and displayed in accordance with one embodiment of the invention.
FIG. 5B shows an exemplary output associated with another type of EM source produced and displayed in accordance with one embodiment of the invention.
FIG. 6 shows an exemplary degrees of freedom table having data which is used in association with blackbody calculations for use in accordance with one embodiment of the invention.
FIG. 7A shows an exemplary method used in producing outputs in accordance with one embodiment of the invention.
FIG. 7B is a continuation of the FIG. 7A description of an exemplary process in accordance with one embodiment of the invention.
FIG. 8 shows an exemplary data plot of temperature data versus frame number of data collected from an EM source using one exemplary embodiment of the invention.



DETAILED DESCRIPTION OF THE DRAWINGS
The embodiments of the invention described herein are not intended to be exhaustive or to limit the invention to precise forms disclosed. Rather, the embodiments selected for description have been chosen to enable one skilled in the art to practice the invention.
Referring to FIG. 1, an exemplary embodiment of an EM field measuring system is shown which includes an EM source 1 (e.g., an infrared source), an EM absorption material 3 (e.g., a sheet of carbon loaded Kapton®), and a temperature sensor array 5 (e.g., an infrared camera) to measure the temperature of the EM absorption material 3. The temperature sensor array 5 is coupled with a processing system 7 which has software (further discussed below) for processing data received and outputting data representations characterizing an EM field which was directed on the EM absorption material 3. These data representations can be displayed on a display system 9 or transmitted to a network via a network hub 11 for further processing, output, or other action. Software on the processing system 7 generally can be operable by a user via a graphical interface which permits interaction between the user and system via an interface system such as a keyboard 13 and a mouse 15 however another interface can be used such as a touch screen (not pictured). Optionally, a gas or air circulation unit 17 can be used to keep the EM absorption material 3 near or at room temperature for the purposes of heat transfer from the absorber to ambient air or atmosphere surrounding the EM absorption material 3.
A variety of types of suitable coolant or temperature control system can be used to provide a temperature control action/effect for the invention particularly where the system is in an environment where environmental factors or heat transfer variables can impact the overall results obtained by the temperature sensor array 5. In other words, temperature sensor array data can be adversely impacted where environmental factors alter the heat transfer characteristics of an EM absorption material 3 e.g., being outside with variable winds, sun, etc which create changes in heat transfer from factors besides a particular EM source. Other types of heat transfer control components/systems can also be used, including multiple cooling mechanisms, heat sinks or other heat transfer or control systems. Also, increased accuracy can be obtained by using such a temperature control mechanism for an EM absorption material in general, including indoor environments, as “hot pockets” or uneven heat transfer areas can form around the surface of a given EM absorption material preventing heat transfer to occur as convection forces alone frequently are insufficient to ensure uniform heat transfer across the applicable surface of a given EM absorption material.
Referring to FIG. 1, EM source 1 is energized and is positioned so an EM field emitted by EM source 1 is incident onto the EM absorption material 3. A temperature sensor array, e.g., an infrared (IR) camera frame capture, which shows the temperature of the absorption material, e.g., Kapton®, before the EM source 1 is turned on, is used along with a frame of when the absorption material 3 has reached thermal equilibrium and the time it took to reach thermal equilibrium. It should be noted that another embodiment can use another temperature sensor array for the EM Absorber besides an IR camera such as an array of thermocouples which could produce simultaneous temperature readings of a region of interest of the EM absorber over a period of time. A variety of different EM absorbers can be used with the invention as well.
Referring to FIG. 2, a functional diagram of processing software associated with an exemplary embodiment of the invention is shown. An IR Camera Output Processing (IRCOP) Module 21 creates a Temperature Data per Pixel (TDPP) array. IRCOP produces the TDPP array based on a series of IR camera image frame captures of an EM Absorber before, during, and, optionally, after exposure to an activated EM source (e.g., directed energy beam or an EM field produced by an antenna or wave guide). In particular, the TDPP array stores a series of IR image frames or snapshots of the EM absorber taken by the temperature sensor array 5 (e.g., IR camera) over time. Each IR image frame comprises temperature data associated with pixels in the IR camera 5 used to capture each IR image frame. The TDPP array includes a file name, a series of IR image frame captures including a temperature data value for each pixel, and time of capture for the temperature data value for each frame/pixel. This exemplary array can be stored as a compressed audio video interleave (AVI) file.
A graphical user interface (GUI) Module 23 is provided for controlling operations in the processing software and providing outputs. The GUI Module 23 creates GUIs such as are shown in FIG. 3.
A Temperature Differential Analysis (TDA) Module 25 includes a section for determining an initial state temperature (e.g., a “cold” frame with a lowest overall temperature value) associated with one of the series of IR image capture frames stored in the TDPP, a final state temperature (e.g., a “hot” frame associated with a highest overall temperature value) associated with one of the series of IR image capture frames acquired over a period of time, and elapsed time between the initial and final state temperatures based on the TDPP array data. For example, the TDA Module 25 can select a “cold” IR image capture frame and a “hot” IR image capture frame based on comparing all frames and finding the frames which have pixels which have lowest overall temperature and the highest overall temperature. Processing in the TDA Module 25 can include data input (e.g., TDPP array), pre-processing (e.g., determination of elapsed time between initial and final states (e.g., “cold” and “hot” IR image capture frames), selection of an initial temperature state IR image capture frame (e.g. “cold”) and a final IR image capture frame (e.g., “hot”)), calculation of differential temperature values per pixel, and temperature differential (TD) array creation. In this example, the initial temperature state is a state prior to EM source activation and final temperature state is a state where an EM Absorber 3 reaches a predetermined temperature such as an approximate peak temperature value. In an embodiment where a pulsed EM (e.g. RF) source is directed onto an EM absorber, the final temperature state can be at any point after the EM source is turned on where a maximum temperature value is reached or another parameter (e.g., time) is met. An exemplary TD array includes an IR image capture frame determined to have the desired initial temperature state (e.g., cold), an IR image capture frame determined to have the desired final temperature state (e.g., “hot”), time values associated with both selected IR image capture frames, and elapsed time associated with each selected IR image capture frame and/or pixels in an array of temperature sensor array 5 e.g., IR Camera.
An exemplary Power Density Calculation (PDC) Module 27 calculates power density associated with all or selected pixels in the selected IR image capture frames (i.e., “cold” and “hot” frames) from a temperature sensor array 5 (e.g., IR Camera) and stores it in a Power Density per Pixel (PDPP) array. The PDC Module 27 receives inputs including inputs from the TDA Module 25 including the TD array 47. Power density in this example is calculated based on at least two temperature sensor array captures selected using temperature data values (e.g., initial IR image capture frame (e.g., “cold) and final IR image capture frame (e.g., “hot”)), elapsed time between at least two temperature sensor array captures, reflection and transmission coefficients of an EM absorber, degrees of freedom data associated with atoms in the EM absorber (e.g., see FIG. 6 table), emissivity of the EM absorber, Stefan-Boltzmann's Constant, specific heat capacity of the EM absorber, density of a medium which is used to cool the EM absorber (e.g., air, heat sink, etc), volume of medium used to cool the EM absorber, and area of interface between the medium used to cool the EM absorber and the EM absorber. An alternate embodiment can include an EM absorber which is sufficiently large enough to act as a heat sink for itself and thereby add to the medium used to cool the EM absorber therefore must be factored into the processing by the PDC Module 27. A general form of a power density calculation performed by the PDC Module 27 is provided in FIG. 7A. Power density calculation performed in this exemplary PDC Module 27 can be affected by selection of different EM absorbers and various factors including cooling of the EM absorber.
An example of processing performed in PDC Module 27 for a specific EM absorber can include a thermal imaging and millimeter wave measurement technique using a carbon loaded polymer sheet called Kapton® and treating it as a diatomic gas. An exemplary calculation for a particular EM absorber, a Kapton® EM absorber, and the exemplary testing configuration shown in FIG. 1, and the software described in FIG. 2 is shown in Equation 1.




  
    
      
        P
        =
        
          
            2
            
              1
              -
              
                
                  
                  Γ
                  
                
                2
              
              -
              
                
                  
                  τ
                  
                
                2
              
            
          
          ⁡
          
            [
            
              
                2
                ⁢
                
                  ɛ
                  ir
                
                ⁢
                
                  
                    σ
                    b
                  
                  ⁡
                  
                    (
                    
                      
                        T
                        0
                        4
                      
                      -
                      
                        T
                        0
                        4
                      
                    
                    )
                  
                
              
              +
              
                
                  
                    C
                    ⁢
                    
                      
                    
                    ⁢
                    
                      ρ
                      K
                    
                  
                  
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    t
                  
                
                ⁢
                
                  (
                  
                    
                      T
                      s
                    
                    -
                    
                      T
                      0
                    
                  
                  )
                
              
            
            ]
          
        
      
    
    
      
        (
        1
        )
      
    
  




where, τ=Coefficient of transmission, Γ=Coefficient of reflection, σb=Stefan-Boltzmann constant (W/cm2-K4), εir=emissivity of Kapton®, Ts=Temperature of sample (Kelvin), T0=Initial temperature of sample (Kelvin), Δt=Time to equilibrium (seconds), ρK=Density of Kapton® (kg/cm2), and C=Specific Heat (J/g-K). Note that the portion of equation 1 which has the coefficient of transmission and coefficient of reflection outside the brackets is an absorption coefficient, in this case for the EM material. A power density formula, such as shown in equation 1, must be adapted to a particular material composition used in the EM absorber. An exemplary Data Output Processing (DOP) Module 29 provides various outputs including a power density image (PDI) array. The PDI array stores data required to produce a digital representation that is displayed in cooperation with the GUI Module 23 and further discussed in connection with FIG. 4.
The Data Array 30 section of this exemplary embodiment stores data associated with processing in an exemplary embodiment of the system. Data arrays in this embodiment include TDPP produced by IRCOP Module 21, TD array produced by the TDA Module 25, a PDPP array produced by the PDC Module 27, and PDI array associated with DOP Module 29. A Configuration Data Array (CDA) can be provided to store data which is used in processing such as the degrees of freedom data. It should be noted that a variety of data structures can be used to perform the function of Data Array 30. For example, data arrays can be local to code modules or global data arrays.
Referring to FIG. 3, a screen view of a main control graphical user interface 31 is shown associated with processing software. A process control 33 is used to commence execution of the processing system 7. After the process control 33 is selected, a file path input dialog window (not shown) is displayed to permit selection of a file path similar to Windows Explorer for use in navigating file structures and selecting a file stored on a computer system running the software. A dialog (or dialogue) box is a special window, used in user interfaces to display information to the user, or to get a response if needed. They are so-called because they form a dialog between the computer and the user—either informing the user of something, or requesting input from the user, or both. It provides controls that allow a user to specify how to carry out an action. A display field 35 is provided to show the file path information selected using the file path input dialog window (not shown). A Mode Selector 41 dialog box permits a user to select an automatic or manual process for identifying an initial IR image capture frame (e.g., “cold” frame) data and final IR image capture frame (e.g., “hot” frame) data associated with a series of IR image captures outputted by a temperature sensor array (e.g., IR camera) and stored in a data array.
Referring to FIG. 4, another organization of the data arrays described in association with the Data Arrays Module of FIG. 2 is provided. In this example, the data arrays are shown in separate format from the FIG. 2 representation but with the same functionality and purpose to include the TDPP Array 45, TD Array 47, TDPP Array 49, PDPP Array 51, PDI Array 53, and CDA Array 55.
Referring to FIG. 5A, an exemplary power density image output 59 associated with Gaussian Optical Lens Antenna (GOLA) antenna produced by the GUI Module 23 is shown. This power density image output 59 can be produced based on data including PDI array data generated by DOP Module 29. This power density image output 59 can be displayed in the image viewer window 37 (FIG. 3). The power density image output 59 in this embodiment is a graphical image showing power density per pixel of an IR Camera (e.g., 5) produced using embodiment of the invention.
The power density image output 59 in this embodiment includes a graphical representation of an EM field produced by an EM source under test where the representation is made up of different colors associated with different ranges of power density values. This output shows one way to represent power density per pixel by characterizing an EM source (e.g., an antenna field pattern) plus pattern density in one representation. Each pixel in this output has a power density data associated with that pixel. In this example, different colors are used for a number of ranges of power density such as red for a peak range and blue for the lowest data range. These ranges are predetermined, selected in advance, and associated with a particular color representation in this embodiment. The output processing software associated a desired color with a given power density data value associated with a particular pixel by determining which range the power density data value falls within.
For example, in FIG. 5A a series of differently colored regions is shown which are roughly surrounding a maximum EM field incident on an EM absorber under test. In particular, in this exemplary representation a background (blue) region 61 represents no EM field incident on the EM absorber. A light blue or hazy area 65 within the background region represents the weakest measured range of EM field(s) incident on the EM absorber. A green region 63 which is within the white or hazy region 65 represents the next higher range of power density measured associated with the EM field incident on the EM absorber under test. A yellow region 67 which is within the green region 63 represents the next higher range of power density measured associated with the EM field incident on the EM absorber under test. A red region 69 which is within the yellow region 67 represents the highest power density measured associated with the EM field incident on the EM absorber under test.
This power density image output 59 permits characterization of an EM field produced by an EM source under test. The exemplary power density image output 59 also provides an ability to display a specific power density value 39 (FIG. 3) for a specific pixel. One embodiment of the invention can include provision of a pointing device or an electronically manipulated cursor which can be moved over each pixel in the power density image output 59. A power density per pixel data value 39 (FIG. 3) will appear when the pointing device or cursor points to each pixel. One embodiment of the invention can include a provision for displaying a power density maximum data value (Pmax 43, FIG. 3) that will appear showing a highest power density data value 43 for a specific pixel contained in the exemplary power density image output 59.
Outputs (e.g., power density image output 59) from the DOP 29 and GUI Module 23 can be used to determine field strength of an EM source in a specific area of space, near field analysis, or far field analysis. Other modeling and analysis outputs can be produced with this invention based on representation of EM fields associated with an EM source generating a field into free space. One example is producing a representation and analysis of 3 db spot sizes of a directed energy beam or an antenna.
FIG. 5B shows a different representation of an EM field produced by a W-band pyramidal standard gain antenna coupled to an EM source under test using an embodiment of the invention. In this example, the same type of power density color representation scheme is used as is shown in FIG. 5A, however the distribution of the power density/colors are different due to the differences in EM field(s) measured using an embodiment of the invention.
Any number of or variety of antennas can be used with the invention. Other EM sources/antennas/coupling devices which can be used with the invention include waveguides, lasers, or any other EM spectrum source which produces a change, e.g., a temperature change, in the EM absorber.
FIG. 6 shows an exemplary degrees of freedom table having data which is used in association with blackbody calculations. In particular, this table can be used for computing the number of degrees of freedom of a system of atoms or molecules which comprise an EM absorber where n is the number of atoms per molecule, e.g., a sheet of carbon loaded Kapton®. The table includes position, rotation, and vibration for monatomic atoms as well as linear and non-linear molecules. These values are predetermined and used to obtain one factor of the co-efficient of a blackbody radiation term used in calculations associated with an embodiment of the invention. In this example, the FIG. 5 table is used in processing which is performed by the PDC Module 27. The degrees of freedom table can include degrees of information data associated with different types of materials used as EM absorbers. For example, in the exemplary table in FIG. 6 different ways to calculate or represent degrees of freedom for different materials is shown. For Kapton®, there is a carbon powder pressed into a polymer sheet. Because the carbon is in powder form, it has three degrees of freedom of position and 3n-5 degrees of freedom in vibration plus 3 in position which gives 3n-2. Based on this position and vibration information, total degrees of freedom for Kapton® is then 4 (i.e., 3n-2 where n=2 given carbon is a diatomic atom with 2 atoms per molecule) for the powder form of carbon. Different table values and potentially different table compositions or data elements may be required for different embodiments of the invention including different EM absorbers.
FIG. 7A and 7B shows steps used in producing outputs in accordance with one embodiment of the invention. In step 1 (73), provide an EM source (e.g., an RF source). In step 2 (75), provide an EM absorber. In step 3 (77), provide an array of temperature sensors (e.g., an infrared camera). In step 4 (79), provide a temperature control mechanism (e.g., circulation fan) (optional). In step 5 (81), provide a processing system, network and input/output system. In step 6 (83), position the array of temperature sensors and the EM source to orient on the EM absorber based on a desired measurement of an EM field produced by the EM source. In step 7 (85), capture a series of temperature image data frames (e.g., infrared image frame captures) and saving each temperature image data frames that includes temperature per pixel data and time associated with each temperature image data frame to a Temperature Data per Pixel (TDPP) data array. In step 8 (87), start processing software. In step 9 (89), initialize system and input data from TDPP Array. In step 10 (91), select an initial temperature image data frame (e.g., coldest/ambient temperature or time prior to activation of EM source) and a final temperature image data frame (e.g., hot or highest temperature) from the series of temperature image data frames and store in a power density per pixel array (PDPP) array. In step 11 (93), calculate power density associated with each pixel in the initial and final frames selected at step 10 (91) based on determining a plurality of characteristics of an EM absorber including absorption coefficient of the EM absorber, degrees of freedom for atoms or molecules associated with the EM absorber, emissivity, Stefan-Boltzmann's Constant, number of sides being cooled, specific heat capacity of coolant, density of coolant, volume of coolant, area that the coolant will be in contact with EM absorber, change in time between initial and final frame capture, and an initial and final temperature data of the EM absorber. Equation 2 provides an exemplary power density calculation formula that can be used in step 11 (93).




  
    
      
        
          1
          
            1
            -
            τ
            -
            Γ
          
        
        ⁡
        
          [
          
            
              D
              ⁢
              
                
              
              ⁢
              
                
                  ɛσ
                  b
                
                ⁡
                
                  (
                  
                    
                      T
                      4
                    
                    -
                    
                      T
                      0
                      4
                    
                  
                  )
                
              
            
            +
            
              
                ∑
                
                  
                
                ⁢
                
                  
                    
                      C
                      ⁢
                      
                        
                      
                      ⁢
                      ρ
                      ⁢
                      
                        
                      
                      ⁢
                      V
                    
                    
                      Δ
                      ⁢
                      
                        
                      
                      ⁢
                      tA
                    
                  
                  ⁢
                  
                    (
                    
                      T
                      -
                      
                        T
                        0
                      
                    
                    )
                  
                
              
              
                All
                ⁢
                
                  
                
                ⁢
                Heat
                ⁢
                
                  
                
                ⁢
                Sinks
              
            
          
          ]
        
      
    
    
      
        (
        2
        )
      
    
  




Where Γ is the coefficient of absorption of the EM absorber, τ is the coefficient of transmission of the EM absorber, T is the final temperature of the EM absorber, T0 is the initial temperature of the EM absorber, ε is the emissivity of the EM absorber, σb is the Stefan-Boltzmann constant, C is the specific heat of the heat sink (coolant), V is the volume of the heat sink (coolant), ρ is the density of the heat sink (coolant), A is the contact area between the absorber and the heat sink (coolant), Δt is the time between the initial frame to the final frame, and D is the degrees of freedom of the atoms or molecules that comprise the EM absorber. In step 12 (95), produce output of an image representation of power density of EM field calculated for every pixel in the selected IR image capture frames (e.g., “cold” and “hot”).
An alternate embodiment of a power density measurement method in accordance with the invention for pulsed power density measurement could include steps 1 (73) through step 10 (91). The alternate embodiment processing for alternate embodiment step 11 would then deviate using an alternate power density calculation using the formula described in equation 3.




  
    
      
        
          ∑
          
            l
            =
            0
          
          L
        
        ⁢
        
          
        
        ⁢
        
          
            1
            
              1
              -
              τ
              -
              Γ
            
          
          ⁡
          
            [
            
              
                4
                ⁢
                D
                ⁢
                
                  
                
                ⁢
                
                  ɛσ
                  b
                
                ⁢
                
                  
                    T
                    l
                    3
                  
                  ⁡
                  
                    (
                    
                      
                        T
                        l
                      
                      -
                      
                        T
                        
                          l
                          -
                          1
                        
                      
                    
                    )
                  
                
              
              +
              
                
                  ∑
                  
                    
                  
                  ⁢
                  
                    
                      
                        C
                        ⁢
                        
                          
                        
                        ⁢
                        ρ
                        ⁢
                        
                          
                        
                        ⁢
                        V
                      
                      
                        
                          (
                          
                            
                              t
                              l
                            
                            -
                            
                              t
                              
                                l
                                -
                                1
                              
                            
                          
                          )
                        
                        ⁢
                        A
                      
                    
                    ⁢
                    
                      (
                      
                        
                          T
                          l
                        
                        -
                        
                          T
                          
                            l
                            -
                            1
                          
                        
                      
                      )
                    
                  
                
                
                  All
                  ⁢
                  
                    
                  
                  ⁢
                  Heat
                  ⁢
                  
                    
                  
                  ⁢
                  Sinks
                
              
            
            ]
          
        
      
    
    
      
        (
        3
        )
      
    
  




Regarding equation 3, the first summation is a sum over data points measuring the temperature of the absorber. L is the total number of temperature data frames being sampled and summed. l is a variable of summation in equation 3 representing the lth temperature data frame or the variable which is being summed. The second summation is over all heat sinks in contact with the absorber. Γ is the coefficient of absorption of the EM absorber, τ is the coefficient of transmission of the EM absorber, Tl is the final temperature of the EM absorber of the lth frame, Tl−1 is the initial temperature of the EM absorber of the l−1 frame, ε is the emissivity of the EM absorber, σb is the Stefan-Boltzmann constant, C is the specific heat of the heat sink (coolant), V is the volume of the heat sink (coolant), ρ is the density of the heat sink (coolant), A is the contact area between the absorber and the heat sink (coolant), Δt is the time between the initial frame (tl−1) to the final frame (tl), D is the degrees of freedom of the atoms or molecules that comprise the EM absorber. In the pulsed power alternative embodiment, a modified step 10 would include a selection of additional said temperature image data frames in between said initial and final temperature image data frames and determines associated time differentials between consecutive said selected additional temperature image data frames. As described above, a modified step 11 further serially calculates power density between each two of each selected temperature image data frames then said processing continues by summing all of the power density data for each two consecutive temperature image data frame from the initial frame to the final frame to produce a final power density data. For pulsed power density calculations, when summing the power density of all power density computations, the final solution is calculated by dividing by the total number of paired frames minus one. For example, where six frames are used to compute power density for pulsed power, the sum of the power density values must be divided by five (i.e., 6 total frames−1).
Next, at step 12, an output would be generated showing an image of an EM field associated with the EM source based on power density calculated for every pixel in an initial or “cold” IR image capture frame and a final or “hot” IR image capture frame.
In this pulsed power alternate embodiment of a power density measurement method, it is possible to make a pulsed power density measurement using as few as two data points which show the temperature rise from initial temperature to final temperature instead of waiting for the temperature to reach equilibrium.
FIG. 8 shows an exemplary data plot of temperature data versus frame number of data collected from an EM source using one exemplary embodiment of the invention. A processing system in accordance with an embodiment of the invention sequentially analyzes each IR image capture frame starting a selected time e.g., t=0, then looks for the highest temperature reading associated with all the pixels in each frame, then plots the temperature value for each frame. Each data plot on this table shows a maximum temperature value for each separate IR image frame capture. The x-values in this table represent an individual IR image frame. In other words, each temperature data on the y-axis represents a pixel having the highest temperature in a frame capture comprising any number of pixels which were sampled during a given integration time per frame capture using an EM sensor, in this example an IR camera. The plot begins with a series of temperature data readings associated with a series of IR image capture frames before an EM source is energized and oriented on an EM absorber. After the EM source is energized, then a temperature data is plotted sequentially for each IR image capture frame on the x axis until the system reaches the end of the IR image capture frames.
An alternative embodiment of the invention can further include a mechanism adapted to move the EM source and/or the EM absorber relation to each other in order to manipulate either the EM field which is being measured or manipulate the EM absorber to move it through the EM field which is being measured. An additional set of components and software would be used to track the position of the EM field and/or the EM absorber and associate positional and time information with each movement of the EM absorber and/or EM source. These additional components, in addition to the components described above, would be used to collect sufficient data required to compute the power density and other characteristics associated which each plane of the EM field which the EM absorber was placed within. A series of software instructions for controlling the movement of the EM absorber and/or the EM source and then collecting separate data or measurements for each position of the EM source and/or EM absorber would be incorporated into the above embodiments of the invention. The output would then be a series of plane images of the EM field under test or a three dimensional image of the EM field under test which could include a three dimensional power density representation using the same or a different color scheme as described above. An output system could use a pointer system or a plane which could then be used to select specific cross sections of the three dimensional field.
Another embodiment could include manipulating the EM absorber and temperature sensor array around an area where the EM field is larger than the EM absorber. In this embodiment, the temperature sensor array and EM absorber would be positioned to take a number of sequential measurements which would be individually processed, then compiled into a mosaic to describe the overall area of interest where the EM field of interest is located.
Another embodiment can include changing the distance between the EM absorber and the EM source in order to take different measurements which would then be used to determine phase and amplitude characteristics of the EM source.
Another alternative embodiment can be used to characterize an EM absorber where the EM source and the temperature sensor array are calibrated and performance characteristics are known. This embodiment can be then be used to analyze physical characteristics of the EM absorber. An example could be determining a sample characteristics based on comparisons with known values such as specific heat capacity, density of material, atomic makeup of the material under test.
Another alternative embodiment of the invention could be used to characterize EM absorption of an atmosphere. The system would be oriented on an atmosphere where a star would be used as the EM source. An initial frame would be when the star is not oriented on atmosphere of the planet of interest (e.g., on its cold side), then take another series of frame captures when the star is oriented on the atmosphere of interest.
Benefits of the invention include being able to make parallel measurements versus serial or piece-meal measurements. Also, the invention permits high resolution measurements with significantly less effort, equipment, time and resources than currently required to produce desired outputs.
Although the invention has been described in detail with reference to certain preferred embodiments, variations and modifications exist within the spirit and scope of the invention as described and defined in the following claims.

## Claims

1. An electromagnetic field characterization apparatus comprising:
an electromagnetic source adapted to emit an electromagnetic field;
an electromagnetic absorption material adapted to receive said electromagnetic field;
a plurality of temperature sensors adapted to orient on said electromagnetic absorption material and output a plurality of temperature sensor data;
a plurality of processing sequences comprising:
a first processing sequence adapted to receive said plurality of temperature sensor data and store a plurality of temperature sensor data frames created based on said plurality of temperature sensor data in a first data structure, said temperature sensor data frames comprise a plurality of temperature data associated with some or all of said temperature sensors stored at different points of time;
a second processing sequence adapted to select at least two of said temperature sensor data frames based at least on a first and second temperature data parameter determined based on comparisons of all said temperature sensor data frames to find an initial state and a final state associated with temperature changes in said electromagnetic absorption material, said at least two selected temperature sensor data frames comprise an initial and final temperature sensor data frame, said second processing sequence is further adapted to determine a time differential data for every two of said selected at least two temperature sensor data frames, said second processing sequence is further adapted to store said selected at least two temperature sensor data frames and said time differential data in a second data structure;
a third processing sequence adapted to determine a plurality of power density data associated with some or all of said plurality of temperature data in said initial and final temperature sensor data frames, said plurality of power density data is determined based on data stored in said second data structure, absorption attributes associated with said electromagnetic absorption material, a plurality of blackbody radiation attributes associated with said electromagnetic absorption material comprising degrees of freedom attributes, and a plurality of cooling attributes associated with a cooling medium in proximity to said electromagnetic absorption material, said second processing sequence further outputs and stores a power density image map in a third data structure, said power density image map comprising a power density data associated with all or a portion of said plurality of temperature sensors; and

a fourth processing sequence adapted to produce a graphical or data output of said power density image map.

2. An electromagnetic field characterization apparatus as in claim 1, further comprising providing a temperature control mechanism which provides cooling to said electromagnetic ab sorption material.

3. An electromagnetic field characterization apparatus as in claim 2, wherein said temperature control mechanism is a circulation fan adapted to ensure approximately uniform heat transfer associated with said electromagnetic absorption material.

4. An electromagnetic field characterization apparatus as in claim 1, wherein said electromagnetic source is a radio frequency source coupled to an electromagnetic field radiating structure.

5. An electromagnetic field characterization apparatus as in claim 1, wherein said plurality of temperature sensors comprises one or more infrared image sensing systems.

6. An electromagnetic field characterization apparatus as in claim 1, wherein said first temperature parameter is a lowest temperature data value and said second temperature parameter is a highest temperature value.

7. An electromagnetic field characterization apparatus as in claim 1, further comprising a mechanism adapted to move said electromagnetic source or the electromagnetic absorption material in relation to each other in order to manipulate either said electromagnetic field which is being measured or manipulate the electromagnetic absorption material to move it through the electromagnetic field which is being measured in order to produce a plurality of power density image maps each representing a different plane measurement of said electromagnetic field or a three-dimensional power density image map.

8. An electromagnetic field characterization apparatus as in claim 1, wherein said first processing sequence further selects additional said temperature sensor data frames in between said initial and final temperature sensor data frames and determines associated time differentials between consecutive said selected additional temperature sensor data frames, wherein said second processing sequence for calculating a power density data further serially calculates power density between each two of said selected temperature sensor data frames then said second processing sequence sums all of the power density data for each two consecutive temperature sensor data frame from said initial frame to said final frame to produce a final power density data.

9. An electromagnetic field characterization apparatus comprising:
an electromagnetic source adapted to emit electromagnetic field;
an electromagnetic absorption material adapted to receive said electromagnetic field;
a plurality of temperature sensors adapted to orient on said electromagnetic absorption material and output a plurality of temperature sensor data frames data, wherein said plurality of temperature sensors comprise a plurality of temperature sensor pixels and said plurality of temperature sensor data frame data each store a temperature data output from each said temperature sensor pixels;
a plurality of processing sequences comprising:
a first processing sequence adapted to select and store an initial and final temperature sensor data frame from said plurality of temperature sensor data frames based on a first and second temperature parameter, said first processing sequence further determines a differential temperature data associated with said initial and final temperature sensor data frames;
a second processing sequence for calculating a power density data associated with all or a portion of said plurality of temperature sensors and said final temperature sensor data frame, said second processing sequence determines said power density data based on a power density per said pixel calculation comprising a peak power density value based on a plurality of calculation elements comprising a coefficient of absorption of the electromagnetic absorption material, a coefficient of transmission of the electromagnetic absorption material, the final temperature of the electromagnetic absorption material, the initial temperature of the electromagnetic absorption material, emissivity of the electromagnetic absorption material, Stefan-Boltzmann constant, a specific heat of at least one heat sink, a volume of the at least one heat sink, a density of the at least one heat sink, a contact area between the electromagnetic absorption material and the at least one heat sink, a time differential between the initial and final temperature sensor data frames, and degrees of freedom of atoms or molecules that comprise the electromagnetic absorption material and a plurality of inputs to said power density per said pixel calculation including said differential temperature data, wherein said calculation is modified based on a particular said electromagnetic absorption material, wherein the modified calculation components are coefficient of absorption of the electromagnetic absorption material, a plurality of degrees of freedom data for atoms or molecules associated with said electromagnetic absorption material, emissivity of blackbody radiation associated with said electromagnetic absorption material, specific heat capacity of the at least one heat sink placed in proximity to said electromagnetic absorption material, density of the at least one heat sink, volume of the at least one heat sink, and area of interface between the at least one heat sink and said electromagnetic absorption material; and

a third processing sequence for producing at least one graphical user interface, said at least one graphical user interface comprises a power density image output produced based on a plurality of ranges of power density values which depicts said power density data, said at least one graphical user interface also is adapted to display one or more said power density data associated with at least one of said plurality of temperature sensors.

10. An electromagnetic field characterization apparatus as in claim 9, further comprising a processing system adapted to process said plurality of processing sequences and produce a plurality of outputs comprising said power density image output.

11. An electromagnetic field characterization apparatus as in claim 9, further comprising providing a temperature control mechanism which provides cooling to said electromagnetic ab sorption material.

12. An electromagnetic field characterization apparatus as in claim 11, wherein said temperature control mechanism is a circulation fan.

13. An electromagnetic field characterization apparatus as in claim 9, wherein said electromagnetic source is a radio frequency source coupled to an electromagnetic field radiating structure.

14. An electromagnetic field characterization apparatus as in claim 9, wherein said plurality of temperature sensors comprises one or more infrared image sensing apparatus.

15. An electromagnetic field characterization apparatus as in claim 9, wherein said plurality of temperature sensors is a video camera adapted to record infrared or near infrared spectrum energy.

16. An electromagnetic field characterization apparatus as in claim 9, wherein said first temperature parameter is a lowest temperature data value and said second temperature parameter is a highest temperature value.

17. An electromagnetic field characterization apparatus as in claim 9, further comprising a mechanism adapted to move said electromagnetic source or the electromagnetic absorption material relation to each other in order to manipulate either said electromagnetic field which is being measured or manipulate the electromagnetic absorption material to move it through the electromagnetic field which is being measured.

18. An electromagnetic field characterization apparatus as in claim 9, wherein said first processing sequence further selects additional temperature sensor data frames in between said initial and final temperature sensor data frames and determines associated time differentials between consecutive said additional temperature sensor data frames, wherein said second processing sequence for calculating a power density data further serially calculates power density between each two selected temperature sensor data frames then said second processing sequence sums all of the power density data for each two consecutive temperature sensor data frame from said initial frame to said final frame to produce a final power density data.

19. An electromagnetic field characterization apparatus comprising:
an electromagnetic source adapted to produce electromagnetic radiation;
an electromagnetic absorption material adapted and positioned to said electromagnetic radiation;
at least one array of temperature sensors adapted to focus on said electromagnetic absorption material and output a plurality of temperature sensor data frames comprising a plurality of temperature sensor data associated with a single temperature sensor in said array, wherein each said temperature sensor data frame is associated with a different time data value than other said temperature sensor data frames;
a plurality of processing sequences comprising:
a first processing sequence adapted to store said plurality of temperature sensor data frames and store said plurality of temperature sensor data frames in a first data structure;
a second processing sequence adapted to analyze data stored in said first data structure and select an initial temperature sensor data frame and a final temperature sensor data frame, wherein said initial and final temperature sensor data frames are selected from said plurality of sensor data frames stored in said first data structure based on determination of an initial and final temperature data associated with at least two of said plurality of temperature sensor data frames, said second processing sequence further determines a differential time data associated with said initial and final temperature data frames based on a first and second temperature parameter, said second processing sequence further stores said initial and final temperature sensor data frames and said differential time data a second data structure;
a third processing sequence comprises a power density per temperature sensor calculation adapted to calculate a plurality of power density data associated with said final temperature sensor data frames and said time differential data stored in said second data structure, said third processing sequence is further adapted to store said plurality of power density data in a third data structure, each power density data is individually associated with one or more temperature sensors in said at least one array of temperature sensors, wherein said plurality of power density data is calculated based on:
a blackbody radiation calculation based on said electromagnetic radiation associated with said initial and final temperature sensor data frames;
a determination of cooling of said electromagnetic absorption material based on said initial and final temperature sensor data frames and said differential time data;
a coefficient of a blackbody radiation term determined based on at least a plurality of degrees of freedom data each associated with at least one category of atoms or molecules making up said electromagnetic absorption material; a plurality of constants associated with said electromagnetic absorption material comprising emissivity, Stefan-Boltzmann Constant, specific heat capacity of a cooling medium placed in proximity to said electromagnetic absorption material, density of said cooling medium, volume of said cooling medium, and an area that said cooling medium makes contact with said electromagnetic absorption material; and
a data value determined based on how much of said electromagnetic radiation incident on said electromagnetic absorption material is absorbed by said electromagnetic absorption material;


a fourth processing sequence adapted to produce a plurality of outputs comprising a power density image representation based on said plurality of power density data.

20. An electromagnetic field characterization apparatus as in claim 19, further comprising a processing system adapted to process said plurality of processing sequences and data structures and produce a plurality of said outputs comprising said power density image representation.

21. An electromagnetic field characterization apparatus as in claim 19, wherein said second processing sequence further selects additional temperature sensor data frames in between said initial and final temperature sensor data frames and determines associated time differentials between consecutive said additional temperature sensor data frames, wherein said third processing sequence for calculating a power density per temperature sensor data further serially calculates power density between each two selected temperature sensor data frames then said second processing sequence sums all of the power density data for each two consecutive temperature sensor data frames from said initial frame to said final frame to produce a final said plurality of power density data.

22. An electromagnetic field characterization apparatus as in claim 19, further comprising providing a temperature control mechanism which provides cooling to said electromagnetic absorption material.

23. An electromagnetic field characterization apparatus as in claim 19, wherein said temperature control mechanism is a circulation fan.

24. An electromagnetic field characterization apparatus as in claim 19, wherein said electromagnetic source is a radio frequency source coupled to an electromagnetic field radiating structure.

25. An electromagnetic field characterization apparatus as in claim 19, wherein said plurality of temperature sensors comprises one or more infrared image sensing apparatus.

26. An electromagnetic field characterization apparatus as in claim 19, wherein said first temperature parameter is a lowest temperature data value and said second temperature parameter is a highest temperature value.

27. An electromagnetic field characterization apparatus as in claim 19, further comprising a mechanism adapted to move said electromagnetic source or the electromagnetic absorption material relation to each other in order to manipulate either said electromagnetic field which is being measured or manipulate the electromagnetic absorption material to move it through the electromagnetic field which is being measured.

28. A method of manufacturing an electromagnetic field characterization apparatus comprising:
providing an electromagnetic source adapted to emit an electromagnetic field;
providing an electromagnetic absorption material adapted to receive said electromagnetic field;
providing a plurality of temperature sensors adapted to orient on said electromagnetic absorption material and output a plurality of temperature sensor data;
providing a plurality of processing sequences comprising:
a first processing sequence adapted to receive said plurality of temperature sensor data and store a plurality of temperature sensor data frames created based on said plurality of temperature sensor data in a first data structure, said temperature sensor data frames comprise a plurality of temperature data associated with some or all of said temperature sensors stored at different points of time;
a second processing sequence adapted to select at least two of said temperature sensor data frames based at least on a first and second temperature data parameter determined based on comparisons of all said temperature sensor data frames to find an initial state and a final state associated with temperature changes in said electromagnetic absorption material, said at least two selected temperature sensor data frames comprise an initial and final temperature sensor data frame, said second processing sequence is further adapted to determine a time differential data for every two of said selected at least two temperature sensor data frames, said second processing sequence is further adapted to store said selected at least two temperature sensor data frames and said time differential data in a second data structure;
a third processing sequence adapted to determine a plurality of power density data associated with some or all of said plurality of temperature data in said initial and final temperature sensor data frames, said plurality of power density data is determined based on data stored in said second data structure, absorption attributes associated with said electromagnetic absorption material, a plurality of blackbody radiation attributes associated with said electromagnetic absorption material comprising degrees of freedom attributes, and a plurality of cooling attributes associated with a cooling medium in proximity to said electromagnetic absorption material, said second processing sequence further outputs and stores a power density image map in a third data structure, said power density image map comprising a power density data associated with all or a portion of said plurality of temperature sensors; and

a fourth processing sequence adapted to produce a graphical or data output of said power density image map.

29. A method of manufacturing an electromagnetic field characterization apparatus as in claim 28, further comprising providing a temperature control mechanism which provides cooling to said electromagnetic absorption material.

30. A method of manufacturing an electromagnetic field characterization apparatus as in claim 29, wherein said temperature control mechanism is a circulation fan adapted to ensure approximately uniform heat transfer associated with said electromagnetic absorption material.

31. A method of manufacturing an electromagnetic field characterization apparatus as in claim 28, wherein said electromagnetic source is a radio frequency source coupled to an electromagnetic field radiating structure.

32. A method of manufacturing an electromagnetic field characterization apparatus as in claim 28, wherein said plurality of temperature sensors comprises one or more infrared image sensing systems.

33. A method of manufacturing an electromagnetic field characterization apparatus as in claim 28, wherein said first temperature parameter is a lowest temperature data value and said second temperature parameter is a highest temperature value.

34. A method of manufacturing an electromagnetic field characterization apparatus as in claim 28, further comprising a mechanism adapted to move said electromagnetic source or the electromagnetic absorption material in relation to each other in order to manipulate either said electromagnetic field which is being measured or manipulate the electromagnetic absorption material to move it through the electromagnetic field which is being measured in order to produce a plurality of power density image maps each representing a different plane measurement of said electromagnetic field or a three-dimensional power density image map.

35. A method of manufacturing an electromagnetic field characterization apparatus as in claim 28, wherein said first processing sequence further selects additional said temperature sensor data frames in between said initial and final temperature sensor data frames and determines associated time differentials between consecutive said selected additional temperature sensor data frames, wherein said second processing sequence for calculating a power density data further serially calculates power density between each two of said selected temperature sensor data frames then said second processing sequence sums all of the power density data for each two consecutive temperature sensor data frame from said initial frame to said final frame to produce a final power density data.

36. A method of manufacturing an electromagnetic field characterization apparatus comprising:
providing an electromagnetic source adapted to emit electromagnetic field;
providing an electromagnetic absorption material adapted to receive said electromagnetic field;
providing a plurality of temperature sensors adapted to orient on said electromagnetic absorption material and output a plurality of temperature sensor data frames data, wherein said plurality of temperature sensors comprise a plurality of temperature sensor pixels and said plurality of temperature sensor data frame data each store a temperature data output from each said temperature sensor pixels;
providing a plurality of processing sequences comprising:
a first processing sequence adapted to select and store an initial and final temperature sensor data frame from said plurality of temperature sensor data frames based on a first and second temperature parameter, said first processing sequence further determines a differential temperature data associated with said initial and final temperature sensor data frames;
a second processing sequence for calculating a power density data associated with all or a portion of said plurality of temperature sensors and said final temperature sensor data frame, said second processing sequence determines said power density data based on a power density per said pixel calculation comprising a peak power density value based on a plurality of calculation elements comprising a coefficient of ab sorption of the electromagnetic ab sorption material, a coefficient of transmission of the electromagnetic absorption material, the final temperature of the electromagnetic ab sorption material, the initial temperature of the electromagnetic ab sorption material, emissivity of the electromagnetic ab sorption material, Stefan-Boltzmann constant, a specific heat of at least one heat sink, a volume of the at least one heat sink, a density of the at least one heat sink, a contact area between the electromagnetic absorption material and the at least one heat sink, a time differential between the initial and final temperature sensor data frames, and degrees of freedom of atoms or molecules that comprise the electromagnetic absorption material and a plurality of inputs to said power density per said pixel calculation including said differential temperature data, wherein said calculation is modified based on a particular said electromagnetic ab sorption material, wherein the modified calculation components are coefficient of ab sorption of the electromagnetic absorption material, a plurality of degrees of freedom data for atoms or molecules associated with said electromagnetic absorption material associated with a blackbody radiation calculation, emissivity of blackbody radiation associated with said electromagnetic absorption material, specific heat capacity of the at least one heat sink placed in proximity to said electromagnetic absorption material, density of the at least one heat sink, volume of the at least one heat sink, and area of interface between the at least one heat sink and said electromagnetic absorption material; and


a third processing sequence for producing at least one graphical user interface, said at least one graphical user interface comprises a power density image output produced based on a plurality of ranges of power density which depicts said power density data, said at least one graphical user interface also is adapted to display one or more said power density data associated with at least one of said plurality of temperature sensors.

37. A method of manufacturing an electromagnetic field characterization apparatus as in claim 36, further comprising providing a processing system adapted to process said plurality of processing sequences and produce a plurality of outputs comprising said power density image output.

38. A method of manufacturing an electromagnetic field characterization apparatus as in claim 36, further comprising providing a temperature control mechanism which provides cooling to said electromagnetic absorption material.

39. A method of manufacturing an electromagnetic field characterization apparatus as in claim 38, wherein said temperature control mechanism is a circulation fan.

40. A method of manufacturing an electromagnetic field characterization apparatus as in claim 36, wherein said electromagnetic source is a radio frequency source coupled to an electromagnetic field radiating structure.

41. A method of manufacturing an electromagnetic field characterization apparatus as in claim 36, wherein said plurality of temperature sensors comprises one or more infrared image sensing apparatus.

42. A method of manufacturing an electromagnetic field characterization apparatus as in claim 36, wherein said first temperature parameter is a lowest temperature data value and said second temperature parameter is a highest temperature value.

43. A method of manufacturing an electromagnetic field characterization apparatus as in claim 36, further comprising a mechanism adapted to move said electromagnetic source or the electromagnetic absorption material relation to each other in order to manipulate either said electromagnetic field which is being measured or manipulate the electromagnetic absorption material to move it through the electromagnetic field which is being measured.