[postgis-users] import mif file to postgis

2012-03-28 Thread 刘超
hello

i have some mif files,i want to import them to postgis directly,is it
possible?

BTW: i usually convert mif to shp,and convert shp to sql,then import sql,it
is too much steps~~

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


Re: [postgis-users] import mif file to postgis

2012-03-28 Thread Nicolas Ribot
 hello

 i have some mif files,i want to import them to postgis directly,is it
 possible?

 BTW: i usually convert mif to shp,and convert shp to sql,then import sql,it
 is too much steps~~



Hi
GDAL/OGR can do that.

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


Re: [postgis-users] More details for errors

2012-03-28 Thread Denis Rouzaud

Hi again!

Well, this is always the same problem, if I do linemerge(union()), I 
have to force multi to prevent errors. Then, simplifying to linestring 
seems to consume resources.


On way, would to do a view which determine if the union would lead to a 
multi or a linestring.

But I don't know how to achieve this.
I have no idea either if this would also be greedy on resources.

Denis

On 03/27/2012 05:43 PM, Nicolas Ribot wrote:

Yes, it is 2.0
http://postgis.refractions.net/documentation/manual-svn/ST_CollectionHomogenize.html

Well, I am turning around the problem...
I have to do a cast, as I used
ST_Multi(ST_LineMerge(ST_Union(wkb_geometry)))::geometry(MultiLineString,xxx)
in the first view to prevent errors of unconnected lines.

Right now the best solution I have is to let the view being unusable, and
another view to report errors.
This is not very satisfying as if a problem occurs, I cannot use my view
until I fix the geometric errors.

Basically, I need to compute a view with a LineMerge returning LineString
(i.e. leaving away features leading to MultiLineString) and another view
returning IDs of features leading to MULTI.

I thought about using ST_Touches but it only uses two geometries and I can
have more in the LineMerge. No idea what to look for...

Thanks for your help. I could go deeper in the problem!

Greetings


I would then create 3 views:
One doing the actual computation (linemerge(Union)), then one view
where st_numgeometries(the_geom) = 1, and one view where
st_numgeometries(the_geom)  1.

nicolas
___
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] adding JTS geometries to a PostGIS database

2012-03-28 Thread nicolas.neubauer
Hello,

I'm using the Java Topology Suite (JTS) and I want to store geometries with 
PostGIS. How do I do that? 

What works:
  statement = conn.prepareStatement(INSERT INTO sometable (geofield) values 
(?));
  org.postgis.PGgeometry m = new org.postgis.PGgeometry(myJtsGeometry.toText());
  statement.setObject(1, m);
  statement.execute();

What I would prefer:
  statement.setObject(1, myJtsGeometry);
  statement.execute(); 

According to 
http://postgis.refractions.net/documentation/javadoc/org/postgis/DriverWrapper.html,
 the DriverWrapper is supposed to take care of the wrapping. However, even when 
I add
  DriverWrapper.addGISTypes80((org.postgresql.PGConnection)conn);
I get the error message
  Can't infer the SQL type to use for an instance of 
com.vividsolutions.jts.geom.MultiPolygon. Use setObject() with an explicit 
Types value to specify the type to use.
I am connecting to 
  jdbc:postgresql_postGIS://localhost:5432/
The driver class is 
  org.postgis.DriverWrapper

What am I missing? Or is the first, working way the right one? Seems like a 
waste of time to me, converting binary - text - binary...
 
Thanks
Nicolas
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] postgis raster export problem in using gdal_translate

2012-03-28 Thread Pierre Racine
Which version of GDAL are you using?

Why do you use mode = 2 if there is only one row in the table?

Now you said the output was not as the original input. how is it different?

Pierre

 -Original Message-
 From: postgis-users-boun...@postgis.refractions.net [mailto:postgis-users-
 boun...@postgis.refractions.net] On Behalf Of mahavir trivedi
 Sent: Wednesday, March 28, 2012 12:25 AM
 To: PostGIS Users Discussion
 Subject: Re: [postgis-users] postgis raster export problem in using 
 gdal_translate
 
 
 
 On Tue, Mar 27, 2012 at 7:10 PM, Pierre Racine pierre.rac...@sbf.ulaval.ca
 wrote:
 
 
   Please send the complete raster2pgsql and gdal_translate command you
 used.
 
   Pierre
 
 
