我正在尝试使用shp文件为georaster(netCDN文件)制作一个掩码。这类似于此线程(How to limit the raster processing extent using a spatial mask?)中描述的内容,但由于我的georaster不是图像,因此我没有R变量,因此在那里提供的解决方案对我不起作用。
我首先将shapefile加载到MATLAB中: S1 = shaperead('Export_Output_3.shp');
数据具有Lambert Conformal Conical投影,单位为米。
然后我加载netCDF文件并提取第一片数据。我得到一个矩阵x,每个单元格有一个值。它还具有相同参数值的Lambert Conformal Conic投影,但数据位于“度北”和“度东”。
那么如何将这些数据放在一起呢?我尝试使用utm2deg函数将shapefile的单位从米转换为度,但我不知道UTM区域,所以我对此感到困惑。
我认为我首先需要将两个数据集放在相同的单元中,然后将shapefile转换为与netCDF具有相同分辨率的矩阵。我无法执行任何这些步骤。
有人可以帮忙吗?