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,500 만 포인트 이상을 원하십니까?
—
Mike T
@MikeT 그렇게 생각합니다. 나는 해결하려는 문제에 접근하는 방법에 대한 더 나은 해결책을 모른다.이 데이터 세트에서 미국 블록의 각 중심까지 가장 가까운 픽셀 좌표를 찾은 다음 검색을 통해 cKDTree 쿼리가 가장 가까운 이웃 검색에 도움이된다는 것을 깨달았습니다. 알고리즘의 Python 함수는 numpy 배열로 쿼리하기 위해 "트리"를 요청합니다. 이 데이터 세트의 모든 픽셀 좌표에서 어떻게 든 저장해야합니다. 더 나은 솔루션이 있으면 알려주십시오!
—
irakhman