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

文章目录


前言

受制于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

猜你喜欢

转载自blog.csdn.net/smartdrvsky/article/details/108752970
今日推荐