gmane.comp.gis.postgis

2.1.0beta3dev upgrade error

Stephen Woodbridge — 2013-05-23T05:01:08

Hi all,

I just tried to upgrade postgis from 2.0.3 -> 2.1.0beta3dev and got the 
following error:

pgr_test=# alter extension postgis update to '2.1.0beta3dev';
ERROR:  could not load library 
"/usr/lib/postgresql/9.2/lib/postgis-2.1.so": 
/usr/lib/postgresql/9.2/lib/postgis-2.1.so: undefined symbol: 
json_tokener_errors

Did I forget a step or do something wrong?

-Steve

I have json-c-0.11 built and installed. In fact is is the ./configure 
report:

   PostGIS is now configured for i686-pc-linux-gnu

  -------------- Compiler Info -------------
   C compiler:           gcc -g -O2
   C++ compiler:         g++ -g -O2
   SQL preprocessor:     /usr/bin/cpp -traditional-cpp -P

  -------------- Dependencies --------------
   GEOS config:          /usr/local/bin/geos-config
   GEOS version:         3.3.8
   GDAL config:          /usr/local/bin/gdal-config
   GDAL version:         1.9.2
   PostgreSQL config:    /usr/lib/postgresql/9.2/bin/pg_config
   PostgreSQL version:   PostgreSQL 9.2.4
   PROJ4 version:

Collada Import

Paragon Corporation — 2013-05-23T03:54:44

This is probably more of a question for Olivier and Hugo.

I think the answer is no but wanted to ask anyway.

The new SFCGAL 3D functionality you have is there already present or plan to
be a collada importer.

Thanks,
Regina
http://www.postgis.us
http://postgis.net

create 10m buffer using input data in EPSG4326

karsten — 2013-05-22T20:18:36

Greetings, I am working on creating 10m buffers around roads:
 
SELECT 10 AS score, ST_Difference(
 ST_Union(ST_Buffer(roads.the_geom, 10)),
 ST_Union(ST_Buffer(roads.the_geom, 0))
 )
AS the_geom
FROM geodata.roads2012_ca roads
GROUP BY score
 
So far I have the roads in dd (EPSG 4326) and was wondering how I could pass
the distances in meters into the buffer function (without necessarily
creating a new column with projected geometries). I guess could use some
geography data type function .. but how ?
 
Thanks
Karsten

Karsten Vennemann
Principal

Terra GIS LTD
2119 Boyer Ave E 
Seattle, WA  98112
USA 
 <http://www.terragis.net/> www.terragis.net

Phone ++1 206 905 1711
Fax      ++1 925 905 1711

 
_______________________________________________
postgis-users mailing list
postgis-users< at >lists.osgeo.org
http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users

Area of non-intersecting parts of layers

Cedric Duprez — 2013-05-21T16:55:26

Hi all,

I have 2 layers composed of multipolygons in Postgis 2.0. Each layer contains polygons with a name that is sometime common to several polygons.
When I try to find area of intersections between the 2 layers, no problem with the following query:
SELECT l1.name1, l2.name2, SUM(ST_Area(ST_Intersection(l1.geom, l2.geom))) / 10000 AS surface_ha
FROM layer1 l1
INNER JOIN layer2 l2 ON ST_Intersects(l1.geom, l2.geom)
GROUP BY 1, 2
ORDER BY 1, 2;

OK.
But I also try to get the area of my layer1 (grouped by name1) that does not contain elements of my layer2.

I tried with LEFT JOIN, with ST_CONTAINS, without success.

Does anyone have an idea on how to write this query?

Thanks in advance,

Cedric Duprez
_______________________________________________
postgis-users mailing list
postgis-users< at >lists.osgeo.org
http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users

permission denied for relation geometry_columns onpostgis 2.x

Pietro Rossin — 2013-05-21T12:13:31

Hello all
how can I solve this problem?
With postgres user no problem (of course), but I made a gisadmin user which
I thought I gave enough grants..
geometry_columns with postgis 2.x is not a table but a view that, if I'm not
wrong, is updated by some trigger/function

How can I make it usable/selectable by users?

Thanks
Pietro




--
View this message in context: http://postgis.17.x6.nabble.com/permission-denied-for-relation-geometry-columns-on-postgis-2-x-tp5003338.html
Sent from the PostGIS - User mailing list archive at Nabble.com.

