图像的最短路径求取（有一点问题和局限性）

最新推荐文章于 2022-10-26 14:39:38 发布
最新推荐文章于 2022-10-26 14:39:38 发布
阅读量979 收藏
点赞数
分类专栏： 图像处理 文章标签： fp delete 存储 include file c
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/zhoulingjie311/article/details/6222012
版权

#include <math.h>
#include <iomanip.h>
#include <stdlib.h>
#include <windows.h>
#include <stdio.h>
#include <stdlib.h>
#include <iostream.h>
#include <fstream.h>
#include "string.h"
//画出任意坐标点到终点最短路径的路线
void roadPath(int i,int j,int m_lineByte,unsigned char *m_pBmpBuf);
//八个方向的最短路
void minDire(int i,int j);
// int bmp[1000][1000];
int N,M;
//char *bmp;
//---------------------------------------------------------------------------------------
//以下该模块是完成BMP图像(彩色图像是24bit RGB各8bit)的像素获取，并存在文件名为xiang_su_zhi.txt中
unsigned char *pBmpBuf;//读入图像数据的指针

int bmpWidth;//图像的宽
int bmpHeight;//图像的高
int lineByte;
RGBQUAD *pColorTable;//颜色表指针

int biBitCount;//图像类型，每像素位数
//-------------------------------------------------------------------------------------------
//读图像的位图数据、宽、高、颜色表及每像素位数等数据进内存，存放在相应的全局变量中

int **d,****road;//d[][]是坐标数据,road[N+2][M+2][3][3]是记录行走方向的数组
double **c,****juli;
int dire[8][2]={{1,0},{1,1},{1,-1},{-1,0},{-1,1},{-1,-1},{0,1},{0,-1}};//八个方向;

bool readBmp(char *bmpName)
{
    FILE *fp=fopen(bmpName,"rb");//二进制读方式打开指定的图像文件

    if(fp==0) return 0;

    //跳过位图文件头结构BITMAPFILEHEADER

    fseek(fp, sizeof(BITMAPFILEHEADER),0);

    //定义位图信息头结构变量，读取位图信息头进内存，存放在变量head中

    BITMAPINFOHEADER head; 

    fread(&head, sizeof(BITMAPINFOHEADER), 1,fp); //获取图像宽、高、每像素所占位数等信息

    bmpWidth = head.biWidth;

    bmpHeight = head.biHeight;

    biBitCount = head.biBitCount;//定义变量，计算图像每行像素所占的字节数（必须是4的倍数）

    lineByte=(bmpWidth * biBitCount/8+3)/4*4;//灰度图像有颜色表，且颜色表表项为256

    if(biBitCount==8)
 {

        //申请颜色表所需要的空间，读颜色表进内存

        pColorTable=new RGBQUAD[256];

        fread(pColorTable,sizeof(RGBQUAD),256,fp);

 }

    //申请位图数据所需要的空间，读位图数据进内存

    pBmpBuf=new unsigned char[lineByte * bmpHeight];
    fread(pBmpBuf,1,lineByte * bmpHeight,fp);

    fclose(fp);//关闭文件

    return 1;//读取文件成功
}
//-----------------------------------------------------------------------------------------
//给定一个图像位图数据、宽、高、颜色表指针及每像素所占的位数等信息,将其写到指定文件中
bool saveBmp(char *bmpName, unsigned char *imgBuf, int width, int height,

             int biBitCount, RGBQUAD *pColorTable)

