gornak писал(а): желательно иметь метаданные на ваши растры
Действительно, растр можно создать программно, но я не владею этой технологией.
Какие метаданные для этого нужны?
Суть моих исследований (если вдруг Вам интересно) заключается в следующем:
Гипотеза:
1. засоление почв в сложной рельефной обстановке происходит при определенных климатических условиях, и только тогда, когда имеет место "выпотной" водный режим почв, т.е. когда почва "потеет" и на ее поверхности кристаллизуются соли. Выпотной водный режим возникает при высоком прогреве поверхности, благодаря восходящему току влаги по почвенным капиллярам из грунтовых вод.
2. Накопление солей в верхних горизонтах почв и на их поверхности в горных условиях происходит при аридном и семиаридном климате. Климат в горных условиях зависит не только от широты, но и в значительной степени от высоты.
3. Таким образом, должен существовать определенный климатический потенциал засоления территории, определяющийся ее широтой и высотой.
Для проверки гипотезы я провел визуальный анализ космоснимков на предмет поиска засоленных участков, в том числе с наибольшими высотными отметками для конкретных широт. В итоге я получил выборку точек с координатами и высотами.
Интересным оказалось то, что в северных широтах засоление проявлялось только на низких высотах (до 600 м), тогда как на юге оно отмечалось значительно выше (до 1700 и более).
Положив эти точки на график, где по X - широта, а Y - высота, я отметил, что верхний предел засоления (Hpz) действительно достаточно хорошо описывается экспоненциальным трендом с высоким значением коэффициента детерминации и решил эту закономерность отобразить картографически.
Формула для бассейна р. Селенги имеет следующий вид: Hpz = 5000000 * exp (-0.17 * широта). Hpz - предельная высота для конкретной широты, где потенциально может произойти засоление. Следует отметить, что на других территориях формула скорее всего будет другой. Это проявляется также в самой западной части бассейна, где климат локально меняется крупными горными сооружениями с усилением континентальности. Некоторые изменения закономерности отмечаются также и за пределами бассейна в южном и восточном направлениях, где меняется характер геоморфологии и расчленения рельефа.
Отклонения от выявленной закономерности нивелируются включением в расчеты суммы эффективных осадков в теплый период года. Но до конца я с этим пока не разобрался. Если кому интересно, можем порассуждать на эту тему и даже сделать совместную научную публикацию.
Далее из Hpz мы вычитаем реальную высоту (Н) и получаем для каждой конкретной точки положительную или отрицательную цифру. Положительная показывает, что климатический потенциал засоления в этой точке присутствует, отрицательная - нет. Фактическое значение указывает на силу проявления выпотного режима, чем выше значение - тем потенциал засоления сильнее. Для отрицательных значений, напротив, увеличение значения, если "-" показывает смену выпотного режима на непромывной и промывной, когда соли уходят в глубь профиля или вообще из него вымываются.
Таким образом мы можем получить карту климатического потенциала засоления на определенный климатически замкнутый однородный участок земной поверхности.
Для этого мне необходим растр (или сетку растров) с широтными координатами. Возможно для каких-то других территорий могут потребоваться координаты долготы. Этот растр можно использовать в дальнейших расчетах с помощью калькулятора. Хорошо было бы, если эти растры по своим метаданным совпадали с растрами SRTM для упрощения расчетов.
Естественно, что засоление зависит от массы прочих факторов, в числе которых можно обозначить наличие источников солей (выклинивание глубинных минерализованных вод, засоленные породы и т.п.), условия стока и уровень залегания грунтовых вод, наличие мерзлотного или иного водоупора, геоморфология поверхности, экспозиция и крутизна склона и т.д. и т.п. Но все эти параметры в первую очередь контролируются климатическим потенциалом. Другими словами, даже при наличии источника засоления, например, высокоминерализованного водного источника (ключа) или близкозалегающих засоленных пород, соли не накопятся на поверхности, если
нет климатического потенциала.
Кстати, я не очень опытный юзер ГИСов, но, полагаю, что для гидрологических и геоморфологических параметров также можно рассчитать их потенциалы засоления и сочленить их с климатическим потенциалом.
Практическим выхлопом, например, может служить оценка потенциала вторичного засоления при проектировании оросительных мелиоративных систем и др. Подобным образом также можно попытаться рассчитать климатические потенциалы распределения других почвенных показателей, а также отдельных видов или сообществ растений и других биологических объектов.