在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
————————————————
版权声明：本文为CSDN博主「smartdrvsky」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/smartdrvsky/article/details/108752970