解决Matlab中reshape报错：栅格影像行列数不一致的问题

2023-03-23 01:00| 来源: 网络整理| 查看: 265

今天在用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)

【本文地址】

公司简介

联系我们