headers
Michael Fuhr | 1 May 2006 01:24
Favicon

Re: Transformation problem using transform()

On Mon, May 01, 2006 at 10:12:47AM +1200, Marc Angelo wrote:
> The following is a random sample of the different srid's I tried to
> transform to:
> 2014,27200,2029,2848,26705,27231
> It seems that apart from trying to transform to 4269 all other
> transformations fail with the error:
> ERROR:  transform: couldn't project point: -45 (geocentric transformation
> missing z or ellps)
> 
> The following is the one and only current record in the spatial table:
> 1;"E004010001CC1B0F";"2006-04-30";"-45.879185";"170.50196";"0101000020E61000
> 001F11532289F046C0CEFC6A0E10506540"

Looks like latitude and longitude are in the wrong order:

postgis=> SELECT AsEWKT('0101000020E61000001F11532289F046C0CEFC6A0E10506540');
                asewkt                 
---------------------------------------
 SRID=4326;POINT(-45.879185 170.50196)
(1 row)

Geometries are (X,Y) so lat/lon coordinates should be (lon,lat).

--

-- 
Michael Fuhr
Permalink | Reply | 5 comments |
headers
Shoaib Burq | 1 May 2006 05:48
Picon

PostGIS in Research/Education

Hi all

I have volunteered to give a presentation on the use of OSGeo Software
for Research/Education. The presentation will be at Melbourne Univ,
Australia. All for the cause :)

I was wondering if there are any people out there who might like to
share details of projects or experiences with using PostGIS or related
technologies for Research/Education purposes. It would give me
something solid to crow about.

Thanks heaps
shoaib
--
http://geospatial.nomad-labs.com
Permalink | Reply | 0 comments |
headers
Bruce Rindahl | 1 May 2006 16:28

RE: Use of PostGIS for Renewable Energy resourceanalysis, mapping, etc

Dana
I am working with an agency here in Denver that displays floodplain
information using PostGIS to a web display.  The display uses WMS
(Mapserver) for raster info and SVG for vector data.  We have the ability
for the user to type in their address and the interface will zoom to that
location (geocoding from Yahoo!).  We are also using a similar interface to
display real-time rainfall.  I am also working on an addition to this
interface to print maps to scale.   If any of this sounds interesting let me
know.
The floodplain interface is:
http://www.lrcwe-data.com/UDFCD.svg
A rainfall interface is at:
http://www.lrcwe-data.com/CoCoRaHS.svg
The interface and coding was based on the work done at:
http://www.carto.net
Bruce

-----Original Message-----
From: postgis-users-bounces <at> postgis.refractions.net
[mailto:postgis-users-bounces <at> postgis.refractions.net] On Behalf Of dnrg
Sent: Sunday, April 30, 2006 10:06 AM
To: postgis postgis
Subject: [postgis-users] Use of PostGIS for Renewable Energy
resourceanalysis, mapping, etc

Anyone here use PostGIS (and/or MapServer) for
Renewable Energy analysis, map generation, etc?

