[R-sig-Geo] Conflicting information between GDALinfo() and gdalinfo

Roger Bivand Roger.Bivand at nhh.no
Thu Jun 14 18:19:27 CEST 2012


On Thu, 14 Jun 2012, Agustin Lobo wrote:

> Taking into consideration that the user expects GDALinfo to report the
> same info as gdalinfo,
> it would make sense that GDALinfo state "lower left origin" instead of
> just "origin"

Committed to R-Forge.

Roger

>
> Thanks
>
> Agus
>
> 2012/6/14 Roger Bivand <Roger.Bivand at nhh.no>:
>> GDALinfo() reports the lower left corner, which agrees with gdalinfo.
>> Whether the origin is upper left or lower left is arbitrary, in sp and hence
>> in rgdal, the origin of rasters is lower left.
>>
>> Roger
>>
>>
>> On Thu, 14 Jun 2012, Agustin Lobo wrote:
>>
>>> Hi!
>>> I'm having different origin according to gdalinfo ( "Origin =
>>> (0.000000000000000,1024.000000000000000)") and GDALinfo
>>> ("TTC0125v2flip.tif"). When displayed in QGIS, the origin is 0,1024
>>> (using ubuntu 10.04, see sessionInfo at the end, tif file here:
>>> http://dl.dropbox.com/u/3180464/TTC0125v2flip.tif
>>>
>>> tif file was saved with
>>> writeRaster(bout,filename="TTC0125v2flip.tif", format="GTiff",
>>> overwrite=TRUE,flip=T)
>>>
>>>> system("gdalinfo TTC0125v2flip.tif",intern=T)
>>>
>>> [1] "Driver: GTiff/GeoTIFF"
>>> "Files: TTC0125v2flip.tif"
>>> [3] "Size is 1280, 1024"
>>> "Coordinate System is `'"
>>> [5] "Origin = (0.000000000000000,1024.000000000000000)"
>>> "Pixel Size = (1.000000000000000,-1.000000000000000)"
>>> [7] "Image Structure Metadata:"                                   "
>>> COMPRESSION=LZW"
>>> [9] "  INTERLEAVE=PIXEL"
>>> "Corner Coordinates:"
>>> [11] "Upper Left  (       0.000,    1024.000) "
>>> "Lower Left  (   0.0000000,   0.0000000) "
>>> [13] "Upper Right (    1280.000,    1024.000) "
>>> "Lower Right (    1280.000,       0.000) "
>>> [15] "Center      (     640.000,     512.000) "
>>> "Band 1 Block=1280x1 Type=Float32, ColorInterp=Gray"
>>> [17] "  Min=12.000 Max=219.000 "                                   "
>>> Minimum=12.000, Maximum=219.000, Mean=0.000, StdDev=0.000"
>>> [19] "  NoData Value=-3.39999999999999996e+38"                     "
>>> Metadata:"
>>> [21] "    STATISTICS_MAXIMUM=219"                                  "
>>> STATISTICS_MEAN=0"
>>> [23] "    STATISTICS_MINIMUM=12"                                   "
>>> STATISTICS_STDDEV=0"
>>> [25] "Band 2 Block=1280x1 Type=Float32, ColorInterp=Undefined"     "
>>> Min=0.000 Max=215.000 "
>>> [27] "  Minimum=0.000, Maximum=215.000, Mean=0.000, StdDev=0.000"  "
>>> NoData Value=-3.39999999999999996e+38"
>>> [29] "  Metadata:"                                                 "
>>> STATISTICS_MAXIMUM=215"
>>> [31] "    STATISTICS_MEAN=0"                                       "
>>> STATISTICS_MINIMUM=0"
>>> [33] "    STATISTICS_STDDEV=0"
>>> "Band 3 Block=1280x1 Type=Float32, ColorInterp=Undefined"
>>> [35] "  Min=0.000 Max=234.000 "                                    "
>>> Minimum=0.000, Maximum=234.000, Mean=0.000, StdDev=0.000"
>>> [37] "  NoData Value=-3.39999999999999996e+38"                     "
>>> Metadata:"
>>> [39] "    STATISTICS_MAXIMUM=234"                                  "
>>> STATISTICS_MEAN=0"
>>> [41] "    STATISTICS_MINIMUM=0"                                    "
>>> STATISTICS_STDDEV=0"
>>> [43] "Band 4 Block=1280x1 Type=Float32, ColorInterp=Undefined"     "
>>> Min=0.000 Max=255.000 "
>>> [45] "  Minimum=0.000, Maximum=255.000, Mean=0.000, StdDev=0.000"  "
>>> NoData Value=-3.39999999999999996e+38"
>>> [47] "  Metadata:"                                                 "
>>> STATISTICS_MAXIMUM=255"
>>> [49] "    STATISTICS_MEAN=0"                                       "
>>> STATISTICS_MINIMUM=0"
>>> [51] "    STATISTICS_STDDEV=0"
>>> "Band 5 Block=1280x1 Type=Float32, ColorInterp=Undefined"
>>> [53] "  Min=0.000 Max=221.000 "                                    "
>>> Minimum=0.000, Maximum=221.000, Mean=0.000, StdDev=0.000"
>>> [55] "  NoData Value=-3.39999999999999996e+38"                     "
>>> Metadata:"
>>> [57] "    STATISTICS_MAXIMUM=221"                                  "
>>> STATISTICS_MEAN=0"
>>> [59] "    STATISTICS_MINIMUM=0"                                    "
>>> STATISTICS_STDDEV=0"
>>> [61] "Band 6 Block=1280x1 Type=Float32, ColorInterp=Undefined"     "
>>> Min=0.000 Max=178.000 "
>>> [63] "  Minimum=0.000, Maximum=178.000, Mean=0.000, StdDev=0.000"  "
>>> NoData Value=-3.39999999999999996e+38"
>>> [65] "  Metadata:"                                                 "
>>> STATISTICS_MAXIMUM=178"
>>> [67] "    STATISTICS_MEAN=0"                                       "
>>> STATISTICS_MINIMUM=0"
>>> [69] "    STATISTICS_STDDEV=0"
>>>>
>>>> GDALinfo("TTC0125v2flip.tif")
>>>
>>> rows        1024
>>> columns     1280
>>> bands       6
>>> origin.x        0
>>> origin.y        0
>>> res.x       1
>>> res.y       1
>>> ysign       -1
>>> oblique.x   0
>>> oblique.y   0
>>> driver      GTiff
>>> projection  NA
>>> file        TTC0125v2flip.tif
>>> apparent band summary:
>>>  GDType Bmin Bmax Bmean Bsd hasNoDataValue NoDataValue
>>> 1 Float32   12  219     0   0           TRUE    -3.4e+38
>>> 2 Float32    0  215     0   0           TRUE    -3.4e+38
>>> 3 Float32    0  234     0   0           TRUE    -3.4e+38
>>> 4 Float32    0  255     0   0           TRUE    -3.4e+38
>>> 5 Float32    0  221     0   0           TRUE    -3.4e+38
>>> 6 Float32    0  178     0   0           TRUE    -3.4e+38
>>>
>>>> sessionInfo()
>>>
>>> R version 2.15.0 (2012-03-30)
>>> Platform: x86_64-pc-linux-gnu (64-bit)
>>>
>>> locale:
>>> [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C
>>> LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8
>>> [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8    LC_PAPER=C
>>>              LC_NAME=C
>>> [9] LC_ADDRESS=C               LC_TELEPHONE=C
>>> LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
>>>
>>> attached base packages:
>>> [1] stats     graphics  grDevices utils     datasets  methods   base
>>>
>>> other attached packages:
>>> [1] raster_2.0-02 rgdal_0.7-5   sp_0.9-99
>>>
>>> loaded via a namespace (and not attached):
>>> [1] grid_2.15.0     lattice_0.19-30 tools_2.15.0
>>>
>>>
>>> Any help?
>>> Agus
>>>
>>> _______________________________________________
>>> R-sig-Geo mailing list
>>> R-sig-Geo at r-project.org
>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>>
>>
>> --
>> Roger Bivand
>> Department of Economics, NHH Norwegian School of Economics,
>> Helleveien 30, N-5045 Bergen, Norway.
>> voice: +47 55 95 93 55; fax +47 55 95 95 43
>> e-mail: Roger.Bivand at nhh.no
>>
>> _______________________________________________
>> R-sig-Geo mailing list
>> R-sig-Geo at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>

-- 
Roger Bivand
Department of Economics, NHH Norwegian School of Economics,
Helleveien 30, N-5045 Bergen, Norway.
voice: +47 55 95 93 55; fax +47 55 95 95 43
e-mail: Roger.Bivand at nhh.no


More information about the R-sig-Geo mailing list