1
|
|
|
<?php |
2
|
|
|
|
3
|
|
|
namespace Smindel\GIS\Model; |
4
|
|
|
|
5
|
|
|
use SilverStripe\Core\Config\Configurable; |
6
|
|
|
use SilverStripe\Core\Injector\Injectable; |
7
|
|
|
use Smindel\GIS\GIS; |
8
|
|
|
|
9
|
|
|
/* |
10
|
|
|
gdaldem hillshade -of PNG public/assets/wellington-lidar-1m-dem-2013.4326.tif hillshade.png |
11
|
|
|
*/ |
12
|
|
|
|
13
|
|
|
class Raster |
14
|
|
|
{ |
15
|
|
|
use Configurable; |
16
|
|
|
|
17
|
|
|
use Injectable; |
18
|
|
|
|
19
|
|
|
private static $tile_renderer = 'raster_renderer'; |
|
|
|
|
20
|
|
|
|
21
|
|
|
protected $filename; |
22
|
|
|
protected $info; |
23
|
|
|
|
24
|
3 |
|
public function __construct($filename = null) |
25
|
|
|
{ |
26
|
3 |
|
$this->filename = $filename; |
27
|
3 |
|
} |
28
|
|
|
|
29
|
3 |
|
public function getFilename() |
30
|
|
|
{ |
31
|
3 |
|
return $this->filename ?: $this->config()->full_path; |
32
|
|
|
} |
33
|
|
|
|
34
|
2 |
|
public function getSrid() |
35
|
|
|
{ |
36
|
2 |
|
if (empty($this->info['srid'])) { |
37
|
|
|
|
38
|
2 |
|
$cmd = sprintf(' |
39
|
|
|
gdalsrsinfo -o wkt %1$s', |
40
|
2 |
|
$this->getFilename() |
41
|
|
|
); |
42
|
|
|
|
43
|
2 |
|
$output = `$cmd`; |
44
|
|
|
|
45
|
2 |
|
if (preg_match('/\WAUTHORITY\["EPSG","([^"]+)"\]\]$/', $output, $matches)) { |
46
|
2 |
|
$this->info['srid'] = $matches[1]; |
47
|
|
|
} |
48
|
|
|
} |
49
|
|
|
|
50
|
2 |
|
return $this->info['srid']; |
51
|
|
|
} |
52
|
|
|
|
53
|
1 |
|
public function getLocationInfo($geo = null, $band = null) |
54
|
|
|
{ |
55
|
1 |
|
$cmd = sprintf(' |
56
|
|
|
gdallocationinfo -wgs84 %1$s %2$s %3$s', |
57
|
1 |
|
$band ? sprintf('-b %d', $band) : '', |
58
|
1 |
|
$this->getFilename(), |
59
|
1 |
|
$geo ? sprintf('%f %f', ...GIS::create($geo)->reproject(4326)->coordinates) : '' |
60
|
|
|
); |
61
|
|
|
|
62
|
1 |
|
$output = `$cmd`; |
63
|
|
|
|
64
|
1 |
|
if (preg_match_all('/\sBand\s*(\d+):\s*Value:\s*([\d\.\-]+)/', $output, $matches)) { |
65
|
1 |
|
$bands = array_combine($matches[1], $matches[2]); |
66
|
1 |
|
array_walk($bands, function(&$item) { |
67
|
1 |
|
$item = (int)$item; |
68
|
1 |
|
}); |
69
|
1 |
|
return $bands; |
70
|
|
|
} |
71
|
|
|
} |
72
|
|
|
|
73
|
1 |
|
public function translateRaster($topLeftGeo, $bottomRightGeo, $width, $height, $destFileName = '/dev/stdout') |
74
|
|
|
{ |
75
|
1 |
|
$topLeftGeo = GIS::create($topLeftGeo)->coordinates; |
76
|
1 |
|
$bottomRightGeo = GIS::create($bottomRightGeo)->coordinates; |
77
|
|
|
|
78
|
1 |
|
$cmd = sprintf(' |
79
|
|
|
gdal_translate -of PNG -q -projwin %1$f, %2$f, %3$f, %4$f -outsize %5$d %6$d %7$s %8$s', |
80
|
1 |
|
$topLeftGeo[0], $topLeftGeo[1], |
81
|
1 |
|
$bottomRightGeo[0], $bottomRightGeo[1], |
82
|
1 |
|
$width, $height, |
83
|
1 |
|
$this->getFilename(), |
84
|
1 |
|
$destFileName |
85
|
|
|
); |
86
|
|
|
|
87
|
1 |
|
return `$cmd`; |
88
|
|
|
} |
89
|
|
|
|
90
|
1 |
|
public function searchableFields() |
91
|
|
|
{ |
92
|
|
|
return [ |
93
|
1 |
|
'Band' => 'Band', |
94
|
|
|
]; |
95
|
|
|
} |
96
|
|
|
} |
97
|
|
|
|