Re: [postgis-users] Recommendations with GRASS v.generalize() on Tiger Data

2011-06-29 Thread Dan Putler

Hi Chris,

Your question on QGIS and GRASS was posted to the PostGIS user list. 
While you might get a response, this seems like a question for the QGIS 
user list since it actually doesn't involve PostGIS.


Dan

On 06/29/2011 03:55 PM, Christian Guirreri wrote:
Narrowed down some test to using the "Hermite" algorithm, but I'm 
having some bizarre issues with it.


In the attached gif of California counties from Tiger data, from left 
to right, I have the following tolerance values:

 - original
 - 1.0
 - 0.08
 - 0.01
 - 0.1

Why do counties disappear entirely as I decrease the tolerance?

Setup is QuantumGIS with GRASS. In the Grass Tools I choose the 
v.generalize function. I choose Boundary as the feature type (though 
I've tried checking others and it doesn't seem to change anything). 
Everything else is default, except for tolerance as notated above.


When I tested this original on only Arkansas and Mississippi, I got 
really nice results. I then tried it on the entire US and had the 
missing counties problem. So I tried only California, and still have 
the same issue.


I've tried other algorithms, but this has so far given me what I want. 
Any thoughts?

 - Chris


On Wed, Jun 29, 2011 at 5:10 PM, Christian Guirreri 
mailto:christ...@guirreri.com>> wrote:


I'm currently going through the GRASS v.generalize() function's
various parameters in QuantumGIS. There's so many options, I'm not
entirely sure what's best. Has anyone tried this on Tiger 2010
counties or district data? Any particular recommendations?

What's most important to me, as mentioned in a previous thread, is
that there are no gaps between counties/districts in a similar
fashion to MapShaper.

Thanks,
 - Chris




___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


[postgis-users] Recommendations with GRASS v.generalize() on Tiger Data

2011-06-29 Thread Christian Guirreri
I'm currently going through the GRASS v.generalize() function's various
parameters in QuantumGIS. There's so many options, I'm not entirely sure
what's best. Has anyone tried this on Tiger 2010 counties or district data?
Any particular recommendations?

What's most important to me, as mentioned in a previous thread, is that
there are no gaps between counties/districts in a similar fashion to
MapShaper.

Thanks,
 - Chris
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] st_collect, st_union

2011-06-29 Thread Bob Pawley
The table I sent you was imported using ogr2ogr. A second table, imported 
the same way, although with fewer rows has had no similar problems.


Bob

-Original Message- 
From: Sandro Santilli

Sent: Wednesday, June 29, 2011 12:33 PM
To: Bob Pawley
Cc: PostGIS Users Discussion
Subject: Re: [postgis-users] st_collect, st_union

On Wed, Jun 29, 2011 at 12:17:07PM -0700, Bob Pawley wrote:

I just retried it using the script of imp_process.

Same empty value returned.


Against an empty database ?
Really weird.
What version of postgis are you running ?

--strk;



Bob

-Original Message- From: Sandro Santilli
Sent: Wednesday, June 29, 2011 12:09 PM
To: PostGIS Users Discussion
Subject: Re: [postgis-users] st_collect, st_union

On Wed, Jun 29, 2011 at 11:15:23AM -0700, Bob Pawley wrote:
>Hi
>
>The table is attached.

psql:E5R.backup:108: ERROR:  function imp_process() does not exist

>Following is the query I used.
>
>select st_collect(import_process.wkb_geometry)
>from import_process;

strk=# select st_summary(st_collect(import_process.wkb_geometry))
from import_process;
st_summary

MultiLineString[B] with 5 elements
  LineString[B] with 46 points
  LineString[B] with 2 points
  LineString[B] with 46 points
  LineString[B] with 46 points
  LineString[B] with 46 points

(1 row)

Did you retry after import ?
Maybe the missing function makes a difference...

--strk;

 ()   Free GIS & Flash consultant/developer
 /\   http://strk.keybit.net/services.html
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users

___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


--

 ()   Free GIS & Flash consultant/developer
 /\   http://strk.keybit.net/services.html 


___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] st_collect, st_union

2011-06-29 Thread Bob Pawley

What do you mean with your reference about an empty database??

I'm running   PostgreSQL v 8.4.1 and  PGAdmin v 1.10.0.

Bob

-Original Message- 
From: Sandro Santilli 
Sent: Wednesday, June 29, 2011 12:33 PM 
To: Bob Pawley 
Cc: PostGIS Users Discussion 
Subject: Re: [postgis-users] st_collect, st_union 


