- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我目前有一张来自 GeoTiff 文件的 171 x 171 图像(尽管在其他情况下,我可能有更大的图像)。我的目标是获取图像中的每个像素并将其转换为纬度/经度对。
我已经能够根据这篇 StackOverflow 帖子将图像的角点转换为纬度/经度对:Obtain Latitude and Longitude from a GeoTIFF File。这篇文章很有帮助,因为我的原始坐标位于 UTM Zone 15。
但是,我现在想将图像的所有像素转换为纬度、经度对,并将结果存储在相同维度的 numpy 数组中。所以输出将是一个 171 x 171 x 2 的 numpy 数组,numpy 数组的每个元素都是(经度,纬度)对的元组。
我看到的最相关的帖子是 https://scriptndebug.wordpress.com/2014/11/24/latitudelongitude-of-each-pixel-using-python-and-gdal/ 。然而,该帖子建议在每个像素上创建一个 for 循环并转换为纬度、经度。有没有更有效的方法?
只是为了提供更多关于我的实际用例的上下文,我的最终目标是我有一堆卫星图像(例如,在这种情况下,每张图像都是 171 x 171)。我正在尝试创建一个建筑分割模型。现在,我正在尝试通过在每个图像上创建一个掩码来生成标记数据点,如果它对应于建筑物,则将像素标记为 1,否则为 0。首先,我使用 Microsoft US Building Footprint 数据:https://github.com/microsoft/USBuildingFootprints 其中他们发布了他们检测到的建筑物的多边形(由纬度、经度定义)的 GeoJSON 文件。我考虑这样做的方式是:
目前正在进行第(1)步,即高效地求出图像中每个像素点的经纬度坐标。
最佳答案
不幸的是,我(目前)找不到比遍历所有像素更好的解决方案。到目前为止,这是我的解决方案:
import glob
import os
import pickle
import sys
import gdal
import geopandas as gpd
import matplotlib
import matplotlib.pyplot as plt
from numba import jit
import numpy as np
from osgeo import osr
import PIL
from PIL import Image, TiffImagePlugin
from shapely.geometry import Point, Polygon, box
import torch
def pixel2coord(img_path, x, y):
"""
Returns latitude/longitude coordinates from pixel x, y coords
Keyword Args:
img_path: Text, path to tif image
x: Pixel x coordinates. For example, if numpy array, this is the column index
y: Pixel y coordinates. For example, if numpy array, this is the row index
"""
# Open tif file
ds = gdal.Open(img_path)
old_cs = osr.SpatialReference()
old_cs.ImportFromWkt(ds.GetProjectionRef())
# create the new coordinate system
# In this case, we'll use WGS 84
# This is necessary becuase Planet Imagery is default in UTM (Zone 15). So we want to convert to latitude/longitude
wgs84_wkt = """
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.01745329251994328,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4326"]]"""
new_cs = osr.SpatialReference()
new_cs.ImportFromWkt(wgs84_wkt)
# create a transform object to convert between coordinate systems
transform = osr.CoordinateTransformation(old_cs,new_cs)
gt = ds.GetGeoTransform()
# GDAL affine transform parameters, According to gdal documentation xoff/yoff are image left corner, a/e are pixel wight/height and b/d is rotation and is zero if image is north up.
xoff, a, b, yoff, d, e = gt
xp = a * x + b * y + xoff
yp = d * x + e * y + yoff
lat_lon = transform.TransformPoint(xp, yp)
xp = lat_lon[0]
yp = lat_lon[1]
return (xp, yp)
def find_img_coordinates(img_array, image_filename):
img_coordinates = np.zeros((img_array.shape[0], img_array.shape[1], 2)).tolist()
for row in range(0, img_array.shape[0]):
for col in range(0, img_array.shape[1]):
img_coordinates[row][col] = Point(pixel2coord(img_path=image_filename, x=col, y=row))
return img_coordinates
def find_image_pixel_lat_lon_coord(image_filenames, output_filename):
"""
Find latitude, longitude coordinates for each pixel in the image
Keyword Args:
image_filenames: A list of paths to tif images
output_filename: A string specifying the output filename of a pickle file to store results
Returns image_coordinates_dict whose keys are filenames and values are an array of the same shape as the image with each element being the latitude/longitude coordinates.
"""
image_coordinates_dict = {}
for image_filename in image_filenames:
print('Processing {}'.format(image_filename))
img = Image.open(image_filename)
img_array = np.array(img)
img_coordinates = find_img_coordinates(img_array=img_array, image_filename=image_filename)
image_coordinates_dict[image_filename] = img_coordinates
with open(os.path.join(DATA_DIR, 'interim', output_filename + '.pkl'), 'wb') as f:
pickle.dump(image_coordinates_dict, f)
return image_coordinates_dict
那些是我的辅助函数。因为这会花费很长时间,所以在 find_image_pixel_lat_lon_coord
中,我将结果保存到字典 image_coordinates_dict
中,我将其写入 pickle 文件以保存结果。
那么我的使用方式是:
# Create a list with all tif imagery
image_filenames = glob.glob(os.path.join(image_path_dir, '*.tif'))
image_coordinates_dict = find_image_pixel_lat_lon_coord(image_filenames, output_filename='image_coordinates')
关于python - 查找 GeoTiff 图像中每个像素的纬度/经度坐标,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/63004971/
我需要将文本放在 中在一个 Div 中,在另一个 Div 中,在另一个 Div 中。所以这是它的样子: #document Change PIN
奇怪的事情发生了。 我有一个基本的 html 代码。 html,头部, body 。(因为我收到了一些反对票,这里是完整的代码) 这是我的CSS: html { backgroun
我正在尝试将 Assets 中的一组图像加载到 UICollectionview 中存在的 ImageView 中,但每当我运行应用程序时它都会显示错误。而且也没有显示图像。 我在ViewDidLoa
我需要根据带参数的 perl 脚本的输出更改一些环境变量。在 tcsh 中,我可以使用别名命令来评估 perl 脚本的输出。 tcsh: alias setsdk 'eval `/localhome/
我使用 Windows 身份验证创建了一个新的 Blazor(服务器端)应用程序,并使用 IIS Express 运行它。它将显示一条消息“Hello Domain\User!”来自右上方的以下 Ra
这是我的方法 void login(Event event);我想知道 Kotlin 中应该如何 最佳答案 在 Kotlin 中通配符运算符是 * 。它指示编译器它是未知的,但一旦知道,就不会有其他类
看下面的代码 for story in book if story.title.length < 140 - var story
我正在尝试用 C 语言学习字符串处理。我写了一个程序,它存储了一些音乐轨道,并帮助用户检查他/她想到的歌曲是否存在于存储的轨道中。这是通过要求用户输入一串字符来完成的。然后程序使用 strstr()
我正在学习 sscanf 并遇到如下格式字符串: sscanf("%[^:]:%[^*=]%*[*=]%n",a,b,&c); 我理解 %[^:] 部分意味着扫描直到遇到 ':' 并将其分配给 a。:
def char_check(x,y): if (str(x) in y or x.find(y) > -1) or (str(y) in x or y.find(x) > -1):
我有一种情况,我想将文本文件中的现有行包含到一个新 block 中。 line 1 line 2 line in block line 3 line 4 应该变成 line 1 line 2 line
我有一个新项目,我正在尝试设置 Django 调试工具栏。首先,我尝试了快速设置,它只涉及将 'debug_toolbar' 添加到我的已安装应用程序列表中。有了这个,当我转到我的根 URL 时,调试
在 Matlab 中,如果我有一个函数 f,例如签名是 f(a,b,c),我可以创建一个只有一个变量 b 的函数,它将使用固定的 a=a1 和 c=c1 调用 f: g = @(b) f(a1, b,
我不明白为什么 ForEach 中的元素之间有多余的垂直间距在 VStack 里面在 ScrollView 里面使用 GeometryReader 时渲染自定义水平分隔线。 Scrol
我想知道,是否有关于何时使用 session 和 cookie 的指南或最佳实践? 什么应该和什么不应该存储在其中?谢谢! 最佳答案 这些文档很好地了解了 session cookie 的安全问题以及
我在 scipy/numpy 中有一个 Nx3 矩阵,我想用它制作一个 3 维条形图,其中 X 轴和 Y 轴由矩阵的第一列和第二列的值、高度确定每个条形的 是矩阵中的第三列,条形的数量由 N 确定。
假设我用两种不同的方式初始化信号量 sem_init(&randomsem,0,1) sem_init(&randomsem,0,0) 现在, sem_wait(&randomsem) 在这两种情况下
我怀疑该值如何存储在“WORD”中,因为 PStr 包含实际输出。? 既然Pstr中存储的是小写到大写的字母,那么在printf中如何将其给出为“WORD”。有人可以吗?解释一下? #include
我有一个 3x3 数组: var my_array = [[0,1,2], [3,4,5], [6,7,8]]; 并想获得它的第一个 2
我意识到您可以使用如下方式轻松检查焦点: var hasFocus = true; $(window).blur(function(){ hasFocus = false; }); $(win
我是一名优秀的程序员,十分优秀!