[postgis-devel] pgraster: st_resample error

[Tom van Tilburg]

Hi Bborie,

What's the best way to check that?

Another thing I noticed:
Using gdalwarp in qgis from 28992 to 900913 gives me this 20 km offset. 
However, when I manually force the warped image to be shown in 3875 
projection, it shows up well. Using gdalinfo on the warped image it 
strikes me that there is no information on the 'authority' of the PROJCS 
(should be 900913). (see below)

Cheers,
  Tom


Driver: GTiff/GeoTIFF
Files: ahn1_900913.asc
Size is 788, 690
Coordinate System is:
PROJCS["Google Maps Global Mercator",
     GEOGCS["WGS 84",
         DATUM["WGS_1984",
             SPHEROID["WGS 84",6378137,298.257223563,
                 AUTHORITY["EPSG","7030"]],
             AUTHORITY["EPSG","6326"]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433],
         AUTHORITY["EPSG","4326"]],
     PROJECTION["Mercator_1SP"],
     PARAMETER["central_meridian",0],
     PARAMETER["scale_factor",1],
     PARAMETER["false_easting",0],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]]]
Origin = (656440.434583136229776,6806151.070541447028518)
Pixel Size = (8.117405487721264,-8.117405487721264)
Metadata:
   AREA_OR_POINT=Area
Image Structure Metadata:
   INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  656440.435, 6806151.071) (  5d53'48.86"E, 52d13'10.86"N)
Lower Left  (  656440.435, 6800550.061) (  5d53'48.86"E, 52d11'19.57"N)
Upper Right (  662836.950, 6806151.071) (  5d57'15.72"E, 52d13'10.86"N)
Lower Right (  662836.950, 6800550.061) (  5d57'15.72"E, 52d11'19.57"N)
Center      (  659638.692, 6803350.566) (  5d55'32.29"E, 52d12'15.23"N)
Band 1 Block=788x2 Type=Float32, ColorInterp=Gray


On 28-12-2011 1:57, Bborie Park wrote:
> Hey Tom,
>
> I'm basically wondering if for a latitude/longitude, the value at that
> coordinate is correct regardless of projection.
>
> -bborie
>
> On Tue, Dec 27, 2011 at 8:06 AM, Tom van Tilburg
>