ArcGIS紧凑型切片读取与应用3-紧凑型批量转分散型(附源码)

系列文章目录

ArcGIS紧凑型切片读取与应用1-解析(附源码)

ArcGIS紧凑型切片读取与应用2-webgis动态加载紧凑型切片(附源码)

ArcGIS紧凑型切片读取与应用3-紧凑型批量转分散型(附源码)

1.前言

上篇介绍了webgis动态加载解析紧凑型切片的例子,现在我们使用逆向思维实现紧凑型切片转分散型切片,在实际工作中很有用处,紧凑型切片易于拷贝,但读取只有部署到Arcgis Server才行。相比分散型切片很通用可以部署在类似Geoewebcache的地图缓存服务器或者也可以直接部署到web服务器下。

软件核心功能:

ArcGIS紧凑型切片读取与应用3-紧凑型批量转分散型(附源码)

(1)支持切片等级范围选择。

(2)支持切图范围的选择,有利于局部数据的更新。

(3)支持多线程解析,充分利用系统资源,加快解析速率。

(4)文件命名格式与Arcgis分散型切片相同。

2.核心代码解析

        1.首先要解析将输入的坐标转为切片对应x、y、z值

/// <summary>
/// 通过经纬度获取切片位置
/// </summary>
/// <param name="lat_deg">纬度</param>
/// <param name="lon_deg">经度</param>
/// <param name="zoom">切片等级</param>
private double[] ConvertTile(double lat_deg, double lon_deg, int zoom)
{
double lat_rad = (Math.PI / ) * lat_deg;
double n = Math.Pow(, zoom);
double xtile = Math.Floor((lon_deg + 180.0) / 360.0 * n);
double ytile = Math.Floor((1.0 - Math.Log(Math.Tan(lat_rad) + ( / Math.Cos(lat_rad))) / Math.PI) / 2.0 * n);
return new double[] { xtile, ytile };
}

       2.通过左上角坐标和右上角坐标,解析出指定切片等级下的所有紧凑型切片文件

/// <summary>
/// 获取某一切片等级下的文件对象
/// </summary>
/// <param name="level"></param>
/// <param name="xy">【x1,y1,x2,y2】</param>
/// <returns></returns>
private List<BundleModel> GetLevelBundle(int level, int[] xy)
{
List<BundleModel> bundleModelList = new List<BundleModel>();
int minx = ((xy[] + ) / ) * ;
int maxx = ((xy[] + ) / ) * ;
int miny = ((xy[] + ) / ) * ;
int maxy = ((xy[] + ) / ) * ;
int xcount = (maxx - minx) / + ;
int ycount = (maxy - miny) / + ;
for (int x = ; x < xcount; x++)
{
for (int y = ; y < ycount; y++)
{
BundleModel bm = new BundleModel();
bm.StartX = minx + (x) * ;
bm.StartY = miny + (y) * ;
var rGroup = Convert.ToInt32( * Convert.ToInt32(bm.StartX / ));
var cGroup = Convert.ToInt32( * Convert.ToInt32(bm.StartY / ));
var bundleBase = getBundlePath(textBox1.Text, level, rGroup, cGroup);
bm.Level = level;
bm.BundlxDire = bundleBase + ".bundlx";
bm.BundleDire = bundleBase + ".bundle";
bm.BundleName = Path.GetFileNameWithoutExtension(bm.BundleDire);
bundleModelList.Add(bm);
}
}
return bundleModelList;
}

    3.多线程切片实现,将文件平均分给不同线程

//将文件平均分给各个线程
int count = bundleModelList.Count() / threadcount;
int yu = bundleModelList.Count() % threadcount; if (count == )
{
for (int i = ; i < bundleModelList.Count; i++)
{
List<BundleModel> model = bundleModelList.Skip(i).Take().ToList();
System.Threading.ThreadPool.QueueUserWorkItem((state) =>
{
foreach (var item in model)
{
ToImg(item);
}
this.BeginInvoke(new Action(() =>
{ }));
}, model);
}
}
else
{
for (int i = ; i < threadcount; i++)
{
List<BundleModel> model = bundleModelList.Skip(i * count).Take(count).ToList();
if (i < yu)
{
model.AddRange(bundleModelList.Skip(threadcount * count + i).Take().ToList());
}
System.Threading.ThreadPool.QueueUserWorkItem((state) =>
{
foreach (var item in model)
{
ToImg(item);
}
this.BeginInvoke(new Action(() =>
{ }));
}, model);
}
}

    4.Arcgis散片文件路径格式的生成

string L="L"+ zeroPad(z, );
string C = "C" + zeroPad(x, ,);
string R = "R" + zeroPad(y, ,); //保存路径
string path = textBox2.Text+"\\"+L+"\\"+R+"\\"+C+ ".png"; //文件命名函数 private string zeroPad(int num, int len,int type=)
{
string str = num.ToString();
if (type==)
{
str = num.ToString("X");
}
while (str.Length < len)
{
str = "" + str;
}
return str;
}

解析结果展示

ArcGIS紧凑型切片读取与应用3-紧凑型批量转分散型(附源码)

3.结束

     开启多线程模式切片的速率比较满意,通过比较简单的代码理解了紧凑型切片的所有的细节,我们现在完全可以实现散片型装紧凑型的文件,有兴趣可以反推一下。所有的源代码已近上传到了GitHub,欢迎大家指教。

百度网盘链接:https://pan.baidu.com/s/1I-Bj3EQSN57pQHvKZ2hBUA   提取码:lliw

github项目地址:https://github.com/HuHongYong/TilerArcgisBundle

作者:ATtuing

出处:http://www.cnblogs.com/ATtuing

本文版权归作者和博客园共有,欢迎转载,但未经作者同意必须保留此段声明,且在文章页面明显位置给出原文链接。

上一篇:基于visual Studio2013解决C语言竞赛题之1047百马问题


下一篇:序列化serialize()与反序列化unserialize()的实例