Changes between Version 84 and Version 85 of PostGIS_Raster_SoC_Idea_2012/Distance_Analysis_Tools

Timestamp:

Aug 17, 2012, 12:45:23 PM (12 years ago)

Author:

qliu

Comment:

—

PostGIS_Raster_SoC_Idea_2012/Distance_Analysis_Tools

@@ -972 +972 @@
             END LOOP;
         END LOOP;
-                
-                -- There is place for optimization here for a more efficient scanning method
-        FOR x IN 1..width LOOP
-            FOR y IN 1..height LOOP
-                                newx := ST_Raster2WorldCoordX(newrast, x, y);
-                                newy := ST_Raster2WorldCoordY(newrast, x, y); 
-                                -- If pixel is the source point then Eculidean distance is zero
-                                IF THEN
-                                        newval := 0;
-                                ELSE
-                                        exesql := "SELECT ST_Distance(ST_GeomFromText('POINT(" || newx || " " || newy || ")'," || newsrid || "),ST_Transform(" || sourcegeomcolumn || "," || newsrid || ")) FROM " || sourcetable || " ORDER BY $1 <-> $2 LIMIT 1;";
-                                        newval := EXECUTE(exesql);
-                                END IF;
-                                newrast := ST_SetValue(newrast, band, x, y, newval);
-            END LOOP;
-        END LOOP;
+
         RETURN newrast;
     END;