आपको स्टोर की जाँच करनी चाहिए , पियरे रेसीन और स्टीव कमिंग द्वारा पोस्टग्रेक्यूएल / पोस्टगिस स्थानिक डेटाबेस दस्तावेज़ के भीतर रेखापुंज डेटा को हेरफेर और विश्लेषण करना चाहिए जो यहां एफओएसएस में प्रस्तुत किया गया था । आपकी समस्या को हल करने के लिए बहुत सारे फ़ंक्शन को रेखापुंज आँकड़ों के रूप में परिभाषित किया गया है। मुझे लगता है कि ST_SummaryStats आपको ज़ोनल स्टैटिक्स के बारे में मदद करेगा, ज़ाहिर है, पर्याप्त नहीं।
ST_SummaryStats(raster)
आपको (min, max, sum, mean, stddev, count (withdata pixels)) रिकॉर्ड का एक सेट चालू कर देगा।
पोस्टगिस 2.0 कुछ एसक्यूएल प्रश्नों के साथ आंचलिक सांख्यिकी विश्लेषण का समर्थन करता है, जिस पर आपके पास काम है। मैं esri बनाम पोस्टगिस 2.0 के साथ आंचलिक सांख्यिकी के बारे में नेट पर एक अच्छे प्रलेखन में आया हूं। यहाँ चलती स्थानिक में आपको कुछ जानकारी के लिए जाँच करनी चाहिए । उन दोनों को कुछ कठोर क्वेरी के साथ समान आँकड़े परिणाम दिए गए थे।
उदाहरण क्वेरी जो एंथोनी लोपेज साइट से ली गई थी, जिसका मैंने पहले उल्लेख किया है ...
CREATE TABLE sum_pop2 AS
WITH
feat AS (SELECT gid, geom FROM perez_grid AS b ),
b_stats AS
(SELECT gid, (stats).*
FROM (
SELECT gid, ST_SummaryStats(ST_Clip(rast,1,geom)) AS stats
FROM ls_den
INNER JOIN feat
ON ST_Intersects(feat.geom,rast) ) AS foo )
SELECT gid, SUM(count) AS cell_count
,SUM(sum) AS population
FROM b_stats
WHERE count > 0
GROUP BY gid
ORDER BY gid;
एसरी जोनल स्टैटिस्टिक्स परिणाम: जनसंख्या: 207,578 सेल गणना: 14,400
पोस्टगिस विधि परिणाम: जनसंख्या: 207,578 सेल गणना: 14,400।
जनसंख्या का अंतर: 0%
मैं यह अभी तक की कोशिश की लेकिन एंथोनी ने उल्लेख किया कि पोस्टगिस के साथ रेखापुंज विश्लेषण के साथ कुछ प्रदर्शन समस्या थी।
मुझे उम्मीद है इससे आपको मदद मिली होगी...