隐藏

根据经纬度绘制坐标点相对位置(分别用php和html5实现)

发布:2020/9/24 11:07:40作者:管理员 来源:本站 浏览次数:954

本文主要讲实现 

1.用php实现根据坐标点经纬度绘制坐标点相对位置。

2.用html5实现根据坐标点经纬度绘制坐标点相对位置。


之前有一个需求,是要把一些坐标点的相对位置画到一张图片上,标注出点,点名,并且点之间两两连线生成一张图片。

最开始考虑用html5实现。 但实际上我们生成的图片要在后台放到word文档中去,所以最后使用的php在后台画图,然后保存到本地的方式 。


这几天正好在研究html5,就用html5的方式再实现了一遍。感觉html5画起来还要顺手一些。 

当然,php也好,html5也好,都只是工具,实现起来主要还是要有算法。 

php跟javascript使用起来真是很像啊。



需求分析:

1.图片不能无限大。有一个最大值maxsize。 我们这里设为300px。也就是生成的图片,坐标区域宽或高最大不超过300px。

2.所有的坐标点都要在我们的图片中。

3.坐标点之间的位置关系尽可能真实反映坐标情况。

4.坐标点。点名。连线显示在地图上。


由需求看出难点在第3条上。经纬度坐标怎么换算成我们可以使用来画图的坐标呢?

经过一番思考,确定了一个粗略的方法:


 大概有这么几个步骤:

首先确认一个前提:就是我们假设经度和纬度在地图上摊开来是1比1的关系,即是说经度1度的距离约等于纬度1度的距离。 当然事实情况并不是这样,我们这里只是求一个粗略的位置关系。


1.求出所有坐标点的最大经度差,最大纬度差。比较,将大的那一个的长度作为maxsize。

2.根据最大经度差最大维度差之间的比例,可以算出我们短的那一边的长度。 

3.我们求出一个Rate。这个Rate值表示单位度数在生成图片上所占的尺寸,生成一个类似比例尺的东西。

4.然后我们现在已经确定了我们的图片的长和宽。同时我们以最小经度,最小纬度为图片的x,y轴原点。然后再根据Rate,以及当前点的经度、纬度与最小经度,最小纬度的差,求出当前点在图片上绘制时候的坐标。

5.求出了绘图坐标之后,就可以使用php。或者html5的方法来进行绘图操作了。

ps:对于经度差或纬度差为0的情况下我也做了单独处理。


