지오 서버에서 컬러 밸런스 및 고해상도 이미지 게시

온라인 응용 프로그램의 기본 맵으로 사용하려는 많은 양의 RapidEye 이미지 (+ 300Gb)가 있습니다. 지금까지 나는 지오 서버에서 "스테로이드의 지오 서버 (Geoserver on steroids)"라는 프리젠 테이션에서 보여준 단계에 따라 지오 서버에서 ImagePyramid로 제공했습니다. 문제는 적절한 색상 균형을 수행 할 수 없었기 때문에 전체 데이터 세트를 볼 때 모자이크가 매우 나빠 보이지만 확대하면 더 좋아집니다 (그림은 아래 참조).

내 프로세스는 다음과 같습니다

gdal을 사용하여 모든 이미지를 8 비트, epsg 4326, 3 차 보간으로 변환
모든 이미지로 gdal 가상 래스터 생성 (gdalbuildvrt)
이미지 피라미드 (압축 및 지오메트리 경작이있는 gdal_retile)를 생성하고 지오 서버에 게시 (레이어 스타일에서 히스토그램 스트레치 사용)

프로세스를 개선하고 색상 균형을 개선하는 방법에 대한 팁이 있습니까?

여기에 이미지 설명을 입력하십시오

— 다니엘
소스

답변:

나는 그것을 믿습니다.

8 비트로 변환하는 동안 각 RapidEye 이미지에서 평균 + -2 표준 편차 히스토그램 스트레치를 강제했습니다.

파이썬 스크립트를 사용하여 이미지 최소, 최대, 평균 및 SD를 식별했습니다. 그런 다음 평균-2SD (또는 더 높은 이미지 최소값) 값을 2로 설정하고 평균 + 2SD (또는 최대) 값을 254로 설정했습니다. 그리고 안전을 위해 이미지의 원래 0 값은 1로 설정되었습니다. 2 ^ 16은 255로 설정되었습니다. NoData는 0으로 설정되었습니다.

전달 함수 벡터는 다음과 같습니다.

original image values = [0, mean-2SD, mean+2SD, 65536]
rescaled values = [1 2, 254, 255]

다음은 일부 화면 캡처입니다. 지금 해결해야 할 문제는

이미지 주위의 검은 색 테두리를 제거하십시오.
확대하여 더 나은 이미지를 얻을 수 있습니다. 이제 잘못된 채널로 조정 된 오래된 TV처럼 보입니다.

여기에 이미지 설명을 입력하십시오

편집
이제 코드를 GitHub에 배치했습니다. 내가 사용한지 오래되었습니다. 코드가 약간 지저분하며 저장소도 마찬가지입니다. 그러나 여전히 작동합니다.
https://github.com/dvictori/rapideye_hist2sd

— 다니엘
소스

브라보 선생님 이것을 위해 10 점 만점에 10 점을 얻습니다. 두 번째 문제 (확대 할 때 더 잘 보임)에 관해서는 LandSat 또는 ASTER와 같이 낮은 이미지를 사용하는 것이 간단 할 수 있습니까? 축소하면 더 거친 이미지가 더 좋아 보인다고 생각합니다. .. 첫 번째 관심사에 관해서는, 나는 여전히 그 너트를 깨려고 노력하고 있습니다. 공유하고자하는 경우 py 스크립트를 포함 하여이 작업을 수행 한 방법에 대한 포괄적 인 글을 쓰고 싶습니다. 어디에서나 접근 방식을 제시 할 계획이 있습니까?

— elrobis 2016 년

나는 이것을 발표하거나 글을 쓰려고 계획하지 않았다. 그러나 이제는 당신이 아는 사람을 언급했습니다. 이를 제시하기위한 적절한 미디어에 대한 제안 사항이 있습니까? py 스크립트는 dropbox.com/s/1hfobfp9ymtku2n/rapideye_hist2sd.py 에서 다운로드 할 수 있지만주의해야합니다. 나는 나 자신이 미친 프로그래머라고 생각하며 내 프로그램이 많은 개선점을 사용할 수 있다고 확신합니다.

— Daniel

위의 의견에서 다운로드 링크를 변경했습니다. 내가 사용한 파이썬 코드를보고 싶다면 goo.gl/ePEc7G

— Daniel