-Original Message-
From: postgis-users-boun...@postgis.refractions.net [mailto:postgis-
 users-
boun...@postgis.refractions.net] On Behalf Of mahavir trivedi
Sent: Tuesday, March 27, 2012 4:28 AM
To: postgis-users@postgis.refractions.net
Subject: [postgis-users] postgis raster export problem in using
 gdal_translate
   
dear
   
   
i have 7053x6078 pixel  georeferrance  tiff image.
   
it is import using raster2pgsql  ok
   
but export it using gdal_translate command output is
 uncleard(Strange) not as
original input.my out put format is same Tiff input
 
   ___
   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] Break lines at intersections

2012-03-28 Thread Pavel Iacovlev
Hello,

Anyone has any ideas/tips of how to break lines at intersections ?
I have a a layer of lines and I want to break them at intersections
with them selfs, also (not that important but would be nice) break
lines forming a collapsed loop, for example 0.0;1.0;0.0 is broken at
1.0.

Currently this can be done in grass with v.clean tool=break,rmdupl but
it's relatively hard to add GRASS to a web app.

So far I came up with this query and out of ideas right now:

SELECT DISTINCT
gid,
ST_AsText(intersect_geom),
ST_Line_Locate_Point(GeometryN(the_geom, 1), intersect_geom) as point_loc1
FROM

(SELECT a.gid, ST_Intersection(a.the_geom, b.the_geom) AS
intersect_geom, a.the_geom
FROM roads as a, roads as b
WHERE ST_Intersects(a.the_geom, b.the_geom) AND (a.gid = 1)
) AS res1

WHERE GeometryType(intersect_geom) = 'POINT' -- Need support for
MULTILINESTRINGS, take start and end point ?
ORDER BY point_loc1

which has the following result: http://jsfiddle.net/9RJ8a/
I have the line, and the point positions on the line, but I can't get
around how to create my NEW broken lines with
ST_Line_Substring(geometry a_linestring, float startfraction, float
endfraction)

Any suggestion of how to accomplished this ? is the approach correct ?
any other ideas ?
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] Break lines at intersections

2012-03-28 Thread Horst Düster
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1

Hi

if you work with PostGIS  2.0 you can use the cleangeometry() function.

http://www.kappasys.ch/pgtools/cleangeometry/cleanGeometry.sql

Regards
Horst

Dr. Horst Düster
Kappasys
Alte Bernstrasse 13
CH-4500 Solothurn

Tel.: +41 (0) 32 514 00 40
Mobil: +41 (0) 76 408 51 06
eMail: horst.dues...@kappasys.ch
skype: horstduester
Twitter: moazagotl
Web: www.kappasys.ch

lat: N 47° 12' 06.408 lon: E 7° 32' 13.740

Am 28.03.2012 12:29, schrieb Pavel Iacovlev:
 Hello,
 
 Anyone has any ideas/tips of how to break lines at intersections ?
 I have a a layer of lines and I want to break them at intersections
 with them selfs, also (not that important but would be nice) break
 lines forming a collapsed loop, for example 0.0;1.0;0.0 is broken at
 1.0.
 
 Currently this can be done in grass with v.clean tool=break,rmdupl but
 it's relatively hard to add GRASS to a web app.
 
 So far I came up with this query and out of ideas right now:
 
 SELECT DISTINCT
 gid,
 ST_AsText(intersect_geom),
 ST_Line_Locate_Point(GeometryN(the_geom, 1), intersect_geom) as point_loc1
 FROM
 
 (SELECT a.gid, ST_Intersection(a.the_geom, b.the_geom) AS
 intersect_geom, a.the_geom
 FROM roads as a, roads as b
 WHERE ST_Intersects(a.the_geom, b.the_geom) AND (a.gid = 1)
 ) AS res1
 
 WHERE GeometryType(intersect_geom) = 'POINT' -- Need support for
 MULTILINESTRINGS, take start and end point ?
 ORDER BY point_loc1
 
 which has the following result: http://jsfiddle.net/9RJ8a/
 I have the line, and the point positions on the line, but I can't get
 around how to create my NEW broken lines with
 ST_Line_Substring(geometry a_linestring, float startfraction, float
 endfraction)
 
 Any suggestion of how to accomplished this ? is the approach correct ?
 any other ideas ?
 ___
 postgis-users mailing list
 postgis-users@postgis.refractions.net
 http://postgis.refractions.net/mailman/listinfo/postgis-users
