<html><head><style>body{font-family:Helvetica,Arial;font-size:13px}</style></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div id="bloop_customfont" style="font-family:Helvetica,Arial;font-size:13px; color: rgba(0,0,0,1.0); margin: 0px; line-height: auto;">So, I have a “solution” but I’m not 100% happy with it. It’s attached to the ticket here,</div><div id="bloop_customfont" style="font-family:Helvetica,Arial;font-size:13px; color: rgba(0,0,0,1.0); margin: 0px; line-height: auto;"><br></div><div id="bloop_customfont" style="font-family:Helvetica,Arial;font-size:13px; color: rgba(0,0,0,1.0); margin: 0px; line-height: auto;">http://trac.osgeo.org/postgis/ticket/2953</div><div id="bloop_customfont" style="font-family:Helvetica,Arial;font-size:13px; color: rgba(0,0,0,1.0); margin: 0px; line-height: auto;"><br></div><div id="bloop_customfont" style="font-family:Helvetica,Arial;font-size:13px; color: rgba(0,0,0,1.0); margin: 0px; line-height: auto;">P.</div> <div id="bloop_sign_1412270689335729152" class="bloop_sign">

        <title></title>

        <div>

            <br>

        </div>

        -- <br>

        Paul Ramsey<br>

        http://cleverelephant.ca<div>http://postgis.net

</div>

<script type="application/ld+json">

{

  "@context": "http://schema.org",

  "@type": "em",

  "name": "John Doe",

  "jobTitle": "Graduate research assistant",

  "affiliation": "University of Dreams",

  "additionalName": "Johnny",

  "url": "http://www.example.com",

  "address": {

    "@type": "PostalAddress",

    "streetAddress": "1234 Peach Drive",

    "addressLocality": "Wonderland",

    "addressRegion": "Georgia"

  }

}

</script>

</div> <br><p style="color:#000;">On October 2, 2014 at 9:58:13 AM, Paul Ramsey (<a href="mailto:pramsey@cleverelephant.ca">pramsey@cleverelephant.ca</a>) wrote:</p> <blockquote type="cite" class="clean_bq"><span><div><div></div><div>OK, with the real test data I can reproduce, and the source of the

<br>problem is pretty "obvious" if not the exact place things are going

<br>wrong yet:

<br>

<br>                        st_astext

<br>---------------------------------------------------------

<br> POINT ZM (-124.22007 41.75967 0 -1.79769313486232e+308)

<br><br>Those are some big-ass "M" values. So the problem is almost certainly

<br>with the higher dimensions. Since you aren't using them, a simple

<br>expedient to fix your problem would be to just alter the type to a 2d

<br>point, and then the stats gathering system will work right. It's worth

<br>doing this, since as the NOTICE says, you aren't currently gathering

any stats.

<br>

<br>P.

<br>

<br>

<br>On Wed, Oct 1, 2014 at 3:26 PM, Paul Ramsey <pramsey@cleverelephant.ca> wrote:

> Well, trying to synthetically reproduce your error is not working.

<br>>

<br>> Butterfly:~/Code/postgis-svn/postgis pramsey$ vacuumdb -v -Z -d postgis22 -t pts

<br>>

<br>> INFO:  analyzing "public.pts"

<br>>

<br>> INFO:  "pts": scanned 14286 of 14286 pages, containing 1000000 live

<br>> rows and 0 dead rows; 30000 rows in sample, 1000000 estimated total

<br>> rows

<br>>

<br>> I used a basic random table,

<br>>

<br>> CREATE TABLE pts AS

<br>> WITH rands AS (

<br>>   SELECT generate_series as id, random() AS u1, random() AS u2 FROM

<br>> generate_series(1,1000000)

<br>> )

<br>> SELECT

<br>>   id,

<br>>   ST_SetSRID(ST_MakePoint(

<br>>     50 * sqrt(-2 * ln(u1)) * cos(2*pi()*u2),

<br>>     50 * sqrt(-2 * ln(u1)) * sin(2*pi()*u2),

<br>>     50 * sqrt(-2 * ln(u1)) * sin(2*pi()*u2),

<br>>     50 * sqrt(-2 * ln(u1)) *

<br>> sin(2*pi()*u2)),4326)::Geometry(PointZM,4326) AS geom

<br>> FROM rands;

<br>>

<br>> P

<br>>

<br>> On Tue, Sep 23, 2014 at 3:54 PM, Paul Ramsey <pramsey@cleverelephant.ca> wrote:

>> I guess in theory I could reproduce it with some random xyz data...

<br>>>

<br>>> P

<br>>>

<br>>> On Tue, Sep 23, 2014 at 1:47 PM, Jason Mathis

<br>>> <jmathis@redzonesoftware.com> wrote:

<br>>>> So unfortunately I did not get the go ahead to share the data, but I solved

<br>>>> the problem. It was a dataset that was merged from several different

<br>>>> datasets which is why I got the “PointZM” from shp2pgsql. I was given a

<br>>>> shape file to load into the database. Although the data just needed to be

<br>>>> two dimensional. I also had the original coordinates, so I drop the column

<br>>>> and re-added with just “POINT” and updated the geometry from the lat/longs

>>> and viola no “notice” message.

<br>>>>

>>> Is it only an issue with three or four dimensional data?

<br>>>>

<br>>>> Thanks,

<br>>>> jason

<br>>>>

<br>>>> On September 22, 2014 at 5:43:43 PM, Paul Ramsey (pramsey@cleverelephant.ca)

<br>>>> wrote:

