Matlab 根据 shp 裁剪矩阵/图像 函数

最新推荐文章于 2023-05-31 16:17:16 发布
最新推荐文章于 2023-05-31 16:17:16 发布
阅读量5.1k 收藏 44
点赞数
分类专栏： # Matlab 文章标签： matlab
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Gou_Hailong/article/details/108209764
版权

关于这个函数，我最近（2020.11.1）又修改了一下，之前的版本第二步通用性不强，这波应该可以通用了。现将几个版本命个名：
2020-8-24 V1.0 有内切圆外切圆，第二步不通用
2020-11-1 V2.0 有内切圆外切圆，第二步通用
2020-11-1 V2.1 无内切圆外切圆，第二步通用

1、全代码

1.1 V1.0

function varargout=tailorsheng(varargin)
%% 此函数用于根据 shp 裁剪矩阵/图像
% 输入：
%       P2file shpfile
%       TDMSPfile tiffile，可以不是tif
%       str  要裁剪谁
%       mode 1省0国家
% 输出：
%       sheng shp
%       photo mat图像
%       ex    经纬度极值
% 调用：
%       TDMSPfile='D:\复习资料\自学\7_科研\夜光遥感\data\DMSP\F121999.v4\F121999.v4b_web.stable_lights.avg_vis.tif';
%       P2file='D:\下载\Useful\shp\国家基础地理数据\bou2_4m\bou2_4p.shp';%省界多边形
%       str='北京市';
%       mode=1;
%       [sheng,photo,ex]=tailorsheng(P2file,TDMSPfile,str);
%       [sheng,photo,ex]=tailorsheng(P2file,TDMSPfile,str,mode);
%-------------------------------------------------------------------
    %%%%    Authors:   Bill O'Hanlon
    %%%%    EMAIL:     ohanlon@qq.com
    %%%%    DATE:      24-08-2020