On Wed, Jun 29, 2011 at 12:17:07PM -0700, Bob Pawley wrote:

I just retried it using the script of imp_process.

Same empty value returned.


Against an empty database ?
Really weird.
What version of postgis are you running ?

--strk;



Bob

-Original Message- From: Sandro Santilli
Sent: Wednesday, June 29, 2011 12:09 PM
To: PostGIS Users Discussion
Subject: Re: [postgis-users] st_collect, st_union

On Wed, Jun 29, 2011 at 11:15:23AM -0700, Bob Pawley wrote:
>Hi
>
>The table is attached.

psql:E5R.backup:108: ERROR:  function imp_process() does not exist

>Following is the query I used.
>
>select st_collect(import_process.wkb_geometry)
>from import_process;

strk=# select st_summary(st_collect(import_process.wkb_geometry))
from import_process;
st_summary

MultiLineString[B] with 5 elements
  LineString[B] with 46 points
  LineString[B] with 2 points
  LineString[B] with 46 points
  LineString[B] with 46 points
  LineString[B] with 46 points

(1 row)

Did you retry after import ?
Maybe the missing function makes a difference...

--strk;

 ()   Free GIS & Flash consultant/developer
 /\   http://strk.keybit.net/services.html
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users

___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


--

 ()   Free GIS & Flash consultant/developer
 /\   http://strk.keybit.net/services.html
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] simple question about views and tables

2011-06-29 Thread David Fawcett
Is is an issue with just using SELECT * and returning multiple columns
named region?  Try naming the columns that you want from your join
explicitly.

On Wed, Jun 29, 2011 at 3:05 PM, Robert Buckley
 wrote:
> Hi,
> I have a table called "windturbines".
> I have created a view called "v_sum_yield".
> I would now like to either create another view which consists of both tables
> joined with the column "bez", or join the view to the original table.
> The table comprises of single wind turbines and their yield per year. I need
> to sum the yields according to the region. This table then needs to be
> served as wfs/wms through geoserver.
> my first step was...
> CREATE VIEW v_sum_yield AS
> select region, sum(yield) AS result
> from windturbines
> group by region
> and then to create another view comprising of an left outer join to the
> original table...
>
> CREATE VIEW v_join AS
> SELECT *
>     FROM windturbines LEFT OUTER JOIN v_sum_yield ON (windturbines.region =
> v_sum_yield.region);
> the result is..
> ERROR: column "region" specified more than once
> SQL Status:42701
> however if I just run the selection with the "CREATE VIEW v_join AS"...it
> works fine.
> can anyone tell me why?
>
> thanks,
> Rob
> ___
> postgis-users mailing list
> postgis-users@postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users
>
>
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


[postgis-users] simple question about views and tables

2011-06-29 Thread Robert Buckley
Hi,

I have a table called "windturbines". 
I have created a view called "v_sum_yield".

I would now like to either create another view which consists of both tables 
joined with the column "bez", or join the view to the original table.

The table comprises of single wind turbines and their yield per year. I need to 
sum the yields according to the region. This table then needs to be served as 
wfs/wms through geoserver.

my first step was...

CREATE VIEW v_sum_yield AS
select region, sum(yield) AS result 
from windturbines
group by region

and then to create another view comprising of an left outer join to the 
original 
table...



CREATE VIEW v_join AS

SELECT *
FROM windturbines LEFT OUTER JOIN v_sum_yield ON (windturbines.region = 
v_sum_yield.region);

the result is..
ERROR: column "region" specified more than onceSQL Status:42701

however if I just run the selection with the "CREATE VIEW v_join AS"...it works 
fine.

can anyone tell me why?


thanks,

Rob___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] UTF8 Error when importing Tiger County shapefile

2011-06-29 Thread Christian Guirreri
That worked - thanks! In the wizard, there's an option to change the
encoding that I didn't realize was directly related ("DBF file character
encoding").

On Tue, Jun 28, 2011 at 1:50 PM, Stephen Woodbridge  wrote:

