遥感图像处理中,有时需要获取正射后影像的有效区域,背景值可以通过GDAL中的
函数获取。如果没有,可以通过设置来给定这个值。一般DOM背景无效值为0,DSM为-9999.00000。而海洋遥感数据背景无效值可能为-32767.00000。
首选,使用无效值将影像二值化,然后对二值数据进行边界追踪。其中最重要的是两点:1初始点的选取;2追踪规则(方向)。代码如下:
struct Point
{
double x,y;
};
template<class TT> vector<Point> GetImageArea(TT *pData, TT nodata, int width, int height)
{
//如果遥感影像过大,先采样,再返回
//首先将影像二值化
unsigned char *pBinaryData= new unsigned char[(width)*(height)];
memset(pBinaryData, sizeof(unsigned char)*(width)*(height), 255);//二值化时,255背景区域,0位有效区域
for (int i = 0; i < height; i++)
{
for (int j = 0; j < width; j++)
{
if (pData[i*width + j] != nodata)
pBin