Many thanks to all who replied to my past 2 weeks ago. I ended up going with Regina's suggested LEFT JOIN, which worked wonders. Question though:<br><br>Given Regin's suggested LEFT JOIN:<br><br>
SELECT A.*<br>
FROM A LEFT JOIN B ON ST_Intersects(A.the_geom, B.the_geom)<br>
LEFT JOIN C ON ST_Intersects(A.the_geom, C.the_geom)<br>
<br>
WHERE B.gid IS NULL AND C.gid IS NULL;<br><br>This finds me all geometries from A that have 0 intersection with B OR C, which is what I was looking for. But what if I want to find all polys in A that have <100% overlap with B or C. Meaning, some overlap is fine so long as A is not contained by B or C. <br>
<br>Cheers,<br>dylan<br><br><br><br><br><br><div class="gmail_quote">On Fri, Apr 4, 2008 at 12:00 PM, <<a href="mailto:postgis-users-request@postgis.refractions.net">postgis-users-request@postgis.refractions.net</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;">Send postgis-users mailing list submissions to<br>
<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:postgis-users-request@postgis.refractions.net">postgis-users-request@postgis.refractions.net</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:postgis-users-owner@postgis.refractions.net">postgis-users-owner@postgis.refractions.net</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of postgis-users digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Re: TIGER/Line Shapefiles released (Stephen Frost)<br>
2. Invoking GIST index on 2 disjoints within a query (Dylan Lorimer)<br>
3. Re: Invoking GIST index on 2 disjoints within a query<br>
(Paul Ramsey)<br>
4. Re: TIGER/Line Shapefiles released (Jeshua Lacock)<br>
5. Re: Invoking GIST index on 2 disjoints within a query<br>
(Kevin Neufeld)<br>
6. RE: TIGER/Line Shapefiles released (Paragon Corporation)<br>
7. TIGER/LINE Shapefiles SMID mapping (Paragon Corporation)<br>
8. Re: Invoking GIST index on 2 disjoints within a query (nw)<br>
9. RE: Invoking GIST index on 2 disjoints within a query<br>
(Paragon Corporation)<br>
10. Re: TIGER/Line Shapefiles released (Jonathan W. Lowe)<br>
11. Re: postgis (Colin Wetherbee)<br>
12. Re: TIGER/Line Shapefiles released (Frank Durstewitz)<br>
13. Re: TIGER/Line Shapefiles released (Jonathan W. Lowe)<br>
14. Re: TIGER/Line Shapefiles released (Stephen Frost)<br>
15. Re: TIGER/Line Shapefiles released (Stephen Frost)<br>
16. RE: TIGER/Line Shapefiles released (Paragon Corporation)<br>
17. Re: TIGER/Line Shapefiles released (Jonathan W. Lowe)<br>
18. Re: TIGER/Line Shapefiles released (Ron M)<br>
19. Re: TIGER/Line Shapefiles released (Jonathan W. Lowe)<br>
20. Re: TIGER/Line Shapefiles released (Stephen Frost)<br>
21. Re: TIGER/Line Shapefiles released (Stephen Woodbridge)<br>
22. Re: TIGER/Line Shapefiles released (Stephen Woodbridge)<br>
23. Re: TIGER/Line Shapefiles released (Stephen Frost)<br>
24. Re: TIGER/Line Shapefiles released (Stephen Frost)<br>
25. RE: TIGER/Line Shapefiles released (Paragon Corporation)<br>
26. Re: TIGER/Line Shapefiles released (Stephen Woodbridge)<br>
27. Re: TIGER/Line Shapefiles released (Stephen Woodbridge)<br>
28. Re: TIGER/Line Shapefiles released (Stephen Frost)<br>
29. Re: TIGER/Line Shapefiles released (Stephen Woodbridge)<br>
30. Cannot display postGIS layers on ms4w 4.6.1<br>
(<a href="mailto:kreshna_iceheart@yahoo.com">kreshna_iceheart@yahoo.com</a>)<br>
31. RE: Cannot display postGIS layers on ms4w 4.6.1<br>
(Paragon Corporation)<br>
32. Re: Problem with the Codification/charset from Postgis to<br>
mapserver (Barend Kobben)<br>
33. Re: TIGER/Line Shapefiles released (Jonathan W. Lowe)<br>
34. Re: TIGER/Line Shapefiles released (Stephen Frost)<br>
35. Re: TIGER/Line Shapefiles released (Nick Black)<br>
36. Re: TIGER/Line Shapefiles released (Stephen Frost)<br>
37. Re: TIGER/Line Shapefiles released (Jonathan W. Lowe)<br>
38. Re: TIGER/Line Shapefiles released (Nick Black)<br>
39. ERROR: Operation on mixed SRID geometries after upgrading to<br>
recent postgres/postgis (Sebastian Reitenbach)<br>
40. RE: ERROR: Operation on mixed SRID geometries afterupgrading<br>
to recent postgres/postgis (Obe, Regina)<br>
41. Re: TIGER/Line Shapefiles released (Stephen Woodbridge)<br>
42. RE: ERROR: Operation on mixed SRID geometries afterupgrading<br>
to recent postgres/postgis (Sebastian Reitenbach)<br>
43. FOSS4G 2008 Call for Workshops / Papers (Paul Ramsey)<br>
44. Loading Tiger2007fe data - Thoughts on using inherited tables<br>
(Stephen Woodbridge)<br>
45. Re: Loading Tiger2007fe data - Thoughts on using inherited<br>
tables (Stephen Frost)<br>
46. Re: Loading Tiger2007fe data - Thoughts on using inherited<br>
tables (Paul Ramsey)<br>
47. Re: Loading Tiger2007fe data - Thoughts on using inherited<br>
tables (Stephen Frost)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Thu, 3 Apr 2008 17:07:46 -0400<br>
From: Stephen Frost <<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:20080403210746.GV4999@tamriel.snowman.net">20080403210746.GV4999@tamriel.snowman.net</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
* Stephen Frost (<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>) wrote:<br>
> I think they may have also upgraded their pipe.. I got about 1.41MB/s<br>
> (11 Mb/s) for the whole transfer. It's about 22G all told. I'll<br>
> probably be trying to load it up into PG on one of our servers tomorrow.<br>
> It was a bit over 4 hours for me to pull down off of their<br>
> <a href="http://ftp2.census.gov" target="_blank">ftp2.census.gov</a> ftp site.<br>
<br>
Just to update those who might be interested- I've finished the data<br>
load into one of our servers at work. It comes to ~60GB on disk in<br>
PostgreSQL/PostGIS with appropriate indexes in most places and whatnot.<br>
Based on what I've seen so far, it looks *very* nice, especially the<br>
hydrogrophy ("areawater"). It also appears to be pretty consistant<br>
across the layers, which is also good.<br>
<br>
If anyone's interested in the scripts used to load the data (they're<br>
pretty simple, really), I'd be happy to provide them.<br>
<br>
Enjoy,<br>
<br>
Stephen<br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: not available<br>
Type: application/pgp-signature<br>
Size: 189 bytes<br>
Desc: Digital signature<br>
Url : <a href="http://lists.refractions.net/pipermail/postgis-users/attachments/20080403/bb74af0c/attachment-0001.bin" target="_blank">http://lists.refractions.net/pipermail/postgis-users/attachments/20080403/bb74af0c/attachment-0001.bin</a><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Thu, 3 Apr 2008 16:12:38 -0500<br>
From: "Dylan Lorimer" <<a href="mailto:edylan@google.com">edylan@google.com</a>><br>
Subject: [postgis-users] Invoking GIST index on 2 disjoints within a<br>
query<br>
To: "PostGIS Users Discussion" <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID:<br>
<<a href="mailto:3ddac38d0804031412h16e7241fu2c3823a0e798edf6@mail.gmail.com">3ddac38d0804031412h16e7241fu2c3823a0e798edf6@mail.gmail.com</a>><br>
Content-Type: text/plain; charset=ISO-8859-1<br>
<br>
Hi,<br>
<br>
So I have probably a simple question but I don't use PostGIS enough to<br>
know the right way to go about this. Hoping someone can help.<br>
<br>
So I've got 3 spatial tables, each containing a bunch of polygon<br>
geometries. Let's call the tables A, B, C. What I want is to find out<br>
which polygons in table A are disjoint from all polygons in tables B<br>
and C. Meaning, which polygons in A have 0 overlap with polygons in B<br>
or C.<br>
<br>
I seem to have no troubles doing this when using only 2 tables, but<br>
adding the 3rd is messing things up. Oh, and I also want to invoke the<br>
GIST indexes to speed things up.<br>
<br>
Here's the (slightly generalized) query I think I should be using:<br>
<br>
SELECT <some stuff> FROM A, B, C WHERE disjoint(A.the_geom,<br>
B.the_geom) AND disjoint(A.the_geom, C.the_geom) AND <some other<br>
attribute filters> AND A.the_geom && B.the_geom AND A.the_geom &&<br>
C.the_geom;<br>
<br>
I've not seen this query complete successfully as it takes so long,<br>
but I think something is wrong with it as running it against only A<br>
and B is really quick and there aren't too many geometries in the<br>
tables.<br>
<br>
Thoughts?<br>
<br>
Cheers,<br>
dylan<br>
<br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Thu, 3 Apr 2008 14:23:27 -0700<br>
From: "Paul Ramsey" <<a href="mailto:pramsey@cleverelephant.ca">pramsey@cleverelephant.ca</a>><br>
Subject: Re: [postgis-users] Invoking GIST index on 2 disjoints within<br>
a query<br>
To: "PostGIS Users Discussion" <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID:<br>
<<a href="mailto:30fe546d0804031423k3169a38at3b02346530e3fe80@mail.gmail.com">30fe546d0804031423k3169a38at3b02346530e3fe80@mail.gmail.com</a>><br>
Content-Type: text/plain; charset=ISO-8859-1<br>
<br>
Drop the &&. The indexes are not much help for testing Disjoint, unfortunately.<br>
<br>
The other thing you could do is re-cast it as a does-not-intersect<br>
test instead...<br>
<br>
SELECT <some stuff> FROM A, B, C WHERE NOT ( ST_Intersects(A.the_geom,<br>
B.the_geom) OR ST_Intersects(A.the_geom, C.the_geom) );<br>
<br>
(The above will be indexed automatically for recent versions of<br>
PostGIS. For older, add in the && clause manually like this.)<br>
<br>
SELECT <some stuff> FROM A, B, C WHERE NOT ( (A.the_geom && B.the_geom<br>
AND Intersects(A.the_geom, B.the_geom)) OR (A.the_geom && c.the_geom<br>
AND ST_Intersects(A.the_geom, C.the_geom) ) );<br>
<br>
YMMV... the B and C tables might be interacting and blowing up the<br>
query, I don't have an intuitive feel for this stuff like Regina :)<br>
<br>
P.<br>
<br>
On Thu, Apr 3, 2008 at 2:12 PM, Dylan Lorimer <<a href="mailto:edylan@google.com">edylan@google.com</a>> wrote:<br>
> Hi,<br>
><br>
> So I have probably a simple question but I don't use PostGIS enough to<br>
> know the right way to go about this. Hoping someone can help.<br>
><br>
> So I've got 3 spatial tables, each containing a bunch of polygon<br>
> geometries. Let's call the tables A, B, C. What I want is to find out<br>
> which polygons in table A are disjoint from all polygons in tables B<br>
> and C. Meaning, which polygons in A have 0 overlap with polygons in B<br>
> or C.<br>
><br>
> I seem to have no troubles doing this when using only 2 tables, but<br>
> adding the 3rd is messing things up. Oh, and I also want to invoke the<br>
> GIST indexes to speed things up.<br>
><br>
> Here's the (slightly generalized) query I think I should be using:<br>
><br>
> SELECT <some stuff> FROM A, B, C WHERE disjoint(A.the_geom,<br>
> B.the_geom) AND disjoint(A.the_geom, C.the_geom) AND <some other<br>
> attribute filters> AND A.the_geom && B.the_geom AND A.the_geom &&<br>
> C.the_geom;<br>
><br>
> I've not seen this query complete successfully as it takes so long,<br>
> but I think something is wrong with it as running it against only A<br>
> and B is really quick and there aren't too many geometries in the<br>
> tables.<br>
><br>
> Thoughts?<br>
><br>
> Cheers,<br>
> dylan<br>
> _______________________________________________<br>
> postgis-users mailing list<br>
> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
><br>
<br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Thu, 3 Apr 2008 15:33:07 -0600<br>
From: Jeshua Lacock <jeshua@3DTOPO.com><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>>,<br>
Stephen Frost <<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>><br>
Message-ID: <C1CC27EC-9092-4FF3-B63E-4C22431211B7@3DTOPO.com><br>
Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes<br>
<br>
<br>
On Apr 3, 2008, at 3:07 PM, Stephen Frost wrote:<br>
<br>
> If anyone's interested in the scripts used to load the data (they're<br>
> pretty simple, really), I'd be happy to provide them<br>
<br>
Hello Stephen,<br>
<br>
Sure; I'll bite. No sense reinventing the wheel if it can be helped.<br>
If it is more appropriate, please feel free to email me off list.<br>
<br>
<br>
Thanks,<br>
<br>
Jeshua Lacock<br>
Founder/Programmer<br>
3DTOPO Incorporated<br>
<<a href="http://3DTOPO.com" target="_blank">http://3DTOPO.com</a>><br>
Phone: 877.240.1364<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Thu, 03 Apr 2008 14:35:47 -0700<br>
From: Kevin Neufeld <<a href="mailto:kneufeld@refractions.net">kneufeld@refractions.net</a>><br>
Subject: Re: [postgis-users] Invoking GIST index on 2 disjoints within<br>
a query<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:47F54DB3.4070802@refractions.net">47F54DB3.4070802@refractions.net</a>><br>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>
<br>
Hi Dylan,<br>
<br>
You could try something like this:<br>
SELECT <some stuff><br>
FROM a<br>
LEFT JOIN b ON (ST_Intersects(a.the_geom, b.the_geom))<br>
LEFT JOIN c ON (ST_Intersects(a.the_geom, c.the_geom))<br>
WHERE b.the_geom IS NULL<br>
AND c.the_geom IS NULL;<br>
<br>
-- Kevin<br>
<br>
Dylan Lorimer wrote:<br>
> Hi,<br>
><br>
> So I have probably a simple question but I don't use PostGIS enough to<br>
> know the right way to go about this. Hoping someone can help.<br>
><br>
> So I've got 3 spatial tables, each containing a bunch of polygon<br>
> geometries. Let's call the tables A, B, C. What I want is to find out<br>
> which polygons in table A are disjoint from all polygons in tables B<br>
> and C. Meaning, which polygons in A have 0 overlap with polygons in B<br>
> or C.<br>
><br>
> I seem to have no troubles doing this when using only 2 tables, but<br>
> adding the 3rd is messing things up. Oh, and I also want to invoke the<br>
> GIST indexes to speed things up.<br>
><br>
> Here's the (slightly generalized) query I think I should be using:<br>
><br>
> SELECT <some stuff> FROM A, B, C WHERE disjoint(A.the_geom,<br>
> B.the_geom) AND disjoint(A.the_geom, C.the_geom) AND <some other<br>
> attribute filters> AND A.the_geom && B.the_geom AND A.the_geom &&<br>
> C.the_geom;<br>
><br>
> I've not seen this query complete successfully as it takes so long,<br>
> but I think something is wrong with it as running it against only A<br>
> and B is really quick and there aren't too many geometries in the<br>
> tables.<br>
><br>
> Thoughts?<br>
><br>
> Cheers,<br>
> dylan<br>
> _______________________________________________<br>
> postgis-users mailing list<br>
> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
><br>
<br>
<br>
------------------------------<br>
<br>
Message: 6<br>
Date: Thu, 3 Apr 2008 17:47:14 -0400<br>
From: "Paragon Corporation" <<a href="mailto:lr@pcorp.us">lr@pcorp.us</a>><br>
Subject: RE: [postgis-users] TIGER/Line Shapefiles released<br>
To: "'PostGIS Users Discussion'"<br>
<<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>>, "'Stephen Frost'"<br>
<<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>><br>
Message-ID: <000e01c895d4$47c29900$4d812e40@l><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
I would be interested too.<br>
<br>
Thanks,<br>
Regina<br>
<br>
-----Original Message-----<br>
From: <a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a><br>
[mailto:<a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a>] On Behalf Of Jeshua<br>
Lacock<br>
Sent: Thursday, April 03, 2008 5:33 PM<br>
To: PostGIS Users Discussion; Stephen Frost<br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
<br>
<br>
On Apr 3, 2008, at 3:07 PM, Stephen Frost wrote:<br>
<br>
> If anyone's interested in the scripts used to load the data (they're<br>
> pretty simple, really), I'd be happy to provide them<br>
<br>
Hello Stephen,<br>
<br>
Sure; I'll bite. No sense reinventing the wheel if it can be helped.<br>
If it is more appropriate, please feel free to email me off list.<br>
<br>
<br>
Thanks,<br>
<br>
Jeshua Lacock<br>
Founder/Programmer<br>
3DTOPO Incorporated<br>
<<a href="http://3DTOPO.com" target="_blank">http://3DTOPO.com</a>><br>
Phone: 877.240.1364<br>
<br>
_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
<a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 7<br>
Date: Thu, 3 Apr 2008 17:51:47 -0400<br>
From: "Paragon Corporation" <<a href="mailto:lr@pcorp.us">lr@pcorp.us</a>><br>
Subject: [postgis-users] TIGER/LINE Shapefiles SMID mapping<br>
To: "'PostGIS Users Discussion'"<br>
<<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <001201c895d4$ea384220$4d812e40@l><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
I've been looking at the TIGER/LINE 2007 shapefiles too and I notice that<br>
they have a field called SMID (Spatial Metadata Identifier) in the edges<br>
table. I thought these were supposed to be defined in the edges.shp.xml<br>
file or in the PDF help, but on closer inspection of that file I see this<br>
<br>
- <attr><br>
<attrlabl>SMID</attrlabl><br>
<attrdef>Spatial metadata identifier</attrdef><br>
<attrdefs>U.S. Census Bureau</attrdefs><br>
- <attrdomv><br>
- <edom><br>
<edomv>1</edomv><br>
<edomvd>Spatial MID is positive, non-zero number</edomvd><br>
<edomvds>U.S. Census Bureau</edomvds><br>
</edom><br>
</attrdomv><br>
</attr><br>
<br>
<br>
Well duh? I loaded all of Florida and there are 69 distinct values in that<br>
and I have no clue what they mean. Does any one by chance know where these<br>
codes are defined?<br>
<br>
Thanks,<br>
Regina<br>
<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 8<br>
Date: Thu, 3 Apr 2008 17:01:43 -0500<br>
From: nw <<a href="mailto:nw@hydaspes.if.org">nw@hydaspes.if.org</a>><br>
Subject: Re: [postgis-users] Invoking GIST index on 2 disjoints within<br>
a query<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:8913B9C0-76F3-467E-8809-EEDAD73A15A3@hydaspes.if.org">8913B9C0-76F3-467E-8809-EEDAD73A15A3@hydaspes.if.org</a>><br>
Content-Type: text/plain; charset=US-ASCII; format=flowed<br>
<br>
On Apr 3, 2008, at 16:12, Dylan Lorimer wrote:<br>
<br>
> What I want is to find out which polygons in table A<br>
> are disjoint from all polygons in tables B and C.<br>
<br>
> Here's the (slightly generalized) query I think I should be using:<br>
><br>
> SELECT <some stuff> FROM A, B, C WHERE disjoint(A.the_geom,<br>
> B.the_geom) AND disjoint(A.the_geom, C.the_geom) AND <some other<br>
> attribute filters> AND A.the_geom && B.the_geom AND A.the_geom &&<br>
> C.the_geom;<br>
><br>
> I've not seen this query complete successfully as it takes so long,<br>
> but I think something is wrong with it as running it against only A<br>
> and B is really quick and there aren't too many geometries in the<br>
> tables.<br>
><br>
> Thoughts?<br>
<br>
Use a union.<br>
<br>
select <some stuff> from <working a vs b query><br>
union<br>
select <some compatible stuff> from <working a vs c query><br>
;<br>
<br>
--<br>
nw<br>
<a href="mailto:nw@hydaspes.if.org">nw@hydaspes.if.org</a><br>
<br>
<br>
<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 9<br>
Date: Thu, 3 Apr 2008 18:14:37 -0400<br>
From: "Paragon Corporation" <<a href="mailto:lr@pcorp.us">lr@pcorp.us</a>><br>
Subject: RE: [postgis-users] Invoking GIST index on 2 disjoints within<br>
a query<br>
To: "'PostGIS Users Discussion'"<br>
<<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <000201c895d8$1b153b70$4d812e40@l><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
If I understand you correctly, then I think the most efficient route is to<br>
figure out which ones overlap with an Intersect as Paul mentioned, and then<br>
throw them out. I presume you mean intersect rather than overlaps - ie.<br>
Overlaps means it can't be contained within where as intersects can be.<br>
Either way flip st_intersect with st_overlaps if you really want overlap<br>
check.<br>
<br>
SELECT A.*<br>
FROM A LEFT JOIN B ON ST_Intersects(A.the_geom, B.the_geom)<br>
LEFT JOIN C ON ST_Intersects(A.the_geom, C.the_geom)<br>
<br>
WHERE B.gid IS NULL AND C.gid IS NULL;<br>
<br>
Hope that helps,<br>
Regina<br>
<br>
<br>
<br>
-----Original Message-----<br>
From: <a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a><br>
[mailto:<a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a>] On Behalf Of Paul<br>
Ramsey<br>
Sent: Thursday, April 03, 2008 5:23 PM<br>
To: PostGIS Users Discussion<br>
Subject: Re: [postgis-users] Invoking GIST index on 2 disjoints within a<br>
query<br>
<br>
Drop the &&. The indexes are not much help for testing Disjoint,<br>
unfortunately.<br>
<br>
The other thing you could do is re-cast it as a does-not-intersect test<br>
instead...<br>
<br>
SELECT <some stuff> FROM A, B, C WHERE NOT ( ST_Intersects(A.the_geom,<br>
B.the_geom) OR ST_Intersects(A.the_geom, C.the_geom) );<br>
<br>
(The above will be indexed automatically for recent versions of PostGIS. For<br>
older, add in the && clause manually like this.)<br>
<br>
SELECT <some stuff> FROM A, B, C WHERE NOT ( (A.the_geom && B.the_geom AND<br>
Intersects(A.the_geom, B.the_geom)) OR (A.the_geom && c.the_geom AND<br>
ST_Intersects(A.the_geom, C.the_geom) ) );<br>
<br>
YMMV... the B and C tables might be interacting and blowing up the query, I<br>
don't have an intuitive feel for this stuff like Regina :)<br>
<br>
P.<br>
<br>
On Thu, Apr 3, 2008 at 2:12 PM, Dylan Lorimer <<a href="mailto:edylan@google.com">edylan@google.com</a>> wrote:<br>
> Hi,<br>
><br>
> So I have probably a simple question but I don't use PostGIS enough<br>
> to know the right way to go about this. Hoping someone can help.<br>
><br>
> So I've got 3 spatial tables, each containing a bunch of polygon<br>
> geometries. Let's call the tables A, B, C. What I want is to find out<br>
> which polygons in table A are disjoint from all polygons in tables B<br>
> and C. Meaning, which polygons in A have 0 overlap with polygons in B<br>
> or C.<br>
><br>
> I seem to have no troubles doing this when using only 2 tables, but<br>
> adding the 3rd is messing things up. Oh, and I also want to invoke the<br>
> GIST indexes to speed things up.<br>
><br>
> Here's the (slightly generalized) query I think I should be using:<br>
><br>
> SELECT <some stuff> FROM A, B, C WHERE disjoint(A.the_geom,<br>
> B.the_geom) AND disjoint(A.the_geom, C.the_geom) AND <some other<br>
> attribute filters> AND A.the_geom && B.the_geom AND A.the_geom &&<br>
> C.the_geom;<br>
><br>
> I've not seen this query complete successfully as it takes so long,<br>
> but I think something is wrong with it as running it against only A<br>
> and B is really quick and there aren't too many geometries in the<br>
> tables.<br>
><br>
> Thoughts?<br>
><br>
> Cheers,<br>
> dylan<br>
> _______________________________________________<br>
> postgis-users mailing list<br>
> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
><br>
_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
<a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 10<br>
Date: Thu, 03 Apr 2008 23:40:57 +0100<br>
From: "Jonathan W. Lowe" <<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <1207262457.4621.6.camel@localhost.localdomain><br>
Content-Type: text/plain<br>
<br>
Stephen,<br>
<br>
Have you yet tried overlaying TIGER 2007 linework or census block/tract<br>
polygons over Google or OpenStreetMap tiles? I'm seeing a good match in<br>
some areas but a significant shift (~50 meters) in others. Thought it<br>
might be a datum conversion issue, but can't seem to find a match.<br>
<br>
Jonathan<br>
<br>
On Thu, 2008-04-03 at 17:07 -0400, Stephen Frost wrote:<br>
> * Stephen Frost (<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>) wrote:<br>
> > I think they may have also upgraded their pipe.. I got about 1.41MB/s<br>
> > (11 Mb/s) for the whole transfer. It's about 22G all told. I'll<br>
> > probably be trying to load it up into PG on one of our servers tomorrow.<br>
> > It was a bit over 4 hours for me to pull down off of their<br>
> > <a href="http://ftp2.census.gov" target="_blank">ftp2.census.gov</a> ftp site.<br>
><br>
> Just to update those who might be interested- I've finished the data<br>
> load into one of our servers at work. It comes to ~60GB on disk in<br>
> PostgreSQL/PostGIS with appropriate indexes in most places and whatnot.<br>
> Based on what I've seen so far, it looks *very* nice, especially the<br>
> hydrogrophy ("areawater"). It also appears to be pretty consistant<br>
> across the layers, which is also good.<br>
><br>
> If anyone's interested in the scripts used to load the data (they're<br>
> pretty simple, really), I'd be happy to provide them.<br>
><br>
> Enjoy,<br>
><br>
> Stephen<br>
> _______________________________________________<br>
> postgis-users mailing list<br>
> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 11<br>
Date: Thu, 03 Apr 2008 18:45:02 -0400<br>
From: Colin Wetherbee <<a href="mailto:cww@denterprises.org">cww@denterprises.org</a>><br>
Subject: Re: [postgis-users] postgis<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:47F55DEE.8050503@denterprises.org">47F55DEE.8050503@denterprises.org</a>><br>
Content-Type: text/plain; charset=UTF-8; format=flowed<br>
<br>
Gregory Williamson wrote:<br>
> It's been a while since I used MMS, but IIRC it does need a reference<br>
> to a unique id; OID is picked as a default if there is none<br>
> specified.<br>
<br>
Yes, this is true. It doesn't have to be any special or meaningful<br>
unique ID, like an OID or a primary key, but MapServer looks for those.<br>
Anything else must be explicitly specified.<br>
<br>
> I think you can put something in the DATA statement like the<br>
> following, assuming that the table has "gid" as a unique id (add a<br>
> serial column and populate it is an easy way to add such a thing).<br>
> Then:<br>
><br>
> DATA "the_geom FROM (SELECT table1.the_geom AS the_geom, table1.gid<br>
> AS gid, table2.data AS data FROM table1 LEFT JOIN table2 ON <a href="http://table1.id" target="_blank">table1.id</a><br>
> = <a href="http://table2.id" target="_blank">table2.id</a>) AS new_table USING UNIQUE gid USING SRID=-1"<br>
><br>
> Which I think will make the interface use "gid" instead of trying for<br>
> OIDs, which are being deprecated in PostgreSQL itself.<br>
<br>
In general, that DATA statement looks correct.<br>
<br>
As another example, here's a DATA statement of mine.<br>
<br>
DATA "line FROM (<br>
SELECT num, id, line<br>
FROM jsview_journeyroutes<br>
WHERE userid = %userid%<br>
) AS foo USING UNIQUE id USING SRID=4326"<br>
<br>
It's always a good idea to set the USING SRID=xxxx to the proper SRID<br>
for the data set.<br>
<br>
> You can also dump the table data, drop the table, recreate it with an<br>
> explicit "WITH OIDS" in the table creation, reload the data and away<br>
> you go without changing the MMS interface. But it's better to get rid<br>
> of OIDs where you can since they are not reliablely unique.<br>
<br>
Yes, it's better to get rid of OIDs. Creating an UNIQUE SERIAL column<br>
is preferable.<br>
<br>
Colin<br>
<br>
<br>
------------------------------<br>
<br>
Message: 12<br>
Date: Fri, 04 Apr 2008 00:58:45 +0200<br>
From: Frank Durstewitz <<a href="mailto:frank.durstewitz@emporis.com">frank.durstewitz@emporis.com</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:47F56125.9050308@emporis.com">47F56125.9050308@emporis.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Skipped content of type multipart/alternative-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: smime.p7s<br>
Type: application/x-pkcs7-signature<br>
Size: 3281 bytes<br>
Desc: S/MIME Cryptographic Signature<br>
Url : <a href="http://lists.refractions.net/pipermail/postgis-users/attachments/20080404/d7c9cbe9/smime-0001.bin" target="_blank">http://lists.refractions.net/pipermail/postgis-users/attachments/20080404/d7c9cbe9/smime-0001.bin</a><br>
<br>
------------------------------<br>
<br>
Message: 13<br>
Date: Fri, 04 Apr 2008 00:09:54 +0100<br>
From: "Jonathan W. Lowe" <<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <1207264194.4621.8.camel@localhost.localdomain><br>
Content-Type: text/plain<br>
<br>
I'm interested as well.<br>
Thanks, Jonathan<br>
<br>
On Fri, 2008-04-04 at 00:58 +0200, Frank Durstewitz wrote:<br>
> Me too, please.<br>
><br>
> Regards, Frank<br>
><br>
> Paragon Corporation am 03.04.2008 23:47:<br>
> > I would be interested too.<br>
> ><br>
> > Thanks,<br>
> > Regina<br>
> ><br>
> > -----Original Message-----<br>
> > From: <a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a><br>
> > [mailto:<a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a>] On Behalf Of Jeshua<br>
> > Lacock<br>
> > Sent: Thursday, April 03, 2008 5:33 PM<br>
> > To: PostGIS Users Discussion; Stephen Frost<br>
> > Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
> ><br>
> ><br>
> > On Apr 3, 2008, at 3:07 PM, Stephen Frost wrote:<br>
> ><br>
> ><br>
> > > If anyone's interested in the scripts used to load the data (they're<br>
> > > pretty simple, really), I'd be happy to provide them<br>
> > ><br>
> ><br>
> > Hello Stephen,<br>
> ><br>
> > Sure; I'll bite. No sense reinventing the wheel if it can be helped.<br>
> > If it is more appropriate, please feel free to email me off list.<br>
> ><br>
> ><br>
> > Thanks,<br>
> ><br>
> > Jeshua Lacock<br>
> > Founder/Programmer<br>
> > 3DTOPO Incorporated<br>
> > <<a href="http://3DTOPO.com" target="_blank">http://3DTOPO.com</a>><br>
> > Phone: 877.240.1364<br>
> ><br>
> > _______________________________________________<br>
> > postgis-users mailing list<br>
> > <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> > <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
> ><br>
> ><br>
> > _______________________________________________<br>
> > postgis-users mailing list<br>
> > <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> > <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
> ><br>
> ><br>
> _______________________________________________<br>
> postgis-users mailing list<br>
> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 14<br>
Date: Thu, 3 Apr 2008 19:15:54 -0400<br>
From: Stephen Frost <<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:20080403231553.GW4999@tamriel.snowman.net">20080403231553.GW4999@tamriel.snowman.net</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
* Stephen Frost (<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>) wrote:<br>
> If anyone's interested in the scripts used to load the data (they're<br>
> pretty simple, really), I'd be happy to provide them.<br>
<br>
Alright, since it's apparently kind of popular, I went ahead and put the<br>
script up here:<br>
<br>
<a href="http://snowman.net/tiger/" target="_blank">http://snowman.net/tiger/</a><br>
<br>
I also included a README.txt that is essentially what I wrote up to the<br>
first person who asked me for the script. Please let me know if you<br>
improve upon it or find problems with it.<br>
<br>
Thanks,<br>
<br>
Stephen<br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: not available<br>
Type: application/pgp-signature<br>
Size: 189 bytes<br>
Desc: Digital signature<br>
Url : <a href="http://lists.refractions.net/pipermail/postgis-users/attachments/20080403/3bc528c7/attachment-0001.bin" target="_blank">http://lists.refractions.net/pipermail/postgis-users/attachments/20080403/3bc528c7/attachment-0001.bin</a><br>
<br>
------------------------------<br>
<br>
Message: 15<br>
Date: Thu, 3 Apr 2008 19:18:00 -0400<br>
From: Stephen Frost <<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: "Jonathan W. Lowe" <<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>>, PostGIS Users<br>
Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:20080403231800.GX4999@tamriel.snowman.net">20080403231800.GX4999@tamriel.snowman.net</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Jonathan,<br>
<br>
* Jonathan W. Lowe (<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>) wrote:<br>
> Have you yet tried overlaying TIGER 2007 linework or census block/tract<br>
> polygons over Google or OpenStreetMap tiles? I'm seeing a good match in<br>
> some areas but a significant shift (~50 meters) in others. Thought it<br>
> might be a datum conversion issue, but can't seem to find a match.<br>
<br>
I hadn't looked at the linework too much yet or tried to overlay it.<br>
I'm curious where you're seeing the differences though because I know<br>
that Census is only about half way through their MAF improvment project<br>
and I actually have some info about what has been done so far and what<br>
hasn't. It'd be interesting to see if it matches up.<br>
<br>
There are a few places (Guam, Hawaii islands) where they actually do use<br>
an SRID other than 4269, but my scripts don't yet handle that and I'm<br>
guessing that's not what you're referring to anyway. :)<br>
<br>
Thanks!<br>
<br>
Stephen<br>
<br>
> On Thu, 2008-04-03 at 17:07 -0400, Stephen Frost wrote:<br>
> > * Stephen Frost (<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>) wrote:<br>
> > > I think they may have also upgraded their pipe.. I got about 1.41MB/s<br>
> > > (11 Mb/s) for the whole transfer. It's about 22G all told. I'll<br>
> > > probably be trying to load it up into PG on one of our servers tomorrow.<br>
> > > It was a bit over 4 hours for me to pull down off of their<br>
> > > <a href="http://ftp2.census.gov" target="_blank">ftp2.census.gov</a> ftp site.<br>
> ><br>
> > Just to update those who might be interested- I've finished the data<br>
> > load into one of our servers at work. It comes to ~60GB on disk in<br>
> > PostgreSQL/PostGIS with appropriate indexes in most places and whatnot.<br>
> > Based on what I've seen so far, it looks *very* nice, especially the<br>
> > hydrogrophy ("areawater"). It also appears to be pretty consistant<br>
> > across the layers, which is also good.<br>
> ><br>
> > If anyone's interested in the scripts used to load the data (they're<br>
> > pretty simple, really), I'd be happy to provide them.<br>
> ><br>
> > Enjoy,<br>
> ><br>
> > Stephen<br>
> > _______________________________________________<br>
> > postgis-users mailing list<br>
> > <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> > <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
><br>
> _______________________________________________<br>
> postgis-users mailing list<br>
> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: not available<br>
Type: application/pgp-signature<br>
Size: 189 bytes<br>
Desc: Digital signature<br>
Url : <a href="http://lists.refractions.net/pipermail/postgis-users/attachments/20080403/d5bf23c9/attachment-0001.bin" target="_blank">http://lists.refractions.net/pipermail/postgis-users/attachments/20080403/d5bf23c9/attachment-0001.bin</a><br>
<br>
------------------------------<br>
<br>
Message: 16<br>
Date: Thu, 3 Apr 2008 19:37:36 -0400<br>
From: "Paragon Corporation" <<a href="mailto:lr@pcorp.us">lr@pcorp.us</a>><br>
Subject: RE: [postgis-users] TIGER/Line Shapefiles released<br>
To: "'PostGIS Users Discussion'"<br>
<<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <006601c895e3$b26a87e0$4d812e40@l><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
As a side note to what you were saying in the .sh file it would be really<br>
nice if the shp2pgsql dealt with standalone dbf files too. I mean the logic<br>
is all there so doesn't seem like it would be that hard to put in a switch<br>
for that.<br>
<br>
For the dbfs I was using Ogr2Ogr which works well except it adds a useless<br>
geometry field.<br>
<br>
<br>
<br>
-----Original Message-----<br>
From: <a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a><br>
[mailto:<a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a>] On Behalf Of Stephen<br>
Frost<br>
Sent: Thursday, April 03, 2008 7:16 PM<br>
To: PostGIS Users Discussion<br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
<br>
* Stephen Frost (<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>) wrote:<br>
> If anyone's interested in the scripts used to load the data (they're<br>
> pretty simple, really), I'd be happy to provide them.<br>
<br>
Alright, since it's apparently kind of popular, I went ahead and put the<br>
script up here:<br>
<br>
<a href="http://snowman.net/tiger/" target="_blank">http://snowman.net/tiger/</a><br>
<br>
I also included a README.txt that is essentially what I wrote up to the<br>
first person who asked me for the script. Please let me know if you improve<br>
upon it or find problems with it.<br>
<br>
Thanks,<br>
<br>
Stephen<br>
<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 17<br>
Date: Fri, 04 Apr 2008 01:07:03 +0100<br>
From: "Jonathan W. Lowe" <<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: Stephen Frost <<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>><br>
Cc: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <1207267623.4932.38.camel@localhost.localdomain><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Stephen,<br>
<br>
My initial testing has been on Alameda County (California) TIGER data.<br>
The two attached image files show an overlay of US Census 2000 Blocks<br>
over an area south of the UC Berkeley campus. The offset is the same<br>
for both Google and OpenStreetMap (OSM). This suggests that I've made a<br>
mistake somewhere, because the OSM tiles in the United States are all<br>
rendered from TIGER linework, so the TIGER census blocks should match<br>
exactly.<br>
<br>
For the same source shapefile (tabblock00.shp), there's a nearly perfect<br>
match between block boundaries and streets in the area just South of<br>
Oakland's Lake Merritt. It smells like a datum conversion issue...<br>
<br>
The conversion path was from shapefile to PostGIS using shp2pgsql. I<br>
used a custom projection of 32767 rather than 4269 because the existing<br>
srtext for 4269 had a degree value as 0.01745329251994328, but the US<br>
Census metadata listed a degree value of 0.017453292519943295. Perhaps<br>
not significant? My spatial_ref_sys entries for 4269 and 32767 are<br>
otherwise pretty similar:<br>
<br>
SRID: 4269<br>
SRTEXT: GEOGCS["NAD83",DATUM["North_American_Datum_1983",<br>
SPHEROID["GRS 1980",6378137,298.257222101,<br>
AUTHORITY["EPSG","7019"]],<br>
AUTHORITY["EPSG","6269"]],<br>
PRIMEM["Greenwich",0,<br>
AUTHORITY["EPSG","8901"]],<br>
UNIT["degree",0.01745329251994328,<br>
AUTHORITY["EPSG","9122"]],<br>
AUTHORITY["EPSG","4269"]]<br>
PROJ4TEXT: +proj=longlat +ellps=GRS80 +datum=NAD83 +no_defs<br>
<br>
SRID: 32767<br>
SRTEXT: GEOGCS["GCS_North_American_1983",<br>
DATUM["D_North_American_1983",<br>
SPHEROID["GRS_1980",6378137,298.257222101]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["Degree",0.017453292519943295]]<br>
PROJ4TEXT: +proj=longlat +ellps=clrk66 +datum=NAD27 +no_defs<br>
<br>
To display census block data in OpenStreetMap, I extract it from PostGIS<br>
with a transform to EPSG 4326, although the coordinates don't seem to<br>
change as a result. (This seems correct, as datum=NAD83 and datum=WGS84<br>
are, for my purposes at least, are essentially identical.)<br>
<br>
Thanks,<br>
Jonathan<br>
<br>
2 attachments: TIGER2007andOSM.png, TIGER2007andGoogle.png<br>
<br>
<br>
On Thu, 2008-04-03 at 19:18 -0400, Stephen Frost wrote:<br>
> Jonathan,<br>
><br>
> * Jonathan W. Lowe (<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>) wrote:<br>
> > Have you yet tried overlaying TIGER 2007 linework or census block/tract<br>
> > polygons over Google or OpenStreetMap tiles? I'm seeing a good match in<br>
> > some areas but a significant shift (~50 meters) in others. Thought it<br>
> > might be a datum conversion issue, but can't seem to find a match.<br>
><br>
> I hadn't looked at the linework too much yet or tried to overlay it.<br>
> I'm curious where you're seeing the differences though because I know<br>
> that Census is only about half way through their MAF improvment project<br>
> and I actually have some info about what has been done so far and what<br>
> hasn't. It'd be interesting to see if it matches up.<br>
><br>
> There are a few places (Guam, Hawaii islands) where they actually do use<br>
> an SRID other than 4269, but my scripts don't yet handle that and I'm<br>
> guessing that's not what you're referring to anyway. :)<br>
><br>
> Thanks!<br>
><br>
> Stephen<br>
><br>
> > On Thu, 2008-04-03 at 17:07 -0400, Stephen Frost wrote:<br>
> > > * Stephen Frost (<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>) wrote:<br>
> > > > I think they may have also upgraded their pipe.. I got about 1.41MB/s<br>
> > > > (11 Mb/s) for the whole transfer. It's about 22G all told. I'll<br>
> > > > probably be trying to load it up into PG on one of our servers tomorrow.<br>
> > > > It was a bit over 4 hours for me to pull down off of their<br>
> > > > <a href="http://ftp2.census.gov" target="_blank">ftp2.census.gov</a> ftp site.<br>
> > ><br>
> > > Just to update those who might be interested- I've finished the data<br>
> > > load into one of our servers at work. It comes to ~60GB on disk in<br>
> > > PostgreSQL/PostGIS with appropriate indexes in most places and whatnot.<br>
> > > Based on what I've seen so far, it looks *very* nice, especially the<br>
> > > hydrogrophy ("areawater"). It also appears to be pretty consistant<br>
> > > across the layers, which is also good.<br>
> > ><br>
> > > If anyone's interested in the scripts used to load the data (they're<br>
> > > pretty simple, really), I'd be happy to provide them.<br>
> > ><br>
> > > Enjoy,<br>
> > ><br>
> > > Stephen<br>
> > > _______________________________________________<br>
> > > postgis-users mailing list<br>
> > > <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> > > <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
> ><br>
> > _______________________________________________<br>
> > postgis-users mailing list<br>
> > <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> > <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: TIGER2007andGoogle.png<br>
Type: image/png<br>
Size: 123343 bytes<br>
Desc: not available<br>
Url : <a href="http://lists.refractions.net/pipermail/postgis-users/attachments/20080404/5c180b3f/TIGER2007andGoogle-0001.png" target="_blank">http://lists.refractions.net/pipermail/postgis-users/attachments/20080404/5c180b3f/TIGER2007andGoogle-0001.png</a><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: TIGER2007andOSM.png<br>
Type: image/png<br>
Size: 252886 bytes<br>
Desc: not available<br>
Url : <a href="http://lists.refractions.net/pipermail/postgis-users/attachments/20080404/5c180b3f/TIGER2007andOSM-0001.png" target="_blank">http://lists.refractions.net/pipermail/postgis-users/attachments/20080404/5c180b3f/TIGER2007andOSM-0001.png</a><br>
<br>
------------------------------<br>
<br>
Message: 18<br>
Date: Thu, 03 Apr 2008 17:08:29 -0700<br>
From: Ron M <<a href="mailto:rm_postgis@cheapcomplexdevices.com">rm_postgis@cheapcomplexdevices.com</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: "Jonathan W. Lowe" <<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>>, PostGIS Users<br>
Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:47F5717D.5080301@cheapcomplexdevices.com">47F5717D.5080301@cheapcomplexdevices.com</a>><br>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>
<br>
<br>
I'm curious if they've fixed some of the self-inconsistencies they<br>
seem to have had (or did I just load the data wrong) in the 2006 (or<br>
was it 2005) data I loaded earlier.<br>
<br>
An example is Interstate 280's odd gap behind Stanford.<br>
<br>
<a href="http://map1.forensiclogic.com/maps/mapcache.pl?userid=1&sessionid=0&features=0&layer=land&layer=roads&layer=userfeatures&map_size=640+480&map=cp.map&mapext=-122.2233+37.38634+-122.1655+37.42966&mode=map" target="_blank">http://map1.forensiclogic.com/maps/mapcache.pl?userid=1&sessionid=0&features=0&layer=land&layer=roads&layer=userfeatures&map_size=640+480&map=cp.map&mapext=-122.2233+37.38634+-122.1655+37.42966&mode=map</a><br>
<br>
Other examples were the Dumbarton and San Mateo bridges across<br>
San Francisco Bay not quite connecting with each other; but I<br>
can't provide a link for that because IIRC we hand-edited those.<br>
<br>
<br>
Stephen Frost wrote:<br>
> Jonathan,<br>
><br>
> * Jonathan W. Lowe (<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>) wrote:<br>
>> Have you yet tried overlaying TIGER 2007 linework or census block/tract<br>
>> polygons over Google or OpenStreetMap tiles? I'm seeing a good match in<br>
>> some areas but a significant shift (~50 meters) in others. Thought it<br>
>> might be a datum conversion issue, but can't seem to find a match.<br>
><br>
> I hadn't looked at the linework too much yet or tried to overlay it.<br>
> I'm curious where you're seeing the differences though because I know<br>
> that Census is only about half way through their MAF improvment project<br>
> and I actually have some info about what has been done so far and what<br>
> hasn't. It'd be interesting to see if it matches up.<br>
><br>
> There are a few places (Guam, Hawaii islands) where they actually do use<br>
> an SRID other than 4269, but my scripts don't yet handle that and I'm<br>
> guessing that's not what you're referring to anyway. :)<br>
><br>
> Thanks!<br>
><br>
> Stephen<br>
><br>
>> On Thu, 2008-04-03 at 17:07 -0400, Stephen Frost wrote:<br>
>>> * Stephen Frost (<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>) wrote:<br>
>>>> I think they may have also upgraded their pipe.. I got about 1.41MB/s<br>
>>>> (11 Mb/s) for the whole transfer. It's about 22G all told. I'll<br>
>>>> probably be trying to load it up into PG on one of our servers tomorrow.<br>
>>>> It was a bit over 4 hours for me to pull down off of their<br>
>>>> <a href="http://ftp2.census.gov" target="_blank">ftp2.census.gov</a> ftp site.<br>
>>> Just to update those who might be interested- I've finished the data<br>
>>> load into one of our servers at work. It comes to ~60GB on disk in<br>
>>> PostgreSQL/PostGIS with appropriate indexes in most places and whatnot.<br>
>>> Based on what I've seen so far, it looks *very* nice, especially the<br>
>>> hydrogrophy ("areawater"). It also appears to be pretty consistant<br>
>>> across the layers, which is also good.<br>
>>><br>
>>> If anyone's interested in the scripts used to load the data (they're<br>
>>> pretty simple, really), I'd be happy to provide them.<br>
>>><br>
>>> Enjoy,<br>
>>><br>
>>> Stephen<br>
>>> _______________________________________________<br>
>>> postgis-users mailing list<br>
>>> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
>> _______________________________________________<br>
>> postgis-users mailing list<br>
>> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
>><br>
>> ------------------------------------------------------------------------<br>
>><br>
>> _______________________________________________<br>
>> postgis-users mailing list<br>
>> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 19<br>
Date: Fri, 04 Apr 2008 01:27:48 +0100<br>
From: "Jonathan W. Lowe" <<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <1207268868.4932.43.camel@localhost.localdomain><br>
Content-Type: text/plain<br>
<br>
...And in case the images don't persist through the mail server, they're<br>
viewable at: <a href="http://www.giswebsite.com/demos/tiger_overlays.html" target="_blank">http://www.giswebsite.com/demos/tiger_overlays.html</a><br>
<br>
On Fri, 2008-04-04 at 01:07 +0100, Jonathan W. Lowe wrote:<br>
> Stephen,<br>
><br>
> My initial testing has been on Alameda County (California) TIGER data.<br>
> The two attached image files show an overlay of US Census 2000 Blocks<br>
> over an area south of the UC Berkeley campus. The offset is the same<br>
> for both Google and OpenStreetMap (OSM). This suggests that I've made a<br>
> mistake somewhere, because the OSM tiles in the United States are all<br>
> rendered from TIGER linework, so the TIGER census blocks should match<br>
> exactly.<br>
><br>
> For the same source shapefile (tabblock00.shp), there's a nearly perfect<br>
> match between block boundaries and streets in the area just South of<br>
> Oakland's Lake Merritt. It smells like a datum conversion issue...<br>
><br>
> The conversion path was from shapefile to PostGIS using shp2pgsql. I<br>
> used a custom projection of 32767 rather than 4269 because the existing<br>
> srtext for 4269 had a degree value as 0.01745329251994328, but the US<br>
> Census metadata listed a degree value of 0.017453292519943295. Perhaps<br>
> not significant? My spatial_ref_sys entries for 4269 and 32767 are<br>
> otherwise pretty similar:<br>
><br>
> SRID: 4269<br>
> SRTEXT: GEOGCS["NAD83",DATUM["North_American_Datum_1983",<br>
> SPHEROID["GRS 1980",6378137,298.257222101,<br>
> AUTHORITY["EPSG","7019"]],<br>
> AUTHORITY["EPSG","6269"]],<br>
> PRIMEM["Greenwich",0,<br>
> AUTHORITY["EPSG","8901"]],<br>
> UNIT["degree",0.01745329251994328,<br>
> AUTHORITY["EPSG","9122"]],<br>
> AUTHORITY["EPSG","4269"]]<br>
> PROJ4TEXT: +proj=longlat +ellps=GRS80 +datum=NAD83 +no_defs<br>
><br>
> SRID: 32767<br>
> SRTEXT: GEOGCS["GCS_North_American_1983",<br>
> DATUM["D_North_American_1983",<br>
> SPHEROID["GRS_1980",6378137,298.257222101]],<br>
> PRIMEM["Greenwich",0],<br>
> UNIT["Degree",0.017453292519943295]]<br>
> PROJ4TEXT: +proj=longlat +ellps=clrk66 +datum=NAD27 +no_defs<br>
><br>
> To display census block data in OpenStreetMap, I extract it from PostGIS<br>
> with a transform to EPSG 4326, although the coordinates don't seem to<br>
> change as a result. (This seems correct, as datum=NAD83 and datum=WGS84<br>
> are, for my purposes at least, are essentially identical.)<br>
><br>
> Thanks,<br>
> Jonathan<br>
><br>
> 2 attachments: TIGER2007andOSM.png, TIGER2007andGoogle.png<br>
><br>
><br>
> On Thu, 2008-04-03 at 19:18 -0400, Stephen Frost wrote:<br>
> > Jonathan,<br>
> ><br>
> > * Jonathan W. Lowe (<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>) wrote:<br>
> > > Have you yet tried overlaying TIGER 2007 linework or census block/tract<br>
> > > polygons over Google or OpenStreetMap tiles? I'm seeing a good match in<br>
> > > some areas but a significant shift (~50 meters) in others. Thought it<br>
> > > might be a datum conversion issue, but can't seem to find a match.<br>
> ><br>
> > I hadn't looked at the linework too much yet or tried to overlay it.<br>
> > I'm curious where you're seeing the differences though because I know<br>
> > that Census is only about half way through their MAF improvment project<br>
> > and I actually have some info about what has been done so far and what<br>
> > hasn't. It'd be interesting to see if it matches up.<br>
> ><br>
> > There are a few places (Guam, Hawaii islands) where they actually do use<br>
> > an SRID other than 4269, but my scripts don't yet handle that and I'm<br>
> > guessing that's not what you're referring to anyway. :)<br>
> ><br>
> > Thanks!<br>
> ><br>
> > Stephen<br>
> ><br>
> > > On Thu, 2008-04-03 at 17:07 -0400, Stephen Frost wrote:<br>
> > > > * Stephen Frost (<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>) wrote:<br>
> > > > > I think they may have also upgraded their pipe.. I got about 1.41MB/s<br>
> > > > > (11 Mb/s) for the whole transfer. It's about 22G all told. I'll<br>
> > > > > probably be trying to load it up into PG on one of our servers tomorrow.<br>
> > > > > It was a bit over 4 hours for me to pull down off of their<br>
> > > > > <a href="http://ftp2.census.gov" target="_blank">ftp2.census.gov</a> ftp site.<br>
> > > ><br>
> > > > Just to update those who might be interested- I've finished the data<br>
> > > > load into one of our servers at work. It comes to ~60GB on disk in<br>
> > > > PostgreSQL/PostGIS with appropriate indexes in most places and whatnot.<br>
> > > > Based on what I've seen so far, it looks *very* nice, especially the<br>
> > > > hydrogrophy ("areawater"). It also appears to be pretty consistant<br>
> > > > across the layers, which is also good.<br>
> > > ><br>
> > > > If anyone's interested in the scripts used to load the data (they're<br>
> > > > pretty simple, really), I'd be happy to provide them.<br>
> > > ><br>
> > > > Enjoy,<br>
> > > ><br>
> > > > Stephen<br>
> > > > _______________________________________________<br>
> > > > postgis-users mailing list<br>
> > > > <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> > > > <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
> > ><br>
> > > _______________________________________________<br>
> > > postgis-users mailing list<br>
> > > <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> > > <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 20<br>
Date: Thu, 3 Apr 2008 20:31:32 -0400<br>
From: Stephen Frost <<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: "Jonathan W. Lowe" <<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>>, PostGIS Users<br>
Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:20080404003132.GA4999@tamriel.snowman.net">20080404003132.GA4999@tamriel.snowman.net</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Jonathan,<br>
<br>
* Jonathan W. Lowe (<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>) wrote:<br>
> ...And in case the images don't persist through the mail server, they're<br>
> viewable at: <a href="http://www.giswebsite.com/demos/tiger_overlays.html" target="_blank">http://www.giswebsite.com/demos/tiger_overlays.html</a><br>
<br>
You know, I just realized that you were talking about the Census 2000<br>
blocks (tabblock00.shp). Is there some reason you're using that<br>
instead of the current data (tabblock.shp)? They might not want to<br>
update the data from 2000 for historical reasons...<br>
(Note: I havn't actually gone and looked, it just occured to me..)<br>
<br>
Thanks,<br>
<br>
Stephen<br>
<br>
> On Fri, 2008-04-04 at 01:07 +0100, Jonathan W. Lowe wrote:<br>
> > Stephen,<br>
> ><br>
> > My initial testing has been on Alameda County (California) TIGER data.<br>
> > The two attached image files show an overlay of US Census 2000 Blocks<br>
> > over an area south of the UC Berkeley campus. The offset is the same<br>
> > for both Google and OpenStreetMap (OSM). This suggests that I've made a<br>
> > mistake somewhere, because the OSM tiles in the United States are all<br>
> > rendered from TIGER linework, so the TIGER census blocks should match<br>
> > exactly.<br>
> ><br>
> > For the same source shapefile (tabblock00.shp), there's a nearly perfect<br>
> > match between block boundaries and streets in the area just South of<br>
> > Oakland's Lake Merritt. It smells like a datum conversion issue...<br>
> ><br>
> > The conversion path was from shapefile to PostGIS using shp2pgsql. I<br>
> > used a custom projection of 32767 rather than 4269 because the existing<br>
> > srtext for 4269 had a degree value as 0.01745329251994328, but the US<br>
> > Census metadata listed a degree value of 0.017453292519943295. Perhaps<br>
> > not significant? My spatial_ref_sys entries for 4269 and 32767 are<br>
> > otherwise pretty similar:<br>
> ><br>
> > SRID: 4269<br>
> > SRTEXT: GEOGCS["NAD83",DATUM["North_American_Datum_1983",<br>
> > SPHEROID["GRS 1980",6378137,298.257222101,<br>
> > AUTHORITY["EPSG","7019"]],<br>
> > AUTHORITY["EPSG","6269"]],<br>
> > PRIMEM["Greenwich",0,<br>
> > AUTHORITY["EPSG","8901"]],<br>
> > UNIT["degree",0.01745329251994328,<br>
> > AUTHORITY["EPSG","9122"]],<br>
> > AUTHORITY["EPSG","4269"]]<br>
> > PROJ4TEXT: +proj=longlat +ellps=GRS80 +datum=NAD83 +no_defs<br>
> ><br>
> > SRID: 32767<br>
> > SRTEXT: GEOGCS["GCS_North_American_1983",<br>
> > DATUM["D_North_American_1983",<br>
> > SPHEROID["GRS_1980",6378137,298.257222101]],<br>
> > PRIMEM["Greenwich",0],<br>
> > UNIT["Degree",0.017453292519943295]]<br>
> > PROJ4TEXT: +proj=longlat +ellps=clrk66 +datum=NAD27 +no_defs<br>
> ><br>
> > To display census block data in OpenStreetMap, I extract it from PostGIS<br>
> > with a transform to EPSG 4326, although the coordinates don't seem to<br>
> > change as a result. (This seems correct, as datum=NAD83 and datum=WGS84<br>
> > are, for my purposes at least, are essentially identical.)<br>
> ><br>
> > Thanks,<br>
> > Jonathan<br>
> ><br>
> > 2 attachments: TIGER2007andOSM.png, TIGER2007andGoogle.png<br>
> ><br>
> ><br>
> > On Thu, 2008-04-03 at 19:18 -0400, Stephen Frost wrote:<br>
> > > Jonathan,<br>
> > ><br>
> > > * Jonathan W. Lowe (<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>) wrote:<br>
> > > > Have you yet tried overlaying TIGER 2007 linework or census block/tract<br>
> > > > polygons over Google or OpenStreetMap tiles? I'm seeing a good match in<br>
> > > > some areas but a significant shift (~50 meters) in others. Thought it<br>
> > > > might be a datum conversion issue, but can't seem to find a match.<br>
> > ><br>
> > > I hadn't looked at the linework too much yet or tried to overlay it.<br>
> > > I'm curious where you're seeing the differences though because I know<br>
> > > that Census is only about half way through their MAF improvment project<br>
> > > and I actually have some info about what has been done so far and what<br>
> > > hasn't. It'd be interesting to see if it matches up.<br>
> > ><br>
> > > There are a few places (Guam, Hawaii islands) where they actually do use<br>
> > > an SRID other than 4269, but my scripts don't yet handle that and I'm<br>
> > > guessing that's not what you're referring to anyway. :)<br>
> > ><br>
> > > Thanks!<br>
> > ><br>
> > > Stephen<br>
> > ><br>
> > > > On Thu, 2008-04-03 at 17:07 -0400, Stephen Frost wrote:<br>
> > > > > * Stephen Frost (<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>) wrote:<br>
> > > > > > I think they may have also upgraded their pipe.. I got about 1.41MB/s<br>
> > > > > > (11 Mb/s) for the whole transfer. It's about 22G all told. I'll<br>
> > > > > > probably be trying to load it up into PG on one of our servers tomorrow.<br>
> > > > > > It was a bit over 4 hours for me to pull down off of their<br>
> > > > > > <a href="http://ftp2.census.gov" target="_blank">ftp2.census.gov</a> ftp site.<br>
> > > > ><br>
> > > > > Just to update those who might be interested- I've finished the data<br>
> > > > > load into one of our servers at work. It comes to ~60GB on disk in<br>
> > > > > PostgreSQL/PostGIS with appropriate indexes in most places and whatnot.<br>
> > > > > Based on what I've seen so far, it looks *very* nice, especially the<br>
> > > > > hydrogrophy ("areawater"). It also appears to be pretty consistant<br>
> > > > > across the layers, which is also good.<br>
> > > > ><br>
> > > > > If anyone's interested in the scripts used to load the data (they're<br>
> > > > > pretty simple, really), I'd be happy to provide them.<br>
> > > > ><br>
> > > > > Enjoy,<br>
> > > > ><br>
> > > > > Stephen<br>
> > > > > _______________________________________________<br>
> > > > > postgis-users mailing list<br>
> > > > > <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> > > > > <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
> > > ><br>
> > > > _______________________________________________<br>
> > > > postgis-users mailing list<br>
> > > > <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> > > > <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
><br>
> _______________________________________________<br>
> postgis-users mailing list<br>
> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: not available<br>
Type: application/pgp-signature<br>
Size: 189 bytes<br>
Desc: Digital signature<br>
Url : <a href="http://lists.refractions.net/pipermail/postgis-users/attachments/20080403/dac89b9a/attachment-0001.bin" target="_blank">http://lists.refractions.net/pipermail/postgis-users/attachments/20080403/dac89b9a/attachment-0001.bin</a><br>
<br>
------------------------------<br>
<br>
Message: 21<br>
Date: Thu, 03 Apr 2008 20:48:49 -0500<br>
From: Stephen Woodbridge <<a href="mailto:woodbri@swoodbridge.com">woodbri@swoodbridge.com</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: "Jonathan W. Lowe" <<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>>, PostGIS Users<br>
Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:47F58901.3050701@swoodbridge.com">47F58901.3050701@swoodbridge.com</a>><br>
Content-Type: text/plain; charset=UTF-8; format=flowed<br>
<br>
Stephen Frost wrote:<br>
> Jonathan,<br>
><br>
> * Jonathan W. Lowe (<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>) wrote:<br>
>> Have you yet tried overlaying TIGER 2007 linework or census block/tract<br>
>> polygons over Google or OpenStreetMap tiles? I'm seeing a good match in<br>
>> some areas but a significant shift (~50 meters) in others. Thought it<br>
>> might be a datum conversion issue, but can't seem to find a match.<br>
><br>
> I hadn't looked at the linework too much yet or tried to overlay it.<br>
> I'm curious where you're seeing the differences though because I know<br>
> that Census is only about half way through their MAF improvment project<br>
> and I actually have some info about what has been done so far and what<br>
> hasn't. It'd be interesting to see if it matches up.<br>
><br>
> There are a few places (Guam, Hawaii islands) where they actually do use<br>
> an SRID other than 4269, but my scripts don't yet handle that and I'm<br>
> guessing that's not what you're referring to anyway. :)<br>
<br>
I had an extensive discussion with some of the Census Geography staff<br>
about the Island provinces and Hawaii because I was trying to align<br>
Navteq routes over Tiger data and the Navteq routes were about a .25-.5<br>
miles east of the Tiger data in Hawaii.<br>
<br>
It turns out that they really have not idea what the "local datums" are<br>
that were originally used. I tried a lot of the local Island datums in<br>
the proj4 epsg file but could not find any close matches.<br>
<br>
It will be nice when they get that fixed up :)<br>
<br>
Stephen - thank you for sharing your scripts.<br>
<br>
Best regards,<br>
-Stephen Woodbridge<br>
<a href="http://imaptools.com/" target="_blank">http://imaptools.com/</a><br>
<br>
> Thanks!<br>
><br>
> Stephen<br>
><br>
>> On Thu, 2008-04-03 at 17:07 -0400, Stephen Frost wrote:<br>
>>> * Stephen Frost (<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>) wrote:<br>
>>>> I think they may have also upgraded their pipe.. I got about 1.41MB/s<br>
>>>> (11 Mb/s) for the whole transfer. It's about 22G all told. I'll<br>
>>>> probably be trying to load it up into PG on one of our servers tomorrow.<br>
>>>> It was a bit over 4 hours for me to pull down off of their<br>
>>>> <a href="http://ftp2.census.gov" target="_blank">ftp2.census.gov</a> ftp site.<br>
>>> Just to update those who might be interested- I've finished the data<br>
>>> load into one of our servers at work. It comes to ~60GB on disk in<br>
>>> PostgreSQL/PostGIS with appropriate indexes in most places and whatnot.<br>
>>> Based on what I've seen so far, it looks *very* nice, especially the<br>
>>> hydrogrophy ("areawater"). It also appears to be pretty consistant<br>
>>> across the layers, which is also good.<br>
>>><br>
>>> If anyone's interested in the scripts used to load the data (they're<br>
>>> pretty simple, really), I'd be happy to provide them.<br>
>>><br>
>>> Enjoy,<br>
>>><br>
>>> Stephen<br>
>>> _______________________________________________<br>
>>> postgis-users mailing list<br>
>>> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
>> _______________________________________________<br>
>> postgis-users mailing list<br>
>> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
>><br>
>> ------------------------------------------------------------------------<br>
>><br>
>> _______________________________________________<br>
>> postgis-users mailing list<br>
>> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 22<br>
Date: Thu, 03 Apr 2008 20:52:50 -0500<br>
From: Stephen Woodbridge <<a href="mailto:woodbri@swoodbridge.com">woodbri@swoodbridge.com</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:47F589F2.4060804@swoodbridge.com">47F589F2.4060804@swoodbridge.com</a>><br>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>
<br>
Yes, it would be ideal if it was handled by shp2pgsql but have you<br>
looked at:<br>
<br>
<a href="http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2pgsql&btnG=Search" target="_blank">http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2pgsql&btnG=Search</a><br>
<a href="http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2psql&btnG=Search" target="_blank">http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2psql&btnG=Search</a><br>
<a href="http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2sql&btnG=Search" target="_blank">http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2sql&btnG=Search</a><br>
<br>
some of these might be helpful.<br>
<br>
Best regards,<br>
-Stephen Woodbridge<br>
<a href="http://imaptools.com/" target="_blank">http://imaptools.com/</a><br>
<br>
Paragon Corporation wrote:<br>
> As a side note to what you were saying in the .sh file it would be really<br>
> nice if the shp2pgsql dealt with standalone dbf files too. I mean the logic<br>
> is all there so doesn't seem like it would be that hard to put in a switch<br>
> for that.<br>
><br>
> For the dbfs I was using Ogr2Ogr which works well except it adds a useless<br>
> geometry field.<br>
><br>
><br>
><br>
> -----Original Message-----<br>
> From: <a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a><br>
> [mailto:<a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a>] On Behalf Of Stephen<br>
> Frost<br>
> Sent: Thursday, April 03, 2008 7:16 PM<br>
> To: PostGIS Users Discussion<br>
> Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
><br>
> * Stephen Frost (<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>) wrote:<br>
>> If anyone's interested in the scripts used to load the data (they're<br>
>> pretty simple, really), I'd be happy to provide them.<br>
><br>
> Alright, since it's apparently kind of popular, I went ahead and put the<br>
> script up here:<br>
><br>
> <a href="http://snowman.net/tiger/" target="_blank">http://snowman.net/tiger/</a><br>
><br>
> I also included a README.txt that is essentially what I wrote up to the<br>
> first person who asked me for the script. Please let me know if you improve<br>
> upon it or find problems with it.<br>
><br>
> Thanks,<br>
><br>
> Stephen<br>
><br>
><br>
> _______________________________________________<br>
> postgis-users mailing list<br>
> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 23<br>
Date: Thu, 3 Apr 2008 20:59:43 -0400<br>
From: Stephen Frost <<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:20080404005943.GB4999@tamriel.snowman.net">20080404005943.GB4999@tamriel.snowman.net</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Stephen,<br>
<br>
* Stephen Woodbridge (<a href="mailto:woodbri@swoodbridge.com">woodbri@swoodbridge.com</a>) wrote:<br>
> I had an extensive discussion with some of the Census Geography staff<br>
> about the Island provinces and Hawaii because I was trying to align<br>
> Navteq routes over Tiger data and the Navteq routes were about a .25-.5<br>
> miles east of the Tiger data in Hawaii.<br>
<br>
Interesting...<br>
<br>
> It turns out that they really have not idea what the "local datums" are<br>
> that were originally used. I tried a lot of the local Island datums in<br>
> the proj4 epsg file but could not find any close matches.<br>
<br>
In their documentation they talk about the datums used for these<br>
regions.. Lemme go look up what they say..<br>
<br>
Here we go:<br>
<br>
North American Datum of 1983 in the 48 contiguous states, the District<br>
of Columbia, Alaska, Hawaii (only Oahu Island within Honolulu County, HI<br>
(15003)), Puerto Rico, and the U.S. Virgin Islands.<br>
<br>
Regional datums are used in the Pacific Island Areas (American Samoa,<br>
Guam, and the Commonwealth of the Northern Mariana Islands) and most of<br>
Hawaii.<br>
<br>
The datums used in the remainder of Hawaii are as follows: Hawaii<br>
County, HI (15001), Old Hawaiian Datum; Honolulu County, HI (15003),<br>
local astronomic datums for all islands northwest of 161 degrees west<br>
longitude; Kalawao County, HI (15005), Old Hawaiian Datum; Kauai County,<br>
HI (15007), Old Hawaiian Datum for Kauai Island and local astronomic<br>
datum for Kaula Rock; Maui County, HI (15009), Old Hawaiian Datum.<br>
<br>
The datums used in American Samoa are as follows: Eastern District, AS<br>
(60010), American Samoa Datum of 1962; Manu'a District, AS (60020),<br>
American Samoa Datum of 1962; Rose Island, AS (60030), local astronomic<br>
datum; Swains Island, AS (60040), local astronomic datum 1939; Western<br>
District, AS (60050), American Samoa Datum of 1962.<br>
<br>
The Guam Datum of 1963 is used in Guam (66010). The datums used in the<br>
Commonwealth of the Northern Mariana Islands are as follows: Northern<br>
Islands Municipality (69085), Guam Datum of 1963 (Agrihan, Alamagan,<br>
Anatahan, Gugan, Medinilla, Pagan, and Sarigan Islands) and local<br>
astronomic datums (Asuncion, Maug, and Farallon De Pajaros Islands);<br>
Rota Municipality (69100), Guam Datum of 1963; Saipan Municipality<br>
(69110), Guam Datum of 1963; Tinian Municipality (69120), Guam Datum of<br>
1963.<br>
<br>
So, there's somewhat more information there than perhaps there was<br>
previously, but I havn't gone hunting for, eg, "Old Hawaiian Datum". If<br>
someone does come up with the correct datums/srids/etc based off of this<br>
(or whatever), I'd be happy to update my script accordingly.<br>
<br>
Thanks!<br>
<br>
Stephen<br>
<br>
> It will be nice when they get that fixed up :)<br>
><br>
> Stephen - thank you for sharing your scripts.<br>
><br>
> Best regards,<br>
> -Stephen Woodbridge<br>
> <a href="http://imaptools.com/" target="_blank">http://imaptools.com/</a><br>
><br>
>> Thanks!<br>
>><br>
>> Stephen<br>
>><br>
>>> On Thu, 2008-04-03 at 17:07 -0400, Stephen Frost wrote:<br>
>>>> * Stephen Frost (<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>) wrote:<br>
>>>>> I think they may have also upgraded their pipe.. I got about 1.41MB/s<br>
>>>>> (11 Mb/s) for the whole transfer. It's about 22G all told. I'll<br>
>>>>> probably be trying to load it up into PG on one of our servers tomorrow.<br>
>>>>> It was a bit over 4 hours for me to pull down off of their<br>
>>>>> <a href="http://ftp2.census.gov" target="_blank">ftp2.census.gov</a> ftp site.<br>
>>>> Just to update those who might be interested- I've finished the data<br>
>>>> load into one of our servers at work. It comes to ~60GB on disk in<br>
>>>> PostgreSQL/PostGIS with appropriate indexes in most places and whatnot.<br>
>>>> Based on what I've seen so far, it looks *very* nice, especially the<br>
>>>> hydrogrophy ("areawater"). It also appears to be pretty consistant<br>
>>>> across the layers, which is also good.<br>
>>>><br>
>>>> If anyone's interested in the scripts used to load the data (they're<br>
>>>> pretty simple, really), I'd be happy to provide them.<br>
>>>><br>
>>>> Enjoy,<br>
>>>><br>
>>>> Stephen<br>
>>>> _______________________________________________<br>
>>>> postgis-users mailing list<br>
>>>> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
>>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
>>> _______________________________________________<br>
>>> postgis-users mailing list<br>
>>> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
>>><br>
>>> ------------------------------------------------------------------------<br>
>>><br>
>>> _______________________________________________<br>
>>> postgis-users mailing list<br>
>>> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
><br>
> _______________________________________________<br>
> postgis-users mailing list<br>
> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: not available<br>
Type: application/pgp-signature<br>
Size: 189 bytes<br>
Desc: Digital signature<br>
Url : <a href="http://lists.refractions.net/pipermail/postgis-users/attachments/20080403/cd560eac/attachment-0001.bin" target="_blank">http://lists.refractions.net/pipermail/postgis-users/attachments/20080403/cd560eac/attachment-0001.bin</a><br>
<br>
------------------------------<br>
<br>
Message: 24<br>
Date: Thu, 3 Apr 2008 21:01:09 -0400<br>
From: Stephen Frost <<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:20080404010108.GC4999@tamriel.snowman.net">20080404010108.GC4999@tamriel.snowman.net</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
* Stephen Woodbridge (<a href="mailto:woodbri@swoodbridge.com">woodbri@swoodbridge.com</a>) wrote:<br>
> Yes, it would be ideal if it was handled by shp2pgsql but have you<br>
> looked at:<br>
><br>
> <a href="http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2pgsql&btnG=Search" target="_blank">http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2pgsql&btnG=Search</a><br>
> <a href="http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2psql&btnG=Search" target="_blank">http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2psql&btnG=Search</a><br>
> <a href="http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2sql&btnG=Search" target="_blank">http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2sql&btnG=Search</a><br>
><br>
> some of these might be helpful.<br>
<br>
Just fyi, my script uses dbview from the, conveniantly named, dbview<br>
Debian package. It works quite well for me. Not that I'd be against<br>
having that functionality in shp2pgsql, it's certainly be useful.<br>
<br>
First I'd like to see an option to add the primary key *after* the data<br>
is loaded though.<br>
<br>
Thanks,<br>
<br>
Stephen<br>
<br>
> Paragon Corporation wrote:<br>
>> As a side note to what you were saying in the .sh file it would be really<br>
>> nice if the shp2pgsql dealt with standalone dbf files too. I mean the logic<br>
>> is all there so doesn't seem like it would be that hard to put in a switch<br>
>> for that.<br>
>><br>
>> For the dbfs I was using Ogr2Ogr which works well except it adds a useless<br>
>> geometry field.<br>
>><br>
>><br>
>><br>
>> -----Original Message-----<br>
>> From: <a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a><br>
>> [mailto:<a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a>] On Behalf Of Stephen<br>
>> Frost<br>
>> Sent: Thursday, April 03, 2008 7:16 PM<br>
>> To: PostGIS Users Discussion<br>
>> Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
>><br>
>> * Stephen Frost (<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>) wrote:<br>
>>> If anyone's interested in the scripts used to load the data (they're<br>
>>> pretty simple, really), I'd be happy to provide them.<br>
>><br>
>> Alright, since it's apparently kind of popular, I went ahead and put the<br>
>> script up here:<br>
>><br>
>> <a href="http://snowman.net/tiger/" target="_blank">http://snowman.net/tiger/</a><br>
>><br>
>> I also included a README.txt that is essentially what I wrote up to the<br>
>> first person who asked me for the script. Please let me know if you improve<br>
>> upon it or find problems with it.<br>
>><br>
>> Thanks,<br>
>><br>
>> Stephen<br>
>><br>
>><br>
>> _______________________________________________<br>
>> postgis-users mailing list<br>
>> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
><br>
> _______________________________________________<br>
> postgis-users mailing list<br>
> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: not available<br>
Type: application/pgp-signature<br>
Size: 189 bytes<br>
Desc: Digital signature<br>
Url : <a href="http://lists.refractions.net/pipermail/postgis-users/attachments/20080403/366a2d47/attachment-0001.bin" target="_blank">http://lists.refractions.net/pipermail/postgis-users/attachments/20080403/366a2d47/attachment-0001.bin</a><br>
<br>
------------------------------<br>
<br>
Message: 25<br>
Date: Thu, 3 Apr 2008 21:10:17 -0400<br>
From: "Paragon Corporation" <<a href="mailto:lr@pcorp.us">lr@pcorp.us</a>><br>
Subject: RE: [postgis-users] TIGER/Line Shapefiles released<br>
To: "'PostGIS Users Discussion'"<br>
<<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <00a801c895f0$a57ef630$4d812e40@l><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Thanks I'll take a look. I was thinking more for completeness and the<br>
marginal effort that I think would be involved from my memory of what the<br>
shp2pgsql C code looks like. I would change it myself, but I'm moderately<br>
afraid of C code.<br>
<br>
-----Original Message-----<br>
From: <a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a><br>
[mailto:<a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a>] On Behalf Of Stephen<br>
Woodbridge<br>
Sent: Thursday, April 03, 2008 9:53 PM<br>
To: PostGIS Users Discussion<br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
<br>
Yes, it would be ideal if it was handled by shp2pgsql but have you looked<br>
at:<br>
<br>
<a href="http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2pgsql&btnG=Search" target="_blank">http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2pgsql&<br>
btnG=Search</a><br>
<a href="http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2psql&btnG=Search" target="_blank">http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2psql&b<br>
tnG=Search</a><br>
<a href="http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2sql&btnG=Search" target="_blank">http://www.google.com/search?num=100&hl=en&newwindow=1&safe=off&q=dbf2sql&bt<br>
nG=Search</a><br>
<br>
some of these might be helpful.<br>
<br>
Best regards,<br>
-Stephen Woodbridge<br>
<a href="http://imaptools.com/" target="_blank">http://imaptools.com/</a><br>
<br>
Paragon Corporation wrote:<br>
> As a side note to what you were saying in the .sh file it would be<br>
> really nice if the shp2pgsql dealt with standalone dbf files too. I<br>
> mean the logic is all there so doesn't seem like it would be that hard<br>
> to put in a switch for that.<br>
><br>
> For the dbfs I was using Ogr2Ogr which works well except it adds a<br>
> useless geometry field.<br>
><br>
><br>
><br>
> -----Original Message-----<br>
> From: <a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a><br>
> [mailto:<a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a>] On Behalf Of<br>
> Stephen Frost<br>
> Sent: Thursday, April 03, 2008 7:16 PM<br>
> To: PostGIS Users Discussion<br>
> Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
><br>
> * Stephen Frost (<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>) wrote:<br>
>> If anyone's interested in the scripts used to load the data (they're<br>
>> pretty simple, really), I'd be happy to provide them.<br>
><br>
> Alright, since it's apparently kind of popular, I went ahead and put<br>
> the script up here:<br>
><br>
> <a href="http://snowman.net/tiger/" target="_blank">http://snowman.net/tiger/</a><br>
><br>
> I also included a README.txt that is essentially what I wrote up to<br>
> the first person who asked me for the script. Please let me know if<br>
> you improve upon it or find problems with it.<br>
><br>
> Thanks,<br>
><br>
> Stephen<br>
><br>
><br>
> _______________________________________________<br>
> postgis-users mailing list<br>
> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br>
_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
<a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 26<br>
Date: Thu, 03 Apr 2008 21:48:38 -0500<br>
From: Stephen Woodbridge <<a href="mailto:woodbri@swoodbridge.com">woodbri@swoodbridge.com</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:47F59706.6090409@swoodbridge.com">47F59706.6090409@swoodbridge.com</a>><br>
Content-Type: text/plain; charset=UTF-8; format=flowed<br>
<br>
Stephen Frost wrote:<br>
> Stephen,<br>
><br>
> * Stephen Woodbridge (<a href="mailto:woodbri@swoodbridge.com">woodbri@swoodbridge.com</a>) wrote:<br>
>> I had an extensive discussion with some of the Census Geography staff<br>
>> about the Island provinces and Hawaii because I was trying to align<br>
>> Navteq routes over Tiger data and the Navteq routes were about a .25-.5<br>
>> miles east of the Tiger data in Hawaii.<br>
><br>
> Interesting...<br>
><br>
>> It turns out that they really have not idea what the "local datums" are<br>
>> that were originally used. I tried a lot of the local Island datums in<br>
>> the proj4 epsg file but could not find any close matches.<br>
><br>
> In their documentation they talk about the datums used for these<br>
> regions.. Lemme go look up what they say..<br>
><br>
> Here we go:<br>
><br>
> North American Datum of 1983 in the 48 contiguous states, the District<br>
> of Columbia, Alaska, Hawaii (only Oahu Island within Honolulu County, HI<br>
> (15003)), Puerto Rico, and the U.S. Virgin Islands.<br>
><br>
> Regional datums are used in the Pacific Island Areas (American Samoa,<br>
> Guam, and the Commonwealth of the Northern Mariana Islands) and most of<br>
> Hawaii.<br>
><br>
> The datums used in the remainder of Hawaii are as follows: Hawaii<br>
> County, HI (15001), Old Hawaiian Datum; Honolulu County, HI (15003),<br>
> local astronomic datums for all islands northwest of 161 degrees west<br>
> longitude; Kalawao County, HI (15005), Old Hawaiian Datum; Kauai County,<br>
> HI (15007), Old Hawaiian Datum for Kauai Island and local astronomic<br>
> datum for Kaula Rock; Maui County, HI (15009), Old Hawaiian Datum.<br>
><br>
> The datums used in American Samoa are as follows: Eastern District, AS<br>
> (60010), American Samoa Datum of 1962; Manu'a District, AS (60020),<br>
> American Samoa Datum of 1962; Rose Island, AS (60030), local astronomic<br>
> datum; Swains Island, AS (60040), local astronomic datum 1939; Western<br>
> District, AS (60050), American Samoa Datum of 1962.<br>
><br>
> The Guam Datum of 1963 is used in Guam (66010). The datums used in the<br>
> Commonwealth of the Northern Mariana Islands are as follows: Northern<br>
> Islands Municipality (69085), Guam Datum of 1963 (Agrihan, Alamagan,<br>
> Anatahan, Gugan, Medinilla, Pagan, and Sarigan Islands) and local<br>
> astronomic datums (Asuncion, Maug, and Farallon De Pajaros Islands);<br>
> Rota Municipality (69100), Guam Datum of 1963; Saipan Municipality<br>
> (69110), Guam Datum of 1963; Tinian Municipality (69120), Guam Datum of<br>
> 1963.<br>
><br>
> So, there's somewhat more information there than perhaps there was<br>
> previously, but I havn't gone hunting for, eg, "Old Hawaiian Datum". If<br>
> someone does come up with the correct datums/srids/etc based off of this<br>
> (or whatever), I'd be happy to update my script accordingly.<br>
<br>
Yeah, that is pretty much what they finally offered to me after a lot of<br>
research. I think some of this is newer info than I originally got. I<br>
would have to go digging for my notes, and this info about is probably<br>
more accurate at this point.<br>
<br>
I was primarily focused on the problem I had with the route alignments<br>
that we resolved by using the Navteq data for the maps when we were<br>
displaying routes.<br>
<br>
From proj4 epsg:<br>
<br>
# Old Hawaiian<br>
<4135> +proj=longlat +ellps=clrk66 +no_defs <><br>
<br>
# American Samoa 1962<br>
<4169> +proj=longlat +ellps=clrk66 +towgs84=-115,118,426,0,0,0,0<br>
+no_defs <><br>
<br>
From proj4 esri<br>
<br>
# GCS Guam 1963<br>
<37220> +proj=longlat +ellps=clrk66 no_defs <><br>
<br>
Best regards,<br>
-Stephen Woodbridge<br>
<a href="http://imaptools.com/" target="_blank">http://imaptools.com/</a><br>
<br>
> Thanks!<br>
><br>
> Stephen<br>
><br>
>> It will be nice when they get that fixed up :)<br>
>><br>
>> Stephen - thank you for sharing your scripts.<br>
>><br>
>> Best regards,<br>
>> -Stephen Woodbridge<br>
>> <a href="http://imaptools.com/" target="_blank">http://imaptools.com/</a><br>
>><br>
>>> Thanks!<br>
>>><br>
>>> Stephen<br>
>>><br>
>>>> On Thu, 2008-04-03 at 17:07 -0400, Stephen Frost wrote:<br>
>>>>> * Stephen Frost (<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>) wrote:<br>
>>>>>> I think they may have also upgraded their pipe.. I got about 1.41MB/s<br>
>>>>>> (11 Mb/s) for the whole transfer. It's about 22G all told. I'll<br>
>>>>>> probably be trying to load it up into PG on one of our servers tomorrow.<br>
>>>>>> It was a bit over 4 hours for me to pull down off of their<br>
>>>>>> <a href="http://ftp2.census.gov" target="_blank">ftp2.census.gov</a> ftp site.<br>
>>>>> Just to update those who might be interested- I've finished the data<br>
>>>>> load into one of our servers at work. It comes to ~60GB on disk in<br>
>>>>> PostgreSQL/PostGIS with appropriate indexes in most places and whatnot.<br>
>>>>> Based on what I've seen so far, it looks *very* nice, especially the<br>
>>>>> hydrogrophy ("areawater"). It also appears to be pretty consistant<br>
>>>>> across the layers, which is also good.<br>
>>>>><br>
>>>>> If anyone's interested in the scripts used to load the data (they're<br>
>>>>> pretty simple, really), I'd be happy to provide them.<br>
>>>>><br>
>>>>> Enjoy,<br>
>>>>><br>
>>>>> Stephen<br>
>>>>> _______________________________________________<br>
>>>>> postgis-users mailing list<br>
>>>>> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
>>>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
>>>> _______________________________________________<br>
>>>> postgis-users mailing list<br>
>>>> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
>>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
>>>><br>
>>>> ------------------------------------------------------------------------<br>
>>>><br>
>>>> _______________________________________________<br>
>>>> postgis-users mailing list<br>
>>>> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
>>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
>> _______________________________________________<br>
>> postgis-users mailing list<br>
>> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
>><br>
>> ------------------------------------------------------------------------<br>
>><br>
>> _______________________________________________<br>
>> postgis-users mailing list<br>
>> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 27<br>
Date: Thu, 03 Apr 2008 22:40:18 -0500<br>
From: Stephen Woodbridge <<a href="mailto:woodbri@swoodbridge.com">woodbri@swoodbridge.com</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:47F5A322.3050802@swoodbridge.com">47F5A322.3050802@swoodbridge.com</a>><br>
Content-Type: text/plain; charset=UTF-8; format=flowed<br>
<br>
Stephen Frost wrote:<br>
> * Stephen Frost (<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>) wrote:<br>
>> If anyone's interested in the scripts used to load the data (they're<br>
>> pretty simple, really), I'd be happy to provide them.<br>
><br>
> Alright, since it's apparently kind of popular, I went ahead and put the<br>
> script up here:<br>
><br>
> <a href="http://snowman.net/tiger/" target="_blank">http://snowman.net/tiger/</a><br>
><br>
> I also included a README.txt that is essentially what I wrote up to the<br>
> first person who asked me for the script. Please let me know if you<br>
> improve upon it or find problems with it.<br>
<br>
Stephen,<br>
<br>
I added to process_tiger.sh:<br>
<br>
# User to connect to database with or use USER=$LOGNAME<br>
USER="postgres"<br>
<br>
and in vi did:<br>
<br>
%s/\<psql\>/psql -U $USER/g<br>
<br>
to update the script.<br>
<br>
Question on the README.txt<br>
<br>
-n Load national-level data<br>
-b Load state-level data<br>
-c Load county-level data<br>
<br>
Do the above options load all the data at the given level or do I also<br>
need to run it for each state and county with the options below?<br>
<br>
-s <2-digit state code> Individual state code to load<br>
-c <3-digit or 5-digit county code> Individual county code to load<br>
<br>
Does:<br>
<br>
cd /u/srcdata/tiger2007fe<br>
./process_tiger.sh -n -b -c<br>
<br>
load all the data? Assuming that that directory is a mirror of the<br>
Census download directory.<br>
<br>
Thanks,<br>
-Stephen Woodbridge<br>
<a href="http://imaptools.com/" target="_blank">http://imaptools.com/</a><br>
<br>
<br>
------------------------------<br>
<br>
Message: 28<br>
Date: Thu, 3 Apr 2008 22:54:08 -0400<br>
From: Stephen Frost <<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:20080404025408.GD4999@tamriel.snowman.net">20080404025408.GD4999@tamriel.snowman.net</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Stephen,<br>
<br>
* Stephen Woodbridge (<a href="mailto:woodbri@swoodbridge.com">woodbri@swoodbridge.com</a>) wrote:<br>
> I added to process_tiger.sh:<br>
><br>
> # User to connect to database with or use USER=$LOGNAME<br>
> USER="postgres"<br>
<br>
Ah, sure, makes sense. We user Kerberos for our authentication, so I<br>
tend to forget about user/pw authentication issues.<br>
<br>
> and in vi did:<br>
><br>
> %s/\<psql\>/psql -U $USER/g<br>
><br>
> to update the script.<br>
<br>
All of the options passed to psql really should be done as options to<br>
the script. I'll look at adding that (and -U support) sometime<br>
tomorrow.<br>
<br>
> Question on the README.txt<br>
><br>
> -n Load national-level data<br>
> -b Load state-level data<br>
> -c Load county-level data<br>
><br>
> Do the above options load all the data at the given level or do I also<br>
> need to run it for each state and county with the options below?<br>
<br>
Yes, by default it'll load all data at that level.<br>
<br>
> -s <2-digit state code> Individual state code to load<br>
> -c <3-digit or 5-digit county code> Individual county code to load<br>
<br>
These are used if you, for example, have to reload a specific state<br>
and/or county. I had the luxury of doing this when someone decided to<br>
restart the database in the middle of my data load.. :)<br>
<br>
> Does:<br>
><br>
> cd /u/srcdata/tiger2007fe<br>
> ./process_tiger.sh -n -b -c<br>
><br>
> load all the data? Assuming that that directory is a mirror of the<br>
> Census download directory.<br>
<br>
Yes, it should, though it expects to be run one level up from<br>
'TIGER2007FE', unless you adjust the "BASE" environment variable at the<br>
top of the script.<br>
<br>
Thanks,<br>
<br>
Stephen<br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: not available<br>
Type: application/pgp-signature<br>
Size: 189 bytes<br>
Desc: Digital signature<br>
Url : <a href="http://lists.refractions.net/pipermail/postgis-users/attachments/20080403/6ba49bb6/attachment-0001.bin" target="_blank">http://lists.refractions.net/pipermail/postgis-users/attachments/20080403/6ba49bb6/attachment-0001.bin</a><br>
<br>
------------------------------<br>
<br>
Message: 29<br>
Date: Fri, 04 Apr 2008 00:05:53 -0500<br>
From: Stephen Woodbridge <<a href="mailto:woodbri@swoodbridge.com">woodbri@swoodbridge.com</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:47F5B731.2080309@swoodbridge.com">47F5B731.2080309@swoodbridge.com</a>><br>
Content-Type: text/plain; charset=UTF-8; format=flowed<br>
<br>
OK, Now I'm feeling kind of dumb.<br>
<br>
./process_tiger.sh -n<br>
Internal error!<br>
./process_tiger.sh -n \*<br>
Internal error!<br>
./process_tiger.sh -n ''<br>
Internal error!<br>
<br>
woodbri@carto:~/work/new-tiger$ sh -x process_tiger.sh -d -n<br>
+ BASE=TIGER2007FE<br>
+ SETBASE=fe_2007<br>
+ PREFIX=tiger<br>
+ SKIP00=y<br>
++ mktemp -d -p .<br>
+ TMPDIR=./tmp.ZKN9vx<br>
+ SRID=4269<br>
+ USER=postgres<br>
+ HOST=carto<br>
+ DB=tiger2007fe<br>
+ ENCODING=LATIN1<br>
+ NATIONAL=false<br>
+ STATELVL=false<br>
+ STATES=<br>
+ COUNTYLVL=false<br>
+ COUNTIES=<br>
+ DROP=false<br>
++ getopt -o n::ls::bc::ed -n process_tiger.sh -- -d -n<br>
+ TEMP= -d -n '' --<br>
+ '[' 0 '!=' 0 ']'<br>
+ eval set -- ' -d -n '\'''\'' --'<br>
++ set -- -d -n '' --<br>
+ true<br>
+ DROP=true<br>
+ shift<br>
+ true<br>
+ NATIONAL=true<br>
+ shift<br>
+ true<br>
+ echo 'Internal error!'<br>
Internal error!<br>
+ exit 1<br>
<br>
What did I miss here.<br>
<br>
-Steve<br>
<br>
Stephen Frost wrote:<br>
> Stephen,<br>
><br>
> * Stephen Woodbridge (<a href="mailto:woodbri@swoodbridge.com">woodbri@swoodbridge.com</a>) wrote:<br>
>> I added to process_tiger.sh:<br>
>><br>
>> # User to connect to database with or use USER=$LOGNAME<br>
>> USER="postgres"<br>
><br>
> Ah, sure, makes sense. We user Kerberos for our authentication, so I<br>
> tend to forget about user/pw authentication issues.<br>
><br>
>> and in vi did:<br>
>><br>
>> %s/\<psql\>/psql -U $USER/g<br>
>><br>
>> to update the script.<br>
><br>
> All of the options passed to psql really should be done as options to<br>
> the script. I'll look at adding that (and -U support) sometime<br>
> tomorrow.<br>
><br>
>> Question on the README.txt<br>
>><br>
>> -n Load national-level data<br>
>> -b Load state-level data<br>
>> -c Load county-level data<br>
>><br>
>> Do the above options load all the data at the given level or do I also<br>
>> need to run it for each state and county with the options below?<br>
><br>
> Yes, by default it'll load all data at that level.<br>
><br>
>> -s <2-digit state code> Individual state code to load<br>
>> -c <3-digit or 5-digit county code> Individual county code to load<br>
><br>
> These are used if you, for example, have to reload a specific state<br>
> and/or county. I had the luxury of doing this when someone decided to<br>
> restart the database in the middle of my data load.. :)<br>
><br>
>> Does:<br>
>><br>
>> cd /u/srcdata/tiger2007fe<br>
>> ./process_tiger.sh -n -b -c<br>
>><br>
>> load all the data? Assuming that that directory is a mirror of the<br>
>> Census download directory.<br>
><br>
> Yes, it should, though it expects to be run one level up from<br>
> 'TIGER2007FE', unless you adjust the "BASE" environment variable at the<br>
> top of the script.<br>
><br>
> Thanks,<br>
><br>
> Stephen<br>
><br>
><br>
> ------------------------------------------------------------------------<br>
><br>
> _______________________________________________<br>
> postgis-users mailing list<br>
> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 30<br>
Date: Thu, 3 Apr 2008 21:34:47 -0700 (PDT)<br>
From: "<a href="mailto:kreshna_iceheart@yahoo.com">kreshna_iceheart@yahoo.com</a>" <<a href="mailto:kreshna_iceheart@yahoo.com">kreshna_iceheart@yahoo.com</a>><br>
Subject: [postgis-users] Cannot display postGIS layers on ms4w 4.6.1<br>
To: <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
Message-ID: <<a href="mailto:168108.36541.qm@web56105.mail.re3.yahoo.com">168108.36541.qm@web56105.mail.re3.yahoo.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
I'm trying to display postGIS layers on MapServer. I<br>
am using the following components:<br>
1) MapServer 4.6.1 for Windows (ms4w)<br>
2) postgreSQL 8.2<br>
3) shp2pgsql release 1.1.6<br>
4) QuantumGIS 0.9.2 to automatically generate my<br>
mapfile<br>
<br>
My mapfile is attached with this email. It only has a<br>
single layer (for testing purpose), and it does not<br>
contain any complex query at all. The layer definition<br>
is as follows:<br>
=======================================================<br>
LAYER<br>
NAME 'batas_administrasi_line2_utm'<br>
TYPE LINE<br>
CONNECTIONTYPE postgis<br>
CONNECTION 'host=localhost user=rtrw<br>
dbname=rtrw-pacitan-trenggalek'<br>
DATA 'the_geom FROM<br>
batas-admin.batas_administrasi_line2_utm'<br>
=======================================================<br>
<br>
<br>
Yet, when I tried to display the mapfile, I got the<br>
following error message:<br>
=======================================================<br>
msDrawMap(): Image handling error. Failed to draw<br>
layer named 'batas_administrasi_line2_utm'.<br>
prepare_database(): Query error. Error executing<br>
POSTGIS DECLARE (the actual query) statement: 'DECLARE<br>
mycursor BINARY CURSOR FOR SELECT<br>
asbinary(force_collection(force_2d(the_geom)),'NDR'),OID::text<br>
from batas-admin.batas_administrasi_line2_utm WHERE ()<br>
and (the_geom && setSRID( 'BOX3D(528656.510416667<br>
9074209,609071.177083333<br>
9134520)'::BOX3D,find_srid('','batas-admin.batas_administrasi_line2_utm','the_geom')<br>
))'<br>
<br>
Postgresql reports the error as 'ERROR: syntax error<br>
at or near "-" at character 116 '<br>
<br>
More Help:<br>
<br>
Error with POSTGIS data variable. You specified<br>
'<check your .map file>'.<br>
Standard ways of specifiying are :<br>
(1) 'geometry_column from geometry_table'<br>
(2) 'geometry_column from (<sub query>) as foo using<br>
unique <column name> using SRID=<srid#>'<br>
<br>
Make sure you put in the 'using unique <column name>'<br>
and 'using SRID=#' clauses in.<br>
<br>
For more help, please see<br>
<a href="http://postgis.refractions.net/documentation.php" target="_blank">http://postgis.refractions.net/documentation.php</a><br>
<br>
Mappostgis.c - version of Jan 23/2004.<br>
=======================================================<br>
<br>
I have tried the following solution:<br>
<br>
(1) changing "FROM" to "from" (uppercase to lowercase)<br>
<br>
(2) enabling oid in my database:<br>
- uncommenting the line "default_with_oids" in my<br>
postgresql.on<br>
- changing the value to "on"<br>
- restarting posgresql<br>
- re-creating all my postGIS tables<br>
<br>
<br>
Yet none of the solution above work. MapServer just<br>
keeps giving me the error message.<br>
<br>
What happens? Anyone know a working solution?<br>
<br>
Thanks,<br>
-Kresh<br>
<br>
<br>
<br>
<br>
____________________________________________________________________________________<br>
You rock. That's why Blockbuster's offering you one<br>
month of Blockbuster Total Access, No Cost.<br>
<a href="http://tc.deals.yahoo.com/tc/blockbuster/text5.com" target="_blank">http://tc.deals.yahoo.com/tc/blockbuster/text5.com</a><br>
<br>
__________________________________________________________________<br>
So be it. If saying "NO" means being alone, then to hell with love, with women, with marriage, with God, religions, bars, nightclubs, computer games, and all the shit life keeps pumping at me. I'll walk alone, but with freedom and a healed pride.<br>
<br>
<br>
____________________________________________________________________________________<br>
You rock. That's why Blockbuster's offering you one month of Blockbuster Total Access, No Cost.<br>
<a href="http://tc.deals.yahoo.com/tc/blockbuster/text5.com" target="_blank">http://tc.deals.yahoo.com/tc/blockbuster/text5.com</a><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: rtrw-pacitan-trenggalek1.map<br>
Type: application/octet-stream<br>
Size: 2994 bytes<br>
Desc: 165038048-rtrw-pacitan-trenggalek1.map<br>
Url : <a href="http://lists.refractions.net/pipermail/postgis-users/attachments/20080403/38e12258/rtrw-pacitan-trenggalek1.dll" target="_blank">http://lists.refractions.net/pipermail/postgis-users/attachments/20080403/38e12258/rtrw-pacitan-trenggalek1.dll</a><br>
<br>
------------------------------<br>
<br>
Message: 31<br>
Date: Fri, 4 Apr 2008 01:37:24 -0400<br>
From: "Paragon Corporation" <<a href="mailto:lr@pcorp.us">lr@pcorp.us</a>><br>
Subject: RE: [postgis-users] Cannot display postGIS layers on ms4w<br>
4.6.1<br>
To: "'PostGIS Users Discussion'"<br>
<<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <001701c89615$f6453730$50812e40@r><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Just a guess I think bata-admin as a schema is a bad choice of schema names.<br>
I think you have to quote it if you are going to use that. I would change<br>
it to rename your schema if I were you. To something like bata_admin<br>
<br>
If you go with what you have, I think you may need to quote it and I'm not<br>
even sure if that is legal in mapserver since I've never needed to quote my<br>
schemas.<br>
<br>
LAYER<br>
NAME 'batas_administrasi_line2_utm'<br>
TYPE LINE<br>
CONNECTIONTYPE postgis<br>
CONNECTION 'host=localhost user=rtrw dbname=rtrw-pacitan-trenggalek'<br>
DATA 'the_geom FROM<br>
"batas-admin".batas_administrasi_line2_utm'<br>
<br>
Hope that helps,<br>
Regina<br>
<br>
<br>
<br>
-----Original Message-----<br>
From: <a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a><br>
[mailto:<a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a>] On Behalf Of<br>
<a href="mailto:kreshna_iceheart@yahoo.com">kreshna_iceheart@yahoo.com</a><br>
Sent: Friday, April 04, 2008 12:35 AM<br>
To: <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
Subject: [postgis-users] Cannot display postGIS layers on ms4w 4.6.1<br>
<br>
I'm trying to display postGIS layers on MapServer. I am using the following<br>
components:<br>
1) MapServer 4.6.1 for Windows (ms4w)<br>
2) postgreSQL 8.2<br>
3) shp2pgsql release 1.1.6<br>
4) QuantumGIS 0.9.2 to automatically generate my mapfile<br>
<br>
My mapfile is attached with this email. It only has a single layer (for<br>
testing purpose), and it does not contain any complex query at all. The<br>
layer definition is as follows:<br>
=======================================================<br>
LAYER<br>
NAME 'batas_administrasi_line2_utm'<br>
TYPE LINE<br>
CONNECTIONTYPE postgis<br>
CONNECTION 'host=localhost user=rtrw dbname=rtrw-pacitan-trenggalek'<br>
DATA 'the_geom FROM<br>
batas-admin.batas_administrasi_line2_utm'<br>
=======================================================<br>
<br>
<br>
Yet, when I tried to display the mapfile, I got the following error message:<br>
=======================================================<br>
msDrawMap(): Image handling error. Failed to draw layer named<br>
'batas_administrasi_line2_utm'.<br>
prepare_database(): Query error. Error executing POSTGIS DECLARE (the actual<br>
query) statement: 'DECLARE mycursor BINARY CURSOR FOR SELECT<br>
asbinary(force_collection(force_2d(the_geom)),'NDR'),OID::text<br>
from batas-admin.batas_administrasi_line2_utm WHERE () and (the_geom &&<br>
setSRID( 'BOX3D(528656.510416667<br>
9074209,609071.177083333<br>
9134520)'::BOX3D,find_srid('','batas-admin.batas_administrasi_line2_utm','th<br>
e_geom')<br>
))'<br>
<br>
Postgresql reports the error as 'ERROR: syntax error at or near "-" at<br>
character 116 '<br>
<br>
More Help:<br>
<br>
Error with POSTGIS data variable. You specified '<check your .map file>'.<br>
Standard ways of specifiying are :<br>
(1) 'geometry_column from geometry_table'<br>
(2) 'geometry_column from (<sub query>) as foo using unique <column name><br>
using SRID=<srid#>'<br>
<br>
Make sure you put in the 'using unique <column name>'<br>
and 'using SRID=#' clauses in.<br>
<br>
For more help, please see<br>
<a href="http://postgis.refractions.net/documentation.php" target="_blank">http://postgis.refractions.net/documentation.php</a><br>
<br>
Mappostgis.c - version of Jan 23/2004.<br>
=======================================================<br>
<br>
I have tried the following solution:<br>
<br>
(1) changing "FROM" to "from" (uppercase to lowercase)<br>
<br>
(2) enabling oid in my database:<br>
- uncommenting the line "default_with_oids" in my postgresql.on<br>
- changing the value to "on"<br>
- restarting posgresql<br>
- re-creating all my postGIS tables<br>
<br>
<br>
Yet none of the solution above work. MapServer just keeps giving me the<br>
error message.<br>
<br>
What happens? Anyone know a working solution?<br>
<br>
Thanks,<br>
-Kresh<br>
<br>
<br>
<br>
<br>
____________________________________________________________________________<br>
________<br>
You rock. That's why Blockbuster's offering you one month of Blockbuster<br>
Total Access, No Cost.<br>
<a href="http://tc.deals.yahoo.com/tc/blockbuster/text5.com" target="_blank">http://tc.deals.yahoo.com/tc/blockbuster/text5.com</a><br>
<br>
__________________________________________________________________<br>
So be it. If saying "NO" means being alone, then to hell with love, with<br>
women, with marriage, with God, religions, bars, nightclubs, computer games,<br>
and all the shit life keeps pumping at me. I'll walk alone, but with freedom<br>
and a healed pride.<br>
<br>
<br>
<br>
____________________________________________________________________________<br>
________<br>
You rock. That's why Blockbuster's offering you one month of Blockbuster<br>
Total Access, No Cost.<br>
<a href="http://tc.deals.yahoo.com/tc/blockbuster/text5.com" target="_blank">http://tc.deals.yahoo.com/tc/blockbuster/text5.com</a><br>
<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 32<br>
Date: Fri, 04 Apr 2008 10:24:21 +0200<br>
From: Barend Kobben <<a href="mailto:kobben@itc.nl">kobben@itc.nl</a>><br>
Subject: Re: [postgis-users] Problem with the Codification/charset<br>
from Postgis to mapserver<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:C41BB255.1BE1%25kobben@itc.nl">C41BB255.1BE1%kobben@itc.nl</a>><br>
Content-Type: text/plain; charset="ISO-8859-1"<br>
<br>
.... also make sure Mapserevre connects as a UTF-8 client to PG by satining<br>
including in the CONNECTION object:<br>
options='-c client_encoding=UTF8'<br>
<br>
<br>
--<br>
Barend Köbben<br>
International Institute for Geo-Information<br>
Sciences and Earth Observation (ITC)<br>
PO Box 6<br>
7500AA Enschede, The Netherlands<br>
+31 (0)53 4874253<br>
<br>
<br>
<br>
On 03-04-2008 18:10, "Stephen Woodbridge" <<a href="mailto:woodbri@swoodbridge.com">woodbri@swoodbridge.com</a>> wrote:<br>
<br>
> Emilio Ponce wrote:<br>
>> Hi everyone!<br>
>><br>
>> Recently I've done a simple mapfile (mapserver) that gets a postgis<br>
>> table and represents it. The problem is about the codification: The<br>
>> Postgis database is coded as UTF8 and when I represent a varchar column<br>
>> with mapserver the accents and especial characters are wrong. What can I<br>
>> do? Can I easily change the mapserver codification?<br>
><br>
> This is probably better posted to the mapserver list. But you might try<br>
> adding to your label block<br>
><br>
> ENCODING "UTF-8"<br>
><br>
> -Steve W<br>
> _______________________________________________<br>
> postgis-users mailing list<br>
> <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br>
International Institute for Geo-Information Science and Earth Observation (ITC)<br>
Chamber of Commerce: 410 27 560<br>
<br>
E-mail disclaimer<br>
The information in this e-mail, including any attachments, is intended for the addressee only. If you are not the intended recipient, you are hereby notified that any disclosure, copying, distribution or action in relation to the content of this information is strictly prohibited. If you have received this e-mail by mistake, please delete the message and any attachment and inform the sender by return e-mail. ITC accepts no liability for any error or omission in the message content or for damage of any kind that may arise as a result of e-mail transmission.<br>
<br>
<br>
------------------------------<br>
<br>
Message: 33<br>
Date: Fri, 04 Apr 2008 08:01:31 +0000<br>
From: "Jonathan W. Lowe" <<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <1207296091.3135.22.camel@localhost.localdomain><br>
Content-Type: text/plain<br>
<br>
On Thu, 2008-04-03 at 20:31 -0400, Stephen Frost wrote:<br>
> Jonathan,<br>
><br>
> * Jonathan W. Lowe (<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>) wrote:<br>
> > ...And in case the images don't persist through the mail server, they're<br>
> > viewable at: <a href="http://www.giswebsite.com/demos/tiger_overlays.html" target="_blank">http://www.giswebsite.com/demos/tiger_overlays.html</a><br>
><br>
> You know, I just realized that you were talking about the Census 2000<br>
> blocks (tabblock00.shp). Is there some reason you're using that<br>
> instead of the current data (tabblock.shp)? They might not want to<br>
> update the data from 2000 for historical reasons...<br>
> (Note: I havn't actually gone and looked, it just occured to me..)<br>
<br>
I haven't yet confirmed whether the 2007 version of Census Blocks<br>
maintains a 1:1 match with the associated statistics in SF1, so was<br>
starting with the 2000 version. However, positionally, the two Census<br>
Block versions (tabblock00 and tabblock) are identical. And they both<br>
positionally match the single 2007 edges data. (I've added a third<br>
screen shot to <a href="http://www.giswebsite.com/demos/tiger_overlays.html" target="_blank">www.giswebsite.com/demos/tiger_overlays.html</a> to<br>
illustrate.)<br>
<br>
Has anyone else seen the same misalignment between TIGER (2007 and/or<br>
2000) when overlaying on other datasets?<br>
<br>
It's still a mystery how OpenStreetMap's tiles for Berkeley, which are<br>
based on TIGER, don't seem to duplicate the zig-zag qualities of the<br>
edges data I've downloaded from TIGER 2007. Something is definitely<br>
missing in the puzzle, but it doesn't sound like the issue is with the<br>
PostGIS conversion part of the process, so, thanks for the help to this<br>
point -- I'll check with the OpenStreetMap community next.<br>
<br>
><br>
> Thanks,<br>
><br>
> Stephen<br>
><br>
> > On Fri, 2008-04-04 at 01:07 +0100, Jonathan W. Lowe wrote:<br>
> > > Stephen,<br>
> > ><br>
> > > My initial testing has been on Alameda County (California) TIGER data.<br>
> > > The two attached image files show an overlay of US Census 2000 Blocks<br>
> > > over an area south of the UC Berkeley campus. The offset is the same<br>
> > > for both Google and OpenStreetMap (OSM). This suggests that I've made a<br>
> > > mistake somewhere, because the OSM tiles in the United States are all<br>
> > > rendered from TIGER linework, so the TIGER census blocks should match<br>
> > > exactly.<br>
> > ><br>
> > > For the same source shapefile (tabblock00.shp), there's a nearly perfect<br>
> > > match between block boundaries and streets in the area just South of<br>
> > > Oakland's Lake Merritt. It smells like a datum conversion issue...<br>
> > ><br>
> > > The conversion path was from shapefile to PostGIS using shp2pgsql. I<br>
> > > used a custom projection of 32767 rather than 4269 because the existing<br>
> > > srtext for 4269 had a degree value as 0.01745329251994328, but the US<br>
> > > Census metadata listed a degree value of 0.017453292519943295. Perhaps<br>
> > > not significant? My spatial_ref_sys entries for 4269 and 32767 are<br>
> > > otherwise pretty similar:<br>
> > ><br>
> > > SRID: 4269<br>
> > > SRTEXT: GEOGCS["NAD83",DATUM["North_American_Datum_1983",<br>
> > > SPHEROID["GRS 1980",6378137,298.257222101,<br>
> > > AUTHORITY["EPSG","7019"]],<br>
> > > AUTHORITY["EPSG","6269"]],<br>
> > > PRIMEM["Greenwich",0,<br>
> > > AUTHORITY["EPSG","8901"]],<br>
> > > UNIT["degree",0.01745329251994328,<br>
> > > AUTHORITY["EPSG","9122"]],<br>
> > > AUTHORITY["EPSG","4269"]]<br>
> > > PROJ4TEXT: +proj=longlat +ellps=GRS80 +datum=NAD83 +no_defs<br>
> > ><br>
> > > SRID: 32767<br>
> > > SRTEXT: GEOGCS["GCS_North_American_1983",<br>
> > > DATUM["D_North_American_1983",<br>
> > > SPHEROID["GRS_1980",6378137,298.257222101]],<br>
> > > PRIMEM["Greenwich",0],<br>
> > > UNIT["Degree",0.017453292519943295]]<br>
> > > PROJ4TEXT: +proj=longlat +ellps=clrk66 +datum=NAD27 +no_defs<br>
> > ><br>
> > > To display census block data in OpenStreetMap, I extract it from PostGIS<br>
> > > with a transform to EPSG 4326, although the coordinates don't seem to<br>
> > > change as a result. (This seems correct, as datum=NAD83 and datum=WGS84<br>
> > > are, for my purposes at least, are essentially identical.)<br>
> > ><br>
> > > Thanks,<br>
> > > Jonathan<br>
> > ><br>
> > > 2 attachments: TIGER2007andOSM.png, TIGER2007andGoogle.png<br>
> > ><br>
> > ><br>
> > > On Thu, 2008-04-03 at 19:18 -0400, Stephen Frost wrote:<br>
> > > > Jonathan,<br>
> > > ><br>
> > > > * Jonathan W. Lowe (<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>) wrote:<br>
> > > > > Have you yet tried overlaying TIGER 2007 linework or census block/tract<br>
> > > > > polygons over Google or OpenStreetMap tiles? I'm seeing a good match in<br>
> > > > > some areas but a significant shift (~50 meters) in others. Thought it<br>
> > > > > might be a datum conversion issue, but can't seem to find a match.<br>
> > > ><br>
> > > > I hadn't looked at the linework too much yet or tried to overlay it.<br>
> > > > I'm curious where you're seeing the differences though because I know<br>
> > > > that Census is only about half way through their MAF improvment project<br>
> > > > and I actually have some info about what has been done so far and what<br>
> > > > hasn't. It'd be interesting to see if it matches up.<br>
> > > ><br>
> > > > There are a few places (Guam, Hawaii islands) where they actually do use<br>
> > > > an SRID other than 4269, but my scripts don't yet handle that and I'm<br>
> > > > guessing that's not what you're referring to anyway. :)<br>
> > > ><br>
> > > > Thanks!<br>
> > > ><br>
> > > > Stephen<br>
> > > ><br>
> > > > > On Thu, 2008-04-03 at 17:07 -0400, Stephen Frost wrote:<br>
> > > > > > * Stephen Frost (<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>) wrote:<br>
> > > > > > > I think they may have also upgraded their pipe.. I got about 1.41MB/s<br>
> > > > > > > (11 Mb/s) for the whole transfer. It's about 22G all told. I'll<br>
> > > > > > > probably be trying to load it up into PG on one of our servers tomorrow.<br>
> > > > > > > It was a bit over 4 hours for me to pull down off of their<br>
> > > > > > > <a href="http://ftp2.census.gov" target="_blank">ftp2.census.gov</a> ftp site.<br>
> > > > > ><br>
> > > > > > Just to update those who might be interested- I've finished the data<br>
> > > > > > load into one of our servers at work. It comes to ~60GB on disk in<br>
> > > > > > PostgreSQL/PostGIS with appropriate indexes in most places and whatnot.<br>
> > > > > > Based on what I've seen so far, it looks *very* nice, especially the<br>
> > > > > > hydrogrophy ("areawater"). It also appears to be pretty consistant<br>
> > > > > > across the layers, which is also good.<br>
> > > > > ><br>
> > > > > > If anyone's interested in the scripts used to load the data (they're<br>
> > > > > > pretty simple, really), I'd be happy to provide them.<br>
> > > > > ><br>
> > > > > > Enjoy,<br>
> > > > > ><br>
> > > > > > Stephen<br>
> > > > > > _______________________________________________<br>
> > > > > > postgis-users mailing list<br>
> > > > > > <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> > > > > > <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
> > > > ><br>
> > > > > _______________________________________________<br>
> > > > > postgis-users mailing list<br>
> > > > > <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> > > > > <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
> ><br>
> > _______________________________________________<br>
> > postgis-users mailing list<br>
> > <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
> > <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 34<br>
Date: Fri, 4 Apr 2008 07:04:50 -0400<br>
From: Stephen Frost <<a href="mailto:sfrost@snowman.net">sfrost@snowman.net</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: PostGIS Users Discussion <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID: <<a href="mailto:20080404110450.GF4999@tamriel.snowman.net">20080404110450.GF4999@tamriel.snowman.net</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Steve,<br>
<br>
* Stephen Woodbridge (<a href="mailto:woodbri@swoodbridge.com">woodbri@swoodbridge.com</a>) wrote:<br>
> OK, Now I'm feeling kind of dumb.<br>
<br>
You shouldn't, that was my goof. :)<br>
<br>
> ++ getopt -o n::ls::bc::ed -n process_tiger.sh -- -d -n<br>
<br>
That should have been 'nl::s::...'. I've fixed that in my script,<br>
please give it another go. I made some changes (like adding parameters,<br>
heh) after I got the data loaded and so havn't completely retested<br>
everything. Actually, I've made some further updates/improvments which<br>
will help (like, I dunno, error checking...).<br>
<br>
Thanks!<br>
<br>
Stephen<br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: not available<br>
Type: application/pgp-signature<br>
Size: 189 bytes<br>
Desc: Digital signature<br>
Url : <a href="http://lists.refractions.net/pipermail/postgis-users/attachments/20080404/e6024002/attachment-0001.bin" target="_blank">http://lists.refractions.net/pipermail/postgis-users/attachments/20080404/e6024002/attachment-0001.bin</a><br>
<br>
------------------------------<br>
<br>
Message: 35<br>
Date: Fri, 4 Apr 2008 14:52:38 +0300<br>
From: "Nick Black" <<a href="mailto:nickblack1@gmail.com">nickblack1@gmail.com</a>><br>
Subject: Re: [postgis-users] TIGER/Line Shapefiles released<br>
To: "Jonathan W. Lowe" <<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>>, "PostGIS Users<br>
Discussion" <<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>><br>
Message-ID:<br>
<<a href="mailto:223020e60804040452i15f6e3b3qd84ec37225edb170@mail.gmail.com">223020e60804040452i15f6e3b3qd84ec37225edb170@mail.gmail.com</a>><br>
Content-Type: text/plain; charset=ISO-8859-1<br>
<br>
On Fri, Apr 4, 2008 at 11:01 AM, Jonathan W. Lowe <<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>> wrote:<br>
> On Thu, 2008-04-03 at 20:31 -0400, Stephen Frost wrote:<br>
> > Jonathan,<br>
> ><br>
> > * Jonathan W. Lowe (<a href="mailto:jlowe@giswebsite.com">jlowe@giswebsite.com</a>) wrote:<br>
><br>
>...<br><br>[Message clipped] </blockquote></div><br><br clear="all"><br>-- <br>Dylan Lorimer | Strategic Partner Management<br>415.573.2909 (Grand Central) | 650.644.0182 (Fax)<br><br>If you received this communication by mistake, please don't forward it to anyone else (it may contain confidential or privileged information), please erase all copies of it, including all attachments, and please let the sender know it went to the wrong person. Thanks.