<?xml version="1.0" encoding="UTF-8"?>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns="http://purl.org/rss/1.0/" xmlns:taxo="http://purl.org/rss/1.0/modules/taxonomy/" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:syn="http://purl.org/rss/1.0/modules/syndication/" xmlns:admin="http://webns.net/mvcb/">
  <channel rdf:about="http://blog.gmane.org/gmane.comp.gis.postgis">
    <title>gmane.comp.gis.postgis</title>
    <link>http://blog.gmane.org/gmane.comp.gis.postgis</link>
    <description/>
    <syn:updatePeriod>hourly</syn:updatePeriod>
    <syn:updateFrequency>1</syn:updateFrequency>
    <syn:updateBase>1901-01-01T00:00+00:00</syn:updateBase>
    <items>
      <rdf:Seq>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31903"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31902"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31891"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31887"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31886"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31877"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31876"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31866"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31858"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31843"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31839"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31830"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31824"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31818"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31813"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31812"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31802"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31801"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31800"/>
        <rdf:li rdf:resource="http://comments.gmane.org/gmane.comp.gis.postgis/31796"/>
      </rdf:Seq>
    </items>
    <image rdf:resource="http://gmane.org/img/gmane-25t.png"/>
    <textinput rdf:resource=""/>
  </channel>
  <image rdf:about="http://gmane.org/img/gmane-25t.png">
    <title>Gmane</title>
    <url>http://gmane.org/img/gmane-25t.png</url>
    <link>http://gmane.org</link>
  </image>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31903">
    <title>loading mostly identical states-based shapefiles</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31903</link>
    <description>&lt;pre&gt;I have US states shape files that are mostly identical, except for those silly legacy columns that have the state code prefix. Is there a way I can import these shapefiles via `shp2pgsql` specifying to *not* load a particular column? As is, the loader croaks because the '-c' switch creates the table on the first state load creating a column starting with that state's code, and then the second state can't find that column.

An additional question -- the '-S' switch promises to "Generate simple geometries instead of MULTI geometries." Why would I want to do that? Please educate me.



--
Puneet Kishor
&lt;/pre&gt;</description>
    <dc:creator>Mr. Puneet Kishor</dc:creator>
    <dc:date>2012-05-26T02:29:32</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31902">
    <title>How to smoothify (rounding the edges)?</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31902</link>
    <description>&lt;pre&gt;Dear list,

I am trying smoothify some linestring.
ST_OffsetCurve looks good, but simply simulates just one side of a
ST_Buffer. 
Chaining of ST_OffsetCurve could be the right way.

SELECT
ST_OffsetCurve(
ST_OffsetCurve(
ST_OffsetCurve(
ST_GeometryFromText(
'LINESTRING(
369172.556 5645414.822
,369083.213632107 5645360.17950829
,369074.985 5645382.472
,368991.798 5645340.312
,368948.814 5645362.793
,368922.301 5645286.628
,368868.668 5645301.318
,368762.112 5645217.804
,368722.516936163 5645220.71672624
,368734.299541879 5645251.65724707
,368696.988917574 5645268.59693863
,368652.612373605 5645243.03860572
,368655.33980816 5645277.83126126
,368622.411 5645234.496)'
,25832)
,-5
,'quad_segs=4 join=round')
,10
,'quad_segs=4 join=round')
,-5
,'quad_segs=4 join=round');

http://postgis.17.n6.nabble.com/file/n4997829/smoothify.png 

How can I get my data edgeless.

Thanks a lot and best regards,
Daniel

--
View this message in context: http://postgis.17.n6.nabble.com/How-to-smoothify-rounding-the-edges-tp4997829.html
Sent from the PostGIS - User mailing list archive at Nabble.com.
&lt;/pre&gt;</description>
    <dc:creator>tee</dc:creator>
    <dc:date>2012-05-25T23:19:02</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31891">
    <title>Topology: Simplifying Coastlines with Islands</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31891</link>
    <description>&lt;pre&gt;Dear list,

I am trying to simplify a dataset of country boundaries by using  
topology functions, as described in Strk's Blog:  
http://strk.keybit.net/blog/2012/04/13/simplifying-a-map-layer-using-postgis-topology/

This has worked quite well, except for coastlines with islands nearby.  
The simplified coastline would probably intersect the island, or pass  
the island at the wrong side, leading to invalid topology. Does that  
mean that countries with coastal islands such as Canada and Norway  
(see attached image) cannot be simplified this way?

