在CAD中加载大影像的一种方法

 

文章目录

• 前言
• 一、开发环境
• 二、具体思路
  • 1.gdal-rasterIO
  • 2.objectarx-自定义实体
  • 3.坐标转换
  • 4.调用方法

 

---

前言

受制于AutoCAD中RasterImage对象的限制，原生CAD环境中加载大型文件的影像时会异常卡顿，本文演示如何在AutoCAD中结合gdal加载几十上百G的tif和WMS/TMS/WTMS发布的影像服务。

全球缩略图

 

 

缩放至帝都

---

 

一、开发工具和环境

本次采用objectarx2020(64位)+vs2017+gdal3.14环境

objectarx2020可在autodesk官网下载

gdal开发环境推荐下载地址

https://www.gisinternals.com/release.php

 

二、具体思路

1.gdal RaterIO

关于rasterIO,gdal大神李民录已经写得很详细，具体参考

https://blog.csdn.net/liminlu0314/article/details/7072224

raterIO可以打开WMS/WMTS/TMS/geogiff中指定位置的RGB数据

 

2.objectarx自定义实体

在cad自定义实体重载subworldDraw时根据当前屏幕范围用gdal中的rasterIO方法去请求获取RGB绘制信息块

在请求数据后根据指定长宽即可在CAD自定义实体中调用mode->rawGeometry()->image(....）方法绘制

CAD在View大小发生改变时即重新根据范围去请求刷新图片即可，最后结果是根据当前显示范围大小加载不同层级的影像，实际坐标与影像坐标1：1匹配。

 

3.坐标转换

CAD屏幕坐标获取方法

void ImageIOThread::getScreenCorners(AcGePoint2d& ptLeftLower, AcGePoint2d& ptRightUpper)
{
	CRect rc;
	CView* pCurView = acedGetAcadDwgView();
	pCurView->GetClientRect(&rc);
	CPoint ptTopLeft = rc.TopLeft();
	CPoint ptBottomRight = rc.BottomRight();
	acedDwgPoint pt1, pt2;
	acedCoordFromPixelToWorld(ptTopLeft, pt1);
	acedCoordFromPixelToWorld(ptBottomRight, pt2);
	ptLeftLower.x = pt1[0];
	ptLeftLower.y = pt2[1];
	ptRightUpper.x = pt2[0];
	ptRightUpper.y = pt1[1];
}

通过当前屏幕坐标去计算对应在Image中对应的行列号，这里需要注意的是WMS等服务采用的是EPSG3857即伪墨卡托托影，实际使用过程中需要将当前坐标转至经纬度后再将经纬度转换成 3857，具体转换方法可参考李民录blog,也可用proj库进行坐标转换，各地不同坐标系定义不同椭球参数或用4参数/7参数进行转换。

4.调用方法

调用示例：

链接：https://pan.baidu.com/s/17xG5Sh1rhndPqq5fqizZqg