<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">POS</journal-id><journal-title-group><journal-title>Positioning</journal-title></journal-title-group><issn pub-type="epub">2150-850X</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/pos.2016.72007</article-id><article-id pub-id-type="publisher-id">POS-66390</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Computer Science&amp;Communications</subject></subj-group></article-categories><title-group><article-title>
 
 
  Generalized Method of Biparametric Sub Pixel Thermal Location
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>.</surname><given-names>Sh. Mehdiyev</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>N.</surname><given-names>A. Abdullayev</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>R.</surname><given-names>N. Abdulov</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>H.</surname><given-names>H. Asadov</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Sevda</surname><given-names>N. Abdullayeva</given-names></name><xref ref-type="aff" rid="aff4"><sup>4</sup></xref></contrib></contrib-group><aff id="aff4"><addr-line>Scientific-Industrial Center, OZONE, Baku, Azerbaijan</addr-line></aff><aff id="aff3"><addr-line>Research Institute of Aerospace Informatics, Baku, Azerbaijan</addr-line></aff><aff id="aff2"><addr-line>Research Institute of Ministry of Defence Industry, Baku, Azerbaijan</addr-line></aff><aff id="aff1"><addr-line>National Aviation Academy, Baku, Azerbaijan</addr-line></aff><pub-date pub-type="epub"><day>12</day><month>05</month><year>2016</year></pub-date><volume>07</volume><issue>02</issue><fpage>75</fpage><lpage>79</lpage><history><date date-type="received"><day>29</day>	<month>January</month>	<year>2016</year></date><date date-type="rev-recd"><day>accepted</day>	<month>9</month>	<year>May</year>	</date><date date-type="accepted"><day>12</day>	<month>May</month>	<year>2016</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  It is well-known that according the Dozier’s method, utilization of integral of Planks function in fusion of signals of two different channels of airborne radiometer makes it possible to compute such components of temperature field within one pixel as temperatures of the object and background. In the paper, the generalization of Dozier method is suggested. The suggested generalization of Dozier’s bispectral method named as biparametric method is applicable for static remote objects. In the suggested biparametric method, the measurements are carried out at the moments 
  t
  <sub>1</sub>
   and 
  t
  <sub>2</sub>
  . It is assumed that the object temperature reaches quantity 
  T
  (
  t
  <sub>1</sub>
  ) and 
  T
  (
  t
  <sub>2</sub>
  ) at these moments. On the bases of operational data of scanning infrared radiometer, the square area of one pixel can be calculated in dependence of distance between object and radiometer. This makes it possible to carry out location of static objects from two basis points using serial single wavelengths measurements of radiation emitted by the sub pixel object.
 
</p></abstract><kwd-group><kwd>Radiometer</kwd><kwd> Pixel</kwd><kwd> Radiation</kwd><kwd> Location</kwd><kwd> Bispectral Measurements</kwd><kwd> Sensor</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>It is well-known that between such spheres of technical cybernetics as location, positioning and remote sensing, the firm interrelation does exist. The information theory based grounding of such an interrelation firstly is described in the work [<xref ref-type="bibr" rid="scirp.66390-ref1">1</xref>] .</p><p>The properties of thermal location are that the pixel type structure of images used for location purposes by scanning airborne radiometers leads to inevitable errors of integrated assessment of signal. Upon thermal scanning of surface of researched object, if the latter contains two different surface materials, all radiation emitted from these materials located within one pixel will be averaged as a single pixel signal depending on wavelength of sensor’s operational channel.</p><p>According to [<xref ref-type="bibr" rid="scirp.66390-ref2">2</xref>] , utilization of integral of Planks function for different channels of airborne radiometer makes it possible to compute following parameters:</p><p>1) Radiation temperature of one of two temperature fields on sub pixel level of resolution;</p><p>2) Share of each component of temperature field within pixel (that is temperatures of the object and background).</p><p>If assume, that effect of atmosphere is lacking, the upward radiation at the input of airborne scanning radiometer will be determined as [<xref ref-type="bibr" rid="scirp.66390-ref1">1</xref>]</p><disp-formula id="scirp.66390-formula351"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-8501123x6.png"  xlink:type="simple"/></disp-formula><p>where:<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-8501123x7.png" xlink:type="simple"/></inline-formula>―emissivity at the wavelength λ;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-8501123x8.png" xlink:type="simple"/></inline-formula>―Plank’s function, W∙m<sup>−3</sup>;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-8501123x9.png" xlink:type="simple"/></inline-formula>―sensor’s spectral instrument function.</p><p>The Planks functions of black body with temperature T at the wavelength λ is determined as</p><disp-formula id="scirp.66390-formula352"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-8501123x10.png"  xlink:type="simple"/></disp-formula><p>where: С<sub>1</sub>―the first constant of Plank, equals to 3.741832 &#215; 10<sup>−16</sup> W∙m<sup>2</sup>;</p><p>С<sub>2</sub>―the second constant of Plank, equals to 1.438786 &#215; 10<sup>−2</sup> m∙k;</p><p>Т―temperature, K;</p><p>λ―wavelength, μm.</p><p>The property of the Planks function is that if in the fixed temperature T<sub>1</sub>, the signal at the wavelength λ<sub>1</sub> surpasses the signal at the wavelength λ<sub>2</sub> in sufficiently higher temperature T<sub>2</sub>, the contrary case does occur.</p><p>This property of the Plank’s function is illustrated in <xref ref-type="fig" rid="fig1">Figure 1</xref>, where the output signals of 3-rd and 4-th channels of AVHRR radiometer installed on the board NOAA-6 are shown. As it seen from shown graphics upon temperature above 460 K, the signal of 4-th channel (102 - 116 μm) surpasses the one of 3-rd channel (35 - 40 μm).</p><p>The above said property is the bases of the Dozier’s method [<xref ref-type="bibr" rid="scirp.66390-ref2">2</xref>] , according which the total radiation emitted from non-homogenous, two-parts structured pixel at the surface of researched site upon lacking of atmospheric effect can be computed as</p><disp-formula id="scirp.66390-formula353"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-8501123x11.png"  xlink:type="simple"/></disp-formula><p>where: p―weight coefficient; 0 &lt; p &lt; 1;</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Dependence of signals of 3-rd and 4-th channels from temperature in spectroradiometer AVHRR of satellite NOAA-6</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-8501123x12.png"/></fig><p>T<sub>t</sub>―objects temperature;</p><p>λ<sub>1</sub>, λ<sub>2</sub>―wavelengths used for measurements.</p><p>According to the Dozier’s method upon carrying out measurements at the wavelengths λ<sub>1</sub> and λ<sub>2</sub>, if the value of T<sub>b</sub> is known, the amount of p and T<sub>t</sub> can be calculated using the system of Equation (3). It should be noted that to remove any possible dynamic errors, all measurements upon realization of this method should be carried out synchronously. This method makes it possible to carry out the sub pixel remote identification of hidden and remote objects.</p></sec><sec id="s2"><title>2. Further Development of Doziers Bispectral Method and Suggested Biparametric Generalization</title><p>The Dozier’s method further was developed and modernized in works [<xref ref-type="bibr" rid="scirp.66390-ref3">3</xref>] - [<xref ref-type="bibr" rid="scirp.66390-ref6">6</xref>] .</p><p>Let us consider in brief the modernized Dozier’s method described in [<xref ref-type="bibr" rid="scirp.66390-ref6">6</xref>] . In this modification the limitation imposed in [<xref ref-type="bibr" rid="scirp.66390-ref2">2</xref>] is characterizing the consideration of only ground component of radiation. In the work [<xref ref-type="bibr" rid="scirp.66390-ref6">6</xref>] the general case of remote assessment of parameters of heated remote sub pixel object is considered taking into account the effect of atmosphere. The Equation (3) in this case should be written as</p><disp-formula id="scirp.66390-formula354"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-8501123x13.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.66390-formula355"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-8501123x14.png"  xlink:type="simple"/></disp-formula><p>where: L<sub>b</sub><sub>,4</sub> and L<sub>b</sub><sub>,11</sub>―atmospheric radiations at the wavelength 4 μm and 11 μm;</p><p>τ<sub>4</sub> and τ<sub>11</sub>―atmospheric transfer at the pertinent wavelengths.</p><p>Operationally, L<sub>b</sub><sub>,4</sub> should be determined<sub> </sub>by averaging the radiation of neighbor pixels, assuming that these pixels are identical on temperature. Solution of Equations (4), (5) relative p an T<sub>i</sub> is carried out as following. From Equations (4) and (5) we get</p><disp-formula id="scirp.66390-formula356"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-8501123x15.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.66390-formula357"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-8501123x16.png"  xlink:type="simple"/></disp-formula><p>Equalling (6) and (7) we get</p><disp-formula id="scirp.66390-formula358"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-8501123x17.png"  xlink:type="simple"/></disp-formula><p>where:<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-8501123x18.png" xlink:type="simple"/></inline-formula>;<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-8501123x19.png" xlink:type="simple"/></inline-formula>.</p><p>If assume, that the radiation parameters of background are known, then obviously both components of the left side of equation (8) upon T<sub>f</sub> &#174; ∞ asymptotically go near to zero. The solution of the task is gradual increase of T<sub>f</sub> till the first component at the left side of (8) would approach zero with acceptable accuracy. Thus, the above said method named as Dozier method make it possible to calculate parameters p and T<sub>t</sub> of sub pixel heated object carrying out radiometric measurements at the wavelength λ<sub>1</sub> = 4 μm; λ<sub>2</sub> = 11 μm. The suggested three- measured interpretation of the Dozier method is illustrated in <xref ref-type="fig" rid="fig2">Figure 2</xref>. The <xref ref-type="fig" rid="fig2">Figure 2</xref> illustrate scheme of radiometric measurements at the wavelengths λ<sub>1</sub> and λ<sub>2</sub> at the moment t<sub>o</sub> carried out for identification of sub pixel object with temperature T<sub>to</sub>. As it is seen from three-measured diagram the line AB determines the graphical interpretation of carried out bispectral measurements.</p><p>The suggested generalization of Dozier’s bispectral method named as biparametric method is applicable for static remote objects. We assume that in the time interval <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-8501123x20.png" xlink:type="simple"/></inline-formula> the temperature of static remove object changes from T<sub>t</sub><sub>1</sub> as far as T<sub>t</sub><sub>2</sub> (<xref ref-type="fig" rid="fig3">Figure 3</xref>). The measurements are carried out at the single wavelength. The characteristics of the searched object is that the function</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-8501123x21.png" xlink:type="simple"/></inline-formula>where: t―time of day, is approximately known.</p><fig-group id="fig2"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Three-measured interpretation of Dozier’s method applicable for remote objects with unchanged temperature.</title></caption><fig id ="fig2_1"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-8501123x22.png"/></fig></fig-group><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Three-measured interpretation of suggested biparametric method applicable for remote objects with changing temperature</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-8501123x23.png"/></fig><p>In the suggested biparametric method the measurements are carried out at the moments t<sub>1</sub> and t<sub>2</sub>. It is assumed that the object temperature reaches quantity T(t<sub>1</sub>) and T(t<sub>2</sub>) at these moments. The Equations (4) and (5) in the suggested biparametric method should be written as</p><disp-formula id="scirp.66390-formula359"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-8501123x24.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.66390-formula360"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-8501123x25.png"  xlink:type="simple"/></disp-formula><p>From the Equation (9) we find</p><disp-formula id="scirp.66390-formula361"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-8501123x26.png"  xlink:type="simple"/></disp-formula><p>From Equation (10) we get</p><disp-formula id="scirp.66390-formula362"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-8501123x27.png"  xlink:type="simple"/></disp-formula><p>From both the Equation (11) and (12) we get</p><disp-formula id="scirp.66390-formula363"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-8501123x28.png"  xlink:type="simple"/></disp-formula><p>If according the initial condition <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-8501123x29.png" xlink:type="simple"/></inline-formula> are known, the parameter <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-8501123x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-8501123x30.png" xlink:type="simple"/></inline-formula> can be calculated using Equation (13). Then using one of Equations (6) or (7) the parameter p can be calculated. In order to determine distance as far as an object the following assumption is used. On the bases of operational data of scanning IR radiometer the square area of one pixel can be calculated in dependence of distance between object and radiometer</p><disp-formula id="scirp.66390-formula364"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-8501123x31.png"  xlink:type="simple"/></disp-formula><p>where: S<sub>pix</sub>―square area of one pixel on the surface of object; l<sub>1</sub>―distance between radiometer and object.</p><p>If</p><disp-formula id="scirp.66390-formula365"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-8501123x32.png"  xlink:type="simple"/></disp-formula><p>where: S<sub>ob</sub>―square area of searched object.</p><p>In view of (14) and (15) we get</p><disp-formula id="scirp.66390-formula366"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-8501123x33.png"  xlink:type="simple"/></disp-formula><p>Having calculated l<sub>1</sub> on the bases of Equation (16) we come to conclusion that the searched object is situated at the distance l<sub>1</sub> from the radiometer. Obviously that in order to locate the object all above procedures should be repeated from another base point. Suppose that the second point is situated in distance l<sub>2</sub> from the radiometer. To locate the object it is quite sufficient to draw circles with radius l<sub>1</sub> and l<sub>2</sub> from these bases points and the crossing point of these circles will determine the point of location of object.</p></sec><sec id="s3"><title>3. Conclusion</title><p>Thus it is shown that the suggested generalization of the Dozier’s method named as biparametric sub pixel method of thermal location makes it possible to carry out such a location of static objects from two basis points using serial single wavelengths measurements of radiation emitted by the sub pixel object. Such a modification of Dozier’s method shows the universal character of two-parametric concept of measurements. Utilization of two wavelengths in known Dozier’s method and two serial time moments in suggested modification of the method prove the big potential of two-parametric concept of remote sub-pixel measurements.</p></sec><sec id="s4"><title>Cite this paper</title><p>A. Sh. Mehdiyev,N. A. Abdullayev,R. N. Abdulov,H. H. Asadov,Sevda N. Abdullayeva, (2016) Generalized Method of Biparametric Sub Pixel Thermal Location. Positioning,07,75-79. doi: 10.4236/pos.2016.72007</p></sec></body><back><ref-list><title>References</title><ref id="scirp.66390-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Asadov, H.H. and Ismaylov, K.K.( 2011) Information Method for Synthesis of Optimal Data Subsystems Designated for Positioning, Location and Remote Sensing Systems. Positioning, 2, 61-64.  
http://dx.doi.org/10.4236/pos.2011.21006</mixed-citation></ref><ref id="scirp.66390-ref2"><label>2</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Dozier</surname><given-names> J. </given-names></name>,<etal>et al</etal>. (<year>1981</year>)<article-title>A Method for Satellite Identification of Surface Temperature Fields of Subpixel Resolution</article-title><source> Remote Sensing of Environment</source><volume> 11</volume>,<fpage> 221</fpage>-<lpage>229</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.66390-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Matson, M. and Dozier, J. (1981) Identification of Subresolution High Temperature Sources Using Termal IR Sensor. Photogrammetric Engineering and Remote Sensing, 47, 1311-1318.</mixed-citation></ref><ref id="scirp.66390-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Matson, M. and Holben, B. (1987) Satellite Detection of Tropical Burning in Brazil. International Journal of Remote Sensing, 8, 509-516. http://dx.doi.org/10.1080/01431168708948657</mixed-citation></ref><ref id="scirp.66390-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Giglio, L. and Justice, C.O. (2003) Effect of Wavelength Selection Characterization of Fire Size and Temperature. International Journal of Remote Sensing, 24, 3515-3520. http://dx.doi.org/10.1080/0143116031000117056</mixed-citation></ref><ref id="scirp.66390-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Giglio, L. and Kendall, J.D. (2001) Application of the Dozier Retrieval to Wildfire Characterization: A Sensitivity Analysis. Remote sensing of Environment, 77, 34-39. http://dx.doi.org/10.1016/S0034-4257(01)00192-4</mixed-citation></ref></ref-list></back></article>