Ideally, islands in the way of a to-be-simplified coastline should be  
merged with the mainland. (I have already removed the smaller islands,  
but I think I should not remove the bigger islands too.)

Best regards,
Michiel.

_______________________________________________
postgis-users mailing list
postgis-users&amp;lt; at &amp;gt;postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users
&lt;/pre&gt;</description>
    <dc:creator>Michiel J. van Heek</dc:creator>
    <dc:date>2012-05-25T16:36:39</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31887">
    <title>Nearest neighbor in 3D</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31887</link>
    <description>&lt;pre&gt;Hello All!

This is my first post on the PostGIS mailing list:-).

I would like to use the &amp;lt;-&amp;gt; to solve the nearest neighbor problem, but in
3D space. I tried to do this in 2D and it works perfectly, but apparently
it doesn't cover the 3D case.

Basically I have a table of 3D points (x, y, z) and would like to do:

    ORDER BY point &amp;lt;-&amp;gt; POINT(0 0 0)

but the distance taken into account should be 3D Euclidean distance. Do I
have to use ST_3DDistance() for it or is there a faster way?

Thanks in advance!

Regards,
&lt;/pre&gt;</description>
    <dc:creator>Bartosz Bekier</dc:creator>
    <dc:date>2012-05-25T15:15:07</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31886">
    <title>ST_Aspect artifacts?</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31886</link>
    <description>&lt;pre&gt;Hi all.
The raster resulting from:
CREATE TABLE aspect as select rid, ST_aspect(rast,1,'8BUI') as rast from
hdr;
does not calculate correctly values along the borders of the tiles.
Any way to avoid this?
All the best.

&lt;/pre&gt;</description>
    <dc:creator>Paolo Cavallini</dc:creator>
    <dc:date>2012-05-25T12:56:47</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31877">
    <title>Radcom GIS</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31877</link>
    <description>&lt;pre&gt;MAny thanks to those who can help. I have attached a word based doc to ease the reading. 
 
Hello, my name is Greg and in years past I use to work as a Radio Communication Engineer. 
 
My background is in Electrical and Electronic engineering, not computer science and I am a little stuck. I was wondering if someone can help answer my Comsci questions. The problem is described below. 
 
As an engineer doing Radio Communications the task is to determine coverage and other Radio Communication attributes such as intermods and other things such as harmonic interference. 
 
In the process of doing the job, you start with a transmitter located at an x-y point. 
Then you give that x-y point a transmitting height and power level. 
 
This Tx power level and height, enables distinct points in a closed geographical rectangle or circle, to have its field/strength (F/S) determined at each point.  
 
To determine field strength (F/S), you use different power levels, along with different Tx heights and different computer propagation models to find the F/S. These F/S results are graphically displayed (over maps) 
 
Once the field strength is determined, the second step is called Carrier-To-interference (C/I). If you have two sites close to each other, if they are on the same frequency, the two signals can interfere, degrading either or both of the signals. But sometimes they don’t. eg. Place your new TX site behind a mountain. The mountain shields the two systems from interfering with each other. 
 
Anyway, for C/I calculations, you can pick an x-y point (in a region) and subtract the field strength (from each TX source). Eg. At an x-y point between the two Tx’s, you subtract F/S1 – F/S2. or Vice Versa. 
 
Again computer propagation models are needed to determine F/S, which then allows for the determination of C/I. 
 
If the C/I is bad, you have to re-site the new Tx site location and start again. 
 
It has been a while since doing radio communications but I remember that the maps had DEM’s associated with it. Ie. Different maps for different DEM’s. 
 
I am not the expert when it comes to GIS systems and geographic systems. 
 
From memory you site a TX point, then the program took the DEM (value) across each point (in a region specified by a map) and then for each point, use it as a distance parameter into Radio propagation models, which determines the field strength (F/S) at each point. 
 
I cant remember whether the smaller DEM gives different field strength results compared to a larger DEM being used. 
 
From memory I think you get different field strengths for different DEM Maps. I am not sure if this is correct. What I do know is that a smaller DEM (for a given area) takes longer to compute than using a larger DEM. The smaller DEM might give the same F/S results but takes ages to compute. I am not sure if this is correct. 
 