-BEGIN PGP SIGNATURE-
Version: GnuPG v1.4.11 (GNU/Linux)
Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org/

iQEcBAEBAgAGBQJPcwSiAAoJENJq1wkKM7pS6NwH/3PzIeSfjG29wfrPx4mLu7gS
6gAYWk6Lsfp4kUEoy6uJdtGz6iFn3vY2blz1cA9diUrjw00O8m5/8IVFuI0jdyLM
j6tzDbyB7jmXyF6F3TLKfnAakVGenhbcXTeHRL3ojAiKcq8UGsRDQYviQ3WZK9sY
EXP9wJIunhawILOfFXEO68AcRDPLZcaLtZAE/qj4yY+p538M/8PNfuK2rBw6fcVO
dPcug4wDewciPYQBcgKdbRpXtpfh9ytq4xYV16aAA8wVg4lFicldIm4cm8mytOrb
0bYYyFmAHtpBQ2Hv6A7eMxywzF8iQmu8oPl1I6ihPwol3VQ/Tws+nAKbfGKs3bs=
=FgwT
-END PGP SIGNATURE-
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] Break lines at intersections

2012-03-28 Thread Pavel Iacovlev
Not quite sure if I understand how ST_Node works, but I don't need all
the line segments of the linestring, I need the linestrings to brake
(divided into multiple linestrings while retaining the id attribute
preferably) where they intersect with other line in the same table.

I have a attached a small image to represent what I need:
Line 1 - will be divided into 3 lines
Line 2 - will be divided into 2 lines
Line 3 - will be divided into 2 lines

