[postgis-users] Get point certain degree of another point

2010-01-22 Thread tommy408

How can I get a point at a certain distance and degree from another point. 
For example 10 meter distance to another point at 45 degrees.

I'm thinking, draw some kind of imaginary circle around point A with radius
R, convert it to a line, and use ST_Line_Interpolate_Point.

Thank you
-- 
View this message in context: 
http://old.nabble.com/Get-point-certain-degree-of-another-point-tp27282672p27282672.html
Sent from the PostGIS - User mailing list archive at Nabble.com.

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


Re: [postgis-users] RE Importing shapefiles to PostGIS

2010-01-22 Thread Shreerang Patwardhan
Hello all,
I have followed steps to import TIGER shapefile to PostGIS . I  got
the following output and errors.

shreer...@ubuntu:~/Shapefiles/06001_Alameda_County/tl_2008_06001_arealm$
shp2pgsql -s 4326 -I -S tl_2008_06001_arealm.shp public.test_db >
tl_2008_06001_arealm.sql
Shapefile type: Polygon
Postgis type: POLYGON[2]

shreer...@ubuntu:~/Shapefiles/06001_Alameda_County/tl_2008_06001_arealm$
psql -d test_db -U postgres -f import.sql
psql: FATAL:  Ident authentication failed for user "postgres"

Need help!!!


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


[postgis-users] Get point certain degree of another point

2010-01-22 Thread tommy408

How can I get a point at a certain distance and degree from another point. 
For example 10 meter distance to another point at 45 degrees.

I'm thinking, draw some kind of imaginary circle around point A with radius
R, convert it to a line, and use ST_Line_Interpolate_Point.

-- 
View this message in context: 
http://old.nabble.com/Get-point-certain-degree-of-another-point-tp27282666p27282666.html
Sent from the PostGIS - User mailing list archive at Nabble.com.

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


Re: [postgis-users] PROBLEM CONVERTING AND MANIPULATING LONGITUDE ANDLATITUDE COORDINATES

2010-01-22 Thread Chris Hermansen

Excellent, congratulations!

MILMIL wrote:

Thank u very much Mr. Hermansen. Few words but very wish. I follow your
suggestion and I solved my problem. I deeply appreciate.
Best regards,

Miltos Miltiadous

Chris Hermansen wrote:
  

Create your circle with radius R metres after you transform your
geographic point to UTM, perhaps?
Chris Hermansenchris.herman...@timberline.ca
tel+1.604.714.2878 · fax+1.604.733.0631 · mob+1.778.840.4625
Timberline Natural Resource Group · www.timberline.ca
401 · 958 West 8th Avenue  · Vancouver BC · Canada · V5Z 1E5


-Original Message-
From: "Miltos Miltiadous" 
Date: Fri, 22 Jan 2010 13:40:18 
To: 

Subject: [postgis-users] PROBLEM CONVERTING AND MANIPULATING LONGITUDE AND
LATITUDE COORDINATES

___
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





  



--
Regards,

Chris Hermansen · chris.herman...@timberline.ca · skype:clhermansen
tel+1.604.714.2878   ·  fax+1.604.733.0631   ·   mob+1.778.840.4625
Timberline Natural Resource Group Ltd   ·  http://www.timberline.ca
401  ·  958 West 8th Avenue  ·  Vancouver BC  ·  Canada  ·  V5Z 1E5

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


Re: [postgis-users] PROBLEM CONVERTING AND MANIPULATING LONGITUDE ANDLATITUDE COORDINATES

2010-01-22 Thread MILMIL

Thank u very much Mr. Hermansen. Few words but very wish. I follow your
suggestion and I solved my problem. I deeply appreciate.
Best regards,

Miltos Miltiadous

Chris Hermansen wrote:
> 
> Create your circle with radius R metres after you transform your
> geographic point to UTM, perhaps?
> Chris Hermansenchris.herman...@timberline.ca
> tel+1.604.714.2878 · fax+1.604.733.0631 · mob+1.778.840.4625
> Timberline Natural Resource Group · www.timberline.ca
> 401 · 958 West 8th Avenue  · Vancouver BC · Canada · V5Z 1E5
> 
> 
> -Original Message-
> From: "Miltos Miltiadous" 
> Date: Fri, 22 Jan 2010 13:40:18 
> To: 
> Subject: [postgis-users] PROBLEM CONVERTING AND MANIPULATING LONGITUDE AND
>   LATITUDE COORDINATES
> 
> ___
> 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
> 
> 

-- 
View this message in context: 
http://old.nabble.com/PROBLEM-CONVERTING-AND-MANIPULATING-LONGITUDE-AND-LATITUDE-COORDINATES-tp27272017p27281826.html
Sent from the PostGIS - User mailing list archive at Nabble.com.

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