{

    //如果位图数据指针为0，则没有数据传入，函数返回
 M=width;N=height;
 int i,j,k;
 //动态分配原始数组d[][],距离数组c[][]
 d=(int **)malloc((N+2)*sizeof(int*));
 c=(double **)malloc((N+2)*sizeof(double *));
 road=(int ****)malloc((N+2)*sizeof(int***));
 juli=(double ****)malloc((N+2)*sizeof(double***));
 for (i=0;i<=N+1;i++)
 {
  d[i]=(int *)malloc((M+2)*sizeof(int));
  c[i]=(double *)malloc((M+2)*sizeof(double));
        road[i]=(int ***)malloc((M+2)*sizeof(int**));
  juli[i]=(double ***)malloc((M+2)*sizeof(double**));
  for (j=0;j<=M+1;j++)
  {
   d[i][j]=0;
   c[i][j]=100000; //初始化
            road[i][j]=(int **)malloc((3)*sizeof(int*));
   juli[i][j]=(double **)malloc((3)*sizeof(double*));
   for (int i1=0;i1<=2;i1++)
   {
                road[i][j][i1]=(int *)malloc((3)*sizeof(int));
    juli[i][j][i1]=(double *)malloc((3)*sizeof(double));
    for(int j1=0;j1<=2;j1++)
    {
     road[i][j][i1][j1]=0;
     //juli[i][j][0][0]代表（-1，-1）这个方向，juli[i][j][1][2]代表（0，1）这个方向，其余类推
     juli[i][j][i1][j1]=0;
    }
   } 
  }
 }

    if(!imgBuf)

        return 0;

    //颜色表大小，以字节为单位，灰度图像颜色表为1024字节，彩色图像颜色表大小为0

    int colorTablesize=0;

    if(biBitCount==8)

        colorTablesize=1024;

    //待存储图像数据每行字节数为4的倍数

    int lineByte=(width * biBitCount/8+3)/4*4;

    //以二进制写的方式打开文件

    FILE *fp=fopen(bmpName,"wb");

    if(fp==0) return 0;

    //申请位图文件头结构变量，填写文件头信息

    BITMAPFILEHEADER fileHead;

    fileHead.bfType = 0x4D42;//bmp类型

    //bfSize是图像文件4个组成部分之和

    fileHead.bfSize= sizeof(BITMAPFILEHEADER) + sizeof(BITMAPINFOHEADER)

        + colorTablesize + lineByte*height;

    fileHead.bfReserved1 = 0;

    fileHead.bfReserved2 = 0;

    //bfOffBits是图像文件前3个部分所需空间之和

    fileHead.bfOffBits=54+colorTablesize;

    //写文件头进文件

    fwrite(&fileHead, sizeof(BITMAPFILEHEADER),1, fp);

    //申请位图信息头结构变量，填写信息头信息

    BITMAPINFOHEADER head;

    head.biBitCount=biBitCount;

    head.biClrImportant=0;

    head.biClrUsed=0;

    head.biCompression=0;

    head.biHeight=height;

    head.biPlanes=1;

    head.biSize=40;

    head.biSizeImage=lineByte*height;

    head.biWidth=width;

    head.biXPelsPerMeter=0;

    head.biYPelsPerMeter=0;

    //写位图信息头进内存

    fwrite(&head, sizeof(BITMAPINFOHEADER),1, fp);

    //如果灰度图像，有颜色表，写入文件

    if(biBitCount==8)

        fwrite(pColorTable, sizeof(RGBQUAD),256, fp);
  
    //写位图数据进文件
 int count=0;
     if(biBitCount==24){//彩色图像

        for(int i=1;i<=N;i++)
  {
          
            for(int j=1;j<=M;j++)
   {

     if (*(imgBuf+(i-1)*lineByte+(j-1)*3)==255&&(*(imgBuf+(i-1)*lineByte+(j-1)*3+1)==255)&&(*(imgBuf+(i-1)*lineByte+(j-1)*3+2)==255))
     {
     // *(imgBuf+(i-1)*lineByte+(j-1)*3)=0;
     // *(imgBuf+(i-1)*lineByte+(j-1)*3+1)=0;
      d[N+1-i][j]=1;
     }
     else
     {
      d[N+1-i][j]=0;
     }

   }
  }
  }
  for (i=1;i<=N;i++)
  {
   for (j=1;j<=M;j++)
   {
    for (k=0;k<8;k++)
    {
     //为每个坐标赋各个方向距离值，只有2点都为1的时候才有距离
     if (d[i][j]!=0)
     {
      //x坐标变大，j变大，y坐标变大，i变小！
      if (d[i-dire[k][1]][j+dire[k][0]]!=0)//保证dire[k][0]+1，dire[k][1]+1非负
      {
       juli[i][j][dire[k][0]+1][dire[k][1]+1]=sqrt(dire[k][0]*dire[k][0]+dire[k][1]*dire[k][1]);
      }
     
     }
    }
   }
  }
  //我的动态规划采取波浪外推得方式，就是令终点周围的点先获取最短路，然后再类似往外层扩散
  i=N;j=M;
  int t=0,m_t;
  if (N<M) m_t=N;
  else m_t=M;  
  //开始按照正方形往外推
  while(t<m_t)
  {
   i=N-t;
   for (j=M;j>=M-t;j--)
   {
    if (i==N&&j==M){
     c[i][j]=0;}
    else
    {
     minDire(i,j);
    }
   }
   j=M-t;
   for (i=N;i>=N-t;i--)
   {
    if (i==N&&j==M){
     c[i][j]=0;}
    else
    {
     minDire(i,j);
    }
   }
  
   t++;
  }
  //然后一排一排的往外推
  if (N<M)
  {
   for(j=M-N;j>=1;j--)
    for (i=N;i>=1;i--)
    {
     minDire(i,j);
    }
   
  }
  else
  {
   for(i=N-M;i>=1;i--)
    for (j=M;j>=1;j--)
    {
    minDire(i,j);
    }
  }
//  i=1;j=1;
 printf("请输入第几行，第几列：/n");
 scanf("%d%d",&i,&j);
 roadPath(i,j,lineByte,imgBuf); 

    fwrite(imgBuf, height*lineByte, 1, fp);
 fclose(fp);

    return 1;

}