On Wed, Mar 28, 2012 at 2:47 PM, Sandro Santilli s...@keybit.net wrote:
 On Wed, Mar 28, 2012 at 01:29:20PM +0300, Pavel Iacovlev wrote:

 Anyone has any ideas/tips of how to break lines at intersections ?

 PostGIS 2.0 has an ST_Node() function for that.

 --strk;

  ,--o-.
  |   __/  |    Delivering high quality PostGIS 2.0 !
  |  / 2.0 |    http://strk.keybit.net - http://vizzuality.com
  `-o--'

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


Re: [postgis-users] Break lines at intersections

2012-03-28 Thread Pavel Iacovlev
As I understand from the code, this will create a node/point if my
line self intersects, but I need to brake my line into N new lines
where it intersects with any other line from the same table.

On Wed, Mar 28, 2012 at 3:31 PM, Horst Düster horst.dues...@kappasys.ch wrote:
 -BEGIN PGP SIGNED MESSAGE-
 Hash: SHA1

 Hi

 if you work with PostGIS  2.0 you can use the cleangeometry() function.

 http://www.kappasys.ch/pgtools/cleangeometry/cleanGeometry.sql

 Regards
 Horst

 Dr. Horst Düster
 Kappasys
 Alte Bernstrasse 13
 CH-4500 Solothurn

 Tel.: +41 (0) 32 514 00 40
 Mobil: +41 (0) 76 408 51 06
 eMail: horst.dues...@kappasys.ch
 skype: horstduester
 Twitter: moazagotl
 Web: www.kappasys.ch

 lat: N 47° 12' 06.408 lon: E 7° 32' 13.740

 Am 28.03.2012 12:29, schrieb Pavel Iacovlev:
 Hello,

 Anyone has any ideas/tips of how to break lines at intersections ?
 I have a a layer of lines and I want to break them at intersections
 with them selfs, also (not that important but would be nice) break
 lines forming a collapsed loop, for example 0.0;1.0;0.0 is broken at
 1.0.

 Currently this can be done in grass with v.clean tool=break,rmdupl but
 it's relatively hard to add GRASS to a web app.

 So far I came up with this query and out of ideas right now:

 SELECT DISTINCT
 gid,
 ST_AsText(intersect_geom),
 ST_Line_Locate_Point(GeometryN(the_geom, 1), intersect_geom) as point_loc1
 FROM

 (SELECT a.gid, ST_Intersection(a.the_geom, b.the_geom) AS
 intersect_geom, a.the_geom
 FROM roads as a, roads as b
 WHERE ST_Intersects(a.the_geom, b.the_geom) AND (a.gid = 1)
 ) AS res1

 WHERE GeometryType(intersect_geom) = 'POINT' -- Need support for
 MULTILINESTRINGS, take start and end point ?
 ORDER BY point_loc1

 which has the following result: http://jsfiddle.net/9RJ8a/
 I have the line, and the point positions on the line, but I can't get
 around how to create my NEW broken lines with
 ST_Line_Substring(geometry a_linestring, float startfraction, float
 endfraction)

 Any suggestion of how to accomplished this ? is the approach correct ?
 any other ideas ?
 ___
 postgis-users mailing list
 postgis-users@postgis.refractions.net
 http://postgis.refractions.net/mailman/listinfo/postgis-users
 -BEGIN PGP SIGNATURE-
 Version: GnuPG v1.4.11 (GNU/Linux)
 Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org/

 iQEcBAEBAgAGBQJPcwSiAAoJENJq1wkKM7pS6NwH/3PzIeSfjG29wfrPx4mLu7gS
 6gAYWk6Lsfp4kUEoy6uJdtGz6iFn3vY2blz1cA9diUrjw00O8m5/8IVFuI0jdyLM
 j6tzDbyB7jmXyF6F3TLKfnAakVGenhbcXTeHRL3ojAiKcq8UGsRDQYviQ3WZK9sY
 EXP9wJIunhawILOfFXEO68AcRDPLZcaLtZAE/qj4yY+p538M/8PNfuK2rBw6fcVO
 dPcug4wDewciPYQBcgKdbRpXtpfh9ytq4xYV16aAA8wVg4lFicldIm4cm8mytOrb
 0bYYyFmAHtpBQ2Hv6A7eMxywzF8iQmu8oPl1I6ihPwol3VQ/Tws+nAKbfGKs3bs=
 =FgwT
 -END PGP SIGNATURE-
 ___
 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] Break lines at intersections

2012-03-28 Thread Sandro Santilli
On Wed, Mar 28, 2012 at 03:46:03PM +0300, Pavel Iacovlev wrote:
 Not quite sure if I understand how ST_Node works, but I don't need all
 the line segments of the linestring, I need the linestrings to brake
 (divided into multiple linestrings while retaining the id attribute
 preferably) where they intersect with other line in the same table.

That's what ST_Node should give you, wrapped into collection.
You could then use ST_Dump to explode the collection into
components.

 I have a attached a small image to represent what I need:
 Line 1 - will be divided into 3 lines
 Line 2 - will be divided into 2 lines
 Line 3 - will be divided into 2 lines

The single-line case is likely bogus unless you have GEOS-3.3.2 or higher.

--strk;

  ,--o-. 
  |   __/  |Delivering high quality PostGIS 2.0 !
  |  / 2.0 |http://strk.keybit.net - http://vizzuality.com
  `-o--'

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


Re: [postgis-users] Break lines at intersections

2012-03-28 Thread Pavel Iacovlev
Thank you for you help Sandro,

As I understand something like this in theory:
SELECT ST_Node(ST_Collect(the_geom)) FROM roads;
Then I generate series from the result and using ST_Contains find to
what initial road geometry the segment belongs and get it's gid or
another attribute.

Also you know (or maybe someone else knows) is it possible to
accomplish this in PostGIS 1.5 ?