%% 输入
disp('-------function tailorsheng------');
tic %开始计时
mode=1;
if nargin==3
    P2file=varargin{1}; %注意，这里是{, (会是元胞
    TDMSPfile=varargin{2};
    str=varargin{3};
elseif nargin==4
    P2file=varargin{1}; %注意，这里是{, (会是元胞
    TDMSPfile=varargin{2};
    str=varargin{3};
    mode=varargin{4};%这里不做变换的话会是元胞
else
    disp('参数不合要求！');
    return;
end
%% 第一步，提取目标省份的 shp，注意看其范围
[henan,ex]=drawsheng(P2file,str,mode);
disp('shp文件搞定！');
%% 第二步，把图像大致范围剪出来，这步根据图像种类不同，代码也不一样，Attention!!
TDMSP=imread(TDMSPfile); %DMSP范围 180°W-180°E，65°S-75°N
Xmax=ex(1);              %X是经度，Y是纬度
Ymax=ex(2);              %这些已经取整了。
Xmin=ex(3);
Ymin=ex(4);
Xmax1=Xmax+180;          %这是裁剪时用的
Ymax1=75-Ymax;
Xmin1=Xmin+180;
Ymin1=75-Ymin;
Henan=TDMSP(Ymax1*120:Ymin1*120,Xmin1*120:Xmax1*120);%根据自己需要裁减
disp('粗裁剪完成！');
%% 第三步，计算不规则边界的内切圆和外接圆
Y=Ymin:0.0083333333:Ymax; %分辨率0.0083333333°  1°=120
X=Xmin:0.0083333333:Xmax; %求逻辑矩阵用到
[a,b]=size(Henan);
Y2=repmat(Y',1,b);
X2=repmat(X,a,1);
Y2=flipud(Y2);                     %这个纬度得上下翻转一下。
xv=henan.X;yv=henan.Y;             %提取边界
xv = xv(1:end-1); yv = yv(1:end-1);%把最后一个NaN去掉
disp('需要计算逻辑的区域为 Figure 2 红色区域');
bianjie=[xv;yv];
[zhongxin1,zhongxin2,smallR,bigR]=getZhongxin(bianjie,X2,Y2);
%% 第四步，根据内切圆和外接圆+边界进行裁剪
disp('开始计算逻辑矩阵..');
photo=zeros(a,b);
photo(:)=nan;
for i=1:a
    for j=1:b
        if mod(i*b+j,10000)==0
            disp([int2str(i*b+j),' / ',int2str(a*b), ' OK!']);
        end
        x=X2(i,j);
        y=Y2(i,j);
        if norm(zhongxin1-[x;y])<=smallR
            photo(i,j)=Henan(i,j);%内接圆里面的都在
            continue;
        elseif norm(zhongxin2-[x;y])>=bigR
            continue;%外接圆外面的都不在
        elseif inpolygon(x,y,xv,yv)
            photo(i,j)=Henan(i,j);
        end
    end
end
disp('细裁剪完成！');
%% 第五步，画图，裁剪后的
[x,x1,y,y1] = getxy(X,Y);
x=(x-Xmin).*120; 
y=(y-Ymin).*120;
photo=flipud(photo);%contourf 上下翻转一下，才变成imshow
figure;
contourf(photo,'LineStyle','none');
colormap(gray);colorbar %jet
set(gca,'XTick',x,'XTicklabel',x1);   %设置x,y轴
set(gca,'YTick',y,'YTicklabel',y1);
title([str,'夜光遥感影像']);
%% 输出
if nargout==3
    varargout{1}=henan; 
    varargout{2}=photo;
    varargout{3}=ex;
elseif nargout==2
    varargout{1}=henan; 
    varargout{2}=photo;
elseif nargout==1
    varargout{1}=henan; 
elseif nargout==0
    return;
else
    disp('参数不合要求！');
    return;
end
disp('--------Finished!--------');
toc  %展示运行时间
end

1.2 V2.0

function varargout=tailorsheng(varargin)
%% 此函数用于根据 shp 裁剪矩阵/图像
% 输入：
%       P2file shpfile
%       TDMSPfile tiffile，可以不是tif
%       str  要裁剪谁
%       mode 1省0国家
% 输出：
%       sheng shp
%       photo mat图像
%       ex    经纬度极值
% 调用：
%       TDMSPfile='D:\复习资料\自学\7_科研\夜光遥感\data\DMSP\F121999.v4\F121999.v4b_web.stable_lights.avg_vis.tif';
%       P2file='D:\下载\Useful\shp\国家基础地理数据\bou2_4m\bou2_4p.shp';%省界多边形
%       str='北京市';
%       mode=1;
%       [sheng,photo,ex]=tailorsheng(P2file,TDMSPfile,str);
%       [sheng,photo,ex]=tailorsheng(P2file,TDMSPfile,str,mode);
%-------------------------------------------------------------------
    %%%%    Authors:   Bill O'Hanlon
    %%%%    EMAIL:     ohanlon@qq.com
    %%%%    DATE:      01-11-2020
%% 输入
disp('-------function tailorsheng------');
tic %开始计时
mode=1;
if nargin==3
    P2file=varargin{1}; %注意，这里是{, (会是元胞
    TDMSPfile=varargin{2};
    str=varargin{3};
elseif nargin==4
    P2file=varargin{1}; %注意，这里是{, (会是元胞
    TDMSPfile=varargin{2};
    str=varargin{3};
    mode=varargin{4};%这里不做变换的话会是元胞
else
    disp('参数不合要求！');
    return;
end
%% 第一步，提取目标省份的 shp，注意看其范围
[henan,ex]=drawsheng(P2file,str,mode);
disp('shp文件搞定！');
%% 第二步，把图像大致范围剪出来，这步根据图像种类不同，代码也不一样，Attention!!
TDMSP=imread(TDMSPfile); %DMSP范围 73.49781527°E-135.0894573°E，3.83544513°N-53.56042524°N
Xmax=ex(1);              %X是经度，Y是纬度
Ymax=ex(2);              %这些已经取整了。
Xmin=ex(3);
Ymin=ex(4);
Info=geotiffinfo(TDMSPfile); %读取图像信息，比如经纬度，分辨率
global PixelScale; %定义了几个全局变量，
global Pixel_;
global XBound;
global YBound;
PixelScale=Info.PixelScale(1,1);
Pixel_=1/PixelScale;
XBound=Info.BoundingBox(1,1);
YBound=Info.BoundingBox(2,2);
[a,b]=size(TDMSP);
if(Xmax>Info.BoundingBox(2,1)) Xmax=Info.BoundingBox(2,1); end
if(Xmin<Info.BoundingBox(1,1)) Xmin=Info.BoundingBox(1,1); end
if(Ymax>Info.BoundingBox(2,2)) Ymax=Info.BoundingBox(2,2); end
if(Ymin<Info.BoundingBox(1,2)) Ymin=Info.BoundingBox(1,2); end
Xmax1=Xmax-XBound;          %这是裁剪时用的
Ymax1=YBound-Ymax;
Xmin1=Xmin-XBound;
Ymin1=YBound-Ymin;
hangs=round(Ymax1*Pixel_); %上界
hangx=round(Ymin1*Pixel_); %下界
liez=round(Xmin1*Pixel_);  %左界
liey=round(Xmax1*Pixel_);  %右
if(hangs<1) hangs=1; end
if(hangx>a) hangx=b;end
if(liez<1) liez=1; end
if(liey>b) liey=b; end

Henan=TDMSP(hangs:hangx,...
    liez:liey);%根据自己需要裁减
sprintf('粗裁剪完成！ %d,%d,%d,%d',hangs,hangx,liez,liey);
%% 第三步，计算不规则边界的内切圆和外接圆
Y=Ymin:PixelScale:Ymax; %分辨率0.0083333333°  1°=120
X=Xmin:PixelScale:Xmax; %求逻辑矩阵用到
%[a,b]=size(Henan);
a=size(Y,2);
b=size(X,2);
Y2=repmat(Y',1,b);
X2=repmat(X,a,1);
Y2=flipud(Y2);                     %这个纬度得上下翻转一下。
xv=henan.X;yv=henan.Y;             %提取边界
xv = xv(1:end-1); yv = yv(1:end-1);%把最后一个NaN去掉
disp('需要计算逻辑的区域为 Figure 2 红色区域');
bianjie=[xv;yv];
[zhongxin1,zhongxin2,smallR,bigR]=getZhongxin(bianjie,X2,Y2);
%% 第四步，根据内切圆和外接圆+边界进行裁剪
disp('开始计算逻辑矩阵..');
photo=zeros(a,b);
photo(:)=nan;
for i=1:a
    for j=1:b
        if mod(i*b+j,10000)==0
            disp([int2str(i*b+j),' / ',int2str(a*b), ' OK!']);
        end
        x=X2(i,j);
        y=Y2(i,j);
        if norm(zhongxin1-[x;y])<=smallR
            photo(i,j)=Henan(i,j);%内接圆里面的都在
            continue;
        elseif norm(zhongxin2-[x;y])>=bigR
            continue;%外接圆外面的都不在
        elseif inpolygon(x,y,xv,yv)
            photo(i,j)=Henan(i,j);
        end
    end
end
disp('细裁剪完成！');
%% 第五步，画图，裁剪后的
[x,x1,y,y1] = getxy(X,Y);
x=(x-Xmin).*Pixel_; 
y=(y-Ymin).*Pixel_;
photo=flipud(photo);%contourf 上下翻转一下，才变成imshow
figure;
contourf(photo,'LineStyle','none');
colormap(gray);colorbar %jet
set(gca,'XTick',x,'XTicklabel',x1);   %设置x,y轴
set(gca,'YTick',y,'YTicklabel',y1);
title([str,'夜光遥感影像']);
%% 输出
photo(photo>=0)=1;
if nargout==3
    varargout{1}=henan; 
    varargout{2}=photo;
    varargout{3}=ex;
elseif nargout==2
    varargout{1}=henan; 
    varargout{2}=photo;
elseif nargout==1
    varargout{1}=henan; 
elseif nargout==0
    return;
else
    disp('参数不合要求！');
    return;
end
disp('--------Finished!--------');
toc  %展示运行时间
end

1.3 V2.1

function varargout=tailorsheng(varargin)
%% 此函数用于根据 shp 裁剪矩阵/图像
% 输入：
%       P2file shpfile
%       TDMSPfile tiffile，可以不是tif
%       str  要裁剪谁
%       mode 1省0国家
% 输出：
%       sheng shp
%       photo mat图像
%       ex    经纬度极值
% 调用：
%       TDMSPfile='D:\复习资料\自学\7_科研\夜光遥感\data\DMSP\F121999.v4\F121999.v4b_web.stable_lights.avg_vis.tif';
%       P2file='D:\下载\Useful\shp\国家基础地理数据\bou2_4m\bou2_4p.shp';%省界多边形
%       str='北京市';
%       mode=1;
%       [sheng,photo,ex]=tailorsheng(P2file,TDMSPfile,str);
%       [sheng,photo,ex]=tailorsheng(P2file,TDMSPfile,str,mode);
%-------------------------------------------------------------------
    %%%%    Authors:   Bill O'Hanlon
    %%%%    EMAIL:     ohanlon@qq.com
    %%%%    DATE:      24-08-2020
%% 输入
disp('-------function tailorsheng------');
tic %开始计时
mode=1;
if nargin==3
    P2file=varargin{1}; %注意，这里是{, (会是元胞
    TDMSPfile=varargin{2};
    str=varargin{3};
elseif nargin==4
    P2file=varargin{1}; %注意，这里是{, (会是元胞
    TDMSPfile=varargin{2};
    str=varargin{3};
    mode=varargin{4};%这里不做变换的话会是元胞
else
    disp('参数不合要求！');
    return;
end
%% 第一步，提取目标省份的 shp，注意看其范围
[henan,ex]=drawsheng(P2file,str,mode);
disp('shp文件搞定！');
%% 第二步，把图像大致范围剪出来，这步根据图像种类不同，代码也不一样，Attention!!
TDMSP=imread(TDMSPfile); %DMSP范围 73.49781527°E-135.0894573°E，3.83544513°N-53.56042524°N
Xmax=ex(1);              %X是经度，Y是纬度
Ymax=ex(2);              %这些已经取整了。
Xmin=ex(3);
Ymin=ex(4);
Info=geotiffinfo(TDMSPfile); %读取图像信息，比如经纬度，分辨率
global PixelScale; %定义了几个全局变量，
global Pixel_;
global XBound;
global YBound;
PixelScale=Info.PixelScale(1,1);
Pixel_=1/PixelScale;
XBound=Info.BoundingBox(1,1);
YBound=Info.BoundingBox(2,2);
[a,b]=size(TDMSP);
if(Xmax>Info.BoundingBox(2,1)) Xmax=Info.BoundingBox(2,1); end
if(Xmin<Info.BoundingBox(1,1)) Xmin=Info.BoundingBox(1,1); end
if(Ymax>Info.BoundingBox(2,2)) Ymax=Info.BoundingBox(2,2); end
if(Ymin<Info.BoundingBox(1,2)) Ymin=Info.BoundingBox(1,2); end
Xmax1=Xmax-XBound;          %这是裁剪时用的
Ymax1=YBound-Ymax;
Xmin1=Xmin-XBound;
Ymin1=YBound-Ymin;
hangs=round(Ymax1*Pixel_); %上界
hangx=round(Ymin1*Pixel_); %下界
liez=round(Xmin1*Pixel_);  %左界
liey=round(Xmax1*Pixel_);  %右
if(hangs<1) hangs=1; end
if(hangx>a) hangx=b;end
if(liez<1) liez=1; end
if(liey>b) liey=b; end

Henan=TDMSP(hangs:hangx,...
    liez:liey);%根据自己需要裁减
sprintf('粗裁剪完成！ %d,%d,%d,%d',hangs,hangx,liez,liey);
%% 第三步，计算不规则边界的内切圆和外接圆
Y=Ymin:PixelScale:Ymax; %分辨率0.0083333333°  1°=120
X=Xmin:PixelScale:Xmax; %求逻辑矩阵用到
%[a,b]=size(Henan);
a=size(Y,2);
b=size(X,2);
Y2=repmat(Y',1,b);
X2=repmat(X,a,1);
Y2=flipud(Y2);                     %这个纬度得上下翻转一下。
xv=henan.X;yv=henan.Y;             %提取边界
xv = xv(1:end-1); yv = yv(1:end-1);%把最后一个NaN去掉
disp('需要计算逻辑的区域为 Figure 2 红色区域');
bianjie=[xv;yv];
%[zhongxin1,zhongxin2,smallR,bigR]=getZhongxin(bianjie,X2,Y2);
%% 第四步，根据内切圆和外接圆+边界进行裁剪
disp('开始计算逻辑矩阵..');
photo=zeros(a,b);
photo(:)=nan;
for i=1:a
    for j=1:b
        if mod(i*b+j,10000)==0
            disp([int2str(i*b+j),' / ',int2str(a*b), ' OK!']);
        end
        x=X2(i,j);
        y=Y2(i,j);
        if inpolygon(x,y,xv,yv)
            photo(i,j)=Henan(i,j);
        end
    end
end
disp('细裁剪完成！');
%% 第五步，画图，裁剪后的
[x,x1,y,y1] = getxy(X,Y);
x=(x-Xmin).*Pixel_; 
y=(y-Ymin).*Pixel_;
photo=flipud(photo);%contourf 上下翻转一下，才变成imshow
figure;
contourf(photo,'LineStyle','none');
colormap(gray);colorbar %jet
set(gca,'XTick',x,'XTicklabel',x1);   %设置x,y轴
set(gca,'YTick',y,'YTicklabel',y1);
title([str,'夜光遥感影像']);
%% 输出
photo(photo>=0)=1;
if nargout==3
    varargout{1}=henan; 
    varargout{2}=photo;
    varargout{3}=ex;
elseif nargout==2
    varargout{1}=henan; 
    varargout{2}=photo;
elseif nargout==1
    varargout{1}=henan; 
elseif nargout==0
    return;
else
    disp('参数不合要求！');
    return;
end
disp('--------Finished!--------');
toc  %展示运行时间
end

依赖：

内切圆和外接圆 :https://blog.csdn.net/Gou_Hailong/article/details/108206335
扣shp：https://blog.csdn.net/Gou_Hailong/article/details/108209395
缩放矩阵或图像：https://blog.csdn.net/Gou_Hailong/article/details/108206521
画地图注释：https://blog.csdn.net/Gou_Hailong/article/details/108208442

注：这个函数裁剪小点的边界还好，如果边界太大或者图像太大，运行时间会超级长的，这酸爽被我记录在了下面博客中：

https://blog.csdn.net/Gou_Hailong/article/details/108147268

2、调用

调用代码：

clc
clear
TDMSPfile='D:\复习资料\自学\7_科研\夜光遥感\data\DMSP\F121999.v4\F121999.v4b_web.stable_lights.avg_vis.tif';
P2file='D:\下载\Useful\shp\国家基础地理数据\bou2_4m\bou2_4p.shp';%省界多边形
str='北京市';
mode=1;
[sheng,photo,ex]=tailorsheng(P2file,TDMSPfile,str,mode);

结果：

在这里插入图片描述

CSDN 脑子秀逗了，说图片违规，dd，
博客园链接：https://www.cnblogs.com/Gou-Hailong/p/13559163.html

流浪猪头拯救地球
关注
  • 0
    点赞
  • 踩
  • 44
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 3
    评论
专栏目录
使用MATLAB对数组进行裁剪
07-10
文档中包含网易云课堂MATLAB教程中 冒号区间裁剪，冒号等差序列，删除整行整列的操作，简单易学，内容详细。
MATLAB中使用shp批量裁剪影像数据（一一对应）
qq_46546227的博客
02-04 1508
实现用一个shp文件中的不同要素一一对应裁剪不同的tiff数据。
Matlab利用shp文件裁剪netCDF（nc）数据及将数据写为nc文件
正在努力生长的小菜鸡
05-31 2274
读取shp文件 shape = shaperead('shapefile.shp');% 读取 netCDF 文件 nc = ncinfo('data.nc');% 确定裁剪区域 lon = ncread('data.nc' , 'lon');% 裁剪数据 data_new = data(ilon , ilat);
Matlab用shp文件对数据进行剪裁
z894730988的博客
05-22 3848
Matlab 剪裁
matlab矩阵剪切
huiyayaya的博客
03-29 528
matlab矩阵剪切。
ArcGIS根据SHP裁剪栅格的Python代码
12-10
编写了一个Python文件，用户可自行修改输入栅格、SHP、字段名、保存路径。代码会根据选择的字段名自动对栅格数据进行裁剪，生成多个栅格文件。
利用shapefile对栅格图像进行裁剪
07-21
利用shapefile对栅格图像进行裁剪. ; :Syntax ; RasterSubsetViaShapefile, Fid, shpFile=string, [pos=array], ; [inside={0|1}], [outFile=｛string|variable}], [r_fid=variable] ; ; :Params: ; Fid -- 输入文件...
Matlab处理shp文件的程序
06-07
很多地学或者海洋学用matlab的朋友，都会碰到shp数据，会遇到在matla中如何快速简单提取或者画出自己想要的图，相信这个程序肯定可以帮到你。shp数据都会有shx格式，最新版的matlab处理shx特别方便可以实现和shp一样...
利用shp矢量数据裁剪遥感影像的Python实现
02-25
以上代码展示了如何使用Python结合`geopandas`和`rasterio`库，根据Shp矢量数据裁剪遥感影像的基本步骤。在实际应用中，可能还需要考虑其他因素，如投影转换、坐标系兼容性、数据预处理等。此外，`Input_Path.txt`和...
matlab读取矢量文件.rar_matlab shp_matlab shp_matlab提取shp_matlab读取shp_
07-14
matlab读取shp文件，显示并裁剪栅格数据
使用SHP数据批量裁剪TIF图像并转为JPG格式
aaaq_的博客
11-08 2688
1.shp数据按照属性分块输出2.使用shp数据批量裁剪遥感影像，将影像裁剪为小块的tif图像3.将tif转为jpg，并调整为同一大小。
仅个人记录：shp裁剪tif；shp裁剪shp；矢量转栅格；多个shp裁剪shp；栅格边界矢量化。汇总：输入shp和影像，输出影像对应的标签（栅格边界矢量化，shp裁剪shp，shp转tif)
weixin_61235989的博客
05-18 972
【代码】仅个人记录：shp裁剪tif；shp裁剪shp；矢量转栅格；多个shp裁剪shp；汇总：输入shp和影像，输出影像对应的标签。
MATLAB 根据经纬坐标生成点的shp文件或流域的边界shp文件
天外来客的花田
08-29 4060
任务描述：有站点坐标，制作.shp文件用于arcgis后续处理 [txt,~] = xlsread('C:\Users\DELL\Desktop\Works\数据整理.xlsx',12); location = txt(:,1:3); for i = 1:length(txt) %前面的部分为必要的字段 Point.Geometry = 'Point';%元素类型，Point，Line，Polygon等 Point.X = location(i,1);%坐标X 对应经度 P
matlab矩阵抽取截取删除及叠加
lzdwlhn的专栏
07-02 7436
查看矩阵维度 函数 size 格式 d = size(X) %查看x的维度(行列)信息 矩阵截取 矩阵X b=X(1:m,:)取1至m行，列取全部。 c=X(:,1:n)取1至n列，行取全部。   删除矩阵中某些行 a='2 3 4 5';%定义字符串 b=str2num(a);% 将字符串转换为矩阵 tr=rand(6,9);%定义随机矩阵 tr(b,:)=[];%删除矩
2019-11-08 使用shp文件批量裁剪GeoTif遥感图像
liu_zonghui的博客
11-08 3026
需求与目的 在批量将MODIS数据hdf转为GeoTif格式之后（具体操作可参考我的另一篇博客MODIS数据hdf转tif），因为是全球数据，而我只需要某一块区域的遥感影像，故而需要进行批量裁剪。（如下图所示） 数据按照月份文件夹存储，每个月份文件夹下有当月的tif文件。如下图所示： 运行环境 Windows 7、10（其它系统版本未尝试）； ArcGIS10.5自带Python 2.7.12，...
Matlab利用已知方阵创建斜对角矩阵
huangzhywin的博客
11-21 1155
Matlab利用已知方阵创建斜对角矩阵
MATLAB 基于流域边界批量裁剪栅格图像
天外来客的花田
12-21 1267
基于流域边界批量裁剪栅格图像
matlab柱状斜线_Matlab小练习：按斜线方向依次赋值矩阵
weixin_36190719的博客
12-30 727
来自知乎问题，觉得挺有意思，留给学生解答之余，我也做了一番思考，得到三种解法。题目如下：以n=80为例，————————————————————一、先要根据 确定矩阵的阶数 如果先生成足够大矩阵，再删掉全为零的列，有点太low了。所以先思考一个算法：观察数据放置的特点，不难发现这样一个关系：变形：观察，若 足够大， 起决定作用， 取整一定是 较小的 (或 ) 对吗？试一试。测试 至 :n=...
matlab取图像矩阵并切割,matlab – 使用图形切割与种子点的图像分割
weixin_42533120的博客
03-17 491
我不熟悉你链接到的FEX的图形实现,但是我会用GCMex matlab wrapper(适当的披露：我实现了这个包装)来显示一个例子.假设你有一个具有n个像素的尺寸大小(S)的图像,并且K是尺寸为n×n的稀疏矩阵,K(ii,jj)表示ii和jj像素连接的程度(对于相邻的ii和jj).此外,您有一个前景像素的掩码和背景像素的掩码v被视为硬约束.首先,使用u和v构建数据项：Dc = 1000*[u(:...
matlab用shp裁剪图像
最新发布
06-24
### 回答1： MATLAB是一种强大的计算机软件，用于数据可视化、数值计算、绘图和编程等多种功能。其中，MATLAB也支持使用SHP（Shapefile）文件来裁剪图像。 首先，我们需要将所需的SHP文件导入MATLAB的工作环境中。这可以通过使用shapefile函数来实现，该函数需要传递SHP文件的完整路径作为参数。 接下来，我们需要加载要裁剪的原始图像。这可以通过使用imread函数来实现，该函数需要传递原始图像的完整路径或者二进制图像数据作为参数。 一旦加载了SHP文件和原始图像，我们可以使用mapshow函数来将其显示在屏幕上。然后，我们可以使用imcrop函数根据SHP文件的边界框来裁剪图像。 具体来说，我们可以使用poly2mask函数将SHP文件中的多边形边界转换为二进制掩模。然后，我们可以将掩模传递给imcrop函数来裁剪原始图像。裁剪后的图像可以使用imshow函数来展示。 最后，我们需要将裁剪后的图像保存为新文件。这可以使用imwrite函数来实现，该函数需要传递裁剪后的图像数据和文件路径作为参数。 总之，通过使用MATLAB中的SHP和图像处理函数，我们可以轻松地裁剪图像。这在地理信息系统（GIS）和遥感领域中非常有用。 ### 回答2： 在MATLAB中，我们可以使用shapefile来裁剪图像。shapefile是一种常见的GIS数据格式，用于存储空间数据。要使用shapefile进行裁剪，需要采取以下步骤： 1. 使用shapefile函数读取要用于裁剪的多边形区域的shp文件，并存储在一个变量中。 2. 使用imread函数读取要进行裁剪的图像，并存储在另一个变量中。 3. 使用imref2d函数获取标准空间参考对象，其代表图像的空间参考，并存储在一个变量中。 4. 使用mapshow函数将多边形区域进行绘制并显示在图像上，以便确定要裁剪哪个区域。 5. 使用roipoly函数选择裁剪区域并存储在一个变量中。 6. 使用imcrop函数对图像进行裁剪并保存。裁剪后的图像是只包含原图像中被选中区域的图像。 这是使用MATLAB进行shp裁剪的基本步骤。使用shapefile进行裁剪可以更精确地剪裁图像，从而提高图像分析的准确性和精度。
评论 3
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

流浪猪头拯救地球

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值