감사합니다 @Daniel. 나는 그런 것들에 대한 방법을 게시하는 블로그를 유지 합니다 . 쉽지 않은 웹을 발견하고 웹에서 많은 도움을 얻을 수없는 경우 몇 가지 이유로 블로그 게시물을 작성합니다. 먼저 언젠가 작업을 반복해야한다고 가정하므로 단계를 게시물로 통합하면 나중에 쉽게 검토 할 수 있습니다. 둘째, 나는 누군가가 비슷한 것을하고 싶을 확률이 항상 있다고 생각하며 게시물에서 이익을 얻을 수 있습니다. 당신이 여기서 한 것은 잠시 동안 유일한 게시물이더라도 블로그에 가치가 있습니다. :)

— elrobis 2016

이것은 놀라운 아이디어이며 8 비트에서 8 비트로 이동하는 더 작은 프레임, 3 밴드 카메라에 맞게 조정하려고하지만 구현에 문제가 있습니다. 입력 이미지 디렉토리를이 스크립트로 어떻게 지시해야합니까? 그들은 어디에서 왔습니까? 이 모든 것이 지오 서버를 통해 실행되어야한다고 생각하지만이를 분리하고 독립 실행 형으로 실행할 수 있습니까?

— Wes

다니엘,이 이미지들은 계절마다 다른가요? 아니면 시간? 계절이 다르면 색상 균형을 잡는 것이 꽤 어려울 수 있습니다. 그러나 시간이 다르면 태양 각도에 대한 수정을 적용하면 눈에 띄게 도움이 될 수 있습니다. 태양 각도에 대한 1 차 근사법은 픽셀에 1.0 / cos (angle_of_sun_off_directly_above)를 곱하는 것입니다. 따라서 태양이 바로 위에 있으면 조정이 없으며 새벽 / 황혼에 가까워 질 때 ... 무한대까지 증가합니다.

과거 에 겹쳐진 영역의 장면 사이에서 히스토그램 일치 를 사용하여 큰 모자이크를 모으는 데 나쁜 결과가 나타났습니다. 모자이크에서 이상한 표류 효과가 있기 때문입니다. 더 유용한 접근 방식은 기본 색상 대상 이미지 (아마도 매력적인 지역 모자이크)와 일치하는 일종의 히스토그램 일 수 있다고 생각합니다. 이 문제를 해결하는 방법에도 관심이 있습니다. MapBox의 Charlie Loyd도 생각이있을 수 있습니다.

— 프랭크 워머 담
소스

프랭크, 이미지는 구름이 적을 때 주로 건기입니다. 그러나 8 월, 11 월 등 다른 달의 이미지가 있습니다. 이제 이미지를 작은 블록으로 나누고 각각에 대해 모자이크를 만들려고합니다. 그것이 어떻게 나오는지 알려줄 것입니다. 나는 또한 CDF 일치를 수행하기 위해 파이썬 스크립트를 작성했으며 적은 수의 이미지가있을 때 잘 작동하지만 (20) 큰 수에 도달하면 이미지 중 하나가 내 외부의 값을 갖도록 바인딩됩니다. 참조 이미지 범위를 설정하면 작동이 중지됩니다. 참조 이미지를 선택하는 더 좋은 방법을 생각해야 할 수도 있습니다.

— Daniel

때로는 색상 균형 잡힌 모자이크를 만들기 위해 히스토그램 일치와 함께 OSSIM 이미지 링커를 사용했습니다. 이미지 링커 튜토리얼은 http://download.osgeo.org/ossim/docs/pdfs/ImageLinker_Tutorial.pdf 예제를 제공합니다. 그러나 이미지 링커는 적극적으로 유지 관리되지 않으며 더 이상 작동하는지 알 수 없습니다. Ossim 지오 셀은 현재 프로그램이지만 이에 대한 사용자 문서는 많지 않습니다. 좋은 결과를 얻기 위해 OpenEV로 Landsat 장면을 하나씩 사전 처리하여 이미지가 시각적으로 동일한 모양으로 보이도록 LUT 스트레치를 적용했습니다. 그런 다음 OSSIM이 최종 매칭을하는 것이 더 쉬웠습니다. 최종 제품은 내가 ER Mapper 모자이크 유틸리티로 만든 것만 큼 좋아 보였습니다.

— 사용자 30184
소스