On Wed, Mar 28, 2012 at 3:51 PM, Sandro Santilli s...@keybit.net wrote:
 On Wed, Mar 28, 2012 at 03:46:03PM +0300, Pavel Iacovlev wrote:
 Not quite sure if I understand how ST_Node works, but I don't need all
 the line segments of the linestring, I need the linestrings to brake
 (divided into multiple linestrings while retaining the id attribute
 preferably) where they intersect with other line in the same table.

 That's what ST_Node should give you, wrapped into collection.
 You could then use ST_Dump to explode the collection into
 components.

 I have a attached a small image to represent what I need:
 Line 1 - will be divided into 3 lines
 Line 2 - will be divided into 2 lines
 Line 3 - will be divided into 2 lines

 The single-line case is likely bogus unless you have GEOS-3.3.2 or higher.

 --strk;

  ,--o-.
  |   __/  |    Delivering high quality PostGIS 2.0 !
  |  / 2.0 |    http://strk.keybit.net - http://vizzuality.com
  `-o--'

 ___
 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] Break lines at intersections

2012-03-28 Thread Nicolas Ribot
On 28 March 2012 15:43, Pavel Iacovlev iacovlev.pa...@gmail.com wrote:
 Thank you for you help Sandro,

 As I understand something like this in theory:
 SELECT ST_Node(ST_Collect(the_geom)) FROM roads;
 Then I generate series from the result and using ST_Contains find to
 what initial road geometry the segment belongs and get it's gid or
 another attribute.

 Also you know (or maybe someone else knows) is it possible to
 accomplish this in PostGIS 1.5 ?


Hi,
Yes it is. Look at st_union() and st_dump.
st_union will cut the lines at their intersection points:

select st_astext((st_dump(st_union(geometry))).geom) from lines;

Nicolas
attachment: Screen shot 2012-03-28 at 5.03.32 PM.png___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


[postgis-users] (no subject)

2012-03-28 Thread mahavir trivedi
hi

output image problem is that  whole image isn't display(first some
upper widthxheight of input image display)




with thanks and regards
mahavir
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] (no subject)

2012-03-28 Thread mahavir trivedi
hi

output image problem is that  whole image isn't display(first some
upper widthxheight of input image display)

gdal 1.9


with thanks and regards
mahavir

On 3/28/12, mahavir trivedi mahavir.triv...@gmail.com wrote:
 hi

 output image problem is that  whole image isn't display(first some
 upper widthxheight of input image display)




 with thanks and regards
 mahavir

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


Re: [postgis-users] LinearReferencing: measure vertex of linestring?

2012-03-28 Thread Sandro Santilli
On Wed, Mar 28, 2012 at 06:22:28PM +0200, Andrés Maneiro wrote:

  LINESTRINGM(1 0 1,2 0 2,4 0 4)
 
 What I'm wondering is if somehow I could calibrate the vertex of the
 linestring, having as a result, for example, this route:
 
 LINESTRINGM(1 0 1,2 0 *1.7*,4 0 4)

See ST_SetPoint

--strk;

  ,--o-. 
  |   __/  |Delivering high quality PostGIS 2.0 !
  |  / 2.0 |http://strk.keybit.net - http://vizzuality.com
  `-o--'

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


Re: [postgis-users] LinearReferencing: measure vertex of linestring?

2012-03-28 Thread Andrés Maneiro

On 28/03/12 18:26, Sandro Santilli wrote:

On Wed, Mar 28, 2012 at 06:22:28PM +0200, Andrés Maneiro wrote:


  LINESTRINGM(1 0 1,2 0 2,4 0 4)


Oh, I get it. And there is a native way in postgis to calibrate the line from a 
given set of measured points instead of using the ends?


--

If I understand it correctly I could:

 1) mylinestringm = ST_AddMeasure(mylinestring, measure_start, measure_end);
 2) Add the given set of points to mylinestringm with ST_AddPoint
 3) Using ST_SetPoint to move the original vertexes.

Is that correct? I'm very interested to know if there is a more direct and 
proper way to create a route from a set of points.


best,
Andrés
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] LinearReferencing: measure vertex of linestring?

2012-03-28 Thread Stephen Woodbridge

On 3/28/2012 1:26 PM, Andrés Maneiro wrote:

On 28/03/12 18:26, Sandro Santilli wrote:

On Wed, Mar 28, 2012 at 06:22:28PM +0200, Andrés Maneiro wrote:


LINESTRINGM(1 0 1,2 0 2,4 0 4)


Oh, I get it. And there is a native way in postgis to calibrate the line
from a given set of measured points instead of using the ends?

