如何将地理图像中目标范围值的经纬度找出来；如何采用基于shp文件对栅格进行分区统计?

文章介绍了使用Matlab进行地理坐标图像中指定值范围经纬度提取的方法，以及如何对不规则shp文件划分的区域进行栅格数据统计。作者分享了处理数据的实用代码，以帮助科研人员解决GIS中的投影和范围问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

文章目录

前言
一、需要环境
二、代码
- 1.提取地理坐标图像中指定的值的经纬度
- 2.对栅格进行分区统计
总结

前言

如何将地理坐标图像中指定的范围值的经纬度找出来？就是根据值来找位置（经纬度），我们只有tif影像，如何来找某个值范围的经纬度？如果我们有多幅影像？
如何采用基于shp文件对栅格进行分区统计，shp中有很多子区域，根据子区域进行分区统计大家都做过，GIS可以完成，但是GIS有个问题因为范围或者投影的问题会很容易报错，因此上代码强制来解决。

一、需要环境

matlab的地理工具箱默认已经安装

二、代码

1.提取地理坐标图像中指定的值的经纬度

比如提取一副影像中大于1000的值的坐标经纬度：

close all
clear 
clc

tif_file = 'E:\LHASADATA\Figurecodes\data\Fig.4\370\DWELexposure.tif';
info = geotiffinfo(tif_file);
% 获取地理坐标信息
[Historicalda, R] = geotiffread(tif_file);
Ls=size(Historicalda);
lat1 = [info.CornerCoords.Lat(1) info.CornerCoords.Lat(2); info.CornerCoords.Lat(4) info.CornerCoords.Lat(3)];
lon1 = [info.CornerCoords.Lon(1) info.CornerCoords.Lon(2); info.CornerCoords.Lon(4) info.CornerCoords.Lon(3)];
[Xi,Yi] = meshgrid(1:2,1:2);
[XI,YI] = meshgrid