给矩阵添加地理信息

最新推荐文章于 2023-08-01 15:06:49 发布
最新推荐文章于 2023-08-01 15:06:49 发布
阅读量948 收藏 9
点赞数 3
分类专栏： 关于工具和软件 文章标签： 矩阵 matlab 线性代数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_44223180/article/details/123835118
版权

使用MATLAB实现

im = cell2mat(struct2cell(load('im.mat')));
%a=cell2mat(struct2cell(load('im.mat')))
%imshow(a,[])
[Data, R]= geotiffread('origin_image.tif');  % origin_image.tif 为带有地理信息的图像
info=geotiffinfo('origin_image.tif');
geotiffwrite('new_image.tif',im, R, 'GeoKeyDirectoryTag', info.GeoTIFFTags.GeoKeyDirectoryTag);   % new_image.tif 为写入地理信息后的新图像，new_image_matrix为待写入地理信息图像的像素值矩阵

需要保证原tif图和待写入的矩阵大小一样

](https://img-blog.csdnimg.cn/c04dc897c256466ead8938f62abaea5a.png)

否则会报错： 错误使用 geotiffwrite (line 257) R.RasterSize is inconsistent with rasterSize

不吃水果的太空人
关注
  • 3
    点赞
  • 踩
  • 9
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
遥感图像格式GeoTIFF解析
11-30
关于tiff文件及geotiff文件的格式的详细说明，并利用geolibtiff库来进行文件操作。
raster | R语言中的空间栅格对象及其基本处理方法（Ⅰ）
R语言学堂
07-19 1万+
前面的系列推文已经完成了对R语言中的两个管理空间矢量数据的工具包（sf和sp）的介绍，以及空间自相关、空间插值等空间分析方法。这里小编再推出一个系列来介绍R语言中管理空间栅格数据的工具包：...
计算两个栅格的相关系数
自度度人 自觉觉他
12-30 8069
有时候我们需要计算两个栅格的相关系数，判断相关性，例如计算NDVI和降水的相关系数，NDVI和温度的相关系数。今天分享一下计算两个栅格相关系数的计算方法。 1 相关系数计算 相关系数的计算公式网上书上有计算公式，这里不再赘述。这里介绍一下Python的numpy库计算相关系数，使用np.corrcoef()函数，示例如下。 import numpy as np x1 = np.array([9.6,17.1,64.8,40.9,136.3,182.5,78.3,3.7,26.,0.4]) x2 = np.a
MATLAB中对tif格式栅格影像读取-保存
热门推荐
stone_tigerLI的博客
07-15 1万+
MATLAB中对tif格式栅格影像读取-保存 欢迎使用Markdown编辑器 你好！ 这是你第一次使用 Markdown编辑器 所展示的欢迎页。如果你想学习如何使用Markdown编辑器, 可以仔细阅读这篇文章，了解一下Markdown的基本语法知识。 新的改变 我们对Markdown编辑器进行了一些功能拓展与语法支持，除了标准的Markdown编辑器功能，我们增加了如下几点新功能，帮助你用它写博...
对两种传感器的栅格数据进行逐像元回归分析，得到的每个像元的转换方程，去处理待校正的栅格图像
shanyao1119的博客
04-13 3157
栅格数据，逐像元线性回归分析
使用matlab计算栅格影像相关性： RasterSize value is not consistent with the raster size vector
wjz8486的博客
11-12 3890
使用matlab计算栅格影像相关性报错错误使用 geotiffwrite (line 257) A map.rasterref.MapCellsReference object was provided to function GEOTIFFWRITE, but its RasterSize value is not consistent with the raster size vector.首先点击geotiffwrite (line 257)查看处理内容错误错误定位结果输出 错误使用 geotiffw
matlab输出带坐标信息的大于4G的tiff影像
qq_33765083的博客
06-07 2019
matlab写入带坐标信息的大于4G的tiff影像时出错，解决办法
matlab将数字矩阵保存为地理空间数据出错，显示下标索引必须为正整数类型或逻辑类型，解决
qq_35794085的博客
05-28 713
matlab将数字矩阵保存为地理空间数据出错，解决！！！ 错误如下： 显示：下标索引必须为正整数类型或逻辑类型（关键是下午6点之前都还能正常运行，吃饭回来之后就报错了）。我反复的检查了我的代码，我的矩阵大小（大小为3678×6067），没有发现错误，我从头到尾都调试了一遍！！晕。。。浪费了一个小时，，然后就尝试了一下，把经纬度直接写进去。代码贴在下面，就这样吧，记录一下，！！！！！ 1、报错的代码（之前一直都能用的，搞不懂我动到什么地方了，就一直报错） lon = [96.99102628183385,
MATLAB读入Landsat8 数据报错处理方法
julming的博客
01-10 1791
MATLAB读入Landsat8 数据提示错误“Multiple images exist in the file and their sizes are different.”。
matlab %,Matlab
weixin_34536454的博客
03-17 605
Matlab读取带有地理信息的图像目前需要读取遥感影像并对像元值进行重新计算后再输出图像。代码如下：[A, R]= geotiffread('pathname:\ filename') %得到文件的像元值矩阵A和地理信息RB= A.* 0.0001 %将A进行处理，得到所需要的值geotiffwr...
长三角空间权重矩阵，01矩阵 地理反距离矩阵 经济距离矩阵 经济地理嵌套矩阵，四种矩阵都有哦，资料都是本人制作，还有分地区矩阵
05-30
长三角空间权重矩阵，01矩阵 地理反距离矩阵 经济距离矩阵 经济地理嵌套矩阵，四种矩阵都有哦，资料都是本人制作，还有分地区矩阵，单独的浙江，单独的安徽，单独的苏沪。一共4*4=16个矩阵，工作量很大的哦，希望...
全国30省经济权重矩阵、地理距离权重矩阵、经济地理嵌套矩阵
04-17
的地理距离权重矩阵、经济权重矩阵、经济地理嵌套矩阵，自用，可在stata中复制黏 贴再保存为dta文件，用于空间计量模型的实证分析，空间杜宾模型实证分析命令可在我 另外一个帖子中下载，个人经验在空间杜宾模型的...
地级市 县域 省级空间权重矩阵 地理距离矩阵 01矩阵
05-16
省级空间01矩阵 地理距离 经济距离矩阵 278个地级市三种空间权重矩阵及其标准化 285个地级市空间权重矩阵（地理距离和反距离） 286个地级市01矩阵 2364个县域空间权重矩阵
地理矩阵、经济矩阵和地理经济嵌套矩阵
最新发布
04-17
地理矩阵、经济矩阵和地理经济嵌套矩阵
stata 空间计量模型中运用的地理距离矩阵
06-18
适合正在准备毕业论文或是相关学术研究的人员
R语言绘制12个相异节点优化的有向和无向网络连通图
临风暖阳的博客
08-01 133
笔者主要使用R语言绘制有向和无向的网络连通图
arcengine-栅格数据详解
GIS_Y_C_B的博客
08-18 1453
打开栅格数据工作空间 //Open a file geodatabase workspace as RasterWorkspace. static IRasterWorkspaceEx OpenFGDB(string FGDBPath) { //FGDBPath string example: c:\data\raster.gdb. IWorkspaceFactory2 workspaceFactory = new FileGDBWorkspaceFactoryClass(); return (IRaste
C++中类模板分文件编写出现问题的原因（涉及函数定义作用、编译过程等）
对的时间点
11-11 3047
在一切之前，我们需要首先了解C++的编译过程 C++编译链接全过程 C++的编译过程分大致为预编译、编译、链接。 预编译过程 在预编译过程中，会进行一些最基本的操作。 将会把#define宏定义进行替换。 #define Max 100 如程序中出现了上面的语句，则预编译结束后，程序中的所有 Max 都已经被替换成了100 执行条件编译：#ifdef，#ifndef，#else，#eli...
C++/MFC-自定义启动窗口
学前班的博客
05-23 2482
在C****App.cpp的InitInstance()函数中，有这样几行代码 CSliderCtlDlg dlg; m_pMainWnd = &dlg; INT_PTR nResponse = dlg.DoModal(); if (nResponse == IDOK) { // TODO: 在此放置处理何时用 // “确定”来关闭对话框的代码 } else if (nResponse =
如何用arcgis做gtwr时空地理加权
11-03
在ArcGIS中，可以使用GTWR（Geographically and Temporally Weighted Regression，时空地理加权回归）进行时空地理加权分析。GTWR是一种空间回归分析方法，它结合了地理和时间信息，并允许回归系数随着空间和时间位置的变化而变化。 下面是使用ArcGIS进行GTWR时空地理加权的步骤： 1. 准备数据：将所需的空间和时间数据导入ArcGIS。例如，导入包含自变量和因变量的点矢量数据。 2. 创建空间邻域：使用ArcGIS中的邻域工具（如"Generate Spatial Weights Matrix"）来创建空间邻域。该工具将根据指定的距离阈值创建一个空间邻域矩阵，以描述每个点与其邻近点之间的空间关系。 3. 创建时间邻域：将时间数据转换为连续的时间变量，在数据中添加一个时间字段。然后，使用时间窗口来定义时间邻域，以描述每个点与其邻近时间点之间的时间关系。 4. 设定权重函数：使用ArcGIS的"GTWR Setup"工具来设定权重函数。在该工具中，可以选择不同的权重函数（如Gaussian、Bisquare等），并设置带宽参数以控制权重函数的形状和范围。 5. 进行GTWR分析：使用ArcGIS中的"GTWR"工具来执行GTWR分析。该工具将使用所选的权重函数和带宽参数来计算每个点的加权回归系数。此步骤将生成一个空间和时间分布的回归系数图。 6. 结果分析：使用ArcGIS的图层属性表和地图显示功能，可以对GTWR结果进行进一步分析和可视化。可以显示加权回归系数的空间和时间分布图，以及其统计显著性和可信区间。 通过上述步骤，可以在ArcGIS中成功执行GTWR时空地理加权分析，从而更准确地理解空间和时间因素对回归模型的影响。

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值