Thanks to the both of you for your help -- that helps me quite a bit to understand that how to reproject a srid=-1 on the fly and how i should go about permanently fixing it in my geometry table.<br><br>Take care,<br>Matt<br>
<br><div class="gmail_quote">On Sat, Jun 28, 2008 at 9:35 AM, Paul Spencer <<a href="mailto:pspencer@dmsolutions.ca">pspencer@dmsolutions.ca</a>> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
There is normally a constraint set on postgis tables that specifies the SRID, if you want to change the SRID then you need to drop the constraint, set the SRID and optionally put back the constraint with the correct SRID.<br>
<br>
I do not believe this is necessary to achieve what you want, though. An SRID of -1 tells postgis to return the geometry as is to mapserver. Mapserver then doesn't know what the projection of the data is, but that is okay because you can provide an PROJECTION object on your layer in mapserver. Then mapserver knows what the projection of the data is without having to consult postgis and will be able to reproject the data for you.<br>
<br>
So you should be able to leave your data in the database with SRID -1 and put<br>
<br>
PROJECTION<br>
"init=epsg:26915"<br>
END<br>
<br>
in your layer. Then if you use a different projection at the MAP level, mapserver knows that your data is in UTM zone 15 and will reproject it on the fly.<br>
<br>
If you are *always* going to reprojecting your data then you should actually transform it in the database so that you don't suffer from reprojection overhead on every request.<br>
<br>
Cheers<br>
<br>
Paul<div><div></div><div class="Wj3C7c"><br>
<br>
On 28-Jun-08, at 2:12 AM, <a href="mailto:pcreso@pcreso.com" target="_blank">pcreso@pcreso.com</a> wrote:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi Matthew,<br>
<br>
If you loaded it from shapefiles using shp2pgsql, then use the -s flag to set the srid to the correct value.<br>
<br>
If you \d your table, I think you'll find a constraint on it setting the srid, as well as the geometrycolumns entry.<br>
<br>
The simplest way to get things working with mapserver may well be to leave it as -1 in the database (& geometrycolumns entry) & in your mapserver mapfile data statement, specify the sql returns setsrid(geom,26915) instead of geom, so postgis returns the values with the appropriate SRID set, irrespective of what the database actually stores.<br>
<br>
That will get things working until you load the data with the correct srid.<br>
<br>
HTH,<br>
<br>
Brent Wood<br>
<br>
<br>
<br>
--- On Sat, 6/28/08, Matthew Pettis <<a href="mailto:matthew.pettis@gmail.com" target="_blank">matthew.pettis@gmail.com</a>> wrote:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
From: Matthew Pettis <<a href="mailto:matthew.pettis@gmail.com" target="_blank">matthew.pettis@gmail.com</a>><br>
Subject: [mapserver-users] Basic Projection Help<br>
To: "Mapserver-users" <<a href="mailto:mapserver-users@lists.osgeo.org" target="_blank">mapserver-users@lists.osgeo.org</a>><br>
Date: Saturday, June 28, 2008, 5:47 PM<br>
Hello,<br>
<br>
I'm trying to get a handle on projections and<br>
reprojections. I have a<br>
dataset in a postgis database table that is in epsg:26915<br>
(it's the state of<br>
MN). I know this from the metadata provided by the<br>
provider (state of<br>
Minnesota itself).<br>
<br>
When I loaded it into postgis, the geometry_columns row for<br>
this<br>
table/column combo set srid=-1. As I know that this is in<br>
epsg:26915, I<br>
manually changed the geometry_columns row entry to it (see:<br>
===<br>
geometry_columns === section).<br>
<br>
Problem is, I get an error of using two different<br>
geometries. Anybody have<br>
any idea how i set the srid to the correct value of<br>
epsg:26915? I'd like to<br>
set this correctly so I can reproject it as necessary. The<br>
map will render<br>
with srid=-1, but then I don't think I can reproject<br>
the image correctly,<br>
right?<br>
<br>
below are what I think are the important sections to this<br>
problem. Thanks<br>
in advance for your help.<br>
<br>
Matt<br>
<br>
=== url ===<br>
======================================================================<br>
<a href="http://localhost/cgi-bin/mapserv.exe?map=C:%5Cms4w%5CApache%5Chtdocs%5Col" target="_blank">http://localhost/cgi-bin/mapserv.exe?map=C:\ms4w\Apache\htdocs\ol</a>\mapreproj.map&mode=map<br>
======================================================================<br>
<br>
=== error message ===<br>
======================================================================<br>
msDrawMap(): Image handling error. Failed to draw layer<br>
named 'state'.<br>
msPOSTGISLayerWhichShapes(): Query error. Error executing<br>
POSTGIS SQL<br>
statement (in FETCH ALL): DECLARE mycursor BINARY CURSOR<br>
FOR SELECT<br>
asbinary(force_collection(force_2d(the_geom)),'NDR'),gid::text<br>
from<br>
shp_state WHERE the_geom &&<br>
setSRID('BOX3D(141095.125 4809728.25,810342.625<br>
5478975.75)'::BOX3D,<br>
find_srid('','shp_state','the_geom')<br>
) -ERROR:<br>
Operation on two geometries with different SRIDs<br>
======================================================================<br>
<br>
=== geometry_columns ===<br>
======================================================================<br>
oid;f_table_catalog;f_table_schema;f_table_name;f_geometry_column;coor_dimension;srid;type<br>
336791;"''";"public";"shp_state";"the_geom";2;26915;"MULTIPOLYGON"<br>
======================================================================<br>
<br>
=== mapreproj.map ===<br>
======================================================================<br>
MAP<br>
<br>
################<br>
# Map attributes<br>
################<br>
IMAGETYPE PNG<br>
EXTENT 184056.505 4809728.25 767381.245 5478975.75<br>
SIZE 600 600<br>
IMAGECOLOR 255 255 255<br>
FONTSET<br>
"C:/ms4w/Apache/htdocs/fonts/fonts.list"<br>
<br>
####################<br>
# Web settings<br>
####################<br>
WEB<br>
IMAGEPATH "C:/ms4w/Apache/htdocs/tmp/"<br>
IMAGEURL "/tmp/"<br>
END<br>
<br>
####################<br>
# State<br>
# Background shows through missing precincts to show<br>
where no votes will<br>
be coming from<br>
####################<br>
LAYER<br>
NAME "state"<br>
<br>
CONNECTIONTYPE postgis<br>
CONNECTION "user=postgres dbname=gisvote<br>
password=postgres"<br>
DATA "the_geom from shp_state"<br>
STATUS default<br>
TYPE polygon<br>
<br>
CLASS<br>
NAME "state"<br>
COLOR 200 200 200<br>
END<br>
<br>
END # Layer state<br>
<br>
END<br>
======================================================================<br>
<br>
<br>
-- <br>
It is from the wellspring of our despair and the places<br>
that we are broken<br>
that we come to repair the world.<br>
-- Murray<br>
Waas_______________________________________________<br>
mapserver-users mailing list<br>
<a href="mailto:mapserver-users@lists.osgeo.org" target="_blank">mapserver-users@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/mapserver-users" target="_blank">http://lists.osgeo.org/mailman/listinfo/mapserver-users</a><br>
</blockquote>
_______________________________________________<br>
mapserver-users mailing list<br>
<a href="mailto:mapserver-users@lists.osgeo.org" target="_blank">mapserver-users@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/mapserver-users" target="_blank">http://lists.osgeo.org/mailman/listinfo/mapserver-users</a><br>
</blockquote>
<br>
<br></div></div>
__________________________________________<br><font color="#888888">
<br>
Paul Spencer<br>
Chief Technology Officer<br>
DM Solutions Group Inc<br>
<a href="http://www.dmsolutions.ca/" target="_blank">http://www.dmsolutions.ca/</a><br>
<br>
</font></blockquote></div><br><br clear="all"><br>-- <br>It is from the wellspring of our despair and the places that we are broken that we come to repair the world.<br>-- Murray Waas