html5绘制经纬度坐标点位置代码:


	
  1. <!DOCTYPE html>
  2. <html>
  3. <head>
  4. <meta charset="utf-8">
  5. <title>用html5绘制坐标点相对位置</title>
  6. </head>
  7. <style>
  8. canvas{
  9. border: 1px solid #ccc;
  10. }
  11. </style>
  12. <script type="text/javascript">
  13. //绘图数据准备
  14. var points=[[],[],[],[]];
  15. points[0]['L'] = 110.83717794;
  16. points[0]['B'] = 13;
  17. points[0]['Name'] = 'pt1';
  18. points[1]['L'] = 113;
  19. points[1]['B'] = 24;
  20. points[1]['Name'] = 'pt2';
  21. points[2]['L'] = 105.794566000;
  22. points[2]['B'] = 20;
  23. points[2]['Name'] = 'pt3';
  24. points[3]['L'] = 120;
  25. points[3]['B'] = 17;
  26. points[3]['Name'] = 'pt4';
  27. function draw (points) {
  28. //坐标区域长或宽,最大为300px。判断标准:坐标经度差与纬度差,大的那个为300px。再根据经度差与纬度差的比例计算出短的一边有多少px。
  29. var MAXSIZE = 300;
  30. var maxL = points[0]['L'];
  31. var maxB=points[0]['B'];
  32. var minL = points[0]['L'];
  33. var minB=points[0]['B'];
  34. var value;
  35. for(var i=0,pointsLen=points.length;i<pointsLen;i++){
  36. value = points[i];
  37. maxL = maxL <value['L']?value['L'] :maxL;
  38. maxB = maxB <value['B']?value['B'] :maxB;
  39. minL = minL > value['L']?value['L'] :minL;
  40. minB = minB >value['B']?value['B'] :minB;
  41. }
  42. var diffL = maxL - minL;//经度差
  43. var diffB = maxB - minB;//纬度差
  44. var width,height,Rate,diff;
  45. //计算坐标区域height width;
  46. if(diffL == 0){
  47. width =MAXSIZE;
  48. height = MAXSIZE;
  49. Rate = MAXSIZE/parseFloat(diffB);
  50. }
  51. else if (diffB == 0) {
  52. width =MAXSIZE;
  53. height = MAXSIZE;
  54. Rate = MAXSIZE/parseFloat(diffL);
  55. }else if(diffL >= diffB){
  56. diff = diffL;
  57. width = MAXSIZE;
  58. Rate = MAXSIZE/parseFloat(diffL);//单位坐标的有多少个px值。
  59. height = diffB/diffL*MAXSIZE;
  60. }else {
  61. diff = diffB;
  62. height = MAXSIZE;
  63. Rate = MAXSIZE/parseFloat(diffB);//单位坐标的有多少个px值。
  64. width = diffL/diffB*MAXSIZE;
  65. }
  66. var ctx = document.getElementById("myCanvas").getContext("2d");
  67. ctx.translate(20,20);//原点往右下方分别移动20px,多出来的是用来防止名字,以及点上的圆点显示不下
  68. // 根据B,L计算像素位置。计算应该有px。
  69. for(var k=0,pointsLen=points.length;k<pointsLen;k++){
  70. value = points[k];
  71. if(diffL == 0){
  72. points[k]['Lpx'] =MAXSIZE/2;
  73. points[k]['Bpx'] = parseInt(height - (value['B'] - minB)*Rate);
  74. }
  75. else if (diffB == 0) {
  76. points[k]['Lpx'] =parseInt((value['L'] - minL)* Rate);
  77. points[k]['Bpx'] = MAXSIZE/2;
  78. } else {
  79. points[k]['Lpx'] = parseInt((value['L'] - minL)* Rate);
  80. points[k]['Bpx'] = parseInt(height - (value['B'] - minB)*Rate);
  81. }
  82. }
  83. ctx.font = "bold 14px Arial";
  84. ctx.textAlign = "left";
  85. ctx.textBaseline = "middle";
  86. for(var x=0,pointsLen=points.length;x<pointsLen;x++){
  87. value = points[x];
  88. i = x +1;
  89. ctx.strokeStyle='black';
  90. //两两相连的处理
  91. while (i<pointsLen) {
  92. ctx.beginPath();
  93. ctx.moveTo(value['Lpx'] , value['Bpx']);
  94. ctx.lineTo(points[i]['Lpx'], points[i]['Bpx']);
  95. ctx.stroke();
  96. i++;
  97. }
  98. //画点名
  99. ctx.beginPath();
  100. ctx.fillStyle ='black';
  101. ctx.fillText(value.Name, value['Lpx']+5, value['Bpx']);
  102. //画点
  103. ctx.fillStyle = 'rgb(29,143,254)';
  104. ctx.beginPath();
  105. ctx.arc(value['Lpx'] , value['Bpx'] , 5, 0, Math.PI*2, false);
  106. ctx.fill();
  107. }
  108. }
  109. </script>
  110. <body οnlοad="draw(points)">
  111. <canvas id="myCanvas" width="350" height="350"></canvas>
  112. </body>
  113. </html>


