前言
受制于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
提取码:bb7g