You are very welcome,
it is the work of Paul Ramsey though ;-)
Cheers,
Rémi-C
2014-12-08 0:40 GMT+01:00 George Merticariu merticar...@gmail.com:
Hello Remi,
Your suggestion with point clouds worked perfectly.
Thank you!
Best regards,
George
On Wed, Dec 3, 2014 at 1:21 PM, Rémi Cura
For the second part,
it depends hugely on how you edit the field.
Are you editing with psql, pgadmin, qgis, ...
Cheers,
Rémi-C
2014-12-07 10:39 GMT+01:00 Hugues François hugues.franc...@irstea.fr:
Hello,
I’m not sure to understand well your problems but it seems you have sql
questions
Hey,
it is not mdb, but I know some people from Oslandia have successfully
developed a postgis -SQLite system with sync so to be able to work in the
field.
Maybe some stuff to re-use there (open source of course)
Cheers,
Rémi-C
2014-12-03 12:30 GMT+01:00 Robert Burgholzer rburg...@vt.edu:
I think this is there:
http://www.oslandia.com/qgis-versioning-plugin-en.html
Cheers,
Rémi-C
2014-12-03 18:09 GMT+01:00 Paolo Cavallini cavall...@faunalia.it:
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Il 03/12/2014 12:37, Rémi Cura ha scritto:
Hey, it is not mdb, but I know some
Hey,
some common answer about precision issue :
1. : offset your coordinate to gain in precision :
723120.23091093975 4379700.195452339 - 120.23091093975XXX ,
700.195452339
2. snap to grid you input data
3. check that your input geoms are valid : st_isvalid : st_makevalid()
Cheers,
Rémi-C
a completly
mess :(
I have tried to clean the layer, but I don't know how exactly should I
proceed.
Can anyone lead me in this task with some links or tips?
Thank you!!
On 12/02/2014 11:49 AM, Rémi Cura wrote:
Hey,
some common answer about precision issue :
1. : offset your coordinate
otherwise, you can use an aggregate function on the id field.
For instance min(emt_paradas.id)
This is not a postgis problem, but a sql problem ;-)
Cheers,
Rémi-C
2014-11-26 19:47 GMT+01:00 Nicolas Ribot nicolas.ri...@gmail.com:
Hi,
You can generate an id on-the-fly:
SELECT 1 as id,
:
vehicleid | status | time_start | time_end
I hope its understandable and not to mixed up...
Thanks!
Am Dienstag, 25. November 2014 16:06:33 UTC+1 schrieb Rémi Cura:
Hey, a small correction :
ST_MakeLine is already an aggregate, and you may want to enforce the order
inside the aggregate (see
Hey,
you may use ST_SNapToGrid
Cheers,
Rémi-C
2014-11-26 11:17 GMT+01:00 toni hernández t...@sigte.udg.edu:
Hello,
I need to round off all the coordinates of some polygon geometries.
So far I am using the st_asgeojson function to use the precision I need,
and then I create a new geometry
Hey, a small correction :
ST_MakeLine is already an aggregate, and you may want to enforce the order
inside the aggregate (see at the end).
Another interesting point is the possiblity to pu somehting in the M value
of each point of the line, for instance the time.
This comes very handy when you
Hey,
you don't even need postgis for this, you could store a cube per line,
with a x,y,z number, indexes on it,then querry like
SELECT *
FROM my_cube_table
WHERE x BETWEEN 100 AND 200
AND y BETWEEN 23 AND 45
AND z BETWEEN 45 AND 67
You could alos use postgres range type.
Of course you could store
. November 2014 16:06:33 UTC+1 schrieb Rémi Cura:
Hey, a small correction :
ST_MakeLine is already an aggregate, and you may want to enforce the
order inside the aggregate (see at the end).
Another interesting point is the possiblity to pu somehting in the M
value of each point of the line
Hey,
first you don't need plpgsql for this, you can do it with windows functions
(http://www.postgresql.org/docs/current/static/functions-window.html)
for the following I assume you have some order available on your point
data, the order is stored in id column :
CREATE TABLE my_azimuth_data AS
-19 12:27 GMT+01:00 Sandro Santilli s...@keybit.net:
On Tue, Nov 18, 2014 at 08:26:29PM +0100, Rémi Cura wrote:
This is suprising indeed.
Im my test using the v.out.postgis was still way faster than constructing
with postgis topology (about x3 I think).
Are you using v.out.postgis -l
Something like this?
SELECT zone,gid,state,name,ST_Asgeojson( geom )
FROM polys
WHERE ST_Overlaps(ST_GeomFromText('POLYGON((-109 42,-108 42, -108 43; -109
43, -109 42))', 4269), geom)=TRUE;
Cheers,
Rémi-C
2014-11-18 9:00 GMT+01:00 Hugues François hugues.franc...@irstea.fr:
Hello,
the edge_data trigger are triggered.
The correct way would be to compute full topology structure, then
deactivate trigger, then fill table , then reactivate trigger.
Cheers,
Rémi-C
2014-11-18 19:48 GMT+01:00 Sandro Santilli s...@keybit.net:
On Wed, Nov 12, 2014 at 10:51:34AM +0100, Rémi Cura wrote:
Note
http://www.postgis.us
--
*From:* postgis-users-boun...@lists.osgeo.org [mailto:
postgis-users-boun...@lists.osgeo.org] *On Behalf Of *Rémi Cura
*Sent:* Wednesday, November 12, 2014 11:41 AM
*To:* PostGIS Users Discussion
*Subject:* [postgis-users] no .bat
-boun...@lists.osgeo.org] *On Behalf Of *Rémi Cura
*Sent:* Thursday, November 13, 2014 4:03 AM
*To:* PostGIS Users Discussion
*Subject:* Re: [postgis-users] no .bat in windows zip installers
Hello Regina,
thanks for the answer.
I'm not a big fan of bat files,
but I do have to use the zip
Note you could have directly exported into postgis topology, which is slow
but very convenient (stil faster than postgis topology conversion).
I'm not a grass user,
from what I understood you have to use the -l switch in the command
*v.out.postgis*
Glad it helped =)
Cheers,
Rémi-C
2014-11-10
Hey,
first qgis can't deal with curve (i think),
so your data may be curved in your postgis but the visualisation would be
incorrect
so if you are going to see some curve in qgis, you have to convert the
curve to approximate theim using http://postgis.net/docs/ST_CurveToLine.html
,
then display
5608212.14184435))
So I'm assuming that its a convert problem, not a display problem.
I'm not sure how to use ST_LineToCurve. This is a road layer and
everything is a MULTILINESTRING, I can't tell the straight from the
curves. Can I?
Thanks,
-Andy
On 11/12/2014 3:58 AM, Rémi Cura wrote:
Hey
So it seems you are stuck with approximating curve with line inside FileGDB
before using ogr.
I guess you can do so using arcgis, but I never used it.
Cheers,
Rémi-C
2014-11-12 17:09 GMT+01:00 Paul Ramsey pram...@cleverelephant.ca:
Correct, the curves were likely ruined on the way through
Hey everybody.
The 2 bat files
makepostgisdb_using_extensions.bat
and
makepostgisdb.bat
seems to have vanished from the zip release file for version newer than
2.1.1 (postgis-bundle-pg93x64-2.1.1)
I'm I missing something?
thanks,
cheers,
Rémi-C
___
Just to answer about Grass.
I'm not a Grass user, and I do use PostGIS topology every day (heart of my
work)
i know only 3 commands :
v.in.ogr dsn=PG:host=localhost port=5433 dbname=... user=... password=...
layer=schema.table out=layer_name
v.clean ...
v.out.postgis -l -2 input=layer_name
Just for completeness sake,
if your input geom is really ambiguous with the precision you choose,
it may be faster in both human and machine time to use grass 7 to import
and clean the topology,
then export from grass 7 to postgis topology with the built-in grass
function.
I ressort sometime to
(emit warnings on downgrades
maybe).
On Wed, Oct 29, 2014 at 4:54 AM, Rémi Cura remi.c...@gmail.com wrote:
I'm sorry I think I was unclear.
I don't want to acces point,
I just want to change the pcid of a patch.
For the moment it seems to imply a complete destruction/reconstruction of
a new
A first raw plpgsql implementation :
Loop on all angle from convex hull edge.
https://github.com/Remi-C/_utilities/blob/master/postgis/rc_oriented_bbox_deom_axis.sql
Cheers,
Rémi-C
2014-10-25 17:06 GMT+02:00 Rémi Cura remi.c...@gmail.com:
Of course I plan to use convex hull,
I didn't
postgis-users-boun...@lists.osgeo.org on behalf of Rémi Cura
remi.c...@gmail.com
*Sent:* Tuesday, October 28, 2014 1:33 PM
*To:* PostGIS Users Discussion
*Subject:* Re: [postgis-users] pointcloud : change pcid ?
Now the only solution I found is to break a patch into point, then break
Hey,
any pointcloud users has an idea about changing the pcid of a patch without
breaking it into pieces?
Cheers,
Rémi-C
___
postgis-users mailing list
postgis-users@lists.osgeo.org
http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users
...@lists.osgeo.org on behalf of Rémi Cura
remi.c...@gmail.com
*Sent:* Tuesday, October 28, 2014 11:01 AM
*To:* PostGIS Users Discussion
*Subject:* [postgis-users] pointcloud : change pcid ?
Hey,
any pointcloud users has an idea about changing the pcid of a patch
without breaking it into pieces
, 2014 at 7:09 PM, Rémi Cura remi.c...@gmail.com wrote:
Ok if I understand right,
implementing it at the sql level would be :
- compute convex envelop
- order edges by angle
- round angle to a given precision, keep distinct values
- for each angle, rotate and compute bbox and area
keep
doing the evaluation.
-Steve
On 10/25/2014 8:52 AM, Rémi Cura wrote:
Found a python version here for minimal area enclosing rectangle
http://docs.opencv.org/modules/imgproc/doc/
structural_analysis_and_shape_descriptors.html#minarearect
It takes a numpy 2D array as point
That means
To: PostGIS Users Discussion
Subject: Re: [postgis-users] Oriented BBox Efficient computing?
On Thu, Oct 23, 2014 at 06:17:56PM +0200, Rémi Cura wrote:
Wov thanks everybody :
Seems it possible to do it fully geometrically :
http://cgm.cs.mcgill.ca/~orm/rotcal.html
Nice, would be a good
Hey list,
somebody has a good idea about how efficiently computing an approximate
oriented bbox for a polygon ?
The simplest approaches I can think of are
- iteratively rotate the geom by few degrees, take the envelop. Keep
ration where envelope is min
- using Principal Components Analysis :
/23/2014 9:33 AM, Rémi Cura wrote:
Hey list,
somebody has a good idea about how efficiently computing an approximate
oriented bbox for a polygon ?
The simplest approaches I can think of are
- iteratively rotate the geom by few degrees, take the envelop. Keep
ration where envelope is min
Hm I don't understand.
The point is to get this bouding rectangle axis.
So why compute convex hull?
Cheers,
Rémi-C
2014-10-23 16:23 GMT+02:00 Sandro Santilli s...@keybit.net:
On Thu, Oct 23, 2014 at 04:15:06PM +0200, Rémi Cura wrote:
_I'm precisely looking for 2 axis (convex hull
, Rémi Cura wrote:
Hey list,
somebody has a good idea about how efficiently computing an approximate
oriented bbox for a polygon ?
The simplest approaches I can think of are
- iteratively rotate the geom by few degrees, take the envelop. Keep
ration where envelope is min
- using Principal
:57 AM, Håvard Tveite havard.tve...@nmbu.no wrote:
Search for Rotating calipers + width
Håvard
On 2014-10-23 15:33, Rémi Cura wrote:
Hey list,
somebody has a good idea about how efficiently computing an approximate
oriented bbox for a polygon ?
The simplest approaches I can think
Your querry doesn't handle properly the case when one line has to be split
by multiple other lines.
a simple list
_are the geom valid (ST_ISValid / ST_MakeValid)
_are the geom in the same referential
_do you have a recent postgis/geos version (postgis_full_version())
_are you geometry at the
Exporting using shapefile is a bad idea.
If anything in your data base is relational, it will be lost.
Same for special indexes, triggers, permission, etc etc.
For the moment it's quit painfull to use dblink or fdw because you have to
do it table by table.
(sur it can be automatized, but it is
Hey,
it is tricky to use qgis value directly, because depending on the version
of qgis (and gdal?) you may have errors with (multi , 3D, M ).
I would advise to strongly enforce a type to any input from qgis.
(for instance ST_SnapToGrid( ST_Force2D(
ST_GeometryN(ST_CollectionExtract(geom,1),1) )
For completness , I think CGAL has a merge functionnality on its half edge.
cheers,
Rémi-C
2014-10-13 18:38 GMT+02:00 Sandro Santilli s...@keybit.net:
On Mon, Oct 13, 2014 at 06:29:49PM +0200, Christopher Mutel wrote:
Dear all-
Is there a way to create a new topology geometry as the union
Hey,
seems like a typical job for postgis topology...
Build your topology using AddLinestring
http://postgis.net/docs/TopoGeo_AddLineString.html
Then you have easy access to shared part, intersection so.
Cheers,
Rémi-C
2014-10-10 23:45 GMT+02:00 Zenon Panoussis ora...@provocation.net:
Hello
Hey,
snap is deceiptively named.
As far as I know It just snaps points to points (or points in geometry to
points in geometry).
For exemple, if you draw 2 square side by side, one being bigger than the
other, and separate theim with small distance, snap will do nothing.
Cheers,
Rémi-C
2014-09-24
Thank you,
it was exactly this!
Cheers,
Rémi-C
2014-09-19 20:03 GMT+02:00 Nicolas Ribot nicolas.ri...@gmail.com:
Hi,
-D option should appear just after shp2pgsql, not after the shape and
table parameters, no ?
Nicolas
On 19 September 2014 19:29, Rémi Cura remi.c...@gmail.com wrote:
Hey
Hey list,
I'm a bit ashamed to be stuck with a seemingly trivial problem :
executing
`
/usr/lib/postgresql/9.3/bin/shp2pgsql 2_poles.shp cmm.temp_import_poles -D
`
output somehting like
`
create table
INSERT INTO cmm.temp_import_poles
Hey,
I think you misunderstood me,
of course array of geom are not a replacement of proper use of aggregate.
It's just that in some case you cannot do with aggregate and array function
are very usefull, or the only solution.
(example :
WITH A AS (
)
for instance when
2014-09-01 14:39 GMT+02:00
sorry was sent too early
2014-09-01 15:18 GMT+02:00 Rémi Cura remi.c...@gmail.com:
Hey,
I think you misunderstood me,
of course array of geom are not a replacement of proper use of aggregate.
It's just that in some case you cannot do with aggregate and array
function are very usefull
. For
this you need a measure of distance, you can simply take euclidian distance
(regular postgis distance) for a start.
Cheers,
Rémi-C
2014-08-27 18:38 GMT+02:00 Lee R sregdor...@gmail.com:
On Monday, August 25, 2014 4:46:11 AM UTC-4, Rémi Cura wrote:
You want to map your dimension
(time
You want to map your dimension
(time, note_pitch) to (X,Y) ?
It seems really difficult to achieve if for instance you would consider a C
and a C up an octave closer than a C and a D.
And if you consider
(time,note_pitch) to (X=time, Y = [A-G], Z = [1-8]) ,
Lot's of function in postgis doesn't
Hey,
you can do it using plpython and python gdal binding.
It could be a first step to evaluate if it is fast enough for you.
Cheers,
Rémi-c
2014-08-22 22:09 GMT+02:00 Rasmus Aveskogh ras...@defero.se:
Hi!
How hard/feasible would it be to implement elevation profile/drape line
functionality
Hey,
I did something akin, I decided then for a 2 step process : find potential
matches, then keep none or the better match.
The problem is: given 2 tables A (old) and B(new), find a mapping between A
and B,
with the following limits: a building could be new, destroyed , same or
changed. (Note
From what I have seen (not much),
finding which tile is in your shape is proven to be very fast given the
right index, table structure, server tuning so.
Then if the image is in the server it is going to be slow to get the data.
(I guess if the image is in filesystem this should be not so slow).
You have several ways.
Essentially,
look here toward the end : https://github.com/pramsey/pointcloud to do it
using PDAL
You can use another way that is
LAS File *--las tools--*
CSV File *---COPY statement--*
table(X,Y,Z,Intensity...) *-- PC_Patch(PC_Point(X,Y,Z ..))--*
patch table
Of course you
oups,
forget you use windows.
My code is bash script, won't be able to use it.
However the sql inside the bash may be usefull.
Cheers,
Rémi-C
2014-08-20 19:37 GMT+02:00 Rémi Cura remi.c...@gmail.com:
You have several ways.
Essentially,
look here toward the end : https://github.com/pramsey
Hey,
small change if you stick with your original querry (which looks fine).
You may want to force the planner to apply
*points.geom 'BOX(425930 7197112, 429605.074 7200582.906)'::BOX2D=TRUE;*
before the filtering (to be tested).
If your polygon are well clustered, you could also simply estimate
11:30 GMT+02:00 Rémi Cura remi.c...@gmail.com:
Hey,
small change if you stick with your original querry (which looks fine).
You may want to force the planner to apply
*points.geom 'BOX(425930 7197112, 429605.074 7200582.906)'::BOX2D=TRUE;
*
before the filtering (to be tested).
If your
Oups, typo :
Something like (not tested)
SELECT id
FROM points
EXCEPT
SELECT DISTINCT pt.id
FROM points AS pt
*LEFT JOIN* polygons AS pol ON (ST_Intersects(pt.geom, pol.geom)=TRUE)
2014-08-12 12:47 GMT+02:00 Rémi Cura remi.c...@gmail.com:
Hm,
also,
if you feel more comfortable,
instead
+0200, Rémi Cura wrote:
Hey,
I noticed that QGIs is introducing slight errors when editing topology
edge_data layer (with the regular edition, no pluggin ).
This is dangerous because function such a s ST_getFaceGeometry won't work
afterward (in some case).
Are you aware of the QGIS built
The answer may depend on the coverage of points you have.
1.) full postgis
If you have a lot of points (meaning few hex without data), you can simply
try to give the 0 hex a value following a function of the value of the
neighbours.
For instance, take the N nearest neighbour (called NN) that
that lead here. Best,
Luís
On 28 July 2014 17:42, Rémi Cura remi.c...@gmail.com wrote:
Hmm, no luck :-(
If you don't want to take any chance,
don't delete the .so, but rename them (like libgeos_c.so.backup).
After any change on your lib, you should ldconfig anyway.
Cheers,
Rémi-C
Hey,
it would be a very bad idea to use
set enable_seqscan to off;
for other things that analyzing planner behaviour.
(in perticular, but not limited to , if you do other think that this simple
querry).
You may also want to check your postgres conf file regarding performances.
(and round a little
Hm could you try the following command :
`sudo ldconfig`
It will refresh your dynamic link to library.
Cheers,
Rémi-C
2014-07-28 9:53 GMT+02:00 Luís de Sousa luis.a.de.so...@gmail.com:
Hi Clifford, thank you for replying. I found an old version of
liblwgeom installed, but after removing it
Hmm, no luck :-(
If you don't want to take any chance,
don't delete the .so, but rename them (like libgeos_c.so.backup).
After any change on your lib, you should ldconfig anyway.
Cheers,
Rémi-C
2014-07-28 17:26 GMT+02:00 António M. Rodrigues amcrgrodrig...@gmail.com:
Hi,
A while back I had
Hey
ST_X
and ST_Y
please take a look here when you look for function :
http://postgis.net/docs/reference.html
Cheers,
Rémi-C
2014-07-17 14:46 GMT+02:00 Antti J. Lind anttij.l...@birdlife.fi:
Hi,
We have previously used postgis 1.3.6 with postgres 8.4.1.
I installed postgres 9.3.4 and
...@mybofy.net:
Comment me désincrire de cette liste ?
J'ai TOUT essayé.
Le 17 juil. 2014 à 15:52, Rémi Cura remi.c...@gmail.com a écrit :
Hey
ST_X
and ST_Y
please take a look here when you look for function :
http://postgis.net/docs/reference.html
Cheers,
Rémi-C
2014-07-17 14:46 GMT+02:00
Hey,
If I take into account your last mail,
you probably have forget to transform your data.
If you have a table in srid 4326
, you can't use a transform in your querry if you want it to uses indexes.
OR you have to buil and index like this :
CREATE INDEX ON sdgis.parcels4326 USING GIST
On Thu, Jul 10, 2014 at 1:11 AM, Rémi Cura remi.c...@gmail.com wrote:
Hey,
If I take into account your last mail,
you probably have forget to transform your data.
If you have a table in srid 4326
, you can't use a transform in your querry if you want it to uses indexes.
OR you have
:
The 2 step query is an interesting idea. Would I need to do that in 2
queries, or can that be accomplished in a single statement?
On Thu, Jul 10, 2014 at 8:58 AM, Rémi Cura remi.c...@gmail.com wrote:
Hey,
I'm glad you found out.
I never tried but you could probably try to index
)
AND (ST_DWITHIN(mip.geom::geography,p.geom::geography,
100)=TRUE)
)
Cheers,
Rémi-C
2014-07-10 18:33 GMT+02:00 Rémi Cura remi.c...@gmail.com:
A single statement
, the 2 step is only metaphorical (or can be enforced if you use CTE).
Somehting liek this (untested)
WITH my_input_polygon
the
results. I really appreciate the suggestions.
On Thu, Jul 10, 2014 at 9:36 AM, Rémi Cura remi.c...@gmail.com wrote:
Oups, sent it to early.
If you trust the planner you can simply try :
WITH my_input_polygon AS (
SELECT geom
FROM sdgis.parcels4326 AS parcels
WHERE apn = 3500600300
Hey,
just a trick I uses to make sure everything will work :
_first check befor that every geom is valid, as a length/area, etc etc
_second if you can node your geom before going ot topology
_when converting to topogeom, uses an order by on your data.
_you can modify to topogeom function to add
Or shapely, still in python.
You can also use Geos I guess, in C.
You can also go the command line way with ogr2ogr and gdal
(there are python and R binding).
I'm pretty sure you can find also web oriented parser.
Cheers,
Rémi-C
2014-07-01 14:38 GMT+02:00 Martijn Meijers
Hey,
your question is to imprecise to get a good answer (better have a schema /
test case, and better describe your problems).
If you mean : given a set of points, find the 5 polygons that have the same
number of points (up to 1 if n_points%5 !=0), adding possibly a constraint
on there shape
The obvious answer here is to use PLR/R (http://www.joeconway.com/plr/).
There are many many ways to perform clustering (
http://cran.r-project.org/web/views/Cluster.html),
be it supervised or un-supervised, or even using machine learning (
that originated 0.001 meters
outside of the polygon!
On Fri, May 23, 2014 at 9:41 AM, Rémi Cura remi.c...@gmail.com wrote:
The idea is ST_Within(A,B) is true if no point (mathematical sense) of A
is outside of B.
Now, when splitting a line with another one (the boundary of the polygon
I can use QGis to see postgis raster if the version of qgis is new enough
(2.0 not working, 2.2 and 2.3 working)
,a s well as a very recent gdal version (1.10).
Still it is very slow.
Cheers,
Rémi-C
2014-05-20 9:44 GMT+02:00 dandrigo lcel...@latitude-geosystems.com:
Hello,
thank you. Would
I'm afraid the only way to visualyze raster fast enough in GQIs I found is
to export the raster as tif :-(
Of course using gdal is not an option because it is very slow,
so I use the Pierre Racine trick with plpythonu (I lost the source but all
credits to him of course).
This is order of magnitude
You can try to analyze a bit to understand better
(for the following I consider your table doesn't contains too much data)
with the_data AS (
SELECT the_geom, ST_ISValid(the_geom) AS is_valid,
ST_IsValidReason(the_geom)
FROM ft_potencial
)
SELECT *, ST_AsText(ST_MakeValid(the_geom)) AS
Sorry it didn't work :-( . Can you check that temp_id is not always 1?
It is more that this approach is flawed
_depending on the raster resolution and the small turn in your line (fixe =
simplify line).
_anyway, you are constructing a new line from centroid, and not putting the
height information
Hey,
the order by has to be on the same statement of your ST_makeline, or the
order could be wrong.
Morevover,
the order by should order the Point3D CTE because you want (points3d.geom)
to be ordered, not line or others.
Cheers,
rémi-C
2014-05-11 13:30 GMT+02:00 Hugues François
You can also cnahge the type of the column of the table afterward to be
more precise than geometry(0).
(see ALTER TABLE ...)
Cheers,
Rémi-C
2014-05-07 16:11 GMT+02:00 Nicolas Ribot nicolas.ri...@gmail.com:
Hi,
You could also cast the call to st_makeline into a specific geometry type,
with
Hmm your problem is not about a function, is about using invalid geometries
!
Why simply don't fix your geometry?
you must do it anyway if you want to use PostGIS
without strange bugs.
Use http://postgis.refractions.net/docs/ST_MakeValid.html
Cheers,
Rémi-C
2014-05-03 15:51 GMT+02:00 Larry
a polygon is a common situation, even with valid
geometries.
Depending of use case, I use ST_PointOnSurface when
ST_Within(St_Centroid(surface),
surface) returns false
steve
*Rémi Cura remi.c...@gmail.com remi.c...@gmail.com@lists.osgeo.org
http://lists.osgeo.org*
Envoyé par
David Smith james.david.sm...@gmail.com:
Hey Remi,
I don't understand what you mean? Why do I need to translate my data?
Thanks
James
On 28 April 2014 10:25, Rémi Cura remi.c...@gmail.com wrote:
I'm obsessed with precision,
so I would say translate your data t = ( -53 ,-178000
Hey,
if you look at the error you get :
LINESTRING (530395 178004, 530396 178004) and
LINESTRING (530396 178004, 530396 178004)
at (530395.5457342, 178004.12613484744)
You see there is no reason to output this error, in fact those two
linestring intersect only on (530396 178004),
yet geos
Hey,
If you just want visualisation it seem sa waste to do a lot of computing,
(i.e if you will do no quantitative use )
simply load the linestring in QGIS (version =2),
set the right width for the line in style and *use the transparency options
selecting the darkening rule for object vs object
Hey,
have a look at ST_CollectionExctract.
http://postgis.refractions.net/docs/ST_CollectionExtract.html
Cheers,
Rémi-C
2014-04-25 19:08 GMT+02:00 isla...@infomed.sld.cu:
Hi people i need help in converting a shape that have a GeometryCollection
data to a multipolygon shape, at least a
strategy are not fully integrated.
beyond that). This works with mounted network disk as well, which is even
nicer.
The whole postgres data folder can be put on network disk (with caveat),
which is very usefull !
GL
On 04/17/14, *Rémi Cura * remi.c...@gmail.com remi.c...@gmail.comwrote
Hey,
It may be unrelated, but have you checked your gdal version?
Another question, may it be encoding related (any special character, space
etc in file name ?)
Cheers,
Rémi-C
2014-04-17 5:24 GMT+02:00 droletguillaume droletguilla...@gmail.com:
Sent from Samsung Mobile
Hi,
The
would be less painful ?
2014-03-10 10:54 GMT+01:00 Rémi Cura remi.c...@gmail.com:
Hey Dear List,
I would appreciate some advice about the best way to store my raster :
1 million tiles,
50*50 pixels each (1 m2 or less in real world), around 24 bands (mostly
doubles)
,in db.
About half
Command options :
http://www.bostongis.com/pgsql2shp_shp2pgsql_quickguide_20.bqg
The schema is specified in the same time as table name.
Regarding the database, you have to specify it when using psql with the
-d option
You can also do this in one step using the shp2pgsql-gui .
Cheers,
Rémi-C
Please note that
_ogr2ogr is faster (5 times?)
_ogr2ogr is less reliable (bugs!)
Cheers,
Rémi-C
2014-04-16 15:12 GMT+02:00 J. Gustavo Rocha j...@di.uminho.pt:
Oi Rogério,
I use ogr2ogr, which is more generic than shp2pgsql. You can choose the
schema (active_schema option) and the name of
2014-04-15 16:58 GMT+02:00 Sandro Santilli s...@keybit.net:
On Tue, Apr 15, 2014 at 04:40:34PM +0200, Reijer Copier wrote:
Thanks!
We are currently considering using PostGIS Topology only to maintain
a topology aware copy of the data (to be used for analyses where
topology matters) and
Hey,
I created the ticket as asked :
http://trac.osgeo.org/postgis/ticket/2700
Cheers,
Rémi-C
2014-04-02 12:06 GMT+02:00 Sandro Santilli s...@keybit.net:
On Wed, Apr 02, 2014 at 11:13:02AM +0200, Rémi Cura wrote:
Ok,
now it works,
I added a filter on geom to filter out null geometries
Hey dear users,
I'm trying to define a new srid using a grid
(french RAF09.gtx).
I put this grid in the
/usr/share/proj folder in my Ubuntu 12.04,
now when I want to transform to this srid I get an error saying that this
grid file can't be found.
I found no direct clues in proj4 doc, as I use
Problem solved :
one has to create a 'nad' dir in the /usr/share/proj,
then put .gtx in nad dir.
Restart of server seems mandatory.
Cheers,
Rémi-C
2014-04-03 10:48 GMT+02:00 Rémi Cura remi.c...@gmail.com:
Hey dear users,
I'm trying to define a new srid using a grid
(french RAF09.gtx).
I
Ok,
now it works,
I added a filter on geom to filter out null geometries.
It seems that pgsql2shp doesn't like NULL geometries.
Cheers,
Rémi-C
2014-03-31 17:53 GMT+02:00 Rémi Cura remi.c...@gmail.com:
I should add that using ogr2ogr straight, it works !
Cheers,
Rémi-C
2014-03-31 17:47
Hey,
my pgsql2shp doesn't output 3D (but 2D).
Am I missing something ?
/usr/lib/postgresql/9.3/bin/pgsql2shp -f registration_surf_3D -h localhost
-p 5433 -u postgres -P uGuess -g geom_surf -r import_xml
public.registration_joined_with_trajectory
I have the latest postgis/geos/postgres/gdal.
Of
,931008)
)
What do you think?
Cheers,
Rémi-C
2014-03-31 17:14 GMT+02:00 Sandro Santilli s...@keybit.net:
On Mon, Mar 31, 2014 at 05:06:57PM +0200, Rémi Cura wrote:
Hey,
my pgsql2shp doesn't output 3D (but 2D).
Am I missing something ?
/usr/lib/postgresql/9.3/bin
101 - 200 of 346 matches
Mail list logo