Геоинструменты создают CRS из полигона и измерений

У меня есть некоторые данные сетки (массив 2d), с которыми связаны некоторые геопространственные метаданные. Границы находятся в многоугольнике jts с координатами долготы WSG84. Многоугольник представляет собой прямоугольник, но не выровнен по оси (не ограничивающий прямоугольник).

Я пытаюсь создать объект Geotools CoordinateReferenceSystem, описывающий это покрытие. Я хочу использовать CRS для повторного проецирования покрытия, но мне также нужно использовать CRS как есть, поэтому пропускать этап создания CRS нехорошо.

Вкратце: как перейти от полигона WSG84, dx, dy, nx, ny к объекту CoordinateReferenceSystem, который описывает сетку с помощью геоинструментов?

Заранее спасибо за помощь


person Brian Clements    schedule 12.05.2011    source источник


Ответы (2)


Итак, вы знаете латлонский прямоугольник, в который вписываются ваши координаты? Похоже, вы могли бы сделать аффинное преобразование, чтобы преобразовать свои координаты в широту/долготу; а дальше нормально?

Поэтому сначала обработайте свою геометрию, чтобы «выровнять» ее с прямоугольником; тогда они будут в WGS84, и вы сможете жить дальше.

Сведения об аффинном преобразовании см. в документах:

Если ваши данные на самом деле представляют собой сетку (как в сетке), вы можете посмотреть на код изображения мозаичности; он поддерживает запись сетки во всех направлениях и сшивает их в цельный результат.

Обратите внимание, что для сеток процесс обратный; а не преобразовывать данные на основе запроса; запрос преобразуется для согласования с данными (чтобы он мог выбрать соответствующие значения из растрового изображения).

person Jody Garnett    schedule 24.05.2011
comment
Мне удалось найти преобразование для перепроецирования сетки в WGS84. Первоначально я хотел рекламировать CRS необработанных данных, но это было невозможно без дополнительной информации. В любом случае для моего приложения WGS84 более полезен. - person Brian Clements; 29.05.2011

Если ваши координаты WGS86, вы должны использовать EPSG:4326 в качестве CRS. Однако, судя по звуку, ваша сетка находится в какой-то другой проекции, и в этом случае вам нужно будет перепроецировать ограничивающий прямоугольник в эту проекцию, а затем создать сетку в правильной проекции.

Если вы не знаете, в какой проекции находятся данные сетки, я подозреваю, что вы облажались.

person Ian Turton    schedule 19.05.2011