Difference between revisions of "Build array"

From GeopsyWiki
Jump to navigation Jump to search
 
(3 intermediate revisions by the same user not shown)
Line 1: Line 1:
 
== Visualization of Array Reponse ==
 
== Visualization of Array Reponse ==
  
build_array ia a graphical interface for simple visualization
+
'''build_array''' ia a graphical interface for simple visualization
 
of the array response function for a specific array geometry.  
 
of the array response function for a specific array geometry.  
 
The array reponse function is computed on a regular cartesian grid
 
The array reponse function is computed on a regular cartesian grid
in the wavenumber domain <math>(k_x, k_y)</math> by computing:
+
in the wavenumber domain <math>\vec{k}</math> by evaluating:
  
<math>AR(\vec{k}-\vec{k}_0) = \sum_{i=1}^{N} \exp(-j\vec{r}_i(\vec{k}-\vec{k}_0))</math>
+
<math>AR(\vec{k}-\vec{k}_0) =  
 +
\| \frac{1}{N} \sum_{i=1}^{N} \exp(-j\vec{r}_i(\vec{k}-\vec{k}_0)) \|</math>
  
though the computation of the array reponse function is simple in principle,
+
Here, <math>\vec{r}_i</math> is the position vector of sensor i in  
there exists the particular problem of being computationally efficient
+
a cartesian coordinate system, <math>j</math> the imaginary unit and
on one hand, i.e. small number of grid points, and not to miss small details
+
<math>\vec{k}_0</math> the ''true'' wavenumber vector of a single
of the peaked surface on the other hand, i.e. sufficiently high wavenumber resolution
+
plane wave.
<math>(\delta k_x, \delta k_y)</math>.
+
As most seismological array applications are restricted to the
 +
deployment of sensors on the earth's surface, the above equation
 +
is evaluated in the horizontal plane, i.e. <math>\vec{r}_i = (r_{ix}, r_{iy})</math>
 +
and <math>\vec{k} = (k_x, k_y)</math> (horizontal wavenumber vector).
 +
 
 +
The evaluation of the array reponse function is simple,
 +
still the limits of the wavenumber range as well as the desired
 +
resolution for the grid computation <math>(\delta k_x, \delta k_y)</math>  
 +
needs to be specified. The resolution needs to be sufficient high in
 +
order to not miss small details of the peaked surface and the
 +
size of the grid should allow to determine the position of
 +
grating lobes or strong secondary side lobes in the 2D surface.
 +
 
 +
--[[User:Mao|Mao]] 20:25, 31 October 2008 (UTC)

Latest revision as of 20:25, 31 October 2008

Visualization of Array Reponse

build_array ia a graphical interface for simple visualization of the array response function for a specific array geometry. The array reponse function is computed on a regular cartesian grid in the wavenumber domain by evaluating:

Here, is the position vector of sensor i in a cartesian coordinate system, the imaginary unit and the true wavenumber vector of a single plane wave. As most seismological array applications are restricted to the deployment of sensors on the earth's surface, the above equation is evaluated in the horizontal plane, i.e. and (horizontal wavenumber vector).

The evaluation of the array reponse function is simple, still the limits of the wavenumber range as well as the desired resolution for the grid computation needs to be specified. The resolution needs to be sufficient high in order to not miss small details of the peaked surface and the size of the grid should allow to determine the position of grating lobes or strong secondary side lobes in the 2D surface.

--Mao 20:25, 31 October 2008 (UTC)