如何从GeoTiff文件中获取投影坐标以及这些坐标处的实际像素值,然后将其保存到numpy数组中?我有arsenci020l.tif文件,其坐标以米为单位。以下是我在上面运行的gdalinfo的摘要输出。
~$ gdalinfo arsenci020l.tif
Driver: GTiff/GeoTIFF
Files: arsenci020l.tif
arsenci020l.tfw
Size is 10366, 7273
Coordinate System is:
PROJCS["Lambert Azimuthal Equal Area projection with arbitrary plane grid; projection center 100.0 degrees W, 45.0 degrees N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Lambert_Azimuthal_Equal_Area"],
PARAMETER["latitude_of_center",45],
PARAMETER["longitude_of_center",-100],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (-6086629.000000000000000,4488761.000000000000000)
Pixel Size = (1000.000000000000000,-1000.000000000000000)
...
关于从tiff(从GeoTIFF文件获取纬度和经度)获取纬度/经度坐标,这里存在类似的问题,答案显示了如何仅获取左上角的x和y像素坐标。我需要获取所有投影的像素坐标,以及获取像素值并将其保存在numpy数组中。我该怎么做?
您想要10366×7273 =超过7千5百万点吗?
—
Mike T
@MikeT我想是的,我真的不知道如何解决我要解决的问题的更好的解决方案:我需要从该数据集中找到最接近的像素坐标到US块的每个质心,然后分配通过搜索,我意识到cKDTree查询将帮助我进行最近邻居搜索。该算法的Python函数要求以“ numpy”数组的形式查询“树”。以生成“树”在此数据集中所有像素坐标中,我需要以某种方式存储所有像素坐标。如果您有更好的解决方案,请告诉我!
—
irakhman 2015年