求采样点像元在影像中第几行，第几列

原创

于 2024-08-30 19:51:11 发布

· 176 阅读

0 ·

版权

文章标签：

#python #pandas

1.导入

import ee
import os
import geemap
import pandas as pd
from pprint import pprint
geemap.set_proxy(****)
ee.Initialize()
import math
Map = geemap.Map()

2.导入经纬度和影像信息

pointTable = pd.read_excel(r"带有经纬度影像信息的表.xlsx")

3.提取行列

def get_pixel_coords(row):
    lon = row['longitude']
    lat = row['latitude']
    image_id = row['ImageId']
    
    # 加载影像并获取投影
    image = ee.Image(image_id)
    projection = image.select(0).projection()
    
    # 创建点并转换为影像投影坐标
    point = ee.Geometry.Point(lon, lat)
    transformed_point = point.transform(projection, 1)
    
    # 获取像元的行列号（浮点型）
    pixel_coords = transformed_point.coordinates().getInfo()
    col, row = pixel_coords
    
    # 获取影像的行数和列数
    image_info = image.getInfo()
    rows = image