1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 | <?php /*======================================================================= // File: JPGRAPH_MESHINTERPOLATE.INC.PHP // Description: Utility class to do mesh linear interpolation of a matrix // Created: 2009-03-09 // Ver: $Id: jpgraph_meshinterpolate.inc.php 1709 2009-07-30 08:00:08Z ljp $ // // Copyright (c) Aditus Consulting. All rights reserved. //======================================================================== */ /** * Utility function to do linear mesh interpolation * @param $aDat Matrix to interpolate * @param $aFactor Interpolation factor */ function doMeshInterpolate( &$aData, $aFactor ) { $m = new MeshInterpolate(); $aData = $m->Linear($aData,$aFactor); } /** * Utility class to interpolate a given data matrix * */ class MeshInterpolate { private $data = array(); /** * Calculate the mid points of the given rectangle which has its top left * corner at $row,$col. The $aFactordecides how many spliots should be done. * i.e. how many more divisions should be done recursively * * @param $row Top left corner of square to work with * @param $col Top left corner of square to work with * $param $aFactor In how many subsquare should we split this square. A value of 1 indicates that no action */ function IntSquare( $aRow, $aCol, $aFactor ) { if ( $aFactor <= 1 ) return; $step = pow( 2, $aFactor-1 ); $v0 = $this->data[$aRow][$aCol]; $v1 = $this->data[$aRow][$aCol + $step]; $v2 = $this->data[$aRow + $step][$aCol]; $v3 = $this->data[$aRow + $step][$aCol + $step]; $this->data[$aRow][$aCol + $step / 2] = ( $v0 + $v1 ) / 2; $this->data[$aRow + $step / 2][$aCol] = ( $v0 + $v2 ) / 2; $this->data[$aRow + $step][$aCol + $step / 2] = ( $v2 + $v3 ) / 2; $this->data[$aRow + $step / 2][$aCol + $step] = ( $v1 + $v3 ) / 2; $this->data[$aRow + $step / 2][$aCol + $step / 2] = ( $v0 + $v1 + $v2 + $v3 ) / 4; $this->IntSquare( $aRow, $aCol, $aFactor-1 ); $this->IntSquare( $aRow, $aCol + $step / 2, $aFactor-1 ); $this->IntSquare( $aRow + $step / 2, $aCol, $aFactor-1 ); $this->IntSquare( $aRow + $step / 2, $aCol + $step / 2, $aFactor-1 ); } /** * Interpolate values in a matrice so that the total number of data points * in vert and horizontal axis are $aIntNbr more. For example $aIntNbr=2 will * make the data matrice have tiwce as many vertical and horizontal dta points. * * Note: This will blow up the matrcide in memory size in the order of $aInNbr^2 * * @param $ &$aData The original data matricde * @param $aInNbr Interpolation factor * @return the interpolated matrice */ function Linear( &$aData, $aIntFactor ) { $step = pow( 2, $aIntFactor-1 ); $orig_cols = count( $aData[0] ); $orig_rows = count( $aData ); // Number of new columns/rows // N = (a-1) * 2^(f-1) + 1 $p = pow( 2, $aIntFactor-1 ); $new_cols = $p * ( $orig_cols - 1 ) + 1; $new_rows = $p * ( $orig_rows - 1 ) + 1; $this->data = array_fill( 0, $new_rows, array_fill( 0, $new_cols, 0 ) ); // Initialize the new matrix with the values that we know for ( $i = 0; $i < $new_rows; $i++ ) { for ( $j = 0; $j < $new_cols; $j++ ) { $v = 0 ; if ( ( $i % $step == 0 ) && ( $j % $step == 0 ) ) { $v = $aData[$i / $step][$j / $step]; } $this->data[$i][$j] = $v; } } for ( $i = 0; $i < $new_rows-1; $i += $step ) { for ( $j = 0; $j < $new_cols-1; $j += $step ) { $this->IntSquare( $i, $j, $aIntFactor ); } } return $this->data; } } ?> |