I will be doing some pro-bono work for my alma mater
this summer, serving up wind resource data (overlaid
(Continue reading)

Permalink | Reply | 0 comments |
headers
Robert Burgholzer | 1 May 2006 18:53

RE: Use of PostGIS for Renewable Energy resourceanalysis, mapping, etc

Dana,
I have not used this, but have done some rudimentary solar with GIS
analysis, and would be quite interested in your project. Keep on
posting, or feel free to contact me offline.

Robert - rburghol <at> vt.edu

-----Original Message-----
From: postgis-users-bounces <at> postgis.refractions.net
[mailto:postgis-users-bounces <at> postgis.refractions.net] On Behalf Of dnrg
Sent: Sunday, April 30, 2006 12:06 PM
To: postgis postgis
Subject: [postgis-users] Use of PostGIS for Renewable Energy
resourceanalysis, mapping, etc

Anyone here use PostGIS (and/or MapServer) for
Renewable Energy analysis, map generation, etc?

I will be doing some pro-bono work for my alma mater
this summer, serving up wind resource data (overlaid
upon parcel data), to citizens of my state. Presently,
the University must field phone calls, fire up ArcMap,
search/pan to a citizen's parcel, muck with map
layout, then mail a paper map. Time consuming and a
waste of resources.

I realize a lot of this could be automated with
ArcMap, but not to the extent that satifies me. I want
to architect a self-service web app, extending it to
microhydro and solar insolation later on.
(Continue reading)

Permalink | Reply | 0 comments |
headers
Randy How | 1 May 2006 21:54

WKT - Exponential Coordinates


I'm executing an intersection in PostGIS 1.1.0 using WKT input multipolygon
geometries.  The WKT result coordinates have exponential values (i.e.
28.110482 -1e-06,28.9210135294118 -1e-06,29.7315450588235 -1e-06).  One side
of the input multipolygons is the equator.

Here's the code -

private com.vividsolutions.jts.geom.Geometry _jtsProductCoverage = null;

string sqlIntersection = "select AsText(Intersection(GeometryFromText('" +
wktAOI + 	"',4326),geometry)) " +
	"from coverage " +
	"where product_id = " + this.Product_id + ";";

// Convert the Product Coverage to JTS geometry string wktProductCoverage =
(new DataAccess()).ExecuteScalar(sqlIntersection);

_jtsProductCoverage = JTSHelper.WKTGeometryToJTSGeometry(wktProductCoverage,
4326);

Presently the results are inputs to the JTS.WKTReader for succeeding
procedures, but the exponential values are throwing an exception in the
WKTReader.

Ideally I'd like to eliminate the WKTReader results conversion.  FYI - I'm
writing this in C#.  I'm wondering what the best option is to transfer the
geometries from PostGIS binary to JTS/GeoToolsNet geometries.

If anyone has any input/feedback/examples, I'd greatly appreciate it.
(Continue reading)

Permalink | Reply | 1 comment |
headers
Randy How | 1 May 2006 22:49

RE: WKT - Exponential Coordinates


Well I think I've basically figured this out.  Unless someone
hollers...."NOOO!!!"

Geotools.IO.GeometryWKBReader wkbReader = new Geotools.IO.GeometryWKBReader(
	new Geotools.Geometries.GeometryFactory(
		new Geotools.Geometries.PrecisionModel(), 
		4326));

Geotools.Geometries.Geometry geom =
wkbReader.Create((byte[])cmd.ExecuteScalar());

Using GeoTools GeometryWKBReader and PostGIS AsBinary() creates the GeoTools
geometry.

Randy

-----Original Message-----
From: postgis-users-bounces <at> postgis.refractions.net
[mailto:postgis-users-bounces <at> postgis.refractions.net] On Behalf Of Randy
How
Sent: Monday, May 01, 2006 1:55 PM
To: 'PostGIS Users Discussion'
Subject: [postgis-users] WKT - Exponential Coordinates

I'm executing an intersection in PostGIS 1.1.0 using WKT input multipolygon
geometries.  The WKT result coordinates have exponential values (i.e.
28.110482 -1e-06,28.9210135294118 -1e-06,29.7315450588235 -1e-06).  One side
of the input multipolygons is the equator.
(Continue reading)

Permalink | Reply | 0 comments |
headers
Martin Davis | 1 May 2006 23:18

RE: WKT - Exponential Coordinates

JTS 1.7 added support for scientific notation in WKTWriter.  My guess is
that this hasn't made its way into PostGIS.

Martin Davis, Senior Technical Architect
Vivid Solutions Inc.      www.vividsolutions.com
Suite #1A-2328 Government Street Victoria, B.C. V8T 5G5
Phone: (250) 385 6040 - Local 308 Fax: (250) 385 6046

> -----Original Message-----
> From: postgis-users-bounces <at> postgis.refractions.net 
> [mailto:postgis-users-bounces <at> postgis.refractions.net] On 
> Behalf Of Randy How
> Sent: May 1, 2006 12:55 PM
> To: 'PostGIS Users Discussion'
> Subject: [postgis-users] WKT - Exponential Coordinates
> 
> 
>  
> 
> I'm executing an intersection in PostGIS 1.1.0 using WKT 
> input multipolygon geometries.  The WKT result coordinates 
> have exponential values (i.e. 28.110482 
> -1e-06,28.9210135294118 -1e-06,29.7315450588235 -1e-06).  One 
> side of the input multipolygons is the equator.
> 
> Here's the code -
> 
> private com.vividsolutions.jts.geom.Geometry 
> _jtsProductCoverage = null;
>
(Continue reading)

Permalink | Reply | 501 comments |
headers
Antti.Roppola | 2 May 2006 01:52
Picon

[UNCLASSIFIED]RE: Use of PostGIS for Renewable En ergy resourceanalysis, mapping, etc

Hi Dana,

We have a bioenergy resources atlas that currently looks
at agricultural residues. We're thinking about extending it
to other renewable energy sources and putting it in PostGIS.
http://adl.brs.gov.au/mapserv/biomass/index.phtml
If you're interested, I'll dig up the project report for the
old site as it included a Decision Support System.

It'll be intersting to see where ESRI take ArcGIS. Have a look
at their new ArcGIS Server product, it makes the entire ArcGIS
object set available as Java or .net. It's not a huge stretch
to add Python bindings to that. I'm hoping that this'll prompt
the development of an ecosystem of analystical GIS objects.

I spoke to some of ESRI's developers, and they are expecting to
extend the Python in ArcGIS, but it'll take a while. I'd love to
be able to do cartography templates in it.

We've done some analysis in PostGIS, but ArcGIS (especially the
model building environment) is much more complete and shiny.

Virtual Linux hosting (Xen instances) are pretty easy to come by
and provide an ultimate level of control over the hosted environemnt.

Hope your project works out well. :o)

Cheers,

Antti
(Continue reading)

Permalink | Reply | 0 comments |
headers
Paul Ramsey | 2 May 2006 02:25
Favicon

Summer of Code : Students

Students interested in the summer of code? Your application process  
starts here...

<http://code.google.com/soc/student_step1.html>
Permalink | Reply | 0 comments |
headers
Daniel | 2 May 2006 11:46
Favicon

Collect problem

Hello,

I've used to proceed queries like:
SELECT collect(lga_geom) as buffer FROM lg_annee
WHERE lga_annee='2005' and lga_nat_cdn=2 and lga_var>-1 GROUP BY lga_annee;

With last versions of Postgis, such queries failed, even if I add tests 
like:

SELECT collect(lga_geom) as buffer FROM lg_annee
WHERE lga_annee='2005' and lga_nat_cdn=2 and lga_var>-1 and lga_geom is 
not null
and isvalid(lga_geom) and dimension(lga_geom)=2 GROUP BY lga_annee;

Errors reported are:

NOTICE:  Self-intersection
[....]
NOTICE:  Self-intersection
NOTICE:  Self-intersection

ERROR:  lwcollection_construct: mixed dimension geometries

Table lg_annee contain 2-dimension polygon geometries. The query worked 
fine with older versions.

What can I do ?

Any idea ?
(Continue reading)

Permalink | Reply | 0 comments |

Gmane