[postgis-users] Stephen Salathiel is out of the office on fire duties.

2010-01-22 Thread Stephen . Salathiel

I will be out of the office starting  21/01/2010 and will not return until
03/02/2010.

For all urgent enquiries, please contact Brett Miller on (03) 9637 9008.


Notice:
This email and any attachments may contain information that is personal, 
confidential, legally privileged and/or copyright.No part of it should be 
reproduced, 
adapted or communicated without the prior written consent of the copyright 
owner. 

It is the responsibility of the recipient to check for and remove viruses.
If you have received this email in error, please notify the sender by return 
email, delete 
it from your system and destroy any copies. You are not authorised to use, 
communicate or rely on the information 
contained in this email.

Please consider the environment before printing this email.
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] ST_distance in meters

2010-01-22 Thread Mike Toews

Oscar Zamudio wrote:
I couldn't find 900913 in the Spatial Reference List. Does anybody 
know which one of the following projections uses Google Maps?


900913 isn't any "official" spatial reference number, in fact it is 
leetspeek for "Google" (where 9=g, etc.)


A good reference I found is here:
http://trac.openlayers.org/wiki/SphericalMercator

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


Re: [postgis-users] ST_distance in meters

2010-01-22 Thread Oscar Zamudio
I couldn't find 900913 in the Spatial Reference List. Does anybody know
which one of the following projections uses Google Maps?

   - SR-ORG:6 : Google Projection
   - SR-ORG:16 : Mercator
   (world)
   - SR-ORG:95 : Google
   Maps/Microsoft Live
   - SR-ORG:98 : Mapy.cz
   Projection
   - SR-ORG:6627 : Google
   Mercator
   - SR-ORG:6722 : google map

I'm also trying to compare with Google Maps my results of measuring
distances reprojecting my original SRID-4326 data to some UTM projection but
I don't know what is the right one.
Oscar

On Fri, Jan 22, 2010 at 8:40 AM, ibrahim saricicek <
ibrahimsarici...@gmail.com> wrote:

> Hi;
>
> For 4326 projection distance unit is in degrees.
> Google maps use 900913 projection code, and units for this projection.
> To get the same distance calculation; firstly add 900913 projection to
> spatial_ref_sys table.
>
> Then use this query
>
> select
> st_distance(st_transform(tableA.colB,900913),st_transform(tableA.colC,900913))
>
> Regards..
> IBO..
>
>
> On Fri, Jan 22, 2010 at 1:17 PM, lomarod wrote:
>
>>
>> Hi,
>>
>> I have a table "tableA" with the column "colB" and "colC" type geometry.
>>
>> tableA.colB = geometryfromtext('POINT('-19.922488 -43.9367')',4326)
>> tableA.colC = geometryfromtext('POINT('-19.926357 -43.934223')',4326)
>>
>> Using google earth or google maps, a have the distance between this two
>> points in 500 meters.
>>
>> When I execute the query
>> select st_distance(tableA.colB, tableA.colC) from tableA my result is
>> 0,004593984...
>>
>> How can I have the same result as google maps ?
>> What am I doing wrong ?
>>
>>
>> --
>> View this message in context:
>> http://old.nabble.com/ST_distance-in-meters-tp27271750p27271750.html
>> Sent from the PostGIS - User mailing list archive at Nabble.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 mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] Questions about using SRID not -1

2010-01-22 Thread Oscar Zamudio
Thanks Kevin but Chris Hermansen gave us the key and the problem is solved
now. Please follow my post (my partner Mario actually answered) with the
subjec 'SOLVED!!' and there's all explained.
Regards,

Oscar


On Fri, Jan 22, 2010 at 1:28 PM, Kevin Neufeld wrote:

> On 1/21/2010 12:14 PM, Oscar Zamudio wrote:
>
>> I already tried this way:
>>
>> UPDATE mypoints SET the_geom = ST_SetSRID(the_geom,4326);
>>
>> But I got:
>>
>> ERROR:  Operation on two GEOMETRIES with different SRIDs
>>
>>
> This doesn't make sense to me.  ST_SetSRID simply updates the internal
> structure of the geometry.  It's not doing any operation on the geometry.
>  Do you maybe have a TRIGGER and/or RULE on the table that is invoking
> another function on an UPDATE operation?
>
> Try this:
> CREATE TABLE mypoints2 AS SELECT the_geom FROM mypoints;
> UPDATE mypoints2 SET the_geom = ST_SetSRID(the_geom,4326);
>
> If that doesn't work, can you please isolate an example?  On what geometry
> exactly is the ERROR thrown?
>
> -- Kevin
>
> ___
> 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] PROBLEM CONVERTING AND MANIPULATING LONGITUDE ANDLATITUDE COORDINATES

