<p>必须用 new T[width height] 申请一维连续内存模拟二维布局,通过 data[y width + x] 访问,配对 delete[];推荐 std::unique_ptr 封装防泄漏,禁用 vector 套娃和 malloc/free 混用。</p>

用 new 和 delete[] 手动管理二维矩阵内存
遥感图像动辄上万×万像素,单波段就超 1GB,不能靠栈分配或 std::vector<:vector>></:vector> 套娃——那会触发大量小块内存分配,缓存不友好,还容易爆 bad_alloc。必须用一维连续内存模拟二维布局。
正确做法是申请一块大数组,用行优先索引计算: data[y * width + x]。注意必须用 new T[width * height],而非 new T*[height] 再循环 new 每行——后者内存不连续,SIMD 和 GPU 传输都会失败。
- 释放时严格配对使用
delete[] data,不能漏掉[],否则只析构首元素 - 别用
malloc/free混搭new,C++ 对象(如含构造函数的类)必须用new - 遥感常用
uint16_t或float,声明指针时明确类型:float* data = new float[width * height];
用 std::unique_ptr 封装裸指针避免内存泄漏
手动 new/delete 在异常路径或提前返回时极易忘删。遥感处理流程长(读取→辐射校正→重采样→写入),中间任意一步抛异常,裸指针就成悬垂指针。
推荐用 std::unique_ptr<float std::default_delete>></float>,它能自动调用 delete[]:
立即学习“C++免费学习笔记(深入)”;
auto data = std::unique_ptr<float[]>(new float[width * height]);
更简洁写法(C++14 起):
auto data = std::make_unique<float[]>(width * height);
- 传参时用
data.get()获取原始指针给底层库(如 GDAL 的GDALRasterBand::RasterIO) - 不要把
unique_ptr转成裸指针后长期持有,所有权语义会丢失 - 若需共享所有权(如多线程并行处理不同区域),改用
std::shared_ptr,但注意原子引用计数开销
处理跨平台字节序和内存对齐问题
遥感数据常来自不同传感器(如 Landsat、Sentinel),TIFF/BSQ/BIL 格式自带字节序标记。直接用指针读二进制文件时,uint16_t* 强转可能在小端机器上错读高位字节。
- 用
htons()/ntohs()显式转换,或 GDAL 的GDALCopyWords()自动处理字节序 - 某些算法(如 FFT、卷积)要求 16 字节对齐,
new不保证对齐;改用std::aligned_alloc(C++17)或_aligned_malloc(Windows) - GDAL 默认按块读取(
GDALRasterBand::ReadBlock),比整图加载更省内存,指针只需指向当前块缓冲区
避免 std::vector 的隐式拷贝陷阱
有人图省事用 std::vector<float> data(width * height)</float>,再取 &data[0] 当指针——这看似安全,但一旦 vector 因插入扩容,原有指针立刻失效。遥感处理中常见「先读一部分做直方图统计,再全量处理」,中间 push_back 或 resize 就崩了。
-
vector只适合生命周期明确、无动态增长的场景;否则宁可用unique_ptr+size_t width, height成员变量封装 - 别依赖
vector的「连续性承诺」做跨函数指针传递,除非你能 100% 控制其容量不变 - 调试时加断点检查
data.data() == &data[0]是否恒成立,不成立说明已被移动或重新分配
真正麻烦的不是分配,而是图像处理中频繁的 ROI 提取、重采样、金字塔构建——这些操作都依赖指针算术的精确偏移,稍有越界就是静默数据污染,得靠 ASan 或 Valgrind 配合单元测试覆盖边界条件。


