Ie. The smaller the DEM, the smaller the “distance parameter” between two points, which takes longer to compute all the F/S at each point. 
 
Further, for a smaller DEM, this means over a given area, there is more Field strength points to calculate. But either model yields good results, only if the propagation model is correct. The size of the DEM determines time to compute. 
 
The larger the DEM model, the less points there are in a given area, meaning it takes less time to compute compared to smaller DEM’s, but the results at a large DEM, may give false information. The resolution of parameters is poor.
 
This is the crux of the problem. 
 
I have been reading about the GRASS GIS systems and PostGIS. I believe I can create a better computerised model than that being used at regulators based on GIS concepts.
 
Unfortunately I cant test these systems (on a computer) at the moment. But to save time I am asking some questions. If you can help, I would really appreciate it. 
 
I don’t understand Vectors and Rasters w.r.t GIS systems at the moment, I will learn it later. 
 
What I intend to do is create a new Radio Communications Tool that is web based. 
 
Ie. get all the Maps and layers and display them in a web-page served by a Web based server. Then write models that will take a specified point eg. A transmitter, along with a DEM to figure out the Field Strength and Carrier-to-Interference all done through a web-browser. 
 
I call this application a Web-RADCOM-GIS tool. (WRG - Web based radio communications)
 
This would be good for groups of Radio Communications engineers. 
 
Instead of paying $100,000 for software crippled by dongles and other means, the idea is that the radio communications engineers (or anyone else interested in radio communications, eg HAM radio people). can see F/S results, displayed in a browser window, without any money laid down.
 
The limit of users is limitless. 
 
Instead, the idea is that they can log into a mapserver using a browser, see different layers and do different vector (and raster) operations to determine radio calculations (eg. F/S and C/I). 
 
This new webserver GIS system, could even be extended to planning mobile Base Station sites. The better the modelling tool, the better the mobile coverage and interference modelling, meaning less RF spectrum is wasted. The better the Mobile Phone models get done, the less of mobile base stations need to be deployed, meaning less RF energy bumping around near buildings. In other words, making more out of less. 
 
The less base stations, the less rent paid for building housing the base station antennas. For Example. 
 
Anyway having a web-based GIS system allows for heaps of group-work. Instead of buying the DEM maps and installing them on each PC (in the group), all they have to do is reference different DEM maps through a MAP server, inside a browser. 
 
I have seen this in action. It takes ages to change the Maps on heaps of PC. And similarly if the computer propagation models change, it is effortless to change. 
 
I intent to make this paid for or free for use, this way anyone on the planet can model and investigate different Radio Communication attributes. 
 
Moreover, more everyday people get a chance to see Radio communications. It gives uses who need a frequency allocation (a channel to operate on) a chance to design them using the tool. Then all the FCC has to do is check their results. 
 
Previously in Australia the regulators do all the calculations and frequency allocation. 
 
Anyway going back to the RADCOM tool. I saw a screen shot showing height profiles (in a GRASS window). This is exactly what is needed to determine Field Strength. 
 
I am not sure if I am right, I figure for web page based system, one can use PHP to generate SQL statements, which puts data into either a postgreSQL system or postGIS systems. Then similarly, use PHP and SQL to extract data from the PostGRE or PostGIS and display it on a web-page in a browser. 
 
Even better save the results (as a dataset) and print the F/S or C/I results. Going a step further overlay the results over a census collection district, for example. 
 
The crux of this, leaves me wondering. To make this Radcom tool faster, you should in theory change the way the system computes. Specifically, I am looking at changing either or both, {PostGIS) or (postgre sql) to be re-done in assembly language. 
 
Assembly is orders of magnitude faster. I am unsure whether the performance increase comes from coding PostGIS or PostgreSQL in assembly. 
 
Your advice on whether to code PostGIS or PostGRE in assembly this is keenly sought. 
 
The second thing is that in theory for a web based app, I would have to redo all the menu (and window forms) in PHP? Could you please advise of any books that are good on these topics. 
 
Thankyou for your time and help on these matters. 
 
 
 _______________________________________________
postgis-users mailing list
postgis-users&amp;lt; at &amp;gt;postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users
&lt;/pre&gt;</description>
    <dc:creator>Greg Evans</dc:creator>
    <dc:date>2012-05-25T08:59:21</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31876">
    <title>Difference between ST_Union and ST_Collect</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31876</link>
    <description>&lt;pre&gt;Folks,

after seeing some polygons disappear during the generalization process, I found 
out that ST_Union and ST_Collect are different beasts when dealing with 
intersecting lines.

I did generalization by: breaking original polygons into lines, generalizing them 
and them assembling back lines into (multi) polygons.

At first I used ST_Collect for the last step:
ST_Multi(ST_BuildArea(ST_Collect(r2.wkb_geometry)))

...but some of the lines intersected, especially over intricate coastlines, 
causing big polygons to be silently dropped.

But when I switched to ST_union:
ST_Multi(ST_BuildArea(ST_Union(r2.wkb_geometry)))

all polygons where built (except for the ones that had to be dropped because their 
area was 0, usually islets).

Is this by design (I wasn't able to find documentation on the different behaviour 
of these two functions) ?

Regards,

P.S.
Sandro, I did not forget to help speeding up the topology construction routines, 
it is just that next week we have a deadline, hence I gave a lower priority to the 
task of setting up a machine with the configuration you requested... stay tuned.

Luca Morandini
Data Architect - AURIN project
Department of Computing and Information Systems
University of Melbourne
&lt;/pre&gt;</description>
    <dc:creator>Luca Morandini</dc:creator>
    <dc:date>2012-05-25T03:52:45</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31866">
    <title>Error browsing database</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31866</link>
    <description>&lt;pre&gt;Hi all.

A simple

CREATE TABLE dtm 
 AS SELECT ST_Transform(rast,4326, 'Bilinear') AS rast
  FROM hdr;

results in

gdalinfo "PG:host=localhost dbname=postgis2 user=corso table=dtm mode=2
password=***"

ERROR 1: Error browsing database for PostGIS Raster properties

Similar tables are read without issues.

Any hint?
All the best.

&lt;/pre&gt;</description>
    <dc:creator>Paolo Cavallini</dc:creator>
    <dc:date>2012-05-24T16:32:04</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31858">
    <title>rt_resample test fails on PostGIS 2.0,Debian 6 amd64</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31858</link>
    <description>&lt;pre&gt;The subject line says everything. I try to compile PostGIS 2.0 on an amd64
Debian 6 (squeeze). ./configure works ok, and so does make, but here's part of
the output of "make test":

rt_resample .. failed (diff expected obtained: /tmp/pgis_reg/test_51_diff)

I attach the test_51_diff file.

The problem occurs both on 2.0.0, and on the latest 2.0-branch svn.

Am I doing something wrong? Could someone help me get past this problem? Thanks!
--- rt_resample_expected2012-05-24 09:21:33.000000000 +0000
+++ /tmp/pgis_reg/test_51_out2012-05-24 10:34:40.562455574 +0000
&amp;lt; at &amp;gt;&amp;lt; at &amp;gt; -11,35 +11,35 &amp;lt; at &amp;gt;&amp;lt; at &amp;gt;
 1.14|992163|11|11|1|1000.000|-1000.000|0.000|0.000|-500100.000|600950.000|t|t|t
 1.15|992163|201|201|1|50.000|50.000|0.000|0.000|-500040.000|589957.000|t|t|t
 1.16|992163|84|84|1|121.000|121.000|0.000|0.000|-500093.000|589875.000|t|t|t
-1.17|993310|243|243|1|50.000|50.000|0.000|0.000|950710.000|1397157.000|t|t|t
-1.18|993309|243|243|1|50.000|50.000|0.000|0.000|950760.000|1396957.000|t|t|t
+1.17|993310|248|235|1|50.000|50.000|0.000|0.000|950710.000|1397557.000|t|t|t
+1.18|993309|248|235|1|50.000|50.000|0.000|0.000|950760.000|1397357.000|t|t|t
 1.19|992163|10|10|1|1000.000|-1000.000|3.000|3.000|-500030.000|600000.000|t|t|t
 1.2|993310|12|12|1|1009.894|-1009.894|0.000|0.000|950732.188|1409281.783|t|t|t
 1.20|993310|12|12|1|1009.894|-1009.894|3.000|3.000|950691.792|1409281.783|t|t|t
 1.21|993309|12|12|1|1009.916|-1009.916|1.000|3.000|950742.107|1409088.896|t|t|t
-1.22|993310|24|24|1|500.000|500.000|3.000|3.000|950657.188|1397356.783|t|t|t
-1.23|993310|26|26|1|500.000|500.000|0.000|6.000|950452.000|1396632.000|t|t|t
+1.22|993310|25|24|1|500.000|500.000|3.000|3.000|950657.188|1397356.783|t|t|t
+1.23|993310|26|25|1|500.000|500.000|0.000|6.000|950452.000|1397132.000|t|t|t
 1.24|992163|207|101|1|50.000|-100.000|3.000|0.000|-500319.000|600056.000|t|t|t
 1.25|992163|207|101|1|50.000|-100.000|3.000|0.000|-500319.000|600056.000|t|t|t
 1.26|992163|150|150|1|66.667|-66.667|0.000|0.000|-500000.000|600000.000|t|t|t
-1.27|993310|150|150|1|80.792|-80.792|0.000|0.000|950732.188|1409281.783|t|t|t
+1.27|993310|150|150|1|82.219|-78.090|0.000|0.000|950732.188|1409281.783|t|t|t
 1.28|992163|5|5|1|2064.200|-2291.200|0.000|0.000|-500321.000|601456.000|t|t|t
 1.29|||||||||||||
 1.3|993309|12|12|1|1009.916|-1009.916|0.000|0.000|950762.305|1409088.896|t|t|t
 1.30|984269|12|8|1|0.012|-0.012|0.000|0.000|-107.029|50.206|t|t|t
 1.31|974269|12|8|1|0.012|-0.012|0.000|0.000|-107.029|50.206|t|t|t
 1.4|994269|12|8|1|0.012|-0.012|0.000|0.000|-107.029|50.206|t|t|t
-1.5|993310|24|24|1|500.000|500.000|0.000|0.000|950732.188|1397281.783|t|t|t
+1.5|993310|25|23|1|500.000|500.000|0.000|0.000|950732.188|1397781.783|t|t|t
 1.6|992163|10|10|1|1000.000|-1000.000|0.000|0.000|-500000.000|600000.000|t|t|t
 1.7|992163|10|10|1|1000.000|-1000.000|0.000|0.000|-500000.000|600000.000|t|t|t
 1.8|992163|20|20|1|500.000|500.000|0.000|0.000|-500000.000|590000.000|t|t|t
 1.9|992163|40|40|1|250.000|250.000|0.000|0.000|-500000.000|590000.000|t|t|t
 2.1|993310|12|12|1|1009.894|-1009.894|0.000|0.000|950732.188|1409281.783|t|t|t
-2.10|993310|24|24|1|500.000|500.000|0.000|0.000|950732.188|1397281.783|t|t|t
-2.11|993309|121|121|1|100.000|100.000|0.000|0.000|950762.305|1396988.896|t|t|t
+2.10|993310|25|23|1|500.000|500.000|0.000|0.000|950732.188|1397781.783|t|t|t
+2.11|993309|123|117|1|100.000|100.000|0.000|0.000|950762.305|1397388.896|t|t|t
 2.12|993310|6|6|1|2000.000|2000.000|0.000|0.000|950732.188|1397281.783|t|t|t
 2.13|993310|8|8|1|1500.000|1500.000|0.000|0.000|950732.188|1397281.783|t|t|t
-2.14|993310|24|24|1|500.000|500.000|0.000|0.000|950732.188|1397281.783|t|t|t
+2.14|993310|25|23|1|500.000|500.000|0.000|0.000|950732.188|1397781.783|t|t|t
 2.15|993310|16|16|1|750.000|750.000|0.000|0.000|950732.188|1397281.783|t|t|t
 2.2|993309|12|12|1|1009.916|-1009.916|0.000|0.000|950762.305|1409088.896|t|t|t
 2.3|994269|12|8|1|0.012|-0.012|0.000|0.000|-107.029|50.206|t|t|t
_______________________________________________
postgis-users mailing list
postgis-users&amp;lt; at &amp;gt;postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users
&lt;/pre&gt;</description>
    <dc:creator>Antonis Christofides</dc:creator>
    <dc:date>2012-05-24T11:44:46</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31843">
    <title>AddToPROJ4SRSCache: couldn't parse proj4 string -major axis or radius = 0 or not given</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31843</link>
    <description>&lt;pre&gt;Hi there,

 

I need to add a projection to support an imported shapefiles dataset with
the following projection:

PROJCS [              "Business Analyst Projection",

                                GEOGCS[             

                                                "GCS_WGS_1984",

                                                DATUM[
"D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],

                                                PRIMEM["Greenwich",0.0],

 
UNIT["Degree",0.0174532925199433]],  

                                PROJECTION["Equidistant_Conic"],

                                PARAMETER["False_Easting",0.0],

                                PARAMETER["False_Northing",0.0],

 
PARAMETER["Central_Meridian",-88.02657853618898],

 
PARAMETER["Standard_Parallel_1",41.82305904132364],

 
PARAMETER["Standard_Parallel_2",41.82305804132364],

                                PARAMETER["Latitude_Of_Origin",0.0],

                                UNIT["Meter",1.0]          ]

 

I used the following SQL to add the projection to the spatial_ref_sys
table:

 

INSERT INTO spatial_ref_sys (srid, auth_name, auth_srid, srtext,
proj4text) VALUES (9998, 'EPSG', 9998, 

                'PROJCS["Business Analyst Projection",

                                GEOGCS["GCS_WGS_1984",

 
DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],

                                                PRIMEM["Greenwich",0.0],

 
UNIT["Degree",0.0174532925199433]],

                                PROJECTION["Equidistant_Conic"],

                                PARAMETER["False_Easting",0.0],

                                PARAMETER["False_Northing",0.0],

 
PARAMETER["Central_Meridian",-88.02657853618898],

 
PARAMETER["Standard_Parallel_1",41.82305904132364],

 
PARAMETER["Standard_Parallel_2",41.82305804132364],

                                PARAMETER["Latitude_Of_Origin",0.0],

                                UNIT["Meter",1.0]]', 

 

                '+proj=eqdc +lat_0=41.82305904132364 

                                +lon_0=-88.02657853618898

                                +lon_1=-88.02657853618898 

                                +lon_2=-88.02657853618898 

                                +x_0=0 

                                +y_0=0 

                                +a=6378137.0

                                +ellps=WGS84 

                                +datum=WGS84 

                                +units=m 

                                +no_defs');

 

When I tried to do a transform test [ update major set
geom=transform(geom, 4326) ], I go the following error:

ERROR:  AddToPROJ4SRSCache: couldn't parse proj4 string: '+proj=eqdc


                                +lat_0=41.82305904132364 

                                +lon_0=-88.02657853618898

                                +lon_1=-88.02657853618898 

                                +lon_2=-88.02657853618898 

                                +x_0=0 

                                +y_0=0 

                                +a=6378137.0

                                +ellps=WGS84 

                                +datum=WGS84 

                                +units=m 

                                +no_defs': major axis or radius = 0 or not
given

 

BTW, I found no documentation that provides a complete listing of
parameters I need to set for eqdc projection.



Thanks for any hint.

DWW


____
This communication may contain information that is confidential, privileged or subject to copyright. If you are not the intended recipient, please advise by return e-mail and delete the message and any attachments immediately without reading, copying or forwarding to others.

_______________________________________________
postgis-users mailing list
postgis-users&amp;lt; at &amp;gt;postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users
&lt;/pre&gt;</description>
    <dc:creator>David Wu</dc:creator>
    <dc:date>2012-05-23T17:00:08</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31839">
    <title>Ubuntu and Postgresql9.1.3-2</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31839</link>
    <description>&lt;pre&gt;Someone know how can linstall postgresql-9.1.3-2-linux.run in Ubuntu 12.04?

I installed it but I can not to launch stackbuilder to configure new tools
like postgis.

I need help, please!!!


Thanks
_______________________________________________
postgis-users mailing list
postgis-users&amp;lt; at &amp;gt;postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users
&lt;/pre&gt;</description>
    <dc:creator>José María Amuedo</dc:creator>
    <dc:date>2012-05-23T15:08:54</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31830">
    <title>optimize query</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31830</link>
    <description>&lt;pre&gt;I'm finding a search query we wrote is running really slow.  Any ideas on how to optimize this better?  Or is this the best I can do?

SELECT ST_Value(rast, ST_Transform(ST_SetSRID(ST_Point(66.12, 37.32), 4326),ST_SRID(rast))) FROM dted_elevations WHERE ST_Intersects(rast, ST_Transform(ST_SetSRID(ST_Point(66.12, 37.32), 4326), ST_SRID(rast)));

--Jack Gold


_______________________________________________
postgis-users mailing list
postgis-users&amp;lt; at &amp;gt;postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users
&lt;/pre&gt;</description>
    <dc:creator>Gold, Jack L  (US SSA</dc:creator>
    <dc:date>2012-05-22T21:16:40</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31824">
    <title>Does -a option create when table does not exist?</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31824</link>
    <description>&lt;pre&gt;Hi,

Does the -a option of shp2pgsql and raster2pgsql creates the table when the table does not exist? If not would not that be useful for batch?

Pierre
&lt;/pre&gt;</description>
    <dc:creator>Pierre Racine</dc:creator>
    <dc:date>2012-05-22T19:02:25</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31818">
    <title>calling tiger geocoder functions from rails?</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31818</link>
    <description>&lt;pre&gt;I'm looking to include my postgis db as an additional db in a rails app and
pass data from the app to the tiger.geocode(), tiger.reverse_geocode(), and
probably tiger.normalize_address() functions.

has anyone done this who could maybe point me in the right direction?

--
View this message in context: http://postgis.17.n6.nabble.com/calling-tiger-geocoder-functions-from-rails-tp4997646.html
Sent from the PostGIS - User mailing list archive at Nabble.com.
&lt;/pre&gt;</description>
    <dc:creator>Jeff Rossi</dc:creator>
    <dc:date>2012-05-22T14:42:01</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31813">
    <title>How to find length of a part of line string...</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31813</link>
    <description>&lt;pre&gt;Hi,

     I want to find the length of the particular part(from starting point
of linestring to a point i given) in line string. Is it possible to do this
in Postgis.If possible how to write query for this.

Also I attached the explanation image with this mail....
&lt;/pre&gt;</description>
    <dc:creator>ram prasad</dc:creator>
    <dc:date>2012-05-22T12:28:05</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31812">
    <title>tiger_geocoder import error</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31812</link>
    <description>&lt;pre&gt;Sorry if this is a double post, had an issue with my mail server and had to
use my personal email.

This is my first project with postgresql/postgis and it is possible and
maybe even probable that I made an error.  I used pgadmin to generate import
scripts by state and started with florida. The errors are below:

####

COMMIT
NOTICE:  INSERT INTO
tiger_data.fl_tabblock(the_geom,intptlon,intptlat,awater,aland,funcstat,uace,ur,mtfcc,name,tabblock_id,blockce,tractce,countyfp,statefp)
SELECT
the_geom,intptlon10,intptlat10,awater10,aland10,funcstat10,uatyp10,uace10,ur10,mtfcc10,name10,tabblock_id,blockce10,tractce10,countyfp10,statefp10
FROM tiger_staging.fl_tabblock10;
CONTEXT:  SQL function "loader_load_staged_data" statement 1
ERROR:  INSERT has more expressions than target columns
LINE 1: ...name10,tabblock_id,blockce10,tractce10,countyfp10,statefp10 ...
                                                            ^
QUERY:  INSERT INTO
tiger_data.fl_tabblock(the_geom,intptlon,intptlat,awater,aland,funcstat,uace,ur,mtfcc,name,tabblock_id,blockce,tractce,countyfp,statefp)
SELECT
the_geom,intptlon10,intptlat10,awater10,aland10,funcstat10,uatyp10,uace10,ur10,mtfcc10,name10,tabblock_id,blockce10,tractce10,countyfp10,statefp10
FROM tiger_staging.fl_tabblock10;
CONTEXT:  PL/pgSQL function "loader_load_staged_data" line 24 at EXECUTE
statement
SQL function "loader_load_staged_data" statement 1

####

It appears to me that the extra column is uatyp10

Many Thanks for your Assistance,

Jeff

--
View this message in context: http://postgis.17.n6.nabble.com/tiger-geocoder-import-error-tp4997634.html
Sent from the PostGIS - User mailing list archive at Nabble.com.
&lt;/pre&gt;</description>
    <dc:creator>Jeff Rossi</dc:creator>
    <dc:date>2012-05-22T11:45:46</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31802">
    <title>PostGIS 2.0 on CentOS 5</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31802</link>
    <description>&lt;pre&gt;Hi,

Anybody knows where CentOS 5 binaries are available (if they exists)?

Thanks,

Pierre
&lt;/pre&gt;</description>
    <dc:creator>Pierre Racine</dc:creator>
    <dc:date>2012-05-21T16:05:34</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31801">
    <title>Installing on Debian</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31801</link>
    <description>&lt;pre&gt;Hi all.
The Debian package for 2.0 is still not done (any volunteer here?). Is there a clean
install howto for Debian unstable/testing? I only found outdated/incomplete instructions.
All the best, and thanks.
&lt;/pre&gt;</description>
    <dc:creator>Paolo Cavallini</dc:creator>
    <dc:date>2012-05-21T11:04:46</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31800">
    <title>GIN index?</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31800</link>
    <description>&lt;pre&gt;Wasn't PostGIS 2.0 supposed to support GIN indexes?

Seems like there was talk of that many moons ago, but I don't see anything
about GIN in the 2.0 documentation.

Aren
_______________________________________________
postgis-users mailing list
postgis-users&amp;lt; at &amp;gt;postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users
&lt;/pre&gt;</description>
    <dc:creator>Aren Cambre</dc:creator>
    <dc:date>2012-05-19T19:49:47</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31796">
    <title>Point calculation based on stationary distance</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31796</link>
    <description>&lt;pre&gt;Hello,

 

A question: I have a linestring (polyline) for which i would like to
calculate a point which is located on the polyline but a specific distance
from the end or beggining of the polyline. Once this point is calculated i
would then like to calculate 2 points on both sides of the linestring, which
are perpendicular to the line like this:

 

Description: cid:image001.jpg&amp;lt; at &amp;gt;01CD3522.C7638240

 

Any idea how to do this? I didnt notice any function within postgis that
deal with stationary distances which is what im looking for probably. Thanks
for any advice!

 

_______________________________________________
postgis-users mailing list
postgis-users&amp;lt; at &amp;gt;postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users
&lt;/pre&gt;</description>
    <dc:creator>Uroš Bernik</dc:creator>
    <dc:date>2012-05-19T07:34:07</dc:date>
  </item>
  <item rdf:about="http://comments.gmane.org/gmane.comp.gis.postgis/31784">
    <title>PostGis Raster and GDAL 1.9</title>
    <link>http://comments.gmane.org/gmane.comp.gis.postgis/31784</link>
    <description>&lt;pre&gt;
Hi,
I stored DEM raster data in PostGIS 2.0 and I used gdal_translate 
(with mode=2) to extract a mosaic of tiles. When I’m using GDAL 1.8 is 
working but when I’m using GDAL 1.9 I can see the values of pixel for 
just one tile. For example I uploaded in the database a tile with an 
attribute filename=056n05_w and an other tile with an filename=056n05_e 
and I executed the command with GDAL 1.8 and GDAL 1.9:
 
gdal_translate" -a_nodata "-32767" -of GTIFF  "PG:host=myhost 
port=14070 dbname='mybd' user='myuser' password='mypassword' 
schema='myschema'table='Elevation' mode=2 where=\"filename LIKE 
\'%056n05%\'\" " 056N05.tif
 
With GDAL 1.8 I obtained a mosaic of 2 tiles (the result expected) 
but with GDAL 1.9, I obtained a raster with the same extent but only the
 values of the tile 056n05_w are present. The pixels of the tile 
056n05_w have a value 0.
 
I tried the stable version 1.9 MSVC2008 (Win32) and development version available today  at  http://www.gisinternals.com/sdk/
 and I have still the same problem. For your information, I'm running 
GDAL on windows7 and PostGIS 2.0 is installed on a Linux server. I tried
 many version of GDAL 1.9 and I always have the same problem. The 
version 1.8 is working very well but I would like to take advantage of 
the improvements in the last version because I'm working on a project 
where processing time is very important.
 
Does anybody have any idea about this problem ?
 
Thanks

David Bélanger
       _______________________________________________
postgis-users mailing list
postgis-users&amp;lt; at &amp;gt;postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users
&lt;/pre&gt;</description>
    <dc:creator>David B�langer</dc:creator>
    <dc:date>2012-05-18T00:13:06</dc:date>
  </item>
  <textinput rdf:about="http://search.gmane.org/?group=$group=gmane.comp.gis.postgis">
    <title>Search Engine</title>
    <description>Search the mailing list at Gmane</description>
    <name>query</name>
    <link>http://search.gmane.org/?group=$group=gmane.comp.gis.postgis</link>
  </textinput>
</rdf:RDF>

