GDAL获取影像信息与ENVI头文件信息

最新推荐文章于 2025-06-15 15:33:28 发布

默而识之者

最新推荐文章于 2025-06-15 15:33:28 发布

阅读量4.6k

点赞数

CC 4.0 BY-SA版权

文章标签： GDAL ENVI

本文链接：https://blog.youkuaiyun.com/u011609113/article/details/46932937

GDAL 专栏收录该内容

6 篇文章

订阅专栏

本文介绍如何使用GDAL库读取并解析遥感图像的基本信息，并通过ENVI头文件获取额外的元数据。包括文件路径、驱动类型、分辨率、投影信息等关键属性。

GDAL读取图像后也会得到许多有用信息，需要显示出来。

同时，若有同名HDR文件，既ENVI的头文件，里面也包含有用信息，需要读取出来。

以下是读取函数与信息实例。

C# Code

        #region 根据GDAL获取信息
        private void Get_InforFromGDAL(string S_File)//根据GDAL获取信息，S_File为影像路径
        {
            Dataset ds = Gdal.Open(S_File, Access.GA_ReadOnly);
            DataManager_listBox.Items.Clear();//DataManager_listBox为一ListBox
            DataManager_listBox.Items.Add("文件: " + S_File);
            DataManager_listBox.Items.Add("驱动: " + ds.GetDriver().GetDescription());
            DataManager_listBox.Items.Add("波段数: " + ds.RasterCount);
            DataManager_listBox.Items.Add("X分辨率: " + ds.RasterXSize);
            DataManager_listBox.Items.Add("Y分辨率: " + ds.RasterYSize);
            string S_projectInfor = ds.GetProjectionRef();//投影信息里包含许多子信息（实例如下），需要逐个提取
            /*
            某个图像所包含的投影信息
            PROJCS["UTM_Zone_46N",GEOGCS["GCS_WGS_1984",DATUM["WGS_1984",SPHEROID["WGS_84",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",93.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1]]
            */
            DataManager_listBox.Items.Add("投影空间参考系: " + Get_InforFromGDAL_AttriInfor(S_projectInfor, "PROJCS"));
            DataManager_listBox.Items.Add("地理参考系: " + Get_InforFromGDAL_AttriInfor(S_projectInfor, "GEOGCS"));
            DataManager_listBox.Items.Add("水平基准面: " + Get_InforFromGDAL_AttriInfor(S_projectInfor, "DATUM"));
            DataManager_listBox.Items.Add("椭球体: " + Get_InforFromGDAL_AttriInfor(S_projectInfor, "SPHEROID"));
            DataManager_listBox.Items.Add("本初子午线: " + Get_InforFromGDAL_AttriInfor(S_projectInfor, "PRIMEM"));
            DataManager_listBox.Items.Add("投影方式: " + Get_InforFromGDAL_AttriInfor(S_projectInfor, "PROJECTION"));
            double[] adfGeoTransform = new double[6];
            ds.GetGeoTransform(adfGeoTransform);
            DataManager_listBox.Items.Add("起始坐标 = (" + adfGeoTransform[0]
                + "," + adfGeoTransform[3] + ")");
            DataManager_listBox.Items.Add("像素大小 = (" + adfGeoTransform[1]
                + "," + adfGeoTransform[5] + ")");
        }

        private string Get_InforFromGDAL_AttriInfor(string S_InPut, string S_AttriName)//获取输入值某个参数的信息
        {
            string S_OutPut;
            string S_AllInfor = Get_InforFromGDAL_AttriInfor_AttriSubInfor(S_InPut, S_AttriName);//获取该类下所有信息
            int I_FirstKeyWord = S_InPut.IndexOf(S_AttriName);//获取该属性在其中的位置
            if (S_AllInfor != "" && I_FirstKeyWord != -1)//如果找到该属性字段
            {
                if (S_AllInfor.IndexOf("[") != -1)//若有下级信息
                {
                    int I_FirstQuota = I_FirstKeyWord + S_AttriName.Length + 1;
                    int I_LastQuota = S_InPut.IndexOf("\"", I_FirstQuota + 1);
                    S_OutPut = S_InPut.Substring(I_FirstQuota, I_LastQuota - I_FirstQuota + 1);//提取其中的属性值
                }
                else
                {
                    S_OutPut = S_AllInfor;//否则整个信息段都为属性值
                }
            }
            else
                S_OutPut = "未知";
            return S_OutPut;
        }

        private string Get_InforFromGDAL_AttriInfor_AttriSubInfor(string S_InPut, string S_AttriName)//获取下级信息
        {
            string S_OutPut = "";
            int I_FirstKeyWord = S_InPut.IndexOf(S_AttriName);
            int I_FirstBracket = I_FirstKeyWord + S_AttriName.Length + 1;
            if (I_FirstKeyWord != -1)
            {
                int I_LastBracket = S_InPut.IndexOf("]", I_FirstKeyWord);
                S_OutPut = S_InPut.Substring(I_FirstBracket, I_LastBracket - I_FirstBracket);
            }
            return S_OutPut;
        }

        #endregion 根据GDAL获取信息

根据ENVI头文件获取信息

某图像头文件如下所示

C# Code

ENVI
description = {
  File Resize Result, x resize factor: 1.000000, y resize factor:1.000000.[Mon
  Nov 17 09:07:59 2014]}
samples = 10355
lines   = 4531
bands   = 24
header offset = 0
file type = ENVI Standard
data type = 12
interleave = bip
sensor type = Unknown
byte order = 0
map info = {UTM, 1.000, 1.000, 604694.250, 3986008.500, 7.5000000000e-001, 7.5000000000e-001, 46, North, WGS-84, units=Meters}

C# Code

#region 根据ENVI头文件获取信息

private string Get_InforFromHDR(string S_File, string S_Attri)
{
    string S_AttriValue = "未知";
    StreamReader sr = new StreamReader(S_File);
    while (sr.Peek() >= 0)//读取所有记录
    {
        string S_Value = sr.ReadLine();
        if (S_Value.Contains(S_Attri))
            S_AttriValue = S_Value.Substring(S_Attri.Length + 2);//删除空格和等号
    }
    return S_AttriValue;
}
#endregion  根据ENVI头文件获取信息