> I have not used the import wizard, but all Tiger data is in LATIN1 so when
> I import that data into postgis using the shp2pgsql utility into a UTF8
> database, I specify -W LATIN1 so the import know to convert the LATIN1 data
> into UTF8.
>
> -Steve
>
>
> On 6/28/2011 1:32 PM, Christian Guirreri wrote:
>
>> I'm trying to import the Tiger 2010 County and equivalent shapefile for
>> the entire US. Its tl_2010_us_county10 found here:
>> http://www.census.gov/cgi-bin/**geo/shapefiles2010/main
>>
>> I create a new database with Encoding UTF8, template_postgis. If I
>> choose anything other than UTF8, it will not allow me to use the postgis
>> template.
>>
>> I then use the Shape to PostGIS Importer wizard, change the SRID to
>> 4269, but otherwise leave everything on default. It starts to import,
>> but stops early and here's there error:
>>
>> .
>> Shapefile type: Polygon
>> Postgis type: MULTIPOLYGON[2]
>> Importing shapefile (3221 records)...
>> Unable to convert data value to UTF-8 (iconv reports "Invalid
>> argument"). Current encoding is "UTF-8". Try "LATIN1" (Western
>> European), or one of the values described at
>> http://www.postgresql.org/**docs/current/static/multibyte.**html
>> .
>>
>> Shapefile import failed.
>>
>>
>>
>>
>>
>> __**_
>> postgis-users mailing list
>> postgis-users@postgis.**refractions.net
>> http://postgis.refractions.**net/mailman/listinfo/postgis-**users
>>
>
> __**_
> postgis-users mailing list
> postgis-users@postgis.**refractions.net
> http://postgis.refractions.**net/mailman/listinfo/postgis-**users
>
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] st_collect, st_union

2011-06-29 Thread Sandro Santilli
On Wed, Jun 29, 2011 at 12:17:07PM -0700, Bob Pawley wrote:
> I just retried it using the script of imp_process.
> 
> Same empty value returned.

Against an empty database ?
Really weird.
What version of postgis are you running ?

--strk;

> 
> Bob
> 
> -Original Message- From: Sandro Santilli
> Sent: Wednesday, June 29, 2011 12:09 PM
> To: PostGIS Users Discussion
> Subject: Re: [postgis-users] st_collect, st_union
> 
> On Wed, Jun 29, 2011 at 11:15:23AM -0700, Bob Pawley wrote:
> >Hi
> >
> >The table is attached.
> 
> psql:E5R.backup:108: ERROR:  function imp_process() does not exist
> 
> >Following is the query I used.
> >
> >select st_collect(import_process.wkb_geometry)
> >from import_process;
> 
> strk=# select st_summary(st_collect(import_process.wkb_geometry))
> from import_process;
> st_summary
> 
> MultiLineString[B] with 5 elements
>   LineString[B] with 46 points
>   LineString[B] with 2 points
>   LineString[B] with 46 points
>   LineString[B] with 46 points
>   LineString[B] with 46 points
> 
> (1 row)
> 
> Did you retry after import ?
> Maybe the missing function makes a difference...
> 
> --strk;
> 
>  ()   Free GIS & Flash consultant/developer
>  /\   http://strk.keybit.net/services.html
> ___
> postgis-users mailing list
> postgis-users@postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users
> 
> ___
> postgis-users mailing list
> postgis-users@postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users

-- 

  ()   Free GIS & Flash consultant/developer
  /\   http://strk.keybit.net/services.html
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] st_collect, st_union

2011-06-29 Thread Bob Pawley

I just retried it using the script of imp_process.

Same empty value returned.

Bob

-Original Message- 
From: Sandro Santilli

Sent: Wednesday, June 29, 2011 12:09 PM
To: PostGIS Users Discussion
Subject: Re: [postgis-users] st_collect, st_union

On Wed, Jun 29, 2011 at 11:15:23AM -0700, Bob Pawley wrote:

Hi

The table is attached.


psql:E5R.backup:108: ERROR:  function imp_process() does not exist


Following is the query I used.

select st_collect(import_process.wkb_geometry)
from import_process;


strk=# select st_summary(st_collect(import_process.wkb_geometry)) from 
import_process;

st_summary

MultiLineString[B] with 5 elements
  LineString[B] with 46 points
  LineString[B] with 2 points
  LineString[B] with 46 points
  LineString[B] with 46 points
  LineString[B] with 46 points

(1 row)

Did you retry after import ?
Maybe the missing function makes a difference...

--strk;

 ()   Free GIS & Flash consultant/developer
 /\   http://strk.keybit.net/services.html
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users 


___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] st_collect, st_union

2011-06-29 Thread Sandro Santilli
On Wed, Jun 29, 2011 at 11:15:23AM -0700, Bob Pawley wrote:
> Hi
> 
> The table is attached.

psql:E5R.backup:108: ERROR:  function imp_process() does not exist

> Following is the query I used.
> 
> select st_collect(import_process.wkb_geometry)
> from import_process;

strk=# select st_summary(st_collect(import_process.wkb_geometry)) from 
import_process;
 st_summary

 MultiLineString[B] with 5 elements
   LineString[B] with 46 points
   LineString[B] with 2 points
   LineString[B] with 46 points
   LineString[B] with 46 points
   LineString[B] with 46 points