Shp2pgsql : Error in fread()

Hugues François — 2013-05-20T12:23:24

Hello,

I have a big amount of data to import ins postgis from shapefile (18.000.000 of squares). As usual, I used shp2pgsql to do that with this command line : shp2pgsql -s 2154 -g geom -W LATIN1 -I /home/myhome/shpname myschema.shpname.

But for each of my file (there are four of them) it throw me the same error :
Error in fread() reading object of size 136 at offset 2147483756 from .shp file
Error reading shape object 15790321

In each of the four file the error happens to the same time / object.

Where did I make a mistake and how can I correct it ?

I found this related info : http://trac.osgeo.org/gdal/ticket/4203

But the answer don't really help me to find any solution... Is anyone there having an idea please ?

Hugues.

toTopoGeom take long time

Eva Linhartová — 2013-05-19T17:49:05

Hi,
I'm trying to insert geometry data to topology structure. I have 7000
polygons and the procedure take very long time now it's about 22 hours.

UPDATE plochy  SET topo = topology.toTopoGeom(geom, 'dp_topo',3,0.005);

I run the postgis database on my laptop:
procesor: Intel(R) Core(TM) i5-2430M CPU 2.40GHz
RAM: 8GB
type: 64bit

I'm newbie and I tryed it only on my sample data.
Is this time still ok or should I stop the query?

Thx
 Eva
_______________________________________________
postgis-users mailing list
postgis-users< at >lists.osgeo.org
http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users

postgis topology

Eva Linhartová — 2013-05-18T10:17:44

Hi,
I'm trying to insert geometry data to topology structure. I tryed it on my
sample data before and it was ok, but now if I use polygon table with 7000
polygons or 20000 lines can't work.

I get error like:

ERROR:  schema "dp_topo2" does not exist
LINE 1: SELECT nextval('DP_TOPO2.node_node_id_seq')

ERROR: schema "dp_topo2" does not exist
Stav SQL: 3F000

SET search_path TO dp, topology, public;
SELECT CreateTopology('DP_TOPO2',2065,0.005);  -- 1491 ms

SELECT AddTopoGeometryColumn('DP_TOPO2', 'dp', 'ku', 'topo', 'POLYGON');
 -- 301ms
SELECT AddTopoGeometryColumn('DP_TOPO2', 'dp', 'body', 'topo', 'POINT');
 -- 701ms
SELECT AddTopoGeometryColumn('DP_TOPO2', 'dp', 'plochy', 'topo',
'POLYGON'); -- 298ms
SELECT AddTopoGeometryColumn('DP_TOPO2', 'dp', 'linie', 'topo',
'LINESTRING'); -- 171ms


