经纬度设置、截图操作

一、实验 目的
巩固图像处理编程的步骤格式，熟悉 matlab 图像 处理基本操作（ colobar 、设置坐标 、截图） 。

二、实验 内容
现有一海表面温度数据，数据格式
mat ，数据空间范围 20 40 º N 110 140 º E ，经纬度等
间隔 0.041666668 º分布。
1、 请画出该海表面温度数据图，标出经纬度，给出颜色表（ colorbar ）。
2、 请截取空间范围 28 34 º N 120 125 º E ，存出截取的海表面温度数据块，并画图，标出经纬度，给出颜色表。并在截取前数据图上标出截取的位置。存储图像要求分辨率不低于 200dpi 。

三、实验效果

 

四、实现代码

clear all
%读取海表面温度数据并绘图
load('A20142132014243.L3m_MO_SST_4.mat');
figure,imagesc(sst);
[m n]=size(sst);
title('The temperature of ocean surface ');
colorbar;
%设置经纬度
set(gca,'xtick',1:(n-1)/3:n,'xticklabel',{'110^oE','120^oE','130^oE','140^oE'});
xlabel('Longitude ^oE','fontsize',16,'fontweight','bold');
set(gca,'ytick',1:(m-1)/2:m,'yticklabel',{'40^oN','30^oN','20^oN'});
ylabel('Longitude ^oN','fontsize',16,'fontweight','bold');
%计算四个顶点坐标并绘制矩形框
xmin=(120-110)/30*n;
xmax=(125-110)/30*n;
ymin=(26-20)/20*m;%北纬34°
ymax=(32-20)/20*m;%北纬28°
width=xmax-xmin;
height=ymax-ymin;
rectangle('Position',[xmin,ymin,width,height]);
%保存图片为origin.jpeg，输出图片dpi 为600，格式为jpeg
print origin -djpeg -r600；
%绘制选区
img=sst(ymin:ymax,xmin:xmax);
figure,imagesc(img)
colorbar
title('A part of ocean surface ');
%设置经纬度
[m n]=size(img);
set(gca,'xtick',1:(n-1)/5:n,'xticklabel',{'120^oE','121^oE','122^oE','123^oE','124^oE','125^oE'});
xlabel('Longitude ^oE','fontsize',16,'fontweight','bold');
set(gca,'ytick',1:(m-1)/3:m,'yticklabel',{'34^oN','32^oN','32^oN','30^oN'});
ylabel('Longitude ^oN','fontsize',16,'fontweight','bold');
%保存图片为img.jpeg，输出图片dpi 为600，格式为jpeg
print img -djpeg -r600；