php绘制经纬度坐标点位置代码:代码看起来要复杂一点,可能是对php画图研究的还不够。。


	
  1. <?php
  2. //测试数据
  3. $point[0]['L'] = 110.83717794;
  4. $point[0]['B'] = 13;
  5. $point[0]['Name'] = 'pt1';
  6. $point[1]['L'] = 113;
  7. $point[1]['B'] = 24;
  8. $point[1]['Name'] = 'pt2';
  9. $point[2]['L'] = 105.794566000;
  10. $point[2]['B'] = 20;
  11. $point[2]['Name'] = 'pt3';
  12. $point[3]['L'] = 120;
  13. $point[3]['B'] = 17;
  14. $point[3]['Name'] = 'pt4';
  15. function draw($point)
  16. {
  17. //1.计算四个范围。
  18. $maxL = $point[0]['L'];
  19. $maxB=$point[0]['B'];
  20. $minL = $point[0]['L'];
  21. $minB=$point[0]['B'];
  22. foreach ($point as $key => $value) {
  23. $maxL = $maxL <$value['L']?$value['L'] :$maxL;
  24. $maxB = $maxB <$value['B']?$value['B'] :$maxB;
  25. $minL = $minL > $value['L']?$value['L'] :$minL;
  26. $minB = $minB >$value['B']?$value['B'] :$minB;
  27. }
  28. $MAXSIZE = 300; //坐标区域长或宽,最大为300px。判断标准:坐标经度差与纬度差,大的那个为300px。再根据经度差与纬度差的比例计算出短的一边有多少px。
  29. $diffL = $maxL - $minL;//经度差
  30. $diffB = $maxB - $minB;//纬度差
  31. //计算坐标区域$height $width;
  32. if($diffL == 0){
  33. $width =$MAXSIZE;
  34. $height = $MAXSIZE;
  35. $Rate = $MAXSIZE/(float)$diffB;
  36. }
  37. elseif ($diffB == 0) {
  38. $width =$MAXSIZE;
  39. $height = $MAXSIZE;
  40. $Rate = $MAXSIZE/(float)$diffL;
  41. }else if($diffL >= $diffB){
  42. $diff = $diffL;
  43. $width = $MAXSIZE;
  44. $Rate = $MAXSIZE/(float)$diffL;//单位坐标的有多少个px值。
  45. $height = $diffB/$diffL*$MAXSIZE;
  46. }else{
  47. $diff = $diffB;
  48. $height = $MAXSIZE;
  49. $Rate = $MAXSIZE/(float)$diffB;//单位坐标的有多少个px值。
  50. $width = $diffL/$diffB*$MAXSIZE;
  51. }
  52. $img_width = $width + 50;//多出来的50是用来防止基站名字,以及点上的圆点显示不下
  53. $img_height = $height + 30;//多出来的30是用来防止基站名字,以及点上的圆点显示不下
  54. $image = imagecreatetruecolor($img_width,$img_height);//生成一个黑色背景的图片。
  55. $back = imagecolorallocate($image, 255, 255, 255);//背景颜色,白色
  56. imagefilledrectangle ($image,0,0,$img_width ,$img_height ,$back);//设置背景,其实是用的白色填充矩形。
  57. $linecolor= imagecolorallocate($image,0,0,0);//连线的颜色
  58. $pointcolor = imagecolorallocate($image,29,143,254);//点的颜色
  59. $textcolor = imagecolorallocate($image,0,0,0);//字体颜色,黑色
  60. // 根据B,L计算像素位置。计算应该有px。多出来的6,和10,表示所有坐标都向右移动6px,向下移动10px,也是避免基站名和圆点显示不下。
  61. foreach ($point as $key => $value) {
  62. if($diffL == 0){
  63. $point[$key]['Lpx'] =$MAXSIZE/2;
  64. $point[$key]['Bpx'] = (int)($height - ($value['B'] - $minB)*$Rate)+10;
  65. }
  66. elseif ($diffB == 0) {
  67. $point[$key]['Lpx'] =(int)(($value['L'] - $minL)* $Rate) + 6;
  68. $point[$key]['Bpx'] = $MAXSIZE/2;
  69. }else{
  70. $point[$key]['Lpx'] = (int)(($value['L'] - $minL)* $Rate) + 6;
  71. $point[$key]['Bpx'] = (int)($height - ($value['B'] - $minB)*$Rate)+10;
  72. }
  73. }
  74. $pointNum = count($point);
  75. foreach ($point as $key => $value) {
  76. $i = $key +1;
  77. while ($i<$pointNum) {
  78. imageline($image , $value['Lpx'] , $value['Bpx'], $point[$i]['Lpx'], $point[$i]['Bpx'] , $linecolor );//点与点之间连线
  79. $i++;
  80. }
  81. imagestring( $image ,10 , $value['Lpx'], $value['Bpx'], $value['Name'] ,$linecolor); //每个点的位置写下点名。10是大小。
  82. imagefilledellipse($image , $value['Lpx'] , $value['Bpx'], 10 , 10, $pointcolor); //每个点的位置,用一个填充的圆点表示。
  83. }
  84. // 输出图像
  85. header("Content-type: image/jpeg");
  86. // imagejpeg($image,'d:\1.jpeg');//如果没有路径就输出图片,有路径就保存图片。
  87. imagejpeg($image);//如果没有路径就输出图片,有路径就保存图片。
  88. }
  89. draw($point);
  90. ?>


接下来,对比图!!!

html5绘制画出来的经纬度坐标位置关系图:

php绘制画出来的经纬度坐标位置关系图:

再接下来,我们到地图上看这几个坐标的真实位置:(数据来自之前用天地图WEB API写的一个坐标展绘的功能,可以把输入的坐标点在地图上显示出来)


如果对精确度要求不是特别高的话,我这个方法应该还是可以用的。    有什么不足的地方,欢迎指出探讨。


最后总结一句,其实会某个技能倒是次要的,技能用的时候可以学,更重要的还是要有想法来用这个技能。。。


演示地址:

html5:http://runningls.com/demos/2015/position/html5.html

php:http://runningls.com/demos/2015/position/php.php

github:https://github.com/liusaint/earth/tree/master/position