UPDATE ku SET topo = topology.toTopoGeom(geom, 'DP_TOPO2',1,0.005);  --
687ms
UPDATE body SET topo = topology.toTopoGeom(geom, 'DP_TOPO2',2,0.005); --
164274ms
UPDATE plochy  SET topo = topology.toTopoGeom(geom, 'DP_TOPO

Strange behavior

Richard Cernava — 2013-05-17T23:40:21

I'm trying to query if a point is present inside a MultiLineString. When I
run the following query I get very strange results in my opinion. If
someone could elaborate as to what's going on it would be very much
appreciated!

From this query I get a false and a true on the identical data. First I
check to see if the_geom intersects with my point and receive false. Second
I take the_geom and convert it to text and then compare it with my point
using ST_GeomFromText and get true. How is this possible? I've outputted
the result from the statement below.

SELECT fullname, ST_X(ST_StartPoint(ST_LineMerge(the_geom))),
ST_Y(ST_StartPoint(ST_LineMerge(the_geom))),
ST_AsText(the_geom) AS text,
ST_Intersects(the_geom, ST_GeomFromText('POINT(-122.837109 45.537746)',
4269)) AS fromdb,
ST_Intersects(ST_GeomFromText(ST_AsText(the_geom), 4269),
ST_GeomFromText('POINT(-122.837109 45.537746)', 4269)) AS fromtxt
from tiger_data.or_edges where   fullname like 'NW Norwich St' limit 1


"NW Norwich St";-122.837109;45.537746;"MUL

Getting argument of AND must not return a set, when trying to calculate the intersection of a end of a multiline and a polygon

Dave Potts — 2013-05-16T23:05:35

Hi,

I am trying to find the results of an intersection between a multi line 
string and a polygon,  I had assumed that saying something
like

select rt.source, 
av.geom,st_astext(st_startpoint((st_dump((rt.the_geom))).geom)) from 
route_table rt ,
avoid_areas av
where av.valid=true and rt.parent_port_id=80360 and
st_within(av.geom,st_startpoint((st_dump((rt.the_geom))).geom)) =true

might work where rt.the_geom is a multiline  and av.geom is a polygon

but I get the error,
ERROR:  argument of AND must not return a set
LINE 4: st_within(av.geom,st_startpoint((st_dump((rt.the_geom))).geo

any idea what I am doing wrong?

Dave

greetings!

manos voyatzakis — 2013-05-16T15:31:31

    http://www.antares.kh.ua/likeit.php?dxigmk782mvm










































































































































mvoyatzakis
manos voyatzakis
____________________
The new storage place was called Sure Lock Homes.
% _______________________________________________
postgis-users mailing list
postgis-users< at >lists.osgeo.org
http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users

Google Directions -> PostGIS

James David Smith — 2013-05-16T15:01:46

Hi there,

This isn't so much a question as such, more a request for opinions and
advice. I was wondering whether anyone had any experience with using the
Google API for getting driving directions/routing - but then storing the
results in a PostgreSQL/PostGIS database. I was thinking that my workflow
might go something like this:

1) Take table of trip_ids, start points, end points, travel modes and start
times from the database
2) Export them as a CSV
3) Import the CSV to R
4) Use R to loop through the data and send queries off to the Google API.
5) The results are a series of polylines each with the geometry as a
polyline, the length of the stage, the time of the stage etc. As a JSON
file. If I could break this apart in it's constiuent parts, could I then
store it back in my PostGIS database?

I guess the hardest bit I can see about this process, is converting the
google encoded polyline to a PostGIS polyline I think. Any thoughts on that?

Once this is all done I then want to use the routes with another l

flow analysis

Jonatan Malaver — 2013-05-14T16:40:05

Hello, is it possible to do any flow analysis? For example, query postgis for what down flow a certain point?

Thanks,

Jonatan Malaver

_______________________________________________
postgis-users mailing list
postgis-users< at >lists.osgeo.org
http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users

line direction

franco base — 2013-05-14T11:03:29

Hi.
I have a question for the list.

Is there a postgis function to have line direction?

For example ehen I use st_offset_curve
I have a correct offset according to direction,
So I have this information inside geometry
and I'd extract.

Thank

regards
fb
_______________________________________________
postgis-users mailing list
postgis-users< at >lists.osgeo.org
http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users

Status of Postgis-related rpms for Fedora 18?

Max Pyziur — 2013-05-13T17:35:41

Or should we wait for Fedora 19?


Max Pyziur
pyz< at >brama.com

PostGIS 2.1.0 beta2 is ready

Paragon Corporation — 2013-05-12T05:03:32

The 2.1.0beta2 release of PostGIS is now available.

The PostGIS development team is proud to release a feature complete beta
version of upcoming PostGIS 2.1.0. As befits a minor release, the focus is
on speed improvements, more features, and bug fixes. While this beta release
is feature complete, we expect some bugs and we'd appreciate it if you test
it before final release and report back with any issues you run into so we
can have a smooth release.
Please refer to our support page for instructions on reporting bugs

http://postgis.net/support 

If you are currently using PostGIS 2.0+, you can go the soft upgrade path:
ALTER EXTENSION postgis UPDATE TO "2.1.0beta2";

Users of 1.5 and below wil need to go the hard-upgrade path.

Best served with a bottle of GEOS 3.4.0dev and PostgreSQL 9.3beta1. (which
will also be released soon)

Source: http://download.osgeo.org/postgis/source/postgis-2.1.0beta2.tar.gz
Html docs: http://download.osgeo.org/postgis/docs/doc-html-2.1.0beta2.tar.gz
PDF: http://download.osgeo.

Generating geometry from a quadtree

Arnaud Vandecasteele — 2013-05-08T19:13:07