2010-01-22 Thread Chris Hermansen
Create your circle with radius R metres after you transform your geographic 
point to UTM, perhaps?
Chris Hermansenchris.herman...@timberline.ca
tel+1.604.714.2878 · fax+1.604.733.0631 · mob+1.778.840.4625
Timberline Natural Resource Group · www.timberline.ca
401 · 958 West 8th Avenue  · Vancouver BC · Canada · V5Z 1E5


-Original Message-
From: "Miltos Miltiadous" 
Date: Fri, 22 Jan 2010 13:40:18 
To: 
Subject: [postgis-users] PROBLEM CONVERTING AND MANIPULATING LONGITUDE AND
LATITUDE COORDINATES

___
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] Questions about using SRID not -1

2010-01-22 Thread Kevin Neufeld

On 1/21/2010 12:14 PM, Oscar Zamudio wrote:

I already tried this way:

UPDATE mypoints SET the_geom = ST_SetSRID(the_geom,4326);

But I got:

ERROR:  Operation on two GEOMETRIES with different SRIDs



This doesn't make sense to me.  ST_SetSRID simply updates the internal structure of the geometry.  It's not doing any 
operation on the geometry.  Do you maybe have a TRIGGER and/or RULE on the table that is invoking another function on an 
UPDATE operation?


Try this:
CREATE TABLE mypoints2 AS SELECT the_geom FROM mypoints;
UPDATE mypoints2 SET the_geom = ST_SetSRID(the_geom,4326);

If that doesn't work, can you please isolate an example?  On what geometry 
exactly is the ERROR thrown?

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


[postgis-users] Running PostgreSQL/Postgis in a virtual environment

2010-01-22 Thread Intengu Technologies
We have a virtual environment powered by vmware, is there anyone with
experience in deploying PostgreSQL/Postgis for production purposes.

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


[postgis-users] PROBLEM CONVERTING AND MANIPULATING LONGITUDE AND LATITUDE COORDINATES

2010-01-22 Thread Miltos Miltiadous
Hello.

I am using postgis 1.4 and following the instructions in PostGis in Action,
I am trying to use real coordinates data (longitude and latitude) that they
fetched from a google android mobile device to construct a polygon (I am
going to use it as a query box).

Particularly, the loaded coordinates (lon,lat), which are translated as
mobile user absolute position will be used as the origin of a circle, having
as radius the positional_accuracy of gps device in meters.

So, having the user position in decimal degrees, the radius of a circle in
meters, the yaw angle of mobile device embeded compass showing the direction
of the user, I want to construct a geometry polygon. Namely, I want to find
the four point cartesian coordinates of the polygon. I have some formulas
but I cant get true results. I tried transforming WGS84 (EPSG:4326
coordinate reference system of lon lat coordinates) to a planar reference
system as Tranverse UTM EPSG:32636 but till now I did not find the solution.

It seems that transformations create a huge polygon. Normally, I must create
a polygon that depends from the distance between user position and a real
world object that the user is observing. Increasing that distance the query
window (polygon) must increase too.

If anyone can help I will really appreciate becuase I am stuck here and I
have to deliver very soon my master of thesis LBS & GIS application. 

Thank u very much.

 

Miltos Miltiadous

Postgraduating Information Systmes student in Open University Cyprus.  

 

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


Re: [postgis-users] ST_distance in meters

2010-01-22 Thread ibrahim saricicek
Hi;

For 4326 projection distance unit is in degrees.
Google maps use 900913 projection code, and units for this projection.
To get the same distance calculation; firstly add 900913 projection to
spatial_ref_sys table.

Then use this query

select
st_distance(st_transform(tableA.colB,900913),st_transform(tableA.colC,900913))

Regards..
IBO..

On Fri, Jan 22, 2010 at 1:17 PM, lomarod  wrote:

>
> Hi,
>
> I have a table "tableA" with the column "colB" and "colC" type geometry.
>
> tableA.colB = geometryfromtext('POINT('-19.922488 -43.9367')',4326)
> tableA.colC = geometryfromtext('POINT('-19.926357 -43.934223')',4326)
>
> Using google earth or google maps, a have the distance between this two
> points in 500 meters.
>
> When I execute the query
> select st_distance(tableA.colB, tableA.colC) from tableA my result is
> 0,004593984...
>
> How can I have the same result as google maps ?
> What am I doing wrong ?
>
>
> --
> View this message in context:
> http://old.nabble.com/ST_distance-in-meters-tp27271750p27271750.html
> Sent from the PostGIS - User mailing list archive at Nabble.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


