Вопрос знатокам Global Mapper

[seno]

Понятно.
Продолжаем: правильно понимаю, в таких случаях для вас самое удобное обвести полигоном нужное и обрезать по полигону?
Полигон на глаз рисуется, или как то по другому?

<noindex></noindex>

[DarkDiver]

Нет, полигон вообще не рисуется поля обрезаются в автоматическом режиме по расставленным при привязке в Ozi углам. Tools -> Control Center -> Options -> Cropping -> Automatically Crop Collar понятное дело, что можно нарисовать полигон и экспортировать растр по границам этого полигона, но это опять дополнительная ручная работа, которой хотелось бы избежать. ---------------------- http://john.bdk.com.ru Пешеходный туризм в Приморском крае. Карты Приморского края для GPS Garmin.</noindex>

[seno]

<<<поля обрезаются в автоматическом режиме по расставленным при привязке в Ozi углам>>>
Или я туплю, или Вы чего-то недоговариваете.
Если есть MAP-ы, в которых прописаны градусные координаты углов поля обрезки, в чём проблема?
Недавно именно в такой ситуации (входная 6 градусный Гаус-Крюгер, СК-42, выход Lat/Lon, WGS-84) обрезал разом 2000 листов.
Правильно обрезал, по градусной сетке. Заняло 5 секунд.

Если почему то надо именно по полигонам,
Экспортируем трансформированные в lat/Lon карты в Geotiff. Потом открыть и обрезать по полигонам.

[DarkDiver]

>Если есть MAP-ы, в которых прописаны градусные координаты углов поля обрезки, в чём проблема? В MAP-ах есть следующая информация: 1)точки привязки: Point01,xy, 409, 289,in, deg, 42, 45,N, 130, 40.5,E, grid, , , ,N Point02,xy, 2358, 285,in, deg, 42, 45,N, 130, 45,E, grid, , , ,N Point03,xy, 4522, 281,in, deg, 42, 45,N, 130, 50,E, grid, , , ,N Point04,xy, 410, 3210,in, deg, 42, 40,N, 130, 40.5,E, grid, , , ,N Point05,xy, 2359, 3205,in, deg, 42, 40,N, 130, 45,E, grid, , , ,N Point06,xy, 4528, 3203,in, deg, 42, 40,N, 130, 50,E, grid, , , ,N Point07,xy, 415, 6124,in, deg, 42, 35,N, 130, 40.5,E, grid, , , ,N Point08,xy, 2363, 6122,in, deg, 42, 35,N, 130, 45,E, grid, , , ,N Point09,xy, 4535, 6118,in, deg, 42, 35,N, 130, 50,E, grid, , , ,N 2) Углы: MMPXY,1,410,289 MMPXY,2,4676,280 MMPXY,3,4699,6993 MMPXY,4,417,7002 MMPLL,1, 130.675019, 42.750003 MMPLL,2, 130.839300, 42.750046 MMPLL,3, 130.839566, 42.558348 MMPLL,4, 130.675174, 42.558297 Так вот, как я уже говорил, программа CropGMMinimax, работает не корректно, т.к. режет не по координатам углов (п.2), а по максимально-минимальным координатам точек привязки(п.1). А эти самые точки привязки не всегда находятся в углах растра, более того в углах растра их может вообще не быть в принципе. В результате растр режется не по координатам углов и обрезается полезная часть. GM наоборот, режет правильно - по координатам углов, но экспортирует при этом со своими воображаемыми полями. >Если почему то надо именно по полигонам, По полигонам как раз не надо :), потому как, чтобы резать по полигонам, их (эти самые полигоны) надо создать, а это лишняя ручная работа (растры ведь не номенклатурные, поэтому именно ручная), которой хотелось бы избежать. >Недавно именно в такой ситуации (входная 6 градусный Гаус-Крюгер, СК-42, выход Lat/Lon, WGS-84) обрезал разом 2000 листов.Правильно обрезал, по градусной сетке. Заняло 5 секунд. Вы, вероятно, опять говорите про номенклатурный генштаб или госгицентр, как я уже говорил, с ним проблем нет. ---------------------- http://john.bdk.com.ru Пешеходный туризм в Приморском крае. Карты Приморского края для GPS Garmin.

[seno]

Уф, вроде дошло. Задачка интересна как кроссворд.
Подумаю на досуге.

<noindex></noindex>

[Moskus]

А с какого лешего это ручная работа?
Используя регулярные выражения (PowerGREP+RegexBuddy для не владеющих программированием) данные из описания обрезки в .map-файле легким движением превращаются в GeoJSON-файлик вот такого вида:

{
"type": "FeatureCollection",
"features": [
{ "type": "Feature", "properties": { }, "geometry": { "type": "Polygon", "coordinates": [
[ [ 130.675019, 42.750003 ], [ 130.839300, 42.750046 ], [ 130.839566, 42.558348 ], [ 130.675174, 42.558297 ], [ 130.675019, 42.750003 ] ]
] } }
]
}
который можно открыть и обрезать данные по нему.

К слову, из этого примера ясно, что обрезка не по углам, а по сторонам (по крайней мере, заданных в градусах) для конкретной карты принципиально невозможна (если углы в примере проставлены идеально точно), т.к. у пар точек, находящихся на одной стороне карты, нет равных пар координат, например, Y-координата точек на верхней стороне различается 42.750003 != 42.750046

[DarkDiver]

>Используя регулярные выражения (PowerGREP+RegexBuddy для не владеющих программированием) данные из описания обрезки в .map-файле легким движением превращаются в GeoJSON-файлик вот такого вида... Да, действительно вполне рабочий (хоть и требующий разовой дополнительной работы) вариант решения - создать тем или иным способом обработчик MAP-файлов, который будет генерить полигоны по координатам углов. Спасибо. >К слову, из этого примера ясно, что обрезка не по углам, а по сторонам (по крайней мере, заданных в градусах) для конкретной карты принципиально невозможна (если углы в примере проставлены идеально точно), т.к. у пар точек, находящихся на одной стороне карты, нет равных пар координат, например, Y-координата точек на верхней стороне различается 42.750003 != 42.750046 Это понятно, но допустим, что такая погрешность меня устраивает, да и ее величина значительно меньше размеров воображаемых полей от GM :) <noindex>---------------------- http://john.bdk.com.ru Пешеходный туризм в Приморском крае. Карты Приморского края для GPS Garmin.</noindex>

[DarkDiver]

>Уф, вроде дошло. Задачка интересна как кроссворд.Подумаю на досуге. Ок, в любом случае спасибо за дискуссию. <noindex>---------------------- http://john.bdk.com.ru Пешеходный туризм в Приморском крае. Карты Приморского края для GPS Garmin.</noindex>

[Werefox]

<A HREF="http://sasgis.ru/forum/viewtopic.php?f= ... ></noindex>

[seno]

Работающий вариант с ручным указанием точек обрезки, немного тягомотный, может можно упростить.

Всё в GM
Карта привязана в Гаусса-Крюгера
Выходная проекция Lat/Lon
Создаем две точки в левом верхнем и правом нижнем углах зоны обрезки (Create New Point/Text Feature
Жмем на User Created Feature правой - (Create Area Feature From The Selected Layer Bounds
Выделяем этот полигон.
Идём в Cropping, Crop to Manualy...
Жмём кнопку Initialize Bouds of Selectyed Area feature
ОК, экспорт, готово