【Python&GIS】判断图片中心点/经纬度点是否在某个面内

Python的exifread库可以获取图片中的源数据信息，包括经纬度、相机厂商、曝光时间、焦距、拍摄时间、拍摄地点等等信息。我们可以通过exifread库从图片中获取图片的经纬度，再通过shapely库判断该经纬度对应的点是否在我们的感兴趣区内。

当然你需要知道图片的经纬度是什么坐标系，大部分是WGS84地理坐标系，具体需要自己了解。如果是GPS定位那么就是WGS84坐标系。

shapely主要是在笛卡尔平面对几何对象进行操作和分析，它是一个BSD授权的Python包。shapely不关心数据格式或坐标系，但可以很容易地与这些文件包集成。

geopandas是用来处理地理空间数据的python第三方库，它是在pandas的基础上建立的，完美地融合了pandas的数据类型，并且提供了操作地理空间数据的高级接口，使得在python中进行GIS操作变成可能。

1.导入所需要的三方库。exifread、geopandas、shapely

import exifread

import geopandas

from shapely.geometry import Point

2.打开本地图片,使用exifread读取图片的源数据信息，读取其中的经纬度信息，获取经纬度后记得关闭图片进程。

f = open(path_image,'rb')中的path_image即为本地图片的路径。

f = open(path_image, 'rb')

# 打开图片

contents = exifread.process_file(f)

# 获取图片属性信息

longitude = contents["GPS GPSLongitude"].values

# 获取Long

longitude_f = longitude[0].num/longitude[0].den + (longitude[1].num/longitude[1].den/60) + (longitude[2].num/longitude[2].den/3600)

latitude = contents["GPS GPSLatitude"].values

# 获取Lat

latitude_f = latitude[0].num/latitude[0].den + (latitude[1].num/latitude[1].den/60) + (latitude[2].num/latitude[2].den/3600)

# print("经度:", longitude_f)  # contents['GPS GPSLatitudeRef']

# print("纬度:", latitude_f)  # contents['GPS GPSLongitudeRef']

f.close()

# 关闭图片

3.创建工作空间，将图片经纬度信息写入矢量图层，使用shapely.geometry的contains函数判断点是否在面内。

根据返回值判断点是否在面内。返回值为True or False,True则点在面内，False则点在面外。

datasource = geopandas.GeoDataFrame.from_file(path_conservation_area, encoding='utf-8')

# 创建工作空间

boundary_shape = datasource['geometry'][0]

# 将shp区域加载到工作空间

determine_contain = boundary_shape.contains(Point([longitude_f, latitude_f]))

# 判断点是否在区间内

4.完整代码，记得改成自己的路径。

1）面矢量数据“path_conservation_area”的路径

2）本地图片“path_image”的路径

# -*- coding: utf-8 -*-

"""

@Time ： 2023/5/15 15:57

@Auth ： RS迷途小书童

@File ：Discriminative point surface.py

@IDE ：PyCharm

"""

import exifread

import geopandas

from shapely.geometry import Point

def Get_LatLon(path_image):

    """

    :param path_image: 输入图片路径

    :return: 返回经纬度

    """

    f = open(path_image, 'rb')

    # 打开图片

    contents = exifread.process_file(f)

    # 获取图片属性信息

    longitude = contents["GPS GPSLongitude"].values

    # 获取Long

    longitude_f = longitude[0].num/longitude[0].den + (longitude[1].num/longitude[1].den/60) + (longitude[2].num/longitude[2].den/3600)

    latitude = contents["GPS GPSLatitude"].values

    # 获取Lat

    latitude_f = latitude[0].num/latitude[0].den + (latitude[1].num/latitude[1].den/60) + (latitude[2].num/latitude[2].den/3600)

    # print("经度:", longitude_f)  # contents['GPS GPSLatitudeRef']

    # print("纬度:", latitude_f)  # contents['GPS GPSLongitudeRef']

    f.close()

    # 关闭图片

    return longitude_f, latitude_f

def Filter_images(path_conservation_area, longitude_f, latitude_f):

    """

    :param path_conservation_area: 面数据路径

    :param longitude_f: 经度

    :param latitude_f: 纬度

    :return: True or False

    """

    # path_conservation_area = r'G:/B60/B60.shp'

    datasource = geopandas.GeoDataFrame.from_file(path_conservation_area, encoding='utf-8')

    # 创建工作空间

    boundary_shape = datasource['geometry'][0]

    # 将shp区域加载到工作空间

    determine_contain = boundary_shape.contains(Point([longitude_f, latitude_f]))

    # 判断点是否在区间内

    return determine_contain

if __name__ == "__main__":

    path_conservation_area = r'G:/B60/B60.shp'

    # 面数据路径

    path_image = "G:/image1.jpg"

    # 图片路径

    longitude_f, latitude_f = Get_LatLon(path_image)

    # 执行获取经纬度函数

    determine_contain = Filter_images(path_conservation_area, longitude_f, latitude_f)

    # 执行判断点是否在区间的函数

    print(determine_contain)