Форумы GIS-Lab.info

Здравствуйте. С наступающим! Подскажите, пжлст, где и как грамотно сгенерировать PROJ.4 - строку по ASCII-файлу, в котором в каждой строке есть и X,Y декартовых картографических и пара географических Longitude/Latitude. Строк "цельный сонм" - все в домейне от 50 с.ш. до 90 с.ш. Есть такой сервис в какой-нибудь free-ГИС, или через GDAL/OGR в Python'е , или в R ?

Нужна проекция, потом определяются параметры. На форуме много раз обсуждали, в основном для проекции Гаусса Крюгера

gamm писал(а): ↑

31 дек 2017, 13:35

Нужна проекция, потом определяются параметры. На форуме много раз обсуждали, в основном для проекции Гаусса Крюгера

есть только этот ASCII... да и "картографические" - какие-то условные... целочисленные индексы i- да-j ("Фортрановские" )... поищу в темах... спасибо, gamm

nickleb писал(а): ↑

31 дек 2017, 13:53

да и "картографические" - какие-то условные... целочисленные индексы i- да-j

это все равно, главное угадать, кто на Север, кто на Восток.

Задается набор преобразований (проекция, или проекция + поворот для местных СК), пишется сумма квадратов невязок, и запускается оптимизатор optim(). Я привык сначала максимизировать корреляцию/минимизировать дисперсию невязок (что не зависит от сдвига), а потом вычислять сдвиг.

https://github.com/Ariki/guessproj
Но тип проекции нужно знать (или перебором действовать). Если в качестве координат целочисленные индексы, начните с проекции Меркатора, подбирая параметр R. Мне недавно попадалась самодельная ГИС, в которой координаты хранились в пикселях карты Web Mercator (масштаб соответствовал одному из уровней стандартной тайловой сетки). Хотя лучше, наверно, сначала нарисовать точки на плоскости: так проще угадать проекцию.

Ariki, спасибо!