(1 row)

Did you retry after import ?
Maybe the missing function makes a difference...

--strk;

  ()   Free GIS & Flash consultant/developer
  /\   http://strk.keybit.net/services.html
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] st_collect, st_union

2011-06-29 Thread Bob Pawley

Hi

The table is attached.

Following is the query I used.

select st_collect(import_process.wkb_geometry)
from import_process;

Bob

-Original Message- 
From: Sandro Santilli

Sent: Wednesday, June 29, 2011 10:11 AM
To: PostGIS Users Discussion
Subject: Re: [postgis-users] st_collect, st_union

On Wed, Jun 29, 2011 at 09:54:49AM -0700, Bob Pawley wrote:

Hi

I have a table with five rows each holding a linestring geometry.

Each of theses geometries indicate true on an st_isvalid query.

All rows are displayed as expected in QGis.

However, when I run an st_union or st_collect query on all five rows the 
return is empty.


If I select only four of the five rows (any four)  it returns a value.

Can anyone suggest what type of problem I should look for? Is there some 
restriction on these functions?


That's a really odd behavior. Can you provide a dump of your table ?

--strk;

 ()   Free GIS & Flash consultant/developer
 /\   http://strk.keybit.net/services.html
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users 


E5R.backup
Description: Binary data
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] Parcel household and dwelling unit data: add column if it doesn't exist

2011-06-29 Thread Paul Ramsey
RTFM on 'ALTER TABLE ADD COLUMN', 'CASE' and 'IS NULL'
P.

On Wed, Jun 29, 2011 at 10:33 AM, Conor Henley  wrote:
> Hi all,
>
>
>
> I’m working with parcel data from several different regional sources.  In
> some cases, the parcel data comes with a households field, in others with a
> dwelling units field (in rare cases, both).  In all cases, I need both.
> Luckily, either can be calculated from the other using census multipliers.
>
>
>
> 1.   What is the best way to add the missing column when it does not
> exist?
>
>
>
> 2.   Then, how should I go about populating the newly created field
> using a calculation based on the other field (keeping in mind that I don’t
> want to overwrite either field in the case that both existed in the input
> data)?
>
>
>
> Ideally, I would be able to add the missing field in cases where it does not
> exist and populate its rows using a calculation based on the existing
> field.  Because this step must nest in a larger set of processes, it is
> impractical to check manually each time a new input goes through.  I’ve seen
> a good deal of discussion of this topic on forums but haven’t seen any clear
> solutions.
>
>
>
> Thanks,
>
>
>
>
>
> Conor Henley
>
> GIS/Regional Planning Intern
>
>
>
> C A L T H O R P E  A S S O C I A T E S
>
> 2095 ROSE STREET, SUITE 201, BERKELEY, CALIFORNIA, 94709 USA
> 510-548-6800 | 510 548-6848 (fax)
>
> co...@calthorpe.com | www.calthorpe.com
>
>
>
> ___
> postgis-users mailing list
> postgis-users@postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users
>
>
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


[postgis-users] Parcel household and dwelling unit data: add column if it doesn't exist

2011-06-29 Thread Conor Henley
Hi all,

 

I'm working with parcel data from several different regional sources.
In some cases, the parcel data comes with a households field, in others
with a dwelling units field (in rare cases, both).  In all cases, I need
both.  Luckily, either can be calculated from the other using census
multipliers.

 

1.   What is the best way to add the missing column when it does not
exist?

 

2.   Then, how should I go about populating the newly created field
using a calculation based on the other field (keeping in mind that I
don't want to overwrite either field in the case that both existed in
the input data)?

 

Ideally, I would be able to add the missing field in cases where it does
not exist and populate its rows using a calculation based on the
existing field.  Because this step must nest in a larger set of
processes, it is impractical to check manually each time a new input
goes through.  I've seen a good deal of discussion of this topic on
forums but haven't seen any clear solutions.

 

Thanks, 

 

 

Conor Henley

GIS/Regional Planning Intern

 

C A L T H O R P E  A S S O C I A T E S

2095 ROSE STREET, SUITE 201, BERKELEY, CALIFORNIA, 94709 USA
510-548-6800 | 510 548-6848 (fax)

co...@calthorpe.com   | www.calthorpe.com

 

___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] st_collect, st_union