Hello Postgis Members,

I spent my day trying to find a nice way to solve my problem but I 
didn't find any.
I hope you will be able t help me.

Let's start by the beginning.
I've got a Gis Table that will be fill by polygon features.
These features can have different sizes and they can overlap them.
These features could be everywhere around the globe.
And some "places" can have a high density of polygons where overs not or 
few.

I'm trying to find a nice way to represent these features (and their 
attributes).
The easiest way will probably be to use a regular grid and to calculate 
the value of each cell using the value of the polygons.
But doing that I will probably have a lot of empty cells.
So I was thinking that something like a QuadTree representation [1] 
could be a better idea.

Do you have any idea how I could do that?

Thanks

Arnaud

[1] http://en.wikipedia.org/wiki/Quadtree

Need a method for "noding" a street network

Stephen Woodbridge — 2013-05-08T14:27:39

Hi all,

This question comes up reasonably often on the pgRouting list and has 
been posted he on occasion under titles like "How to break streets at 
intersections?"

It seems to me that this would be a good function to create in either 
postgis or pgrouting.

THE PROBLEM:

I have a table of 10's of thousands of street segments to 10's of 
millions of street segments. These street segments are LINSTRING or 
MULTILINESTRING geometries with some arbitrary number of attribute 
columns. The geometries may cross one another and are not noded 
correctly for use with pgRouting.

THE RESULTS:

We want to process the table and create a new table with the same 
structure (see comment about primary key below), and in the new table 
all the geometries are broken at intersections and all the new pieces of 
the original segment that have been broken have the original attributes 
propagated to them. So if the original segment has column foo='abc' and 
was split into 3 new segments, each of the three new segments would als

ST_CoveredBy() query run very very very slow

Pietro Rossin — 2013-05-07T14:30:32

Hello everybody
I'm trying to select features from two spatial tables, geometry 2D.

The first table is polylines from my whole area (984000 records).
The second table is a polygon table, 970000 records. These polygons were
created from the first table with an external tool (FME).

Lots of polylines weren't converted to polygons and now I want to extract
(into a new table) from the first table all the lines not converted to
polygon to perform additional operations.

The two table have indexes (gist) on the geom column

So, my query is:

SELECT a.id, a.layer, a.the_geom into newtable
  FROM my_first_table as a, my_second_sable as b
  where
  not ST_CoveredBy(a.the_geom, b.the_geom);

Yesterdary I let it run for 72.000.000ms and then I had to stop it.
I thought it was because of the Z values of all the geometries and the query
was:

SELECT a.id, a.layer, a.the_geom into newtable
  FROM my_first_table as a, my_second_sable as b
  where
  not ST_CoveredBy(st_force_2d(a.the_geom), st_force_2d(b.the_geom));

Today

ST_CoveredBy() query run very very very slow

Pietro Rossin — 2013-05-07T14:29:57

Hello everybody
I'm trying to select features from two spatial tables, geometry 2D.

The first table is polylines from my whole area (984000 records).
The second table is a polygon table, 970000 records. These polygons were
created from the first table with an external tool (FME).

Lots of polylines weren't converted to polygons and now I want to extract
(into a new table) from the first table all the lines not converted to
polygon to perform additional operations.

The two table have indexes (gist) on the geom column

So, my query is:

SELECT a.id, a.layer, a.the_geom into newtable
  FROM my_first_table as a, my_second_sable as b
  where
  not ST_CoveredBy(a.the_geom, b.the_geom);

Yesterdary I let it run for 72.000.000ms and then I had to stop it.
I thought it was because of the Z values of all the geometries and the query
was:

SELECT a.id, a.layer, a.the_geom into newtable
  FROM my_first_table as a, my_second_sable as b
  where
  not ST_CoveredBy(st_force_2d(a.the_geom), st_force_2d(b.the_geom));

Today

Clipping A portion of a raster file , in postgis

manish sharma — 2013-05-05T07:13:24

Dear all
Can anybody help me how to clip a raster in postgis, 

as me new to this ,

Regards
Manish Sharma



--
View this message in context: http://postgis.17.x6.nabble.com/Clipping-A-portion-of-a-raster-file-in-postgis-tp5003181.html
Sent from the PostGIS - User mailing list archive at Nabble.com.