Re: [postgis-users] ST_distance in meters

2010-01-22 Thread Suhr, Ralf
ST_Distance_Spheroid( pointa, pointb, 'SPHEROID["WGS 
84",6378137,298.257223563]');

Gr Ralf

-Ursprüngliche Nachricht-
Von: postgis-users-boun...@postgis.refractions.net 
[mailto:postgis-users-boun...@postgis.refractions.net] Im Auftrag von lomarod
Gesendet: Freitag, 22. Januar 2010 12:18
An: postgis-users@postgis.refractions.net
Betreff: [postgis-users] ST_distance in meters


Hi,

I have a table "tableA" with the column "colB" and "colC" type geometry.

tableA.colB = geometryfromtext('POINT('-19.922488 -43.9367')',4326)
tableA.colC = geometryfromtext('POINT('-19.926357 -43.934223')',4326)

Using google earth or google maps, a have the distance between this two
points in 500 meters.

When I execute the query
select st_distance(tableA.colB, tableA.colC) from tableA my result is
0,004593984...

How can I have the same result as google maps ?
What am I doing wrong ?


-- 
View this message in context: 
http://old.nabble.com/ST_distance-in-meters-tp27271750p27271750.html
Sent from the PostGIS - User mailing list archive at Nabble.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] ST_distance in meters

2010-01-22 Thread lomarod

Hi,

I have a table "tableA" with the column "colB" and "colC" type geometry.

tableA.colB = geometryfromtext('POINT('-19.922488 -43.9367')',4326)
tableA.colC = geometryfromtext('POINT('-19.926357 -43.934223')',4326)

Using google earth or google maps, a have the distance between this two
points in 500 meters.

When I execute the query
select st_distance(tableA.colB, tableA.colC) from tableA my result is
0,004593984...

How can I have the same result as google maps ?
What am I doing wrong ?


-- 
View this message in context: 
http://old.nabble.com/ST_distance-in-meters-tp27271750p27271750.html
Sent from the PostGIS - User mailing list archive at Nabble.com.

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


Re: [postgis-users] Postgis on AIX

2010-01-22 Thread Mark Cave-Ayland

Derek Jones wrote:


Hi folks,

I have a need to run postgis on AIX.


(lots cut)

It doesn't appear to complete the postgis shared lib load. Compared with 
the same o/p from strace on x86, where it goes on to load up geos and 
then proj just fine.


The only time I've seen the postgis-1.4.so module fail to load on 
startup is when PostgreSQL has been unable to locate the GEOS and/or 
PROJ.4 libraries.


The first thing I would check is that if you've installed everything 
with --prefix=/usr1/local to ensure that /usr1/local/lib is in the 
currently library search path. I know you can set this in Linux with 
ldconfig and Solaris with crle, but I'm not sure about the AIX 
equivalent. Once this is set, you should then find that restarting 
PostgreSQL will then allow the PostGIS module to load.


Any thoughts appreciated. Stuck with an 18CPU partition and 64G of 
memory and can't use it until I fix this :-)


Nice :)


HTH,

Mark.

--
Mark Cave-Ayland - Senior Technical Architect
PostgreSQL - PostGIS
Sirius Corporation plc - control through freedom
http://www.siriusit.co.uk
t: +44 870 608 0063

Sirius Labs: http://www.siriusit.co.uk/labs
___
postgis-users mailing list
postgis-users@postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users


Re: [postgis-users] database transformation

2010-01-22 Thread Paolo Rizzi

If your data is already correct, you should _not_ use transform().
You should simply force the new SRID, using something like:
UPDATE  table_name SET geom=ST_SetSRID(geom,3763)

So if you have a coordinate like 1234567.890123, that is correct
for Portugal in the old 27492 SRID and also in the new 3763,
you must not transform it, but simply tell PostGIS that now
the SRID is 3763.
If you use transform, you may end up with a changed coordinate,
that is no longer valid.
If, instead, the coordinate 1234567.890123 is _not_ valid
in the new SRID 3763, you have to transform it, but beware that
the trasformation may not go well.
If the old SRID was deprecated, there may be errors in its definition,
so the transform may give you odd results...

Bye
Paolo



Hello everyone,

I have a database created with the program pgadmin3. The geo-referencing 
system is the Datum 73 / Modified Portuguese


Grid (deprecated) "EPSG: 27492" and wanted to change the whole database 
to the coordinate system ETRS89 / Portugal TM06 (EPSG: 3763).


How can I make this transformation?

Thank you for your help

-- bernardo --




___
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