Microsoft MVP성태의 닷넷 이야기
글쓴 사람
정성태 (techsharer at outlook.com)
홈페이지
첨부 파일
 

Python - 공개 패키지를 이용한 위성 이미지 생성 (pystac_client, odc.stac)

아래와 같은 트윗이 있군요. ^^



쉽게 테스트하기 위해 Jupyer Notebook에서 3개의 패키지를 설치한 후,

!pip install pystac_client
!pip install matplotlib
!pip install odc.stac

이렇게 코딩하면 됩니다.

from pystac_client import Client
from odc.stac import load

client = Client.open("http://earth-search.aws.element84.com/v1")
collection = "sentinel-2-l2a"
tas_box = [146.5, -43.6, 146.7, -43.4]
search = client.search(collections=[collection], bbox=tas_box, datetime="2023-12")

data = load(search.items(), bbox=tas_box, groupby="solar_day", chunks={})
data[["red", "green", "blue"]].isel(time=2).to_array().plot.imshow(robust=True)

그 시점의 호출 부하에 따라 다르겠지만 제 경우에 8초 정도 걸려서야 결과가 나왔습니다.




실습을 위해 순천만 지역을 Google Earth로 방문한 후,

https://earth.google.com/web/@34.93009968,127.50141682,3.33634054a,3009.47610511d,35y,0.92573167h,5.8516359t,0r/data=OgMKATA

대충 위치를 잡았더니 위/경도 데이터가 이렇게 나옵니다.

[좌-하] 34°55'21"N 127°29'30"E
[우-상] 34°56'08"N 127°31'10"E

하지만 이 좌표를 이전의 파이썬 소스코드에 적용하면 안 됩니다. 왜냐하면, 저 구글 어스가 보여주는 좌표계는 DMS (degrees, minutes, seconds) 체계이고, pystac_client 호출 시 넘겨주는 좌표계는 DD (decimal degrees)이기 때문입니다. 따라서, 다음과 같은 페이지를 이용해, (또는 코드로),

GPS Coordinates
; https://www.gps-coordinates.net/

구글 어스의 DMS 값을 DD 값으로 변경해,

[좌-하] 127.49, 34.92
[우-상] 127.51, 34.93

다음과 같이 바꿔서 호출하면,

tas_box = [127.49, 34.92, 127.51, 34.93]

아래와 같은 실행 결과가 나옵니다. ^^

pystac_client_1.png

보는 바와 같이, 이런 유의 코드 실행은 멋있긴 하지만 정작 써먹으려고 하면 (속도도 무척 느려서) 현실성 있게 사용하기 위한 시나리오가 딱히 떠오르진 않습니다. (사업적 마인드가 없어서 그럴 수도 있지만. ^^;)




[이 글에 대해서 여러분들과 의견을 공유하고 싶습니다. 틀리거나 미흡한 부분 또는 의문 사항이 있으시면 언제든 댓글 남겨주십시오.]







[최초 등록일: ]
[최종 수정일: 1/8/2024]

Creative Commons License
이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-변경금지 2.0 대한민국 라이센스에 따라 이용하실 수 있습니다.
by SeongTae Jeong, mailto:techsharer at outlook.com

비밀번호

댓글 작성자
 




1  2  3  4  [5]  6  7  8  9  10  11  12  13  14  15  ...
NoWriterDateCnt.TitleFile(s)
13552정성태2/13/20242703닷넷: 2214. windbg - Monitor.Enter의 thin lock과 fat lock
13551정성태2/12/20243016닷넷: 2213. ASP.NET/Core 웹 응용 프로그램 - 2차 스레드의 예외로 인한 비정상 종료
13550정성태2/11/20243237Windows: 256. C# - Server socket이 닫히면 Accept 시켰던 자식 소켓이 닫힐까요?
13549정성태2/3/20243527개발 환경 구성: 706. C# - 컨테이너에서 실행하기 위한 (소켓) 콘솔 프로젝트 구성
13548정성태2/1/20243229개발 환경 구성: 705. "Docker Desktop for Windows" - ASP.NET Core 응용 프로그램의 소켓 주소 바인딩(IPv4/IPv6 loopback, Any)
13547정성태1/31/20242986개발 환경 구성: 704. Visual Studio - .NET 8 프로젝트부터 dockerfile에 추가된 "USER app" 설정
13546정성태1/30/20242817Windows: 255. (디버거의 영향 등으로) 대상 프로세스가 멈추면 Socket KeepAlive로 연결이 끊길까요?
13545정성태1/30/20242632닷넷: 2212. ASP.NET Core - 우선순위에 따른 HTTP/HTTPS 호스트:포트 바인딩 방법
13544정성태1/30/20242720오류 유형: 894. Microsoft.Data.SqlClient - Could not load file or assembly 'System.Security.Permissions, ...'
13543정성태1/30/20242762Windows: 254. Windows - 기본 사용 중인 5357 포트 비활성화는 방법
13542정성태1/30/20242651오류 유형: 893. Visual Studio - Web Application을 실행하지 못하는 IISExpress - 두 번째 이야기
13541정성태1/29/20242877VS.NET IDE: 188. launchSettings.json의 useSSL 옵션
13540정성태1/29/20243060Linux: 69. 리눅스 - "Docker Desktop for Windows" Container 환경에서 IPv6 Loopback Address 바인딩 오류
13539정성태1/26/20242986개발 환경 구성: 703. Visual Studio - launchSettings.json을 이용한 HTTP/HTTPS 포트 바인딩
13538정성태1/25/20243332닷넷: 2211. C# - NonGC(FOH) 영역에 .NET 개체를 생성파일 다운로드1
13537정성태1/24/20243421닷넷: 2210. C# - Native 메모리에 .NET 개체를 생성파일 다운로드1
13536정성태1/23/20243405닷넷: 2209. .NET 8 - NonGC Heap / FOH (Frozen Object Heap) [1]
13535정성태1/22/20243371닷넷: 2208. C# - GCHandle 구조체의 메모리 분석
13534정성태1/21/20243152닷넷: 2207. C# - SQL Server DB를 bacpac으로 Export/Import파일 다운로드1
13533정성태1/18/20243183닷넷: 2206. C# - TCP KeepAlive의 서버 측 구현파일 다운로드1
13532정성태1/17/20243212닷넷: 2205. C# - SuperSimpleTcp 사용 시 주의할 점파일 다운로드1
13531정성태1/16/20243218닷넷: 2204. C# - TCP KeepAlive에 새로 추가된 Retry 옵션파일 다운로드1
13530정성태1/15/20242842닷넷: 2203. C# - Python과의 AES 암호화 연동파일 다운로드1
13529정성태1/15/20242863닷넷: 2202. C# - PublishAot의 glibc에 대한 정적 링킹하는 방법
13528정성태1/14/20243031Linux: 68. busybox 컨테이너에서 실행 가능한 C++, Go 프로그램 빌드
13527정성태1/14/20242953오류 유형: 892. Visual Studio - Failed to launch debug adapter. Additional information may be available in the output window.
1  2  3  4  [5]  6  7  8  9  10  11  12  13  14  15  ...