Геоинформационные системы (ГИС) и Дистанционное зондирование Земли
https://gis-lab.info/forum/
gimran, если бы вы написали статью по этому поводу — было бы просто отлично.gimran писал(а):Делал pan sharpening в Грассе. Если надо — могу кратко описать.
Как видите, тема актуальна.
Вечером проделаю что Вы советуете. У меня стоит версия GRASS GIS 6.4.3RC2. Что касается вылета программ, то они просто вылетают без сообщений. Молча закрываются.gimran писал(а):Скачайте и установите OSGeo4W http://gis-lab.info/qa/qgis-osgeo4w.html
В поставке точно есть Grass 7.0.
Как рушатся именно? Программа что-то пишет?
По этому поводу, пож-та, создайте отдельную тему на форуме GRASS с описанием, что и как вылетает (скриншоты не помешают).kobo35 писал(а):У меня стоит версия GRASS GIS 6.4.3RC2. Что касается вылета программ, то они просто вылетают без сообщений. Молча закрываются.
Попробуйте файл подкачки увеличить, возможно, что из-за переполнения памяти.Что касается вылета программ, то они просто вылетают без сообщений. Молча закрываются.
А кто сказал, что он снимкозависимый? Там просто формула зашита, в которую можно подставлять любые растры (даже и не снимки вообще). В чем выражается, что снимки не обрабатываются?nadiopt писал(а):товарищи, а почему пан-шарпенинг вообще снимкозависимый? мы вот тоже недавно столкнулись, что в спецмодуле грассе почему-то не все снимки обрабатываются.
Код: Выделить всё
DN_b1
DN_fused = ------------------------------- * DN_pan
DN_b1 + DN_b2 + DN_b3
У меня вопрос. Не сочтите за труд подскажите мне, несведущему, в приведенной формуле DN это разрешение в пикселях? И если не сложно направьте, где можно почитать по растровой арифметике, чтобы не задавать глупых вопросов, а самому прочитать и разобраться как это все работает. Я pansharpening выполнял до сих пор используя photoshop,преобразуя rgb в lab, увеличивая размер изображения цветовых каналов до параметров панхроматического, а затем подменял панхроматическим яркостный канал, объединял и преобразовывал в RGB. Но способ этот, мне думается не самый правильный. Я не пытаюсь получить от вас готовых решений, прошу лишь подсказать, где я могу почерпнуть нужные знания.KolesovDmitry писал(а):По существу вопроса. Если посмотреть справку по модулю i.fusion.brovey, то видно, что для того, чтобы воспользоваться идеей и произвести транформацию, не нужно каких-то специальных модулей. Достаточно простого калькулятора растров, например, для первого канала формула будет выглядеть так:Аналогично для других каналов.Код: Выделить всё
DN_b1 DN_fused = ------------------------------- * DN_pan DN_b1 + DN_b2 + DN_b3
DN -- это значения яркостей пикселей. Так, под DN_b1 понимается первый канал снимка (нумерация условна и зависит от того, какие снимки вы используете для получения композита).kobo35 писал(а): У меня вопрос. Не сочтите за труд подскажите мне, несведущему, в приведенной формуле DN это разрешение в пикселях? И если не сложно направьте, где можно почитать по растровой арифметике, чтобы не задавать глупых вопросов, а самому прочитать и разобраться как это все работает.
Код: Выделить всё
g.region rast=Bpan
r.mapcalc "B1.f = 1.0 * Bpan * B10/(B10 + B20 + B30)"
r.mapcalc "B2.f = 1.0 * Bpan * B20/(B10 + B20 + B30)"
r.mapcalc "B3.f = 1.0 * Bpan * B30/(B10 + B20 + B30)"