2011-06-29 Thread Sandro Santilli
On Wed, Jun 29, 2011 at 09:54:49AM -0700, Bob Pawley wrote:
> Hi
> 
> I have a table with five rows each holding a linestring geometry.
> 
> Each of theses geometries indicate true on an st_isvalid query.
> 
> All rows are displayed as expected in QGis.
> 
> However, when I run an st_union or st_collect query on all five rows the 
> return is empty.
> 
> If I select only four of the five rows (any four)  it returns a value.
> 
> Can anyone suggest what type of problem I should look for? Is there some 
> restriction on these functions?

That's a really odd behavior. Can you provide a dump of your table ?

--strk;

  ()   Free GIS & Flash consultant/developer
  /\   http://strk.keybit.net/services.html
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


[postgis-users] st_collect, st_union

2011-06-29 Thread Bob Pawley
Hi

I have a table with five rows each holding a linestring geometry.

Each of theses geometries indicate true on an st_isvalid query.

All rows are displayed as expected in QGis.

However, when I run an st_union or st_collect query on all five rows the return 
is empty.

If I select only four of the five rows (any four)  it returns a value.

Can anyone suggest what type of problem I should look for? Is there some 
restriction on these functions?

Bob

___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] How to Store Metadata?

2011-06-29 Thread Brett Antonides
Thank you all for the help.  I'm going to try storing it in
a separate metadata table with an XML column.  I don't foresee any issues
there, but I'll let you know if I run into any.

Thanks again.

Brett
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] ST_OffsetCurve

2011-06-29 Thread Felipe Provenzale
So, I think that we can have this function in Postgis in the next version.
But now, do you know if I can use a combination of functions with the
existent ST_Buffer function to generate parallel lines from a table withi
more of 600 lines?
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] Tiger Geocoder 2010 setup issues

2011-06-29 Thread Daniel Weitzenfeld
Thank you for all your help - I got it to work!

For posterity, the postgis-related fixes I made:
-unmarked the ALTER DATABASE line in create_geocode.sh
-changed pg_hba.conf to trust, trust, trust
-tweaked a file with a name like "postgis conf" to allow access on port
5432
-added the projections.sql file to my postgis-enabled db

I also had a bunch of silly problems related to my mediocre linux skills,
e.g.
-wget wasn't working because I didn't have permissions to create a directory
-I hadn't installed unzip

Basically I was in over my head but I surfaced.  The precision from the
geocodes I've run so far has been really amazing, good stuff Regina and Leo!



On Wed, Jun 29, 2011 at 12:13 PM, Paragon Corporation  wrote:

