Re: GDAL package error message at install

[PhilipNienhuis]

Lester Anderson wrote
> Hi Philip,
> 
> Thanks for the information about 4.0.1. Have installed the updated
> version and output of rasterinfo.m seems to be ok, mainly looking for
> the projection information:
> 
> Just a test output, does this look like the correct output for a
> geotiff, with mapping 1.2.1 installed:
> 
>>> info = rasterinfo('Topo1.tif')
> warning: your version of GraphicsMagick limits images to 8 bits per pixel
> warning: called from
>     imformats>default_formats at line 256 column 11
>     imformats at line 79 column 3
>     imageIO at line 106 column 11
>     imfinfo at line 159 column 8
>     rasterinfo at line 58 column 13
> info =
> 
>   scalar structure containing the fields:
> 
>     Filename = Q:\geophys\Lester\Topo1.tif
>     FileModDate = 20-Apr-2016 09:31:08
>     FileSize =  4822319
>     Format = TIFF
>     FormatVersion =
>     Width =  1559
>     Height =  917
>     BitDepth =  8
>     ColorType = grayscale
>     DelayTime = 0
>     DisposalMethod =
>     LoopCount = 0
>     ByteOrder = undefined
>     Gamma = 0
>     Chromaticities = [](1x0)
>     Comment =
>     Quality =  75
>     Compression = undefined
>     Colormap = [](0x0)
>     Orientation =  1
>     ResolutionUnit = Inch
>     XResolution = 0
>     YResolution = 0
>     Software =
>     Make =
>     Model =
>     DateTime =
>     ImageDescription =
>     Artist =
>     Copyright =
>     DigitalCamera =
> 
>       scalar structure containing the fields:
> 
> 
>     GPSInfo =
> 
>       scalar structure containing the fields:
> 
> 
>     FileType = GeoTIFF
>     datatype_name = Float32
>     GeoTransformation =
> 
>        47.72605    0.00439    0.00000   27.31861    0.00000   -0.00398
> 
>     Projection = 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"]]

Looks good :-)

Behind the scenes rasterinfo invokes the core Octave function imfinfo() and
that one sometime returns unexpected info for indexed geotiff files.  I saw
a bug report coming by that seems to address this but haven't looked
further.


> Is there a way to just limit the output info to the projection
> information?

No not yet.
But as info is a struct, you could do:

info = rasterinfo (....);
projection = info.Projection

I haven't implemented parsing that Projection field yet into separate fields
like Matlab's geotiffinfo does. Reasons: 
- no time yet, 
- I think someone somewhere might already have done that so I don't need to
reinvent the wheel; but again no time to dive into it. 

Would you be able to search for that and point us here/me to it if you find
something useful?

Philip




--
View this message in context: 
http://octave.1599824.n4.nabble.com/GDAL-package-error-message-at-install-tp4676356p4676394.html
Sent from the Octave - General mailing list archive at Nabble.com.