GRASS에서 LANDSAT 8을 어떻게 선명하게 하는가?


9

QGIS를 매핑하는 데 꽤 오랜 시간 사용했지만, pansharpening이 QGIS에서 직접 실행 가능하지 않기 때문에 GRASS를 파고 들었습니다. 구조를 다소 이해 한 후에 i.fusion.brovey로 몇 가지 시험을 수행했습니다. 15m 해상도의 래스터를 생성 할 수 있었지만 QGIS에로드 할 때 원래 30m ​​래스터에있는 멋진 색상 대비를 달성 / 복구 할 수 없습니다.

에서 이전 포스트 사람이 확인하고 탐사를위한 다른 기능을 사용하는 제안 : 'i.his.rgb'(i.landsat.toar 사용 후, 내가 사용하지 않는). 그는 i.his.rgb를 사용하기 전에 min / max를 1-255로 변환하기 시작합니다 ... 원래 래스터의 범위는 0-65535입니다. 그것은 나를 위해 작동하지 않았다 ...

며칠 동안 팬 샤프를 수행하기 위해 고군분투하고 갇힌 느낌이 들었습니다 ... 왜 다른 사람들이 제안한 것처럼 i.pansharp를 Landsat에 적용해야합니까? Landsat 8 이미지의 팬 샤프닝을 실행하는보다 간단한 워크 플로우를 제안 할 수 있습니까? 조언을 기대하고 미리 감사드립니다 ...

답변:


7

i.landsat.rgb 살펴보기 -팬 샤프닝을 실행하기 전에 LANDSAT 이미지의 색상 자동 균형 조정을 수행합니다. i.landsat.toar 을 사용하여 개별 채널의 디지털 숫자를 대기의 방사 또는 반사로 변환하는 것을 고려할 수도 있습니다 . 참조 http://grasswiki.osgeo.org/wiki/LANDSAT

BTW : 원본 LANDSAT 데이터의 범위는 0-65535이며 8 비트 채널 (0-255)로 전달되기 때문에 이상합니다. 먼저 입력 데이터를 확인하십시오.

업데이트 : 2013 년 10 월

현재 GRASS GIS 7을 사용한 단순화 된 절차는 GRASS GIS 7의 Landsat 8 데이터 처리 페이지 : RGB 컴포지트 및 팬 샤프닝 페이지를 참조하십시오.여기에 이미지 설명을 입력하십시오


나는 원래 neteler.org에서 링크하는 기사를 보았지만 지시 사항은 전혀 작동하지 않았습니다 (로 생성 된 래스터의 부동 소수점 형식으로 인해 gis.stackexchange.comi.landsat.toar 에서 다음과 같은 조언을 얻었습니다. / a / 39239 / 9857
Juan

3
조언은 수년 전에 제공되었으며 최신 GRASS GIS 코드를 반영하지 않습니다. BTW : 정말 멋진 팬 샤프닝 모듈은 새로운 애드온입니다. i.fusion.hpf-고역 통과 필터 추가 기술을 기반으로 한 고해상도 팬 크로마 틱 및 저해상도 다중 스펙트럼 데이터 융합, grass.osgeo.org/grass70/manuals/ addons / i.fusion.hpf.html
markusN

4

Michal과 Markus에게 감사드립니다. 마지막으로 표시된 기능 brov로 pansharpening을 할 수있었습니다. GRASS에서 수행 한 후 RGB 래스터를 데이터 유형 Uint16으로 gdal GTiff로 내보냈습니다. Markus, 저는 여러 개의 Landsat 8 이미지 (11 개 밴드 각각의 강성)를 다운로드하여 사용하여 자연스러운 색과 거짓 색으로 멋진 시각적 결과를 얻었습니다. 내 지식으로, 원래 landsat tiffs는 0에서 65535 사이의 값을 가진 Uint16으로 제공됩니다 (또는 무엇을 놓치고 있습니까?) ... GRASS LANDSAT 링크에 감사드립니다 ...


3
예, Landsat 8 이미지는 Uint16 유형으로 저장됩니다.
webrian

3

다음과 같은 기사의 공식을 따르는 경우 FIHS 또는 Brovery Transform의 팬 샤프닝은 매우 간단합니다. http://www.sensorsportal.com/HTML/DIGEST/august_2014/Vol_177/P_RP_0183.pdf

GRASS GIS를 사용하고 Landsat 8을 pansharpening한다고 가정합니다.

팬 8 tif를 사용하여 지역을 설정하십시오.이 경우 15m 해상도를 가져야합니다.

가중치가없는 Brovery의 경우 각각 파랑, 녹색 및 빨강에 대해 맵 계산기에 입력하십시오.

band8 / (band2 + band3 + band4) / 3 * band2

band8 / (band2 + band3 + band4) / 3 * band3

band8 / (band2 + band3 + band4) / 3 * band4

위의 각 출력에 대해 r.colors를 실행하고 막대 그래프 이퀄라이제이션 및 색상 표를 회색으로 설정하십시오.

r.composite를 실행하고 적절한 밴드를 선택하십시오 (빨간색 = 4, 녹색 = 3, 파란색 = 2)

이미지를 확인하십시오.

분명히이 과정은 위 기사에서 논의 된 것처럼 대기 보정, 색상 스트레치 및 밴드 가중치를 더 고려하여 수행 할 수 있습니다.


2

Monteverdi 를 사용하여 팬 샤프닝을 수행 할 수도 있습니다 . 이 도구는 일반 GeoTIFF를 입력으로 사용할 수 있지만 (별도의 다중 스펙트럼 대역을 하나의 다중 대역 래스터로 병합해야 함) 성능이 좋습니다.


2

이 문제로 어려움을 겪고 나서 SCP를 사용하여 QGIS (2.18.21은 내 버전 임)에서 팬 샤프 할 수있는 간단한 솔루션을 찾았습니다. 이 방법을 통합 한 SCP 업데이트는 https://fromgistors.blogspot.com/2015/07/major-update-semi-automatic-44.html에 설명되어 있습니다 .

이것은 쉬운 일이며 매력처럼 달렸습니다! 나는 최종 출력을 확인했으며 픽셀은 15m에서 해상도를 측정했습니다-panchromatic 8 밴드와 일치합니다. 궁금한 점이 있으면 매뉴얼에 설명 된대로 Brovey 변환 방법을 사용하십시오 ( https://media.readthedocs.org/pdf/semiautomaticclassificationmanual/latest/semiautomaticclassificationmanual.pdf 참조 ).

SCP pansharpen 도구의 1 단계 SCP pansharpen의 3-4 단계



당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.