+<?php declare(strict_types=1);
+
+/*
+ * This file is part of the Rapsys PackBundle package.
+ *
+ * (c) Raphaël Gertz <symfony@rapsys.eu>
+ *
+ * For the full copyright and license information, please view the LICENSE
+ * file that was distributed with this source code.
+ */
+
+namespace Rapsys\PackBundle\Util;
+
+use Symfony\Component\Filesystem\Exception\IOExceptionInterface;
+use Symfony\Component\Filesystem\Filesystem;
+
+/**
+ * Helps manage osm images
+ */
+class OsmUtil {
+ /**
+ * The tile size
+ */
+ const tz = 256;
+
+ /**
+ * The cache directory
+ */
+ protected $cache;
+
+ /**
+ * The public path
+ */
+ protected $path;
+
+ /**
+ * The tile server
+ */
+ protected $server;
+
+ /**
+ * The tile servers
+ *
+ * @see https://wiki.openstreetmap.org/wiki/Slippy_map_tilenames#Tile_servers
+ */
+ protected $servers = [
+ 'osm' => 'https://tile.openstreetmap.org/{Z}/{X}/{Y}.png',
+ 'cycle' => 'http://a.tile.thunderforest.com/cycle/{Z}/{X}/{Y}.png',
+ 'transport' => 'http://a.tile.thunderforest.com/transport/{Z}/{X}/{Y}.png'
+ ];
+
+ /**
+ * The public url
+ */
+ protected $url;
+
+ /**
+ * Creates a new osm util
+ *
+ * @param string $cache The cache directory
+ * @param string $path The public path
+ * @param string $url The public url
+ * @param string $server The server key
+ */
+ function __construct(string $cache, string $path, string $url, string $server = 'osm') {
+ //Set cache
+ $this->cache = $cache.'/'.$server;
+
+ //Set path
+ $this->path = $path.'/'.$server;
+
+ //Set url
+ $this->url = $url.'/'.$server;
+
+ //Set server key
+ $this->server = $server;
+ }
+
+ /**
+ * Return the image
+ *
+ * @desc Generate image in jpeg format or load it from cache
+ *
+ * @param string $pathInfo The path info
+ * @param string $alt The image alt
+ * @param int $updated The updated timestamp
+ * @param float $latitude The latitude
+ * @param float $longitude The longitude
+ * @param int $zoom The zoom
+ * @param int $width The width
+ * @param int $height The height
+ * @return array The image array
+ */
+ public function getImage(string $pathInfo, string $alt, int $updated, float $latitude, float $longitude, int $zoom = 18, int $width = 1280, int $height = 1280): array {
+ //Set path file
+ $path = $this->path.$pathInfo.'.jpeg';
+
+ //Without existing path
+ if (!is_dir($dir = dirname($path))) {
+ //Create filesystem object
+ $filesystem = new Filesystem();
+
+ try {
+ //Create path
+ //XXX: set as 0775, symfony umask (0022) will reduce rights (0755)
+ $filesystem->mkdir($dir, 0775);
+ } catch (IOExceptionInterface $e) {
+ //Throw error
+ throw new \Exception(sprintf('Output path "%s" do not exists and unable to create it', $dir), 0, $e);
+ }
+ }
+
+ //With path file
+ if (
+ false &&
+ is_file($path) &&
+ ($stat = stat($path)) &&
+ $stat['mtime'] >= $updated
+ ) {
+ //Return image data
+ return [
+ 'src' => $this->url.'/'.$stat['mtime'].$pathInfo.'.jpeg',
+ 'alt' => $alt,
+ 'height' => $height,
+ 'width' => $width
+ ];
+ }
+
+ //Create image instance
+ $image = new \Imagick();
+
+ //Add new image
+ $image->newImage($width, $height, new \ImagickPixel('transparent'), 'jpeg');
+
+ //Create tile instance
+ $tile = new \Imagick();
+
+ //Init context
+ $ctx = stream_context_create(
+ [
+ 'http' => [
+ #'header' => ['Referer: https://www.openstreetmap.org/'],
+ 'max_redirects' => 5,
+ 'timeout' => (int)ini_get('default_socket_timeout'),
+ #'user_agent' => 'Mozilla/5.0 (X11; Linux x86_64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/93.0.4577.63 Safari/537.36',
+ 'user_agent' => (string)ini_get('user_agent')?:'rapsys_air/2.0.0',
+ ]
+ ]
+ );
+
+ //Get tile xy
+ $tileX = floor($centerX = $this->longitudeToX($longitude, $zoom));
+ $tileY = floor($centerY = $this->latitudeToY($latitude, $zoom));
+
+ //Calculate start xy
+ $startX = floor(($tileX * self::tz - $width) / self::tz);
+ $startY = floor(($tileY * self::tz - $height) / self::tz);
+
+ //Calculate end xy
+ $endX = ceil(($tileX * self::tz + $width) / self::tz);
+ $endY = ceil(($tileY * self::tz + $height) / self::tz);
+
+ for($x = $startX; $x <= $endX; $x++) {
+ for($y = $startY; $y <= $endY; $y++) {
+ //Set cache path
+ $cache = $this->cache.'/'.$zoom.'/'.$x.'/'.$y.'.png';
+
+ //Without cache path
+ if (!is_dir($dir = dirname($cache))) {
+ //Create filesystem object
+ $filesystem = new Filesystem();
+
+ try {
+ //Create path
+ //XXX: set as 0775, symfony umask (0022) will reduce rights (0755)
+ $filesystem->mkdir($dir, 0775);
+ } catch (IOExceptionInterface $e) {
+ //Throw error
+ throw new \Exception(sprintf('Output directory "%s" do not exists and unable to create it', $dir), 0, $e);
+ }
+ }
+
+ //Without cache image
+ if (!is_file($cache)) {
+ //Set tile url
+ $tileUri = str_replace(['{Z}', '{X}', '{Y}'], [$zoom, $x, $y], $this->servers[$this->server]);
+
+ //Store tile in cache
+ file_put_contents($cache, file_get_contents($tileUri, false, $ctx));
+ }
+
+ //Set dest x
+ $destX = intval(floor(($width / 2) - self::tz * ($centerX - $x)));
+
+ //Set dest y
+ $destY = intval(floor(($height / 2) - self::tz * ($centerY - $y)));
+
+ //Read tile from cache
+ $tile->readImage($cache);
+
+ //Compose image
+ $image->compositeImage($tile, \Imagick::COMPOSITE_OVER, $destX, $destY);
+
+ //Clear tile
+ $tile->clear();
+ }
+ }
+
+ //Add circle
+ //XXX: see https://www.php.net/manual/fr/imagick.examples-1.php#example-3916
+ $draw = new \ImagickDraw();
+
+ //Set text antialias
+ $draw->setTextAntialias(true);
+
+ //Set fill color
+ $draw->setFillColor('#cff');
+
+ //Set stroke color
+ $draw->setStrokeColor('#00c3f9');
+
+ //Set stroke width
+ $draw->setStrokeWidth(2);
+
+ //Draw circle
+ $draw->circle($width/2 - 5, $height/2 - 5, $width/2 + 5, $height/2 + 5);
+ #$draw->circle($tileX/self::tz - 5, $tileY/self::tz - 5, $tileX/self::tz + 5, $tileY/self::tz + 5);
+
+ //Draw on image
+ $image->drawImage($draw);
+
+ //Strip image exif data and properties
+ $image->stripImage();
+
+ //Add latitude
+ //XXX: not supported by imagick :'(
+ $image->setImageProperty('exif:GPSLatitude', $this->latitudeToSexagesimal($latitude));
+
+ //Add longitude
+ //XXX: not supported by imagick :'(
+ $image->setImageProperty('exif:GPSLongitude', $this->longitudeToSexagesimal($longitude));
+
+ //Add description
+ //XXX: not supported by imagick :'(
+ $image->setImageProperty('exif:Description', $alt);
+
+ //Save image
+ if (!$image->writeImage($path)) {
+ //Throw error
+ throw new \Exception(sprintf('Unable to write image "%s"', $dest));
+ }
+
+ //Get dest stat
+ $stat = stat($path);
+
+ //Return image data
+ return [
+ 'src' => $this->url.'/'.$stat['mtime'].$pathInfo.'.jpeg',
+ 'alt' => $alt,
+ 'height' => $height,
+ 'width' => $width
+ ];
+ }
+
+ /**
+ * Convert longitude to tile x number
+ *
+ * @see https://wiki.openstreetmap.org/wiki/Slippy_map_tilenames#Lon..2Flat._to_tile_numbers_5
+ *
+ * @param float $longitude The longitude
+ * @param int $zoom The zoom
+ *
+ * @return float The tile x
+ */
+ public function longitudeToX(float $longitude, int $zoom): float {
+ return (($longitude + 180) / 360) * pow(2, $zoom);
+ }
+
+ /**
+ * Convert latitude to tile y number
+ *
+ * @see https://wiki.openstreetmap.org/wiki/Slippy_map_tilenames#Lon..2Flat._to_tile_numbers_5
+ *
+ * @param $latitude The latitude
+ * @param $zoom The zoom
+ *
+ * @return float The tile y
+ */
+ public function latitudeToY(float $latitude, int $zoom): float {
+ return (1 - log(tan(deg2rad($latitude)) + 1 / cos(deg2rad($latitude))) / pi()) / 2 * pow(2, $zoom);
+ }
+
+ /**
+ * Convert tile x to longitude
+ *
+ * @param float $x The tile x
+ * @param int $zoom The zoom
+ *
+ * @return float The longitude
+ */
+ public function xToLongitude(float $x, int $zoom): float {
+ return $x / pow(2, $zoom) * 360.0 - 180.0;
+ }
+
+ /**
+ * Convert tile y to latitude
+ *
+ * @param float $y The tile y
+ * @param int $zoom The zoom
+ *
+ * @return float The latitude
+ */
+ public function yToLatitude(float $y, int $zoom): float {
+ return rad2deg(atan(sinh(pi() * (1 - 2 * $y / pow(2, $zoom)))));
+ }
+
+ /**
+ * Convert decimal latitude to sexagesimal
+ *
+ * @param float $latitude The decimal latitude
+ *
+ * @return string The sexagesimal longitude
+ */
+ public function latitudeToSexagesimal(float $latitude): string {
+ //Set degree
+ $degree = $latitude % 60;
+
+ //Set minute
+ $minute = ($latitude - $degree) * 60 % 60;
+
+ //Set second
+ $second = ($latitude - $degree - $minute / 60) * 3600 % 3600;
+
+ //Return sexagesimal longitude
+ return $degree.'°'.$minute.'\''.$second.'"'.($latitude >= 0 ? 'N' : 'S');
+ }
+
+ /**
+ * Convert decimal longitude to sexagesimal
+ *
+ * @param float $longitude The decimal longitude
+ *
+ * @return string The sexagesimal longitude
+ */
+ public function longitudeToSexagesimal(float $longitude): string {
+ //Set degree
+ $degree = $longitude % 60;
+
+ //Set minute
+ $minute = ($longitude - $degree) * 60 % 60;
+
+ //Set second
+ $second = ($longitude - $degree - $minute / 60) * 3600 % 3600;
+
+ //Return sexagesimal longitude
+ return $degree.'°'.$minute.'\''.$second.'"'.($longitude >= 0 ? 'E' : 'W');
+ }
+}
Rapsys Git / packbundle / commitdiff