まず、 http://www.nws.noaa.gov/geodata/catalog/national/html/us_state.htmのように、それぞれが米国の州に対応する一連のポリゴンを用意する必要があります。
次に、 http://www.ecse.rpi.edu/Homepages/wrf/Research/Short_Notes/pnpoly.html#The C Codeに記載されている次のアルゴリズムを使用できます。
int pnpoly(int nvert, float *vertx, float *verty, float testx, float testy)
{
int i, j, c = 0;
for (i = 0, j = nvert-1; i < nvert; j = i++) {
if ( ((verty[i]>testy) != (verty[j]>testy)) &&
(testx < (vertx[j]-vertx[i]) * (testy-verty[i]) / (verty[j]-verty[i]) + vertx[i]) )
c = !c;
}
return c;
}
nvert: ポリゴンの頂点の数。最後に最初の頂点を繰り返すかどうかについては、以下で説明します。
vertx, verty: ポリゴンの頂点の x 座標と y 座標を含む配列。
testx, testy: テスト ポイントの X 座標と Y 座標。