Selección de eventos dentro de un polígono y en la superficie de la Tierra con MATLAB

Resumen

Sea un catálogo de puntos en la superficie de la Tierra del cual necesitamos seleccionar un subconjunto con el criterio de que los puntos estén dentro de un polígono. El problema computacional tiene solución conocida a distancias locales y regionales (<~500 km) ya que, entonces, la superficie de la Tierra se puede aproximar como plana. El propósito de este trabajo es entregar un programa de cómputo que aprovecha funciones de MATLAB para resolver compactamente el problema. Al mismo tiempo, se divulgan técnicas de programación; en particular, de código vectorizado. El procedimiento se comprueba con un catálogo de 76545 hipocentros, cada uno descrito por 23 parámetros (columnas, en el catálogo). El catálogo contiene relocalizaciones de sismos que ocurren en la región fronteriza de ambas Californias (Shearer et al., 2003). Un resultado lateral es que un nido con 17 epicentros ubicados en la falla Agua Blanca tiene un error sistemático fuerte y no señala una actividad de la falla.