Перетворення між wgs84 та nad83

Це питання для новачків, що стосується перетворення між wgs84 та nad83. У мене є аерофотозйомки, які я намагаюся обчислити довгі координати якомога ближче. Існує кілька факторів, які роблять складні типові виправлення, наприклад, сумнівні параметри, ... тому виникають деякі спроби та помилки.

Таким чином, мені потрібен спосіб перевірити свій прогрес, коли я вношу виправлення. Спочатку просто перевірити наявність Google Earth було добре, але в цей момент різниця між моїми даними (nad83) і Google Earth (wgs84) стала обмежувальним фактором. Тому мені потрібно перетворити nad83 в wgs84.

Я використовую rgdal в R та dotspatial в C # для того, щоб робити прогнози, але я думаю, що зміна між nad83 та wgs84 не вважається "проекцією" (я прав?) Спроба викликати проект із "+ proj = NAD83 ellps = WGS84" просто видає помилку.

Я добре вписую зміни координат 7 параметрів, за винятком того, що я не впевнений, як перетворити будь-яку в (X, Y, Z) координати (що, як я вважаю, означає перетворення з довгої лати на поверхню еліпсоїдів, на якій вони засновані ). Наявні питання щодо nad83 та wgs84, здається, стосуються того, який метод використовувати, а не відповідати на таке основне питання, як у мене, тобто як це взагалі зробити?

Зазвичай NAD83 і WGS84 знаходяться в межах одного метра один від одного. Ваша стурбованість різницею в 2,5 футів, яка менше метра, вказує на те, що вам потрібно здійснити цю трансформацію даних .

Коротко цей розрахунок вимагає знання того, коли були зібрані координати, щоб можна було врахувати їх рух у часі (в першу чергу за рахунок руху підстилаючих тектонічних плит і частково через "епізодичного руху, пов'язаного з землетрусами"). Ви спочатку перетворюєте необроблені дані в еквівалентні координати, на які посилаються конкретні стандартні дати, такі як 1997.0 або 2002.0. Потім ви шукаєте 14 параметрів, що описують (а) зміщення, обертання та масштабування (тобто семипараметричне перетворення), необхідні для зміни даних на задану дату разом з (b) першими похідними цих параметрів щодо часу. Помноження (b) на пройдений час і додавання до (a) дає наближення першого порядку до потрібного семипараметричного перетворення, яке ви потім застосуєте.

Незважаючи на те, що процедура не є математично складною, вона досить задіяна і залежить від деталей того, як спочатку ваші дані були зібрані та оброблені. Документ NOAA про трансформації між NAD83 та WGS84 дає чітке, зрозуміле пояснення. Він рекомендує використовувати інтерактивне програмне забезпечення "Горизонтальна залежна від часу позиція", HTDP, доступне на веб-сайті NOAA за адресою http://www.ngs.noaa.gov/TOOLS/Htdp/Htdp.shtml . Вихідний код Fortran складається головним чином з коефіцієнтів тектонічної моделі пластини плюс деякого інтерфейсного коду, але я підрахував, що близько 3000 рядків насправді роблять розрахунки. Це може бути значним болем для порту.

Я б не був впевнений у точності повітряних фотографій Google Планета Земля, якщо вона знижена до 2,5 футів. Це може бути лише питання не зовсім ортогональної фотографії.

Можливо, вам знадобиться лише афінна трансформація, щоб вирівняли свої фотографії з їхніми.