zoukankan      html  css  js  c++  java
  • ArcGIS Engine栅格数据使用总结

    jojojojo2002 原文 ArcGIS Engine栅格数据使用总结

    简介:ArcGIS Engine栅格数据使用总结,一个栅格数据集由一个或者多个波段(RasterBand)的数据组成,一个波段就是一个数据矩阵。对于格网数据(DEM数据)和单波段的影像数据,表现为仅仅只有一个波段数据的栅格数据集,而对于多光谱影像数据则表现为具有多个波段的栅格数据集。栅格编目(RasterCatalog)用于显示某个研究区域内各种相邻的栅格数据,这些相邻的栅格数据没有经过拼接处理合成一副大的影像图…
     
    最近做栅格数据,对AO的栅格部分应该有了一定的理解,希望高手指教:
     
    1、栅格数据的存储类型
     
    栅格数据一般可以存储为ESRI GRID(由一系列文件组成),TIFF格式(包括一个TIF文件和一个AUX文件),IMAGINE Image格式 在AE中一般调用ISaveAs接口来保存栅格数据
     
    2、栅格数据集和栅格编目的区别
    一个栅格数据集由一个或者多个波段(RasterBand)的数据组成,一个波段就是一个数据矩阵。对于格网数据(DEM数据)和单波段的影像数据,表现为仅仅只有一个波段数据的栅格数据集,而对于多光谱影像数据则表现为具有多个波段的栅格数据集
    栅格编目(RasterCatalog)用于显示某个研究区域内各种相邻的栅格数据,这些相邻的栅格数据没有经过拼接处理合成一副大的影像图
     
    3、IRasterWorkspaceEx与IRasterWorkspace ,IRsterWorkspace2的区别
    1).IRasteWorkspaceEx接口主要是用来读取GeoDatabase中的栅格数据集和栅格编目
    2) . IRasterWorkspace ,IRsterWorkspace2主要是用来读取以文件格式存储在本地的栅格数据
     
    4、加载栅格数据
    以存储在本地的栅格数据文件为例
    1.直接用IRasterLayer接口打开一个栅格文件并加载到地图控件
    IRasterLayer rasterLayer = new RasterLayerClass();
    rasterLayer.CreateFromFilePath(fileName); // fileName指存本地的栅格文件路径
    axMapControl1.AddLayer(rasterLayer, 0);
    2. 用IRasterDataset接口打开一个栅格数据集
    IWorkspaceFactory workspaceFactory = new RasterWorkspaceFactory();
    IWorkspace workspace;
    workspace = workspaceFactory.OpenFromFile(inPath, 0); //inPath栅格数据存储路径
    if (workspace == null)
    {
    Console.WriteLine(“Could not open the workspace.”);
    return;
    }
    IRasterWorkspace rastWork = (IRasterWorkspace)workspace;
    IRasterDataset rastDataset;
    rastDataset= rastWork.OpenRasterDataset(inName);//inName栅格文件名
    if (rastDataset == null)
    {
    Console.WriteLine(“Could not open the raster dataset.”);
    return;
    }
     
    5、如何读取栅格数据的属性和遍历栅格数据
    栅格数据的属性包括栅格大小,行数,列数,投影信息,栅格范围等等,见下面代码
    (假设当前加载的栅格文件栅格值存储方式为:UShort类型)
    IRasterProps rasterProps = (IRasterProps)clipRaster;
    int dHeight = rasterProps.Height;//当前栅格数据集的行数
    int dWidth = rasterProps.Width; //当前栅格数据集的列数
    double dX = rasterProps.MeanCellSize().X; //栅格的宽度
    double dY = rasterProps.MeanCellSize().Y; //栅格的高度
    IEnvelope extent=rasterProps.Extent; //当前栅格数据集的范围
    rstPixelType pixelType=rasterProps.PixelType; //当前栅格像素类型
    IPnt pntSize = new PntClass();
    pntSize.SetCoords(dX, dY);
    IPixelBlock pixelBlock = clipRaster.CreatePixelBlock(pntSize);
    IPnt pnt = new PntClass();
    for (int i = 0; i < dHeight; i++)
    for (int j = 0; j < dWidth; j++)
    {
    pnt.SetCoords(i, j);
    clipRaster.Read(pnt, pixelBlock);
    if (pixelBlock != null)
    {
    object bj = pixelBlock.GetVal(0, 0, 0);
    MessageBox.Show( Convert.ToUInt32(obj).ToString());
    }
    }
     
    6、如何提取指定的范围的栅格数据
    提取指定范围内的栅格数据通常用两种方法IRasterLayerExport(esriCarto), IExtractionOp, IExtractionOp2 (esriSpatialAnalyst),IRasterLayerExport接口提供的栅格数据提取功能有限,只能以矩形范围作为提取范围,而IExtractionOp接口提供了多边形,圆,属性,矩形等几种形式作为提取栅格数据.
    1).IRasterLayerExport接口
    IRasterLayerExport rLayerExport = new RasterLayerExportClass();
    rLayerExport.RasterLayer = rasterLayer;// rasterLayer指当前加载的栅格图层
    rLayerExport.Extent = clipExtent;//clipExtent指提取栅格数据的范围
    if (proSpatialRef != null)
    rLayerExport.SpatialReference = proSpatialRef;// proSpatialRef当前栅格数据的投影信息
    IWorkspaceFactory pWF = new RasterWorkspaceFactoryClass();
    try
    {
    IWorkspace pRasterWorkspace = pWF.OpenFromFile(_folder, 0);// _folder指栅格文件保存路径
    IRasterDataset utGeoDataset = rLayerExport.Export(pRasterWorkspace, code, strRasterType);
    //调用ISaveAs接口将导出的数据集保存
    ……………………..
    }
    Catch(Exception ex)
    {
    Throw new Argumention(ex.Message);
    }
    2.IExtractionOp接口(调用此接口前,应该先检查空间许可)
    IExtractionOp extraction = new RasterExtractionOpClass();
    try
    {
    IGeoDataset geoDataset = extraction.Rectangle((IGeoDataset)clipRaster, clipExtent, true);
    IRaster raster = geoDataset as IRaster;
    if (raster != null)
    {
    IWorkspaceFactory WF = new RasterWorkspaceFactoryClass();
    IWorkspace rasterWorkspace = WF.OpenFromFile(_folder, 0);
    ISaveAs saveAs = (ISaveAs)raster;
    saveAs.SaveAs(“Result.tif”, rasterWorkspace, “TIFF”);
    }
    }
    catch (Exception ex)
    {
    MessageBox..Show(Ex.message);
    }
     
    7.栅格数据重采样
    栅格数据的重采样主要基于三种方法:最邻近采样(NEAREST),双线性
    ILINEAR)和三次卷积采样(CUBIC)。
    (1).最邻近采样:它用输入栅格数据中最临近栅格值作为输出值。因此,在重采
    样后的输出栅格中的每个栅格值, 都是输入栅格数据中真实存在而未加任何改变的值。这种方法简单易用,计算量小,重采样的速度最快。
    (2).双线性采样:此重采样法取待采样点(x,y)点周围四个邻点,在y方向(或X方向)内插两次,再在x方向(或y方向)内插一次,得到(x,y)点的栅格值。
    (3).三次卷积采样:这是进一步提高内插精度的一种方法。它的基本思想是增加邻点来获
    得最佳插值函数。取待计算点周围相邻的16个点,与双线性采样类似,可先在某一方向上内插,如先在x方向上,每四个值依次内插四次,再根据四次的计算结果在y方上内插,最终得到内插结果
    代码示例:采用双线性采样
    IRasterGeometryProc rasterGeometryProc = new RasterGeometryProcClass();
    rasterGeometryProc.Resample(rstResamplingTypes.RSP_CubicConvolution, newCellSize, clipRaster);
    MessageBox..Show(Ex.message);
    }
     
  • 相关阅读:
    嗨,我在京东发现个好东东,还挺便宜
    电车男_百度百科
    罗生门 百度百科
    丹朱AI的日记豆瓣社区(手机版)
    关于美国Amazon的工作情况 JiansNet
    仓廪实而知礼节,衣食足而知荣辱_百度百科
    分享:对linux的IO的再认识
    怎样征服美丽少女的评论 (4)
    代入感 百度百科
    高富帅们 ucla的衣服怎么样啊
  • 原文地址:https://www.cnblogs.com/arxive/p/6017392.html
Copyright © 2011-2022 走看看