> **
> Daniel,
>
> Yap.  what we had done is move the declaration up above the CD call since
> we rna into same issue when testing in Linux.  As I recall I think when
> people complained about the issue, the CD call came before the STATEDIR.
>
> Regarding 1.5.3.  You should be fine to run with that.  For our production
> work involving tiger geocoder, we are using PostGIS 1.5 and I think most
> people are actually and just testing on PostGIS 2.0.
>
> However we did make changes that now require PostgreSQL 8.4 or above.  So
> you need to be running at least PostGIS 1.5 and PostgreSQL 8.4+.
>
>
> Did you verify you have a tiger schema in your database and that it does
> have tables -- addr,  county, edges, faces etc?
>
> From the addr fatal error you described we guessed 1 of 3 things was
> happening
>
> a) Your script was pointing at the wrong database and not the one you
> loaded tiger in
> b) The tiger schema wasn't in your search path --- If you look at the
> create_geocode.sh (you'll see a remarked out line to alter db to add tiger
> to your database search path), this is a common mistake people make not
> having that line
> which causes the error you described. We left it remarked out since some
> people have custom search paths and we didn't want to overwrite them.
>
> c) Your database just doesn't have the skeleton and lookup tables or is
> missing one of them
>
> Like I said your next issue is your pg_hba.conf or you have another
> instance of PostgreSQL running and you are pointing at the wrong instance.
>
> Regina and Leo,
> http://www.postgis.us
>
>
>  --
>  *From:* postgis-users-boun...@postgis.refractions.net [mailto:
> postgis-users-boun...@postgis.refractions.net] *On Behalf Of *Daniel
> Weitzenfeld
> *Sent:* Tuesday, June 28, 2011 6:52 PM
> *To:* PostGIS Users Discussion
> *Subject:* Re: [postgis-users] Tiger Geocoder 2010 setup issues
>
> Thank you Regina and Leo. It sounds like I might have a version issue.
>
> I'm using PostGIS 1.5.3, as per the installation instructions here:
> http://postgis.refractions.net/docs/ch02.html
>
> When you suggested redownloading, did you mean all of PostGIS, or just the
> geocoder?
>
> Is the code here the up to date version?  here =
> http://svn.osgeo.org/postgis/trunk/extras/tiger_geocoder/tiger_2010/
> Because in tiger_loader.sql at that location, there's the line:
>
> E'STATEDIR="${staging_fold}/${website_root}/${state_fold}"
>
>
> Is that the fixed version?
>
> In the meantime, I'll work on my ident issues.
>
>
>
>
>
>
>
> On Tue, Jun 28, 2011 at 11:46 PM, Paragon Corporation  wrote:
>
>> Daniel,
>>
>> Haven't had a chance to recheck our script to confirm its not something
>> amiss.
>>
>> Which version of PostGIS  are you working with?
>>
>>  You can ignore the does not exist, skipping errors -- those are trying to
>> drop tables which wouldn't exist for new installs.
>>
>> The addr does not exist often happens if you don't have tiger in your
>> database search_path.  Make sure its in your search_path.  Also you can
>> verify the table exists
>> By making sure you have a tiger schema and it has tables like addr, edges
>> etc.
>>
>> Regarding your ident problem.
>>
>> PostgreSQL can listen on localhost or local port and depending on which
>> that
>> has to be registered in pg_hba.conf
>>
>> You can try taking out the PGHOST line since you may be running only on a
>> non-IP port and see if that helps or change your pg_hba.conf to allow.
>>
>> Regarding this issue:
>> STATEDIR="/gisdata/www2.census.gov/geo/pvs/tiger2010st/36_New_York"
>>
>> especially because an early command is
>>
>>  cd $STATEDIR
>>
>> I'm pretty sure we fixed that so I suspect your version is probably out of
>> data and you may want to just redownload from PostGIS site (either the SVN
>> tar ball or directly from svn)
>>
>> Regina and Leo
>> http://www.postgis.us
>>
>>
>> -Original Message-
>> From: postgis-users-boun...@postgis.refractions.net
>> [mailto:postgis-users-boun...@postgis.refractions.net] On Behalf Of
>> Daniel
>> Weitzenfeld
>> Sent: Tuesday, June 28, 2011 6:51 AM
>> To: postgis-users@postgis.refractions.net
>> Subject: [postgis-users] Tiger Geocoder 2010 setup issues
>>
>> Hi,
>> I'm trying to setup the TIGER Geocoder, 

Re: [postgis-users] ST_OffsetCurve

2011-06-29 Thread Sandro Santilli
On Wed, Jun 29, 2011 at 11:11:33AM +0200, Andreas Neumann wrote:
> 
> >>The single sided buffers are already available in GEOS and Postgis,
> >>right?
> >
> >They are available in GEOS, not yet in PostGIS I'm afraid.
> 
> ok - can you do an estimate on that one as well?

I think this should come as an additional parameter to the existing
ST_Buffer() function. Will require GEOS-3.3.0 to work and some care
in order to retain PostGIS support for older versions of GEOS (with
which you wouldn't have such support).

I can implement, document and regress-test that within another 10 hours.

--strk;

  ()   Free GIS & Flash consultant/developer
  /\   http://strk.keybit.net/services.html
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] ST_OffsetCurve

2011-06-29 Thread Sandro Santilli
On Wed, Jun 29, 2011 at 11:17:39AM +0200, Stephan Holl wrote:
> Andreas Neumann , [20110629 - 11:11:33]

> >  I hope that someone else can pay for the single-sided buffer to be 
> >  available in Postgis? Any takers?
> 
> I thought that this was included in the estimation from strk as well?
> The initallay mentioned patch provided this already, but needs to get
> reworked as well as added tests and documentation.

The problem arised from a confusion about the term "single-sided buffer".
Rafalmag, author of the initial patch, called "single-sided buffer"
a function that would return a line, not an area. 

Such a lineal output was (by Martin Davis of JTS, Paul and me) seen as
best called "offset curve" instead, whereas "single-sided buffer" would
be a function returning an area instead. See the first comments exchange
on http://trac.osgeo.org/postgis/ticket/413

Later on Martin implemented proper "single-sided buffer" (areal) in
JTS, which I ported to GEOS-3.3.0 as far as to the C-API, where I also
exposed a proper "GEOSOffsetCurve" with a good signature to replace
the badly named "GEOSSingleSidedBuffer" on which rafalmag patch is based.

--strk; 

  ()   Free GIS & Flash consultant/developer
  /\   http://strk.keybit.net/services.html
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] ST_Transform troubles

2011-06-29 Thread Frans Knibbe


On 2011-06-27 12:44, Mike Toews wrote:

