My own error was causing the problem. I was editing a mapfile that
wasn't actually being used.
On Jan 29, 2008, at 5:53 PM, Paul Moen wrote:
This is my layer.
LAYER
NAME "Montana"
PROJECTION
"init=epsg:1038007"
END
CONNECTION "user=postgres dbname=1038007 host=127.0.0.1
port=5432"
CONNECTIONTYPE postgis
DATA "the_geom from mt_counties using unique gid using
SRID=1038007"
TYPE POLYGON
STATUS ON
CLASS
NAME "MT counties"
STYLE
COLOR 245 245 245
OUTLINECOLOR 200 200 200
END
END
END
I am using the user postgres for the sql statement .
This is the actual command that I used from the command line.
psql -d 1038007 -U postgres -c "SELECT
asbinary(force_collection(force_2d(the_geom)),'NDR'),gid::text from
mt_counties WHERE the_geom && setSRID('BOX3D(1792876.6493994
343910,2160707.3506006 612229)'::BOX3D, 1038007 )"
These are all performed on the same server that PostgreSQL is running.
I do have the database named 1038007. Is there a problem with
database names being all numerals within PostgreSQL? Everything
seems to work fine except when mapserver is making the call.
On Jan 29, 2008, at 5:33 PM, Paul Ramsey wrote:
There's something wrong, because the same SQL should return the
same result... there should not be a non-deterministic difference
depending on from where you enter it. Are you sure you're entering
it in the same database, using the same user, as mapserver is?
P
Paul Moen wrote:
1 row equal to 1038007
On Jan 29, 2008, at 3:46 PM, Paul Ramsey wrote:
what does
SELECT DISTINCT SRID(the_geom) from mt_counties;
return?
what does
SELECT * FROM geometry_columns WHERE f_table_name = 'mt_counties';
return?
Paul Moen wrote:
I get the following error when trying to draw a layer that
exists in postgis.
[MapServer Error]: msPOSTGISLayerWhichShapes(): Error executing
POSTGIS SQL statement (in FETCH ALL): DECLARE mycursor BINARY
CURSOR FOR SELECT
asbinary(force_collection(force_2d(the_geom)),'NDR'),gid::text
from mt_counties WHERE the_geom &&
setSRID('BOX3D(1792876.6493994 343910,2160707.3506006
612229)'::BOX3D, 1038007 )
-ERROR: Operation on two geometries with different SRIDs
If I issue the following sql command in postgis two rows are
returned rather than the error.
SELECT
asbinary(force_collection(force_2d(the_geom)),'NDR'),gid::text
from mt_counties WHERE the_geom &&
setSRID('BOX3D(1792876.6493994 343910,2160707.3506006
612229)'::BOX3D, 1038007 )
The custom SRID has been added to the spatial_ref_sys table in
postgis and the epsg file in the proj resources. These layers
display exactly the way I want when viewing with Qgis connected
as postGIS layers.
Does anyone have any ideas on what I can do to further
troubleshoot this? Is this a problem with mapserver, postgis,
or something I am doing? Version info.
Mac OS 10.4.11
"POSTGIS="1.2.1" GEOS="3.0.0rc4-CAPI-1.3.3" PROJ="Rel. 4.5.0, 22
Oct 2006" USE_STATS"
MapServer Version MapServer version 4.10.1 OUTPUT=GIF
OUTPUT=PNG OUTPUT=JPEG OUTPUT=WBMP OUTPUT=PDF OUTPUT=SWF
OUTPUT=SVG SUPPORTS=PROJ SUPPORTS=FREETYPE SUPPORTS=WMS_SERVER
SUPPORTS=WMS_CLIENT SUPPORTS=WFS_SERVER SUPPORTS=WFS_CLIENT
SUPPORTS=WCS_SERVER SUPPORTS=SOS_SERVER SUPPORTS=GEOS
INPUT=EPPL7 INPUT=POSTGIS INPUT=OGR INPUT=GDAL INPUT=SHAPEFILE
PHP MapScript Version ($Revision: 1.260.2.1 $ $Date:
2006/11/22 19:34:45 $)
Thanks,
Paul
--
Paul Ramsey
Refractions Research
http://www.refractions.net
[EMAIL PROTECTED]
Phone: 250-383-3022
Cell: 250-885-0632
--
Paul Ramsey
Refractions Research
http://www.refractions.net
[EMAIL PROTECTED]
Phone: 250-383-3022
Cell: 250-885-0632
