概述
今天在用Matlab做栅格数据偏相关分析的时候,遇到了如下问题:
参考此回答,可以了解reshape的用法
打开ArcGIS,查看两数据的行列数,果然不一样:
解决方法:
将列数多的影像按列数少的位边界提取,借助批量掩膜工具(工具箱来自ZHOU-LONG)
下载地址:
https://download.csdn.net/download/qq_33356563/11020148
再次检查,发现列数已经变为1031了。
问题解决。
不过批量转换出来的结果是.img文件,再批量转换为.tif文件
在Arcpy中运行批量转换代码即可(代码来自mygisforum):
# Import system modules
import sys, string, os
dir = 'F:\gzNDVI500'
# Import arcpy module
import arcpy
files = os.listdir(dir)
for f in files:
if os.path.splitext(f)[1] == '.img':
# Script arguments...
Input_raster_file = dir + os.sep + f
# Local variables...
Output_data_type = "FLOAT"
Raster_Format = "TIFF"
Output_Workspace = "f:\gzNDVI500tf"
# =============== file name process ======================
basename = os.path.splitext(f)[0];
Output_raster = Output_Workspace + os.sep + basename + ".tif";
if os.path.exists(Output_raster) == False:
print (Input_raster_file)
# Process: Raster To Other Format (multiple)...
arcpy.RasterToOtherFormat_conversion(Input_raster_file,
Output_Workspace, Raster_Format)
print(Output_raster)
最后
以上就是爱笑哑铃为你收集整理的解决Matlab中reshape报错:栅格影像行列数不一致的问题的全部内容,希望文章能够帮你解决解决Matlab中reshape报错:栅格影像行列数不一致的问题所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。
本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
发表评论 取消回复