--

If I understand it correctly I could:

1) mylinestringm = ST_AddMeasure(mylinestring, measure_start, measure_end);
2) Add the given set of points to mylinestringm with ST_AddPoint
3) Using ST_SetPoint to move the original vertexes.

Is that correct? I'm very interested to know if there is a more direct
and proper way to create a route from a set of points.


ST_MakeLine [1] can be used to aggregate a sequence of points into a line.

[1] http://www.postgis.org/documentation/manual-1.5/ST_MakeLine.html
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] LinearReferencing: measure vertex of linestring?

2012-03-28 Thread George Silva
I guess the best way to do this is to combine the actual line and a
virtual line, probably a view, with measured points.

This involves a bit of thinking but it might be possible.

George

On Wed, Mar 28, 2012 at 3:03 PM, Stephen Woodbridge wood...@swoodbridge.com
 wrote:

 On 3/28/2012 1:26 PM, Andrés Maneiro wrote:

 On 28/03/12 18:26, Sandro Santilli wrote:

 On Wed, Mar 28, 2012 at 06:22:28PM +0200, Andrés Maneiro wrote:

  LINESTRINGM(1 0 1,2 0 2,4 0 4)


 Oh, I get it. And there is a native way in postgis to calibrate the line
 from a given set of measured points instead of using the ends?

 --

 If I understand it correctly I could:

 1) mylinestringm = ST_AddMeasure(mylinestring, measure_start,
 measure_end);
 2) Add the given set of points to mylinestringm with ST_AddPoint
 3) Using ST_SetPoint to move the original vertexes.

 Is that correct? I'm very interested to know if there is a more direct
 and proper way to create a route from a set of points.


 ST_MakeLine [1] can be used to aggregate a sequence of points into a line.

 [1] 
 http://www.postgis.org/**documentation/manual-1.5/ST_**MakeLine.htmlhttp://www.postgis.org/documentation/manual-1.5/ST_MakeLine.html

 __**_
 postgis-users mailing list
 postgis-users@postgis.**refractions.netpostgis-users@postgis.refractions.net
 http://postgis.refractions.**net/mailman/listinfo/postgis-**usershttp://postgis.refractions.net/mailman/listinfo/postgis-users




-- 
George R. C. Silva

Desenvolvimento em GIS
http://geoprocessamento.net
http://blog.geoprocessamento.net
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] LinearReferencing: measure vertex of linestring?

2012-03-28 Thread Paul Ramsey
Create a ticket describing the calibration, that's useful info and
would be a good utility for LRS users for a future version.

P.

2012/3/28 Andrés Maneiro amane...@icarto.es:
 Hello people,

 I'm reviewing what postgis is able to do in terms of Linear Referencing. As
 fas as I could see, to calibrate a linestring there is only the method
 ST_AddMeasure:
 http://www.postgis.org/docs/ST_AddMeasure.html

 So, as stated in the example, the operation will create a route like:

    SELECT ST_AsEWKT(ST_AddMeasure(
    ST_GeomFromEWKT('LINESTRING(1 0, 2 0, 4 0)'),1,4)) As ewelev;
               ewelev
    
     LINESTRINGM(1 0 1,2 0 2,4 0 4)

 What I'm wondering is if somehow I could calibrate the vertex of the
 linestring, having as a result, for example, this route:

    LINESTRINGM(1 0 1,2 0 *1.7*,4 0 4)

 The project I'm working on is related to road management, and it's very
 common that the ratio between a route's length and its measures is not
 constant. So, my doubt goes: is it somehow possible to create a route in
 postgis where the measure and the distance are not directly related?

 (For those who know how ArcGIS dynamic segmentation works, what I'd like
 to achive is something similar to calibrate a route using points [1]).

 Let me know if I explain myself properly.

 best,
 Andrés

 [1]
 http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/About_calibrating_route_measures_using_points/0039000t00/
 ___
 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] POSTGIS 2.0 on WIn

2012-03-28 Thread Piotr Pachół

Hello,
I tried all hints from this thread and that: 
http://postgis.refractions.net/pipermail/postgis-users/2012-March/033111.html