//----------------------------------------------------------------------------------------
//以下为像素的读取函数
void xiang_su_du_qu()

{

    //读入指定BMP文件进内存

    char readPath[]="map.bmp";

    readBmp(readPath);
    //输出图像的信息

    cout<<"width="<<bmpWidth<<" height="<<bmpHeight<<" biBitCount="<<biBitCount<<endl;
   
    //循环变量，图像的坐标

    //每行字节数

    int lineByte=(bmpWidth*biBitCount/8+3)/4*4;

    //循环变量，针对彩色图像，遍历每像素的三个分量

    int m=0,n=0,count_xiang_su=0;

    //将图像左下角1/4部分置成黑色
 
    ofstream outfile("图像像素.txt",ios::in|ios::trunc);

    if(biBitCount==24){//彩色图像

        for(int i=bmpHeight-1;i>=0;i--)
  {
            for(int j=0;j<bmpWidth;j++)
   {
                for(int k=0;k<3;k++)//每像素RGB三个分量分别置0才变成黑色

    {
                    m=*(pBmpBuf+i*lineByte+j*3+k);
     outfile<<m<<" ";
     count_xiang_su++;
     if(count_xiang_su%bmpWidth==0)
     {
      outfile<<endl;
     }
    }
              n++;
            }

        }
  cout<<"总的像素个素为:"<<n<<endl;
  cout<<"----------------------------------------------------"<<endl;

    }
   
    //将图像数据存盘
   
    char writePath[]="nvcpy.BMP";//图片处理后再存储

    saveBmp(writePath, pBmpBuf, bmpWidth, bmpHeight, biBitCount, pColorTable);

    //清除缓冲区，pBmpBuf和pColorTable是全局变量，在文件读入时申请的空间
 delete []pBmpBuf;
    if(biBitCount==8)
        delete []pColorTable;

}


//画出任意坐标点到终点最短路径的路线
void roadPath(int i,int j,int m_lineByte,unsigned char *m_pBmpBuf)
{
 int flag=0,i1,j1;
 for (i1=0;i1<=2;i1++)
 {
  for(j1=0;j1<=2;j1++)
  {
   if (road[i][j][i1][j1]==1)
   {
    *(m_pBmpBuf+(N-1-(i-1))*m_lineByte+(j-1)*3+2)=255;
    *(m_pBmpBuf+(N-1-(i-1))*m_lineByte+(j-1)*3+1)=0;
    *(m_pBmpBuf+(N-1-(i-1))*m_lineByte+(j-1)*3)=0;

    *(m_pBmpBuf+(N-1-(i-(j1-1)-1))*m_lineByte+((j+(i1-1)-1))*3+2)=255;
    *(m_pBmpBuf+(N-1-(i-(j1-1)-1))*m_lineByte+(j+(i1-1)-1)*3+1)=0;
    *(m_pBmpBuf+(N-1-(i-(j1-1)-1))*m_lineByte+(j+(i1-1)-1)*3+0)=0;

    if ((i-(j1-1))!=N||(j+(i1-1))!=M)
    {
     roadPath(i-(j1-1),j+(i1-1),m_lineByte,m_pBmpBuf);
    }
   }
  
  }
 }           

}

//八个方向的最短路
void minDire(int i,int j)
{
 double min=1000;
 for (int k=0;k<8;k++)
 {
  if (juli[i][j][dire[k][0]+1][dire[k][1]+1])
  {
   if( (c[i-dire[k][1]][j+dire[k][0]]+juli[i][j][dire[k][0]+1][dire[k][1]+1]) <min)
   {
    min=c[i-dire[k][1]][j+dire[k][0]]+juli[i][j][dire[k][0]+1][dire[k][1]+1];
    c[i][j]=min;
    //先对原来存储的清0
    for (int i1=0;i1<=2;i1++)
    {
     for(int j1=0;j1<=2;j1++)
     {
      road[i][j][i1][j1]=0;
     }
    }
    
    road[i][j][dire[k][0]+1][dire[k][1]+1]=1;
   }   
  }
 }
}

void main()
{
    xiang_su_du_qu();
}

