1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
|
## Copyright (C) 2015-2020 Philip Nienhuis
##
## This program is free software; you can redistribute it and/or modify it
## under the terms of the GNU General Public License as published by
## the Free Software Foundation; either version 3 of the License, or
## (at your option) any later version.
##
## This program is distributed in the hope that it will be useful,
## but WITHOUT ANY WARRANTY; without even the implied warranty of
## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
## GNU General Public License for more details.
##
## You should have received a copy of the GNU General Public License
## along with this program. If not, see <http://www.gnu.org/licenses/>.
## -*- texinfo -*-
## @deftypefn {} {@var{rinfo} =} rasterinfo (@var{fname})
## Return various info about a GIS raster file: a.o., file type, bit
## depth, raster size, projection and geotransformation. If the raster
## file is a geotiff file, additional info is returned.
##
## rasterinfo needs the GDAL library.
##
## @seealso{rasterread}
## @end deftypefn
## Author: Philip Nienhuis <prnienhuis@users.sf.net>
## Created: 2015-12-23
function [rinfo] = rasterinfo (fname)
persistent HAVE_GDAL = [];
## Checks
if (! ischar (fname))
print_usage ();
elseif (isempty (fname))
print_usage ();
elseif (exist (fname) != 2)
error ("rasterinfo: file '%s' not found\n", fname);
endif
bands = info = {};
## Check if we have gdalread.oct installed (depends on GDAL)
if (isempty (HAVE_GDAL))
HAVE_GDAL = ! isempty (which ("gdalread"));
endif
if (! HAVE_GDAL)
error ("rasterinfo: GDAL library is apparently not installed");
endif
## Try to read data file
[rstatus, binfo] = gdalread (fname, 0);
if (! rstatus)
## Check if we have a geotiff file. If so, get extra info
if (strcmpi (binfo.FileType, "geotiff"))
rinfo = imfinfo (fname);
rinfo.FileType = binfo.FileType;
rinfo.datatype_name = binfo.datatype_name;
rinfo.GeoTransformation = binfo.GeoTransformation;
rinfo.Projection = binfo.Projection;
rinfo.nbands = binfo.nbands;
if (isfield (binfo, "bbox"))
rinfo.bbox = binfo.bbox;
else
gt = binfo.GeoTransformation;
rinfo.bbox = zeros (2);
rinfo.bbox(1, 1) = gt(1);
rinfo.bbox(2, 1) = gt(1) + (binfo.Width * gt(2)) + (binfo.Height * gt(3));
if (gt(2) < 0)
rinfo.bbox(:, 1) = flipud (rinfo.bbox(:, 1));
endif
rinfo.bbox(1, 2) = gt(4);
rinfo.bbox(2, 2) = gt(4) + (binfo.Height * gt(6)) + (binfo.Width * gt(5));
if (gt(6) < 0)
rinfo.bbox(:, 2) = flipud (rinfo.bbox(:, 2));
endif
endif
else
rinfo = binfo;
rinfo.Filename = canonicalize_file_name (fname);
rinfo.FileSize = stat (fname).size;
rinfo.FileModDate = strftime ("%d-%b-%y %H:%M:%S", localtime (stat (fname).mtime));
endif
endif
endfunction
|