but it is still impossible to install PostGIS 2.0 on Windows XP.
I also copied libproj-0.dll in bin folder - it didn't help.

Do you have any other suggestion ?

Piotr.

W dniu 2012-03-27 04:26, Paragon Corporation pisze:

PROJSO  settings shouldn't be necessary anymore because libproj.dll is not
used.

I upgraded the install to use proj 4.8.0 so the dll is now libproj-0.dll
which is what gdal assumes by default.


It's possible I missed a dll and my tests systems already had them -- make
sure you have libproj-0.dll in your bin folder.


I think if postgis is in use it won't be able to copy the dll so you will
probably need to shut down your service / restart and then run the copy
script.

Hope that helps,
Regina
http://www.postgis.us




-Original Message-
From: postgis-users-boun...@postgis.refractions.net
[mailto:postgis-users-boun...@postgis.refractions.net] On
Behalf Of Pierre Racine
Sent: Monday, March 26, 2012 2:23 PM
To: PostGIS Users Discussion
Subject: Re: [postgis-users] POSTGIS 2.0 on WIn


Hi at all,
I wished test POSTGIS 2.0 on Win, so after installed on WinXP32   PostgreSQL 
9.1.1,  and following instruction for  postgis-pg91-binaries-2.0.
0beta4w32  but I get the following error
...C:\postgis-pg91-binaries-
2.0.0beta4w32\postgis-pg91-binaries-
2.0.0beta4w32C:\Programmi\PostgreSQL\9.1\bin\\psql  -d testpstg -c CREATE 
EXTENSION postgis; ERROR:  could not load library
C:/Programmi/PostgreSQL/9.1/lib/postgis-2.0.dll
: unknown error 126 .

Of course 1) postgis-2.0.dll is on the lib dir; 2) I stopped and
restared Postgres service; 3) I've tried also another installation dir for Postgers in 
order to prevent eventually space in name trouble, but...
I get the same problem


Any hints?

Try setting the PROJSO to the path of libproj.dll

Pierre


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


Re: [postgis-users] Break lines at intersections

2012-03-28 Thread Pavel Iacovlev
Thank you for your suggestion Nicolas Ribot, I will give it a go.

On Wed, Mar 28, 2012 at 6:04 PM, Nicolas Ribot nicolas.ri...@gmail.com wrote:
 On 28 March 2012 15:43, Pavel Iacovlev iacovlev.pa...@gmail.com wrote:
 Thank you for you help Sandro,

 As I understand something like this in theory:
 SELECT ST_Node(ST_Collect(the_geom)) FROM roads;
 Then I generate series from the result and using ST_Contains find to
 what initial road geometry the segment belongs and get it's gid or
 another attribute.

 Also you know (or maybe someone else knows) is it possible to
 accomplish this in PostGIS 1.5 ?


 Hi,
 Yes it is. Look at st_union() and st_dump.
 st_union will cut the lines at their intersection points:

 select st_astext((st_dump(st_union(geometry))).geom) from lines;

 Nicolas

 ___
 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] Making Tiger shapefiles routable

2012-03-28 Thread Hilton Long


Has a set of standard queries been developed to make Tiger shapefiles
routable?  It would need to eliminate gaps and ensure that all end nodes and
intersecting nodes corresponded with each other.

 

 

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


[postgis-users] PostGIS 2.0.0rc1

2012-03-28 Thread Paul Ramsey
The first release candidate of PostGIS 2.0.0 is now available:

  http://postgis.org/download/postgis-2.0.0rc1.tar.gz

We hope you'll test and let us know (http://trac.osgeo.org/postgis) if
you find problems!

Yours,

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


Re: [postgis-users] Making Tiger shapefiles routable

2012-03-28 Thread Stephen Woodbridge

On 3/28/2012 6:58 PM, Hilton Long wrote:

Has a set of standard queries been developed to make Tiger shapefiles
routable? It would need to eliminate gaps and ensure that all end nodes
and intersecting nodes corresponded with each other.


This might be better asked in the pgRouting list, but I'll address it 
here and cc that list.


I have loaded tiger data into pgrouting and you do not need to mess with 
the end points, just run the just process it like any shapefile and it 
will work for want it is worth.