I think it is well established among most (all?) spatial DBs that
comply with OGC Standards that the coordinate order is Cartesian: "X
Y" or "Longitude Latitude". If there is any issues with this, then it
is with the OGC Standards you disagree with.
As far as I know, I am not disagreeing with any OGC standard. This is 
about PostGIS disagreeing with the EPSG standard, which it uses.

  Mind you, this is only an
internal ordering in WKB, that you may occasionally see in WKT. Both
of these are computer markup languages, not designed to look pretty in
a report. However, if you do want to make the coordinate look more
human readable, it is rather simple to write your own function to
report coordinates in any style you like.
My primary concern is not with human readability, but with automated 
processing by computers. A computer only knows what it is being told. 
How is it to know that the coordinates are switched if it is being told 
that the SRS used is defined by an EPSG specification that uses 
latitude-longitude?

An interesting case study of following OGC Standards for coordinate
ordering is with MS SQL Server 2008 Spatial's GEOGRAPHY type. When the
preview was released, coordinate order of WKT was "Lat Lon", but was
reversed to "Lon Lat" for the subsequent preview and final release.
See the following links for more:
http://social.msdn.microsoft.com/forums/en-US/sqlspatial/thread/41250c42-25e6-4de7-953e-a6c41ada383f/
http://blogs.msdn.com/b/isaac/archive/2007/12/27/latitude-longitude-ordering.aspx
http://www.spatiallyadjusted.com/2007/12/27/microsoft-turns-the-world-right-side-up-again/

Thanks, it is interesting read.

I think it is perfectly all right to always use the same axis order for 
all coordinates. But in that case you can not state that you are using 
EPSG definitions, you will have to use some other system of describing 
coordinate reference systems.

-Mike
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users



___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] ST_OffsetCurve

2011-06-29 Thread Stephan Holl
Hello all,

Andreas Neumann , [20110629 - 11:11:33]

> > I guess you could just collect the results of two ST_OffsetCurve
> > calls to obtain what you call "open buffer".
> 
>  yes, that would work.
> 
> >> The single sided buffers are already available in GEOS and Postgis,
> >> right?
> >
> > They are available in GEOS, not yet in PostGIS I'm afraid.
> 
>  ok - can you do an estimate on that one as well?
> 
> >> Any estimates how much work (hours/money) it would be to improve
> >> the patch and also do the tests and documentation?
> >
> > I belive it'll take around 10 hours, which would cost me 600 EUR.
> > It'd take 5 funders (120 EUR each) or 2 funders (300 EUR each).
> 
>  ok - that is reasonable. We can pay for that - no need to split it
> up. Let's get in contact off-list.
> 
>  I hope that someone else can pay for the single-sided buffer to be 
>  available in Postgis? Any takers?

I thought that this was included in the estimation from strk as well?
The initallay mentioned patch provided this already, but needs to get
reworked as well as added tests and documentation.

But I might me wrong.

Best
Stephan

-- 
Stephan Holl  | Tel.: +49 (0)541-33 508 3663
Intevation GmbH, Neuer Graben 17, 49074 OS  |  AG Osnabrück - HR B 18998
Geschäftsführer:  Frank Koormann, Bernhard Reiter, Dr. Jan-Oliver Wagner


signature.asc
Description: PGP signature
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


[postgis-users] ST_OffsetCurve

2011-06-29 Thread Andrea Peri
>
> I would be interested in having this feature in GEOS and Postgis. If it
> is available in GEOS, we can also use it in QGIS. I would also be
> interested in helping to finance it.

As I know (using the 1.7.0 or the trunk version with windows),
actually qgis don't use the geos 3.3.0,
instead it use the geos 3.2.2

I opened a ticket asking this update.
http://trac.osgeo.org/osgeo4w/ticket/232

but actually it is in waiting state.

Andrea.



-- 
-
Andrea Peri
. . . . . . . . .
qwerty àèìòù
-
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] ST_OffsetCurve

2011-06-29 Thread Andreas Neumann

I guess you could just collect the results of two ST_OffsetCurve
calls to obtain what you call "open buffer".


yes, that would work.


The single sided buffers are already available in GEOS and Postgis,
right?


They are available in GEOS, not yet in PostGIS I'm afraid.


ok - can you do an estimate on that one as well?


Any estimates how much work (hours/money) it would be to improve the
patch and also do the tests and documentation?


I belive it'll take around 10 hours, which would cost me 600 EUR.
It'd take 5 funders (120 EUR each) or 2 funders (300 EUR each).


