Thank you Raul that was the issue! Need to find a different algorithm or approach for generating my dot density map tiles.
Sent from my iPad > On Jun 13, 2018, at 3:00 PM, [email protected] wrote: > > Send postgis-users mailing list submissions to > [email protected] > > To subscribe or unsubscribe via the World Wide Web, visit > https://lists.osgeo.org/mailman/listinfo/postgis-users > or, via email, send a message with subject or body 'help' to > [email protected] > > You can reach the person managing the list at > [email protected] > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of postgis-users digest..." > > > Today's Topics: > > 1. Clarification ST_AsMVTGEOM on Multipoint Geometry (David Krett) > 2. Re: Clarification ST_AsMVTGEOM on Multipoint Geometry > (Raúl Marín Rodríguez) > 3. raster operations (David Haynes) > 4. Re: raster operations (Regina Obe) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Tue, 12 Jun 2018 19:06:40 -0400 > From: David Krett <[email protected]> > To: [email protected] > Subject: [postgis-users] Clarification ST_AsMVTGEOM on Multipoint > Geometry > Message-ID: <[email protected]> > Content-Type: text/plain; charset=us-ascii > > Hopefully someone can point me in the direction, > > I am currently using Postgresql 10 with the latest stable version of postgis > and serving Tiles from a table in my DB. The geometry column is comprised of > Multipoint Features, and not withstanding rows contain multiple points only > 1 point is being plotted on the Vector Tiles generated for each geographic > area. Tried finding documentation re this but really couldnt find much.. So > could a fellow user confirm whether there are any restrictions working with > multipoint geometries are whether I am encountering a bug. > > Many Thanks, > > David Krett > > ------------------------------ > > Message: 2 > Date: Wed, 13 Jun 2018 09:51:56 +0200 > From: Raúl Marín Rodríguez <[email protected]> > To: PostGIS Users Discussion <[email protected]> > Subject: Re: [postgis-users] Clarification ST_AsMVTGEOM on Multipoint > Geometry > Message-ID: > <cam6_um4db39j8qkuo6i585kkvc7vz0hovny2f5o+pxmq3fw...@mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Hi David, > > ST_AsMVTGeom will remove any duplicated points if they have the same > coordinates after switching to tile coordinates. Can you check the output > of using it on one of your rows? > > SELECT ST_ASText(ST_AsMVTGeom(ST_GeomFromText('MULTIPOINT(25 17, 26 18)'), > > ST_MakeBox2D(ST_Point(0, 0), ST_Point(4096, 4096)), 4096, 0, > false)); > st_astext > ----------------------------- > MULTIPOINT(25 4079,26 4078) > (1 row) > > >> On Wed, Jun 13, 2018 at 1:06 AM, David Krett <[email protected]> wrote: >> >> Hopefully someone can point me in the direction, >> >> I am currently using Postgresql 10 with the latest stable version of >> postgis and serving Tiles from a table in my DB. The geometry column is >> comprised of Multipoint Features, and not withstanding rows contain >> multiple points only 1 point is being plotted on the Vector Tiles >> generated for each geographic area. Tried finding documentation re this but >> really couldnt find much.. So could a fellow user confirm whether there are >> any restrictions working with multipoint geometries are whether I am >> encountering a bug. >> >> Many Thanks, >> >> David Krett >> _______________________________________________ >> postgis-users mailing list >> [email protected] >> https://lists.osgeo.org/mailman/listinfo/postgis-users > > > > > -- > > *Raúl Marín Rodríguez *carto.com > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: > <http://lists.osgeo.org/pipermail/postgis-users/attachments/20180613/0f657e92/attachment-0001.html> > > ------------------------------ > > Message: 3 > Date: Wed, 13 Jun 2018 11:56:29 -0500 > From: David Haynes <[email protected]> > To: [email protected] > Subject: [postgis-users] raster operations > Message-ID: > <CAM4C04qZoOaM2b3Z8ZjwBqP=nuxqfmjlc_tvk1s1bj++bkp...@mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > I have a question about how to focal operations for raster statistics are > working. I started off comparing the old depreciated > function st_mapalgebrafctngb with the current version st_mapalgebra. The > difference I have determined is that the st_mapalgebrafctngb will pass the > call back function to the raster dataset similar to the st_mapalgebra, but > as the function approaches the edge of the raster tile things get wonky. > What are values that are assigned to pixels that do not have sufficient > neighbors? When I run this bit of code. The ST_Value function reports that > no values have been assigned. I assume this can be remedied by adding > information through userargs. > > with smoothed_raster as > ( > SELECT st_mapalgebrafctngb(rast, 1, NULL, 1, 1, 'st_mean4ma(double > precision[][][],text,text[])'::regprocedure, 'ignore', NULL) as rast > FROM glc2000_clipped_250 r > WHERE r.rid = 25 > ), original_raster as > ( > SELECT r.rast > FROM glc2000_clipped_250 r > WHERE r.rid = 25 > ) > SELECT x, y, ST_Value(s.rast, x,y) as smoothed, ST_Value(o.rast, x,y) as > original > FROM smoothed_raster s, original_raster o > CROSS JOIN generate_series(1,10) as x > CROSS JOIN generate_series(1,10) as y > > This leads me to a larger question regarding focal analysis functions for > raster datasets. Do these functions cross the tiles? Assume you have 9 > tiles and each tile is composed of 3x3 pixels. Tiles 1-3 in row 1, 4-6 row > 2 etc. If you have a neighborhood analysis function performing on tile 5. > Will it grab adjacent pixels in tile 6? Unfortunately the examples I > created returned with the same result. > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: > <http://lists.osgeo.org/pipermail/postgis-users/attachments/20180613/336cb6c4/attachment-0001.html> > > ------------------------------ > > Message: 4 > Date: Wed, 13 Jun 2018 13:29:52 -0400 > From: "Regina Obe" <[email protected]> > To: "'PostGIS Users Discussion'" <[email protected]> > Subject: Re: [postgis-users] raster operations > Message-ID: <[email protected]> > Content-Type: text/plain; charset="utf-8" > > It doesn't grab adjacent tiles. To consider adjacent tiles you'd need to use > ST_Union to stitch together multiple tiles. > > > > From: postgis-users [mailto:[email protected]] On Behalf > Of David Haynes > Sent: Wednesday, June 13, 2018 12:56 PM > To: [email protected] > Subject: [postgis-users] raster operations > > > > I have a question about how to focal operations for raster statistics are > working. I started off comparing the old depreciated function > st_mapalgebrafctngb with the current version st_mapalgebra. The difference I > have determined is that the st_mapalgebrafctngb will pass the call back > function to the raster dataset similar to the st_mapalgebra, but as the > function approaches the edge of the raster tile things get wonky. What are > values that are assigned to pixels that do not have sufficient neighbors? > When I run this bit of code. The ST_Value function reports that no values > have been assigned. I assume this can be remedied by adding information > through userargs. > > > > with smoothed_raster as > > ( > > SELECT st_mapalgebrafctngb(rast, 1, NULL, 1, 1, 'st_mean4ma(double > precision[][][],text,text[])'::regprocedure, 'ignore', NULL) as rast > > FROM glc2000_clipped_250 r > > WHERE r.rid = 25 > > ), original_raster as > > ( > > SELECT r.rast > > FROM glc2000_clipped_250 r > > WHERE r.rid = 25 > > ) > > SELECT x, y, ST_Value(s.rast, x,y) as smoothed, ST_Value(o.rast, x,y) as > original > > FROM smoothed_raster s, original_raster o > > CROSS JOIN generate_series(1,10) as x > > CROSS JOIN generate_series(1,10) as y > > > > This leads me to a larger question regarding focal analysis functions for > raster datasets. Do these functions cross the tiles? Assume you have 9 tiles > and each tile is composed of 3x3 pixels. Tiles 1-3 in row 1, 4-6 row 2 etc. > If you have a neighborhood analysis function performing on tile 5. Will it > grab adjacent pixels in tile 6? Unfortunately the examples I created returned > with the same result. > > > > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: > <http://lists.osgeo.org/pipermail/postgis-users/attachments/20180613/07b39323/attachment-0001.html> > > ------------------------------ > > Subject: Digest Footer > > _______________________________________________ > postgis-users mailing list > [email protected] > https://lists.osgeo.org/mailman/listinfo/postgis-users > > ------------------------------ > > End of postgis-users Digest, Vol 196, Issue 8 > ********************************************* _______________________________________________ postgis-users mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/postgis-users