<br>>>>

>>> If you can share data, I'm happy to try and reproduce it here.

<br>>>>

<br>>>> P

<br>>>>

<br>>>> On Mon, Sep 22, 2014 at 2:01 PM, Jason Mathis

<br>>>> <jmathis@redzonesoftware.com> wrote:

<br>>>>> I am stumped on it. I think i remember seeing something like this before,

>>>> but can’t quite recall.

<br>>>>>

>>>> Queries both return zero.

<br>>>>>

<br>>>>> 'POSTGIS="2.1.3 r12547" GEOS="3.4.2-CAPI-1.8.2 r3921" PROJ="Rel. 4.8.0, 6

<br>>>>> March 2012" GDAL="GDAL 1.9.2, released 2012/10/08" LIBXML="2.7.6"

<br>>>>> LIBJSON="UNKNOWN" RASTER'

<br>>>>>

<br>>>>> i even tried putting a bogus null row in there but I still got the

<br>>>>> message.

<br>>>>>

<br>>>>> thanks!

<br>>>>>

<br>>>>>

<br>>>>> On September 22, 2014 at 2:52:00 PM, Paul Ramsey

<br>>>>> (pramsey@cleverelephant.ca)

<br>>>>> wrote:

<br>>>>>

<br>>>>> It means it sampled 30000 rows and didn't find any geometries that

<br>>>>> were non-null or non-empty. I suppose it could mean a bug, though that

>>>> code's been around for a while.

<br>>>>>

<br>>>>> select count(*) from yourtable where geom is null;

<br>>>>> select count(*) from yourtable where st_isempty(geom);

<br>>>>> select postgis_full_version();

<br>>>>>

<br>>>>> P.

<br>>>>>

<br>>>>> On Mon, Sep 22, 2014 at 12:11 PM, Jason Mathis

<br>>>>> <jmathis@redzonesoftware.com> wrote:

<br>>>>>> Does anyone know what this “notice” means? I loaded some new data into a

>>>>> table from a shape file.Visually it looks good and loaded without error.

<br>>>>>> I

<br>>>>>> have a nightly analyze job that runs and since loading the data it spits

<br>>>>>> out

<br>>>>>> this “notice” I don’t know what it means. If I drop the geom column I do

<br>>>>>> not

<br>>>>>> get this message. I have tried drop/create the index, and reloading the

>>>>> data without success.

<br>>>>>>

<br>>>>>> Thanks!

<br>>>>>>

<br>>>>>>

<br>>>>>>

<br>>>>>> $ /usr/bin/vacuumdb -v -Z -d redzone -t 'demo.homes'

<br>>>>>> INFO: analyzing "demo.homes"

<br>>>>>> INFO: "homes": scanned 2237 of 2237 pages, containing 101937 live rows

<br>>>>>> and

<br>>>>>> 0 dead rows; 30000 rows in sample, 101937 estimated total rows

<br>>>>>> NOTICE: no non-null/empty features, unable to compute statistics

<br>>>>>> NOTICE: no non-null/empty features, unable to compute statistics

<br>>>>>>

<br>>>>>> CREATE TABLE demo.homes

<br>>>>>> (

<br>>>>>> gid integer NOT NULL DEFAULT nextval('demo.homes_gid_seq'::regclass),

<br>>>>>> geolat numeric,

<br>>>>>> geolong numeric,

<br>>>>>> geolevel text,

<br>>>>>> geonumber integer,

<br>>>>>> geostreet text,

<br>>>>>> geocity text,

<br>>>>>> geocounty text,

<br>>>>>> geostate text,

<br>>>>>> zip text,

<br>>>>>> firstname text,

<br>>>>>> lastname text,

<br>>>>>> tiv integer,

<br>>>>>> geom geometry(PointZM,4326),

<br>>>>>> CONSTRAINT homes_pkey PRIMARY KEY (gid)

<br>>>>>> )

<br>>>>>>

<br>>>>>>

<br>>>>>>

<br>>>>>>

<br>>>>>>

<br>>>>>>

<br>>>>>> This transmission contains confidential and privileged information

<br>>>>>> intended

<br>>>>>> solely for the party identified above. If you receive this message in

<br>>>>>> error,

<br>>>>>> you must not use it or convey it to others. Please destroy it immediately

<br>>>>>> and contact the sender at (303) 386-3955 or by return e-mail to the

<br>>>>>> sender.

<br>>>>>>

<br>>>>>>

<br>>>>>> _______________________________________________

<br>>>>>> postgis-users mailing list

<br>>>>>> postgis-users@lists.osgeo.org

<br>>>>>> http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users

<br>>>>> _______________________________________________

<br>>>>> postgis-users mailing list

<br>>>>> postgis-users@lists.osgeo.org

<br>>>>> http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users

<br>>>>>

<br>>>>>

<br>>>>> This transmission contains confidential and privileged information

<br>>>>> intended

<br>>>>> solely for the party identified above. If you receive this message in

<br>>>>> error,

<br>>>>> you must not use it or convey it to others. Please destroy it immediately

<br>>>>> and contact the sender at (303) 386-3955 or by return e-mail to the

<br>>>>> sender.

<br>>>>

<br>>>>

<br>>>> This transmission contains confidential and privileged information intended

<br>>>> solely for the party identified above. If you receive this message in error,

<br>>>> you must not use it or convey it to others. Please destroy it immediately

>>> and contact the sender at (303) 386-3955 or by return e-mail to the sender.

<br></div></div></span></blockquote></body></html>