ok - that is reasonable. We can pay for that - no need to split it up. 
Let's get in contact off-list.


I hope that someone else can pay for the single-sided buffer to be 
available in Postgis? Any takers?


Andreas

--
Andreas Neumann
Böschacherstrasse 10A
8624 Grüt (Gossau ZH)
Switzerland
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] ST_OffsetCurve

2011-06-29 Thread Sandro Santilli
On Wed, Jun 29, 2011 at 10:22:31AM +0200, Andreas Neumann wrote:

> I would also be interested in having "open buffers" - which means
> two parallel lines, but no line ends. I would assume that "open
> buffers" and ST_OffsetCurve would be quite similar to implement?

I guess you could just collect the results of two ST_OffsetCurve
calls to obtain what you call "open buffer".

> The single sided buffers are already available in GEOS and Postgis,
> right?

They are available in GEOS, not yet in PostGIS I'm afraid.

> Any estimates how much work (hours/money) it would be to improve the
> patch and also do the tests and documentation?

I belive it'll take around 10 hours, which would cost me 600 EUR.
It'd take 5 funders (120 EUR each) or 2 funders (300 EUR each).

--strk; 

  ()   Free GIS & Flash consultant/developer
  /\   http://strk.keybit.net/services.html
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] ST_OffsetCurve

2011-06-29 Thread Andreas Neumann

Hi Strk,

I would be interested in having this feature in GEOS and Postgis. If it 
is available in GEOS, we can also use it in QGIS. I would also be 
interested in helping to finance it.


I would also be interested in having "open buffers" - which means two 
parallel lines, but no line ends. I would assume that "open buffers" and 
ST_OffsetCurve would be quite similar to implement?


The single sided buffers are already available in GEOS and Postgis, 
right?


Any estimates how much work (hours/money) it would be to improve the 
patch and also do the tests and documentation?


Thanks,
Andreas

On Wed, 29 Jun 2011 09:21:08 +0200, Sandro Santilli wrote:

On Tue, Jun 28, 2011 at 07:05:04PM -0300, Felipe Provenzale wrote:

I need to add a .patch file in Postgis, ST_OffsetCurve developed by 
Rafal
Magda. I'm using Windows XP, Postgre 9.0 and Postgis 1.5, and I 
couldn't
find a way to add this .patch file. I need to create parallel lines 
in a

MULTILINESTRING table, with more of 600 lines, and I think that this
function will be helpfull.

Somebody could tell me how can I include this function on Postgis?


Hi Felipe,
the ST_OffsetCurve patch needs some more work before getting
included in the core. One thing is a change in the interface
so that sign of distance expresses which side of the curve
you want, another is rewriting to use the new GEOS API, as
the old one was deprecated. Finally, regress testing should
be added.

For more informations, refer to:
http://trac.osgeo.org/postgis/ticket/413

I'll be glad to work on such cleanup and integration but It'll
take a few hours and although a few people expressed interest
in this feature noone came up with fundings.

Maybe this could be an opportunity to try some crowdfunding, for 
example

trough http://pledgebank.com or similar.


--strk;

  ()   Free GIS & Flash consultant/developer
  /\   http://strk.keybit.net/services.html
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


--
--
Andreas Neumann
Böschacherstrasse 10A
8624 Grüt (Gossau ZH)
Switzerland
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] ST_OffsetCurve

2011-06-29 Thread Sandro Santilli
On Tue, Jun 28, 2011 at 07:05:04PM -0300, Felipe Provenzale wrote:

> I need to add a .patch file in Postgis, ST_OffsetCurve developed by Rafal
> Magda. I'm using Windows XP, Postgre 9.0 and Postgis 1.5, and I couldn't
> find a way to add this .patch file. I need to create parallel lines in a
> MULTILINESTRING table, with more of 600 lines, and I think that this
> function will be helpfull.
> 
> Somebody could tell me how can I include this function on Postgis?

Hi Felipe,
the ST_OffsetCurve patch needs some more work before getting
included in the core. One thing is a change in the interface
so that sign of distance expresses which side of the curve
you want, another is rewriting to use the new GEOS API, as
the old one was deprecated. Finally, regress testing should
be added.

For more informations, refer to:
http://trac.osgeo.org/postgis/ticket/413

I'll be glad to work on such cleanup and integration but It'll
take a few hours and although a few people expressed interest
in this feature noone came up with fundings.

Maybe this could be an opportunity to try some crowdfunding, for example
trough http://pledgebank.com or similar.


--strk;

  ()   Free GIS & Flash consultant/developer
  /\   http://strk.keybit.net/services.html
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users