The issue with Tiger is that the data is not really suitable for 
routing, I was using it for drivetime analysis where the results are 
approximate and that worked fine for my use cases.


Here are some of the problems with Tiger data:

* segments are not flagged for oneway streets
* there is no zlevel information at intersections ( so it will generate 
routes like drive to an overpass and turn right onto the underpass )
* the road classifications do not allow you to set the average speed for 
a road class very accurately, although this is minor compared to the 
first two.


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


Re: [postgis-users] POSTGIS 2.0 on WIn

2012-03-28 Thread Paragon Corporation
We'll be releasing a 2.0.0rc1 soon probably in the next day or so.

Regarding why it's still not working
1) Did you run one of the batch scripts.
2) Can you verify you have the following files in your bin folder
a) libgdal-1.dll
b) libgeos_c-1.dll
c) libgeos-3-3-3dev.dll
d) libproj-0.dll (that I assume you already checked)

3) and you have postgis-2.0.dll and rtpostgis-2.0.dll in your lib folder
(Note these must go in the lib folder not the bin folder)

Let us know if it still doesn't work with those.  We haven't recently tested
on windows xp so there
might be some dependency that exists on windows 2008 and windows 7 (which
are the two platforms we've been testing mostly with) that we don't have in
place that is needed by xp.
We'll try that tomorrow on an xp box if you are still having issues.

Thanks,
Regina and Leo
http://www.postgis.us




 -Original Message-
 From: Piotr Pachól [mailto:piotrpac...@gmail.com] 
 Sent: Wednesday, March 28, 2012 3:57 PM
 To: PostGIS Users Discussion
 Cc: Paragon Corporation
 Subject: Re: [postgis-users] POSTGIS 2.0 on WIn
 
 Hello,
 I tried all hints from this thread and that: 
 http://postgis.refractions.net/pipermail/postgis-users/2012-Ma
rch/033111.html
 but it is still impossible to install PostGIS 2.0 on Windows XP.
 I also copied libproj-0.dll in bin folder - it didn't help.
 
 Do you have any other suggestion ?
 
 Piotr.
 
 W dniu 2012-03-27 04:26, Paragon Corporation pisze:
  PROJSO  settings shouldn't be necessary anymore because 
 libproj.dll is 
  not used.
 
  I upgraded the install to use proj 4.8.0 so the dll is now 
  libproj-0.dll which is what gdal assumes by default.
 
 
  It's possible I missed a dll and my tests systems already 
 had them -- 
  make sure you have libproj-0.dll in your bin folder.
 
 
  I think if postgis is in use it won't be able to copy the 
 dll so you 
  will probably need to shut down your service / restart and then run 
  the copy script.
 
  Hope that helps,
  Regina
  http://www.postgis.us
 
 
 
  -Original Message-
  From: postgis-users-boun...@postgis.refractions.net
  [mailto:postgis-users-boun...@postgis.refractions.net] On 
 Behalf Of 
  Pierre Racine
  Sent: Monday, March 26, 2012 2:23 PM
  To: PostGIS Users Discussion
  Subject: Re: [postgis-users] POSTGIS 2.0 on WIn
 
  Hi at all,
  I wished test POSTGIS 2.0 on Win, so after installed on 
 WinXP32   PostgreSQL 9.1.1,  and following instruction for  
 postgis-pg91-binaries-2.0.
  0beta4w32  but I get the following error
  ...C:\postgis-pg91-binaries-
  2.0.0beta4w32\postgis-pg91-binaries-
  2.0.0beta4w32C:\Programmi\PostgreSQL\9.1\bin\\psql  -d 
 testpstg 
  -c CREATE EXTENSION postgis; ERROR:  could not load 
 library C:/Programmi/PostgreSQL/9.1/lib/postgis-2.0.dll
  : unknown error 126 .
 
  Of course 1) postgis-2.0.dll is on the lib dir; 2) I stopped and 
  restared Postgres service; 3) I've tried also another 
 installation dir for Postgers in order to prevent eventually 
 space in name trouble, but...
  I get the same problem
 
 
  Any hints?
  Try setting the PROJSO to the path of libproj.dll
 
  Pierre
 
 


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