现有一个全国的tiff其记录了全国的建成区其数值范围为0和11代表非建成区0代表建成区然后现有一个全国分省市的shp如何用matlab将全国的建成区tiff裁剪为各个省市的建成区tiff裁剪后的tiff包含三个值0120代表该区域在上海市的shp内但不属于建成区1代表不在上海市的shp内2代表在上海市shp内且为建成区

您可以使用MATLAB中的Spatial Statistics Toolbox来实现这个任务。以下是一个可能的解决方案：

1. 首先，加载全国的建成区tiff，可以使用imread函数：

building_tiff = imread('全国建成区.tif');

2. 加载省市的shp文件，可以使用shaperead函数：

provinces = shaperead('全国省市.shp');

3. 创建一个与建成区tiff相同大小的矩阵，用于存储裁剪后的结果：

cropped_tiff = zeros(size(building_tiff));

4. 对每个省市进行循环，裁剪建成区tiff并将结果存储到cropped_tiff中：

for i = 1:numel(provinces)
    % 获取当前省市的边界框
    bbox = provinces(i).BoundingBox;
    
    % 裁剪建成区tiff
    cropped_building = imcrop(building_tiff, bbox);
    
    % 将裁剪的结果存储到cropped_tiff中
    cropped_tiff(bbox(2):bbox(2)+bbox(4), bbox(1):bbox(1)+bbox(3)) = cropped_building;
end

5. 对于上海市的shp，可以使用同样的方法裁剪建成区tiff，并将结果标记为2：

% 假设上海市的shp在provinces中的索引为sh_index
sh_bbox = provinces(sh_index).BoundingBox;
sh_cropped_building = imcrop(building_tiff, sh_bbox);
cropped_tiff(sh_bbox(2):sh_bbox(2)+sh_bbox(4), sh_bbox(1):sh_bbox(1)+sh_bbox(3)) = 2;

6. 最后，根据裁剪后的结果，将0值标记为不在上海市shp内的建成区，将1值标记为不在上海市shp内的非建成区：

cropped_tiff(cropped_tiff == 0) = 1;

现在，cropped_tiff中的值为0表示在上海市shp内但不属于建成区，1表示不在上海市shp内，2表示在上海市shp内且为建成区。

请注意，上述代码仅为示例，您需要根据实际情况进行调整和修改