小杰博客
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
遗传算法解决最短路径问题matlab程序加注释
12-06
然而，遗传算法也有其局限性，如可能需要较长的计算时间和较高的计算资源，且参数设置（如种群大小、交叉概率、变异概率等）对结果有很大影响，需要通过实验调整。 总的来说，这个MATLAB程序提供了一个使用遗传算法...
求坐标系(x,y)中各点间的最短路径，并使用matplotlib作图
weixin_43616995的博客
08-05 1843
由于新冠肺炎疫情影响，医用口罩的需求量日益增加。某中学班主任了解到班里 20 名学生的家庭所在小区每周均发放到家庭一定数量的免费医用口罩，于是班主任想到了一个口罩共享方案。通过联系各个家庭，了解到每个学生家庭每周拥有的口罩数量以及最低需求口罩数量，有些家庭是口罩数量不足，而有些家庭是有充足的口罩，如表 1 所示。假设这 20 个家庭每周拥有固定口罩总量及固定需求总量。 1.表 1 显示了班主任和各个学生的家庭地址地理坐标（单位 km）。班主任想针对这个方案到那些口罩充足的学生家里做一次方案说明的家访。
几种最短路径算法
mrzhuang
10-12 1934
Dijisktra 不能处理负边，时间复杂度O（n^2),堆优化O（nlogE） 基本思路：首先它是求单源点到其他边的最短路。方法就是创建一个集合来存放一些点，一开始放源点进去。然后每次取到集合外哪个点的距离最短，将其加入集合。（这个时候要更新，即通过该点再到其他点会不会比现在到其他点的距离短，会就更新） 在一开始输入图的信息后，肯定至少有一条边已经是最短边了，如果说能通过其他点再到目标点距离更短...
Dijkstra算法求解三维坐标下的最短路径（附代码）
qq_44776064的博客
04-02 5333
一、迪杰斯特拉算法 Dijkstra算法是由E.W.Dijkstra于1959年提出，又叫迪杰斯特拉算法，它应用了贪心算法模式，是目前公认的最好的求解最短路径的方法。算法解决的是有向图中单个源点到其他顶点的最短路径问题，其主要特点是每次迭代时选择的下一个顶点是标记点之外距离源点最近的顶点。该算法复杂度为n^2,但是对此可以考虑用堆这种数据结构进行优化，可以将复杂度降为O((m+n)logn)。 D...
实验2-动态规划编程题3. 所有点对的最短路径问题
最新发布
qq_44651930的博客
10-26 426
【输入形式】输入的第1行包含2个整数n和m，表示图G中包含n个节点和m条边。接下来m行，每行中有3个整数i,j,w，表示从节点i到节点j存在一条边(节点编号从1开始)，该边的权重为w。输出最短路径矩阵，共包含n行，每行包含n个数字(数字之间使用空格分割)，表示该节点到其他节点的最短距离。特别地，节点到自身的距离定义为0；如果节点之间无法到达，使用1e9+7表示他们之间的距离。【问题描述】给定一个非负的加权有向图G，求其中任意两个节点之间的最短路径。【思考】如果出现负边，如何进行改进。
最短路径问题——Floyd算法详解（所有顶点之间的最短路径）
weixin_51472673的博客
01-25 5288
任意两点间最短路径 与求单源最短路径不同，可以利用Floyd算法求任意两点之间的最短路径，这是源点就从单个确定的点变成了所有的点，即求每一个点到各个点之间的最短路径。 当然，我们也可以用dijkstra算法对每个顶点都循环一次，也可以求得任意两点间最短路。 弗洛伊德（Floyd）算法介绍： 设顶点vi 到顶点vk 的最短路径已知为D[i,k]，顶点vk 到vj 的最短路径已知为D[k,j]，顶点vi 到vj 的路径为D[i,j]，则vi 到vj 的最短路径为：min((D[i,k],D[k,j]
A星搜索算法教程确定目标最短路径的A*搜索算法教程
01-03
A星(A*)搜索算法是一种广泛应用在...总之，A*搜索算法以其高效的性能和寻找最优路径的能力，在解决寻路问题时表现出色。通过理解其工作原理并选择合适的启发式函数，我们可以有效地应用A*算法解决各种路径规划挑战。
最短路径算法
11-30
"ShortestPathByBlueForece"文件提供的Java实现展示了如何用Brute Force方法解决最短路径问题，虽然它对于小规模问题是个可行的解决方案，但在处理大规模图时，应考虑使用其他更优化的算法，以提高性能和效率。...
人工智能专题报告-蚁群算法最短路径规划.rar
09-30
报告最后总结了蚁群算法在最短路径规划中的优势和局限性，并对未来的研究方向提出了展望，包括如何进一步提高算法效率、降低计算复杂度，以及将蚁群算法与其他智能算法融合，以解决更复杂的问题。 该压缩包中的源码...
The-ant-colony-algorithm.zip_routing algorithm_最短路径 蚁群_最短路径优化
07-13
在实际应用中，蚁群算法已被广泛应用于各种最短路径问题，如旅行商问题（Traveling Salesman Problem, TSP）、车辆路径问题（Vehicle Routing Problem, VRP）等。通过对这些问题的优化，蚁群算法可以提供高效、近似...
bmp数字图像处理，C++
06-16
8bit和24bit的BMP图像的一些基本处理，VC6编译
vs2010 BMP图像处理
04-14
对BMP图像的简单的处理，有读写和像素的简单处理，适用于初学者。
C++ 读取bmp图片示例程序
09-01
功能: 读取bmp 代码示例: typedef struct _bmpsize { long bmpWidth;//位图宽，单位为象素LONG=long（有符号） long bmpHeight;//位图高，单位为象素LONG=long（有符号） unsigned short BitCount;//每个象素所占用的比特位WORD=unsigned short long BytesPerLine;//位图每行数据所占的字节数LONG=long（有符号） bool operator ==(_bmpsize bs) { if(bmpWidth==bs.bmpWidth && bmpHeight==bs.bmpHeight && BitCount==bs.BitCount /*&&BytesPerLine==bs.BytesPerLine*/) return true; else return false; } bool operator !=(_bmpsize bs) { return !(*this==bs); } }BMPSIZE, *LPBMPSIZE;//位图大小信息 class scBmp { public: scBmp(); virtual ~scBmp(); BITMAPFILEHEADER GetBmpFileHeader() const {return m_bmfHeader;}; const char * GetFileName() const {return m_pFileName.GetDataPointer();} BITMAPINFO * GetBmpInfo() const {return (BITMAPINFO *)m_pbi.GetDataPointer();} BMPSIZE GetBmpSize() const {return m_bmpSize;} // scFdmtlType * GetPixels() {return &m_pixels;} scBYTE * GetPixels() {return &m_pixels;}
打开BMP文件C++代码
06-14
用C++写的打开BMP文件,比网上的其他很多资源简单很多,很容易读懂
c++读取bmp图片详解
热门推荐
u013509299的专栏
12-01 2万+
先介绍几个会用到的函数。 1、FILE * fopen(const char * path,const char * mode)； path是字符串类型的bmp图片路径；mode读取方式，等下回用到"rb"，读写打开一个二进制文件，允许读写数据，文件必须存在。 2、int fseek(FILE *stream, long offset, int fromwhere); 函数设置文件
C++实现bmp文件的读写
ting
09-21 1781
C++实现BMP文件读写
图像处理学习笔记（一）：C++正确打开真彩色bmp格式图像文件
qq_42181309的博客
03-11 669
参考知识：https://www.runoob.com/cplusplus/cpp-files-streams.html（文件操作） https://www.cnblogs.com/wainiwann/p/7086844.html（BMP文件详解） 数据集上网找bmp真彩色的有很多。 编译环境：VS2017 最近觉得学信号与系统，数字信号处理学的太枯燥了，想写一些数字图像处理的算法。作为一...
纯C++实现24位bmp格式图片的读取和修饰
m0_62505136的博客
11-05 7812
问题：现有一张bmp图片，要求将它读取到程序中并实现灰度化、水平翻转、模糊、棕色滤镜四种效果，如下所示： 原图 棕色滤镜效果 灰度化效果 模糊效果（不太明显） 当我第一次接触到这个问题时，是无从下手的。但在查阅了不少资料之后，整整一天，我成功地只用C++实现了打开、修饰、保存bmp文件的功能！ 1.bmp文件的基本信息 (1).bmp文件的种类 打开Windows自带的画图软件，发现bmp的存储格式有好几种。 单色位图：只有黑白两种颜色，每个像素占1位...
用C++处理BMP图像
海月的博客
03-27 4148
用户界面 显示图像数据 计算直方图（8位图像） 仅仅使用一个通道：rgbBlue 蓝 灰度192 红 灰度64 调整图像亮度 图像反色 代码： #include<iostream> #include<windows.h> #include <fstream> using nam...
迪杰斯特拉算法计算最短路径.doc
05-30
2. **最短路径问题**：福格环游地球的问题被定义为一个最短路径问题，需要找到从起点到终点的旅行时间最短的路线。 3. **图的二维表示**：为了应用Dijkstra算法，需要将三维的交通网络图转化为二维图，这通常涉及图...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值