src/jp/nyatla/nyartoolkit/core/match/NyARMatchPattDeviationColorData.java

   1 /* \r
   2  * PROJECT: NyARToolkit\r
   3  * --------------------------------------------------------------------------------\r
   4  * This work is based on the original ARToolKit developed by\r
   5  *   Hirokazu Kato\r
   6  *   Mark Billinghurst\r
   7  *   HITLab, University of Washington, Seattle\r
   8  * http://www.hitl.washington.edu/artoolkit/\r
   9  *\r
  10  * The NyARToolkit is Java version ARToolkit class library.\r
  11  * Copyright (C)2008 R.Iizuka\r
  12  *\r
  13  * This program is free software; you can redistribute it and/or\r
  14  * modify it under the terms of the GNU General Public License\r
  15  * as published by the Free Software Foundation; either version 2\r
  16  * of the License, or (at your option) any later version.\r
  17  * \r
  18  * This program is distributed in the hope that it will be useful,\r
  19  * but WITHOUT ANY WARRANTY; without even the implied warranty of\r
  20  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the\r
  21  * GNU General Public License for more details.\r
  22  * \r
  23  * You should have received a copy of the GNU General Public License\r
  24  * along with this framework; if not, write to the Free Software\r
  25  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA\r
  26  * \r
  27  * For further information please contact.\r
  28  *      http://nyatla.jp/nyatoolkit/\r
  29  *      <airmail(at)ebony.plala.or.jp>\r
  30  * \r
  31  */\r
  32 package jp.nyatla.nyartoolkit.core.match;\r
  33 \r
  34 \r
  35 import jp.nyatla.nyartoolkit.core.raster.*;\r
  36 import jp.nyatla.nyartoolkit.core.rasterreader.*;\r
  37 \r
  38 /**\r
  39  * INyARMatchPattのRGBColor差分データを格納するクラスです。\r
  40  *\r
  41  */\r
  42 public class NyARMatchPattDeviationColorData\r
  43 {\r
  44         private int _data[];\r
  45         private double _pow;\r
  46         //\r
  47         private int _number_of_pixels;\r
  48         private int _optimize_for_mod;\r
  49         public int[] refData()\r
  50         {\r
  51                 return this._data;\r
  52         }\r
  53         public double getPow()\r
  54         {\r
  55                 return this._pow;\r
  56         }\r
  57                           \r
  58         public NyARMatchPattDeviationColorData(int i_width,int i_height)\r
  59         {\r
  60                 this._number_of_pixels=i_height*i_width;\r
  61                 this._data=new int[this._number_of_pixels*3];\r
  62                 this._optimize_for_mod=this._number_of_pixels-(this._number_of_pixels%8);       \r
  63                 return;\r
  64         }\r
  65 \r
  66         \r
  67         /**\r
  68          * NyARRasterからパターンデータをセットします。\r
  69          * この関数は、データを元に所有するデータ領域を更新します。\r
  70          * @param i_buffer\r
  71          */\r
  72         public void setRaster(INyARRaster i_raster)\r
  73         {\r
  74                 //画素フォーマット、サイズ制限\r
  75                 assert i_raster.getBufferReader().isEqualBufferType(INyARBufferReader.BUFFERFORMAT_INT1D_X8R8G8B8_32);\r
  76                 assert i_raster.getSize().isEqualSize(i_raster.getSize());\r
  77 \r
  78                 final int[] buf=(int[])i_raster.getBufferReader().getBuffer();\r
  79                 //i_buffer[XRGB]→差分[R,G,B]変換                  \r
  80                 int i;\r
  81                 int ave;//<PV/>\r
  82                 int rgb;//<PV/>\r
  83                 final int[] linput=this._data;//<PV/>\r
  84 \r
  85                 // input配列のサイズとwhも更新// input=new int[height][width][3];\r
  86                 final int number_of_pixels=this._number_of_pixels;\r
  87                 final int for_mod=this._optimize_for_mod;\r
  88 \r
  89                 //<平均値計算(FORの1/8展開)>\r
  90                 ave = 0;\r
  91                 for(i=number_of_pixels-1;i>=for_mod;i--){\r
  92                         rgb = buf[i];ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff);\r
  93                 }\r
  94                 for (;i>=0;) {\r
  95                         rgb = buf[i];ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff);i--;\r
  96                         rgb = buf[i];ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff);i--;\r
  97                         rgb = buf[i];ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff);i--;\r
  98                         rgb = buf[i];ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff);i--;\r
  99                         rgb = buf[i];ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff);i--;\r
 100                         rgb = buf[i];ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff);i--;\r
 101                         rgb = buf[i];ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff);i--;\r
 102                         rgb = buf[i];ave += ((rgb >> 16) & 0xff) + ((rgb >> 8) & 0xff) + (rgb & 0xff);i--;\r
 103                 }\r
 104                 //<平均値計算(FORの1/8展開)/>\r
 105                 ave=number_of_pixels*255*3-ave;\r
 106                 ave =255-(ave/ (number_of_pixels * 3));//(255-R)-ave を分解するための事前計算\r
 107 \r
 108                 int sum = 0,w_sum;\r
 109                 int input_ptr=number_of_pixels*3-1;\r
 110                 //<差分値計算(FORの1/8展開)>\r
 111                 for (i = number_of_pixels-1; i >= for_mod;i--) {\r
 112                         rgb = buf[i];\r
 113                         w_sum = (ave - (rgb & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//B\r
 114                         w_sum = (ave - ((rgb >> 8) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//G\r
 115                         w_sum = (ave - ((rgb >> 16) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//R\r
 116                 }\r
 117                 for (; i >=0;) {\r
 118                         rgb = buf[i];i--;\r
 119                         w_sum = (ave - (rgb & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//B\r
 120                         w_sum = (ave - ((rgb >> 8) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//G\r
 121                         w_sum = (ave - ((rgb >> 16) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//R\r
 122                         rgb = buf[i];i--;\r
 123                         w_sum = (ave - (rgb & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//B\r
 124                         w_sum = (ave - ((rgb >> 8) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//G\r
 125                         w_sum = (ave - ((rgb >> 16) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//R\r
 126                         rgb = buf[i];i--;\r
 127                         w_sum = (ave - (rgb & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//B\r
 128                         w_sum = (ave - ((rgb >> 8) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//G\r
 129                         w_sum = (ave - ((rgb >> 16) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//R\r
 130                         rgb = buf[i];i--;\r
 131                         w_sum = (ave - (rgb & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//B\r
 132                         w_sum = (ave - ((rgb >> 8) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//G\r
 133                         w_sum = (ave - ((rgb >> 16) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//R\r
 134                         rgb = buf[i];i--;\r
 135                         w_sum = (ave - (rgb & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//B\r
 136                         w_sum = (ave - ((rgb >> 8) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//G\r
 137                         w_sum = (ave - ((rgb >> 16) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//R\r
 138                         rgb = buf[i];i--;\r
 139                         w_sum = (ave - (rgb & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//B\r
 140                         w_sum = (ave - ((rgb >> 8) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//G\r
 141                         w_sum = (ave - ((rgb >> 16) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//R\r
 142                         rgb = buf[i];i--;\r
 143                         w_sum = (ave - (rgb & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//B\r
 144                         w_sum = (ave - ((rgb >> 8) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//G\r
 145                         w_sum = (ave - ((rgb >> 16) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//R\r
 146                         rgb = buf[i];i--;\r
 147                         w_sum = (ave - (rgb & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//B\r
 148                         w_sum = (ave - ((rgb >> 8) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//G\r
 149                         w_sum = (ave - ((rgb >> 16) & 0xff)) ;linput[input_ptr--] = w_sum;sum += w_sum * w_sum;//R\r
 150                 }\r
 151                 //<差分値計算(FORの1/8展開)/>\r
 152                 final double p=Math.sqrt((double) sum);\r
 153                 this._pow=p!=0.0?p:0.0000001;\r
 154                 return;\r
 155         }\r
 156 \r
 157 }\r