<?php
$polygons = LoadMIF("1.mif");

$image_sx = 400;
$image_sy = 400;
if(empty($scale))$scale    = 35000;
$my_long  = -63.10774861954596;
$my_lat   = 46.2899306519141;
$sx 	  = 2 * $scale;
$sy       = $scale;

$center   = getlocationcoords($my_lat, $my_long, $sx,$sy) ;

$min_x    = $center["x"] - ($image_sx / 2);
$min_y    = $center["y"] - ($image_sy / 2);

$im        = imagecreate($image_sx,$image_sy);

$land      = imagecolorallocate ($im, 0xF7,0xEF,0xDE);
$sea       = imagecolorallocate ($im, 0xB5,0xC7,0xD6);
$red       = imagecolorallocate ($im, 0xff,0x00,0x00);

imagefilledrectangle($im,0,0,$image_sx,$image_sy,$sea);

foreach($polygons as $poly)
{
	$converted_points = array();
	$points        = explode(" ",$poly["poly_string"]);
	$number_points = count($points);
	$i = 0;
	while($i<$number_points)
	{
		$lon  = $points[$i];
		$lat  = $points[$i+1];
		$pt = getlocationcoords($lat, $lon, $sx, $sy);
		$converted_points[] = $pt["x"] - $min_x;
		$converted_points[] = $pt["y"] - $min_y;
		$i+=2;
	}
	imagefilledpolygon($im,$converted_points,$number_points/2,$land);
}

$pt["x"] = $center["x"] - $min_x;
$pt["y"] = $center["y"] - $min_y;

imagefilledrectangle($im,$pt["x"]-2,$pt["y"]-2,$pt["x"]+2,$pt["y"]+2,$red);


header("Content-type: image/png");
imagepng($im);
imagedestroy($im);


function LoadMIF($file)
{
	error_reporting(0xFFFFFFFF);

	$hfile   = fopen("$file", "r");
			
	$polygons = array();
	
	$in_data = false;
		
	while (!feof($hfile))
	{
 		$line = strtoupper(fgets($hfile, 1024));

		if(substr($line,0,4)=="DATA")
		{
			$in_data=true;
		}
		else if($in_data)
		{
			if(substr($line,0,4)=="LINE")
			{
				$array = explode(" ",$line);
				$poly_info = array();
				$poly_info["min_long"] = $long    = trim($array[1]);
				$poly_info["min_lat"]  = $lat     = trim($array[2]);
				$poly_info["max_long"] = $long_to = trim($array[3]);
				$poly_info["max_lat"]  = $lat_to  = trim($array[4]);
				$poly_info["vector_type"] = 1;
				$poly_info["poly_count"]  = 2;
				$poly_info["poly_string"] = "$long $lat $long_to $lat_to";
			}
			else if(substr($line,0,5)=="PLINE")
			{
				$array = explode(" ",$line);
				$poly_info = GetPolyString($hfile,$array[1]);
				$poly_info["vector_type"] = 2;
			}
			else if(substr($line,0,6)=="REGION")
			{
				$count = substr($line,7);
				while($count>0)
				{
					$line = fgets($hfile, 1024);
					$poly_info = GetPolyString($hfile,$line);
					$poly_info["vector_type"] = 3;
				
					$polygons[] = $poly_info;
					unset($poly_info);
					$count-=1;
				}
			}

			if(isset($poly_info))
			{
				$polygons[] = $poly_info;
				unset($poly_info);
			}
		}
	}
   	fclose($hfile);
	return $polygons;
}
	
function GetPolyString($hfile,$poly_count)
{
	$ret_vector = array();
	$ret_vector["min_long"]    = 9999999;
	$ret_vector["min_lat"]     = 9999999;
	$ret_vector["max_long"]    = -9999999;
	$ret_vector["max_lat"]     = -9999999;
	$ret_vector["poly_string"] = "";
	$ret_vector["poly_count"] = $poly_count;

	for($i=0;$i<$ret_vector["poly_count"];$i++)
	{
		$line  = fgets($hfile, 1024);
		$array = explode(" ",$line);
		$long = $array[0];
		$lat  = $array[1];
					
		$ret_vector["min_long"]  = min($long,$ret_vector["min_long"]);
		$ret_vector["min_lat"]   = min($lat ,$ret_vector["min_lat"]);
		$ret_vector["max_long"]  = max($long,$ret_vector["max_long"]);
		$ret_vector["max_lat"]   = max($lat ,$ret_vector["max_lat"]);
		if(!empty($ret_vector["poly_string"]))$ret_vector["poly_string"] .= " ";
		$ret_vector["poly_string"] .= "$long $lat";
	}

	return $ret_vector;
}

function getlocationcoords($lat, $lon, $width, $height)
{  
    $x = (($lon + 180) * ($width / 360));
    $y = ((($lat * -1) + 90) * ($height / 180));
    return array("x"=>round($x),"y"=>round($y));
}
?>