http://blog.gmane.org/gmane.comp.gis.postgis hourly 1 1901-01-01T00:00+00:00 http://gmane.org/img/gmane-25t.png http://gmane.org http://comments.gmane.org/gmane.comp.gis.postgis/31903 <pre>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 </pre> Mr. Puneet Kishor 2012-05-26T02:29:32 http://comments.gmane.org/gmane.comp.gis.postgis/31902 <pre>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,</pre> tee 2012-05-25T23:19:02 http://comments.gmane.org/gmane.comp.gis.postgis/31891 <pre>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&lt; at &gt;postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users </pre> Michiel J. van Heek 2012-05-25T16:36:39 http://comments.gmane.org/gmane.comp.gis.postgis/31887 <pre>Hello All! This is my first post on the PostGIS mailing list:-). I would like to use the &lt;-&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 &lt;-&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, </pre> Bartosz Bekier 2012-05-25T15:15:07 http://comments.gmane.org/gmane.comp.gis.postgis/31886 <pre>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. </pre> Paolo Cavallini 2012-05-25T12:56:47 http://comments.gmane.org/gmane.comp.gis.postgis/31877 <pre>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 comp</pre> Greg Evans 2012-05-25T08:59:21 http://comments.gmane.org/gmane.comp.gis.postgis/31876 <pre>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</pre> Luca Morandini 2012-05-25T03:52:45 http://comments.gmane.org/gmane.comp.gis.postgis/31866 <pre>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. </pre> Paolo Cavallini 2012-05-24T16:32:04 http://comments.gmane.org/gmane.comp.gis.postgis/31858 <pre>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 &lt; at &gt;&lt; at &gt; -11,35 +11,35 &lt; at &gt;&lt; at &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</pre> Antonis Christofides 2012-05-24T11:44:46 http://comments.gmane.org/gmane.comp.gis.postgis/31843 <pre>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 th</pre> David Wu 2012-05-23T17:00:08 http://comments.gmane.org/gmane.comp.gis.postgis/31839 <pre>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&lt; at &gt;postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users </pre> José María Amuedo 2012-05-23T15:08:54 http://comments.gmane.org/gmane.comp.gis.postgis/31830 <pre>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&lt; at &gt;postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users </pre> Gold, Jack L (US SSA 2012-05-22T21:16:40 http://comments.gmane.org/gmane.comp.gis.postgis/31824 <pre>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 </pre> Pierre Racine 2012-05-22T19:02:25 http://comments.gmane.org/gmane.comp.gis.postgis/31818 <pre>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. </pre> Jeff Rossi 2012-05-22T14:42:01 http://comments.gmane.org/gmane.comp.gis.postgis/31813 <pre>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.... </pre> ram prasad 2012-05-22T12:28:05 http://comments.gmane.org/gmane.comp.gis.postgis/31812 <pre>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,uac</pre> Jeff Rossi 2012-05-22T11:45:46 http://comments.gmane.org/gmane.comp.gis.postgis/31802 <pre>Hi, Anybody knows where CentOS 5 binaries are available (if they exists)? Thanks, Pierre </pre> Pierre Racine 2012-05-21T16:05:34 http://comments.gmane.org/gmane.comp.gis.postgis/31801 <pre>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. </pre> Paolo Cavallini 2012-05-21T11:04:46 http://comments.gmane.org/gmane.comp.gis.postgis/31800 <pre>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&lt; at &gt;postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users </pre> Aren Cambre 2012-05-19T19:49:47 http://comments.gmane.org/gmane.comp.gis.postgis/31796 <pre>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&lt; at &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&lt; at &gt;postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users </pre> Uroš Bernik 2012-05-19T07:34:07 http://comments.gmane.org/gmane.comp.gis.postgis/31784 <pre> 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</pre> David B�langer 2012-05-18T00:13:06 Search the mailing list at Gmane query http://search.gmane.org/?group=$group=gmane.comp.gis.postgis