PostGIS recuperar polígonos dentro de un rectángulo

Question

Estoy construyendo un sistema para llevar en el shapefile de datos de una base de datos PostGIS y render en la parte superior de OpenStreetMap baldosas con el Leaflet.js biblioteca.

Mi tabla (sierra_county_parcels) tiene un polígono, geometría de campo (geom) en SRID 4326, y un campo llamado "upc", que contiene la IDENTIFICACIÓN de la parcela. Hay cerca de 16.000 filas en sierra_county_parcels. Lo que tengo que hacer es recuperar la totalidad de la upc, los valores de las filas donde la geometría de polígono es en la actualidad parte visible de la hoja de control del mapa. Ya he escrito el código para enviar el noroeste y el sureste de lat/lon pares y generar una consulta SQL, que tiene este aspecto:

SELECT upc AS pt_parcel_key
FROM  sierra_county_parcels
WHERE  ST_Within(geom, ST_Envelope(ST_GeomFromText('MULTIPOINT(33.167444534375925 -107.14931488037108, 33.08909661229697 -107.35513687133788)')))

El MULTIPUNTO WKT es lo que estoy utilizando actualmente para especificar las esquinas (nw_latitude, nw_longitude, se_latitude, se_longitude).

El problema es que la consulta devuelve siempre 0 resultados. Para esta consulta, el mapa fue ampliada para mostrar la mayor parte de la región donde las geometrías debe ser.

Alguna idea?

Gracias de antemano,

John W.

Answer 1

Su x y y se invierten. Longitud es X y es lo primera, latitud es y viene en segundo lugar.

Answer 2

Trate de usar ST_MakeEnvelope para hacer un cuadro:

SELECT upc AS pt_parcel_key
FROM sierra_county_parcels
WHERE geom && ST_MakeEnvelope($long_min, $lat_min, $long_max, $lat_max, 4326);

(Nota: reemplazar el $ parámetros con el valor apropiado)