Más

Transformación de datos mediante la transformación de costum 2D en Arcmap 10.2 para escritorio


Tengo capas de ArcMap de carreteras, edificios, etc. en un sistema de referencia (Israel_TM). Necesito transformar todos estos datos (manteniendo la topología) al "sistema de referencia" de espacio de píxeles. Ya he calculado los parámetros de transformación (usando Matlab) de una Transformación Conformal (2D) entre los dos sistemas de referencia: Israel_TM a Pixeles. Es una transformación de 4 parámetros que tengo: 1xRotation, 1xScale, 2xTranslations. ¿Cómo puedo crear una transformación 2D definida por el usuario y usarla para transformar todos mis datos en Pixel-space? Todo lo que pude encontrar es Helmert 7-Par. transformación u otros métodos de transformación 3D en Propiedades de capa -> Sistema de coordenadas -> Transformación. Intenté usar la Caja de herramientas para esto, pero no pude encontrar nada relevante.


No se puede transformar el sistema de coordenadas de British National Grid a WGS_OSGB_1984 Mercat

Comenzaré diciendo que soy muy nuevo en esto, y aunque he tomado el curso introductorio a ArcGIS Desktop 10.0, realmente no me está ayudando con mi problema.

He estado recopilando información de varias partes interesadas y marcándola en un mapa (ráster trazado). Para analizar la información, he importado los mismos mapas de fondo y los datos recopilados en un nuevo documento de mapa, sin embargo, no se alinean. Parece que la razón de esto es que en el software que estamos usando (un cuestionario especialmente desarrollado, el CS se configuró como British National Grid, pero el ráster trazado es WGS_OSGB_1984 Mercator; cuando importé por primera vez los mapas al paquete de cuestionario, me informó que era necesario transformar las coordenadas que hice.

He investigado esto tanto como puedo, y aparentemente hay una transformación de BNG a WGS84 disponible para 9.1, pero parece que no puedo encontrar nada para 10.0. Todas las transformaciones proporcionadas en Arc 10 son inútiles, y he pasado días tratando de calcular una transformación mislef, pero está yendo por todos lados sin un patrón lógico real.

En este punto, agregaría que trabajo para una organización sin fines de lucro y, por lo tanto, no tenemos el presupuesto para pagarle a un desarrollador para que solucione esto, pero si puedes ayudarme, ¡serás un héroe absoluto!

por MelitaKennedy

¿Está el ráster en WGS_1984_Web_Mercator o WGS_1984_Web_Mercator_Auxiliary_Sphere? Eso hará una diferencia. El primero usa un sistema de coordenadas geográficas personalizado (WGS 1984 modificado, "GCS_WGS_1984_Major_Auxiliary_Sphere) y no hay transformación directa entre este y OSGB 1936. Tiene que usar dos transformaciones, lo cual es. Difícil. En ArcMap, pero se puede hacer convirtiendo el datos en ArcToolbox usando las herramientas Project o Project Raster, el segundo usa el estándar WGS 1984 y debe haber disponibles varias transformaciones.

Gracias por su aporte. El sistema de coordenadas se describe como WGS_1984_Mercator, no se menciona web_mercator o web_mercator_auxillary_sphere.

Intentaré utilizar las transformaciones que ha mencionado, pero cualquier información adicional que pueda proporcionar será recibida con gratitud.

por MelitaKennedy

Oh, sí, esa es una definición genérica de Mercator (centrada en 0,0) que agregamos como un caso de prueba. Si agrego datos WGS84 a ArcMap y configuro el sistema de coordenadas del marco de datos en British National Grid, obtengo esta lista de transformaciones con las áreas de uso y precisiones agregadas:

OSGB_1936_To_WGS_1984_1 "Reino Unido - Gran Bretaña Isla de Man" 21.0
OSGB_1936_To_WGS_1984_2 "Reino Unido - Inglaterra" 10.0
OSGB_1936_To_WGS_1984_3 "Reino Unido - Inglaterra Gales Isla de Man" 21.0
OSGB_1936_To_WGS_1984_4 "Reino Unido - Escocia" 18.0
OSGB_1936_To_WGS_1984_5 "Reino Unido - Gales" 35,0
OSGB_1936_To_WGS_1984_NGA_7PAR "Reino Unido - Gran Bretaña Isla de Man" 21.0
OSGB_1936_To_WGS_1984_Petroleum "Reino Unido - Gran Bretaña Isla de Man" 2.0

Al verificar posibles transformaciones, solo importan los sistemas de coordenadas geográficas.

Probé las transformaciones enumeradas tanto para transformar los datos de mis partes interesadas (recopilados en GCS 1936 / British National Grid) a 1984 (utilizados para el mapa de fondo) como para transformar los datos del mapa en 1936 / BNG GCS. He puesto un marcador en un punto distintivo para intentar alinearlos sin mucho éxito, a continuación están las coordenadas de cada conjunto de transformaciones para mostrar la escala del problema. No sé si esto le permitirá ayudarme más, pero es todo lo que puedo pensar en este momento (la lógica no me lleva muy lejos).

Transformación de datos de 1936 a 1984 (todos los valores en metros)

Punto de referencia -475663.757, 6844471.977
Transformación 1 -283557.103, 4080120.052
Transformación 2 -283551.527, 4080094.476
Transformación 3 -283602.731, 4080103.425
Transformación 4 -283547.172, 4080068.931
Transformación 5 -283547.172, 4080107.315
Transformación_NGA_7PAR -283568.237, 4080103.425
Transformation_Petroleum -283568.237, 4080068.931

Transformación de datos de 1984 a 1936 (todos los valores en metros)

Punto de referencia 245573.383, 283857.525
Transformación 1 215822.76, 2396583.867
Transformación 2 215825.84, 2396594.223
Transformación 3 215825.84, 2396594.223
Transformación 4 215825.84, 2396610.354
Transformación 5 215809.709, 2396578.091
Transformación_NGA_7PAR 215825.84, 2396642.617
Transformation_Petroleum 215825.24, 2396642.617

Anteriormente había intentado resolver mi propia transformación usando este punto de referencia y podía acercarme bastante (en un kilómetro más o menos), pero creo que hubo un problema de escala o rotación (estaba usando una transformación xy simple). Cuando he intentado desarrollar una transformación más compleja, los datos parecen saltar por todos lados sin ningún patrón lógico.


Crear un archivo de transformación geográfica personalizado

Se crea un archivo de transformación geográfica personalizado mediante la herramienta Crear transformación geográfica personalizada. Este archivo .gtf se crea en una de dos carpetas, según los permisos de escritura de su cuenta de usuario.

  • En Windows, el archivo normalmente se guarda en C: Users & ltUSERNAME & gt AppData Roaming ESRI Desktop10.X ArcToolbox CustomTransformations. Se puede acceder directamente a esta ruta ingresando% appdata% ESRI Desktop10.X ArcToolbox CustomTransformations en el Explorador de Windows. (Asegúrese de actualizar la X en la ruta para que coincida con su versión de ArcGIS for Desktop).
  • Si no tiene permiso de escritura en la ruta anterior (se ha configurado como de solo lectura), el archivo se creará en la carpeta temporal de su cuenta de usuario. Se puede acceder a esto ingresando% temp% ArcToolbox CustomTransformations en el Explorador de Windows.
  • En un sistema Unix, el .gtf se guardará en el directorio de inicio del usuario dentro de ArcToolbox CustomTransformations, ya sea en la carpeta Application Data o en la carpeta tmp, nuevamente, dependiendo de sus permisos de escritura.

Consulte las notas de uso de la herramienta Crear transformación geográfica personalizada para obtener información más detallada sobre dónde se guarda este archivo.


Nordpil

El mejor conjunto de datos disponible gratuitamente sobre placas tectónicas, orógenos y límites de placas fue publicado en 2003 por Peter Bird (anteriormente Universidad de California)..

Los datos se publicaron originalmente en el artículo. Un modelo digital actualizado de límites de placa (Geoquímica Geofísica Geosistemas, 4 (3), 1027, doi: 10.1029 / 2001GC000252, 2003). El conjunto de datos está disponible para su descarga gratuita sin restricciones, actualmente en el sitio web privado del investigador. Sin embargo, solo está disponible en un formato de texto complicado, lo que significa que no es fácilmente accesible para su uso en aplicaciones SIG.

Para llevar este conjunto de datos a la era moderna, los datos originales se analizaron, limpiaron y verificaron con ArcGIS 10.2 y se convirtieron en archivos de formas.

Todos los datos originales se recuperaron de http://peterbird.name/oldFTP/PB2002/ en junio de 2014 y se analizaron con Python y Global Mapper 11. Se realizó una mayor manipulación de datos en ArcGIS 10.2. Las ediciones principales se referían a los segmentos que abarcaban el límite de -180/180, que tenían que dividirse y moverse manualmente.

Esta colección está disponible bajo el Licencia de atribución de Open Data Commons, que impone pocas restricciones al uso de los datos.

El conjunto de datos presenta placas tectónicas y sus límites, así como orógenos e información sobre los límites. Los datos son útiles para aplicaciones geológicas, análisis y educación, y deben ser fáciles de usar en cualquier aplicación de software GIS moderna. Para obtener información sobre los campos y los valores, consulte la documentación original y el artículo científico.


1.3 Tutorial interactivo

El resto del taller se dedicó a un ejercicio práctico en el que vivimos codificados juntos. Los datos que usamos fueron los datos de los límites de los distritos de concejales de Chicago de 1986, que se pueden encontrar en la página de Colección de mapas de UChicago aquí.

El guión de este taller se puede encontrar aquí.

1.3.1 Funciones aprendidas

1.3.2 Importar datos espaciales

Lo primero que haremos es importar los datos espaciales. Nuestros datos espaciales se almacenan como un shapefile (.shp, pero en realidad 4 archivos). Estos se están poniendo menos de moda (consulte GeoJSON, etc.), pero muchos datos espaciales todavía se almacenan de esta manera.

Primero cargue la biblioteca sf. Si no lo tiene, instálelo en su consola o en el panel Packages RStudio.

Usaremos la función st_read. Esto también lee GeoJSON, bases de datos PostGIS y más.

Puede usar el atajo de teclado Ctrl-Enter para ejecutar una línea de código en R.

Presione el tabulador después de escribir st_read (& quot & quot) con el cursor entre las comillas y obtendrá una buena función de autocompletar.

Compruebe en qué proyección se encuentran los datos:

¡Sike, esto no está proyectado! Puede saberlo porque proj4string comienza con + proj = longlat. Aún puede trazar esto, pero las cosas comenzarán a ponerse dudosas si intenta hacer cálculos de distancia o área, o trazar estos límites de barrio con otras capas.

Aún puedo trazar el mapa, pero como buen analista de datos geográficos, tendré que proyectarlo.

1.3.3 Proyecte sus datos

Necesitamos proyectar los datos: pero ¿cómo elegimos qué proyección necesitamos? Esta es una pregunta filosófica profunda, pero afortunadamente la biblioteca UChicago ya nos ha dicho las mejores proyecciones para usar para nuestros datos: Avión estatal del este de Illinois o Zona UTM 16.

Por lo general, hago una búsqueda en línea para buscar los códigos EPSG para la proyección que quiero usar, o uso Spatialreference.org, que tiene una base de datos para todos los códigos EPSG para proyecciones.

Después de investigar un poco, encuentro que el código EPSG que quiero para UTM Zone 16 es 32616. Luego uso st_transform () para proyectar los datos y los guardo como un nuevo marco de datos sf.

Reviso el CRS, ¡sí, se ve bien!

Pregunta: ¿en qué dato está esto? ¿Cuáles son las unidades para esta proyección?

  • Proyecte los datos del distrito de 1986 en la proyección del plano del estado del este de Illinois. Siga los pasos anteriores para hacerlo. Sugerencia: busque el código EPSG en línea.

Tenga en cuenta que la parte + unidades = us-ft de proj4string significa que cualquier cálculo de distancia que haga en los datos del plano de estado del ward86 se hará en pies.

Tenga en cuenta que existen proyecciones "malas" para sus datos. Por ejemplo, si elegí accidentalmente la proyección Alaska Albers para mis datos ...

¡Claramente esta proyección no es muy buena para Chicago! Pero algo a tener en cuenta si, por ejemplo, está haciendo mapas de lugares más cercanos a los polos. ''


Editar los parámetros del transformador

Usando un conjunto de opciones de menú, los parámetros del transformador se pueden asignar haciendo referencia a otros elementos en el espacio de trabajo. Algunas funciones más avanzadas, como un editor avanzado y un editor aritmético, también están disponibles en algunos transformadores. Para acceder a un menú de estas opciones, haga clic junto al parámetro aplicable. Para obtener más información, consulte Opciones del menú de parámetros del transformador.

Definiendo valores

Hay varias formas de definir un valor para su uso en un transformador. La más simple es simplemente escribir un valor o una cadena, que puede incluir funciones de varios tipos, como referencias de atributos, funciones matemáticas y de cadena, y parámetros del espacio de trabajo. Hay una serie de herramientas y atajos que pueden ayudar a construir valores, generalmente disponibles en el menú contextual desplegable adyacente al campo de valor.

Usando el editor de texto

El Editor de texto proporciona una forma conveniente de construir cadenas de texto (incluidas las expresiones regulares) a partir de varias fuentes de datos, como atributos, parámetros y constantes, donde el resultado se usa directamente dentro de un parámetro.

Usando el editor aritmético

El editor aritmético proporciona una forma conveniente de construir expresiones matemáticas a partir de varias fuentes de datos, como atributos, parámetros y funciones de características, donde el resultado se usa directamente dentro de un parámetro.

Valores condicionales

Establezca valores en función de una o más condiciones de prueba que pasen o no.

Contenido

Las expresiones y cadenas pueden incluir varias funciones, caracteres, parámetros y más.

Al configurar valores, ya sea que se ingresen directamente en un parámetro o se construyan con uno de los editores, las cadenas y expresiones que contengan funciones de cadena, matemáticas, fecha / hora o características FME tendrán esas funciones evaluadas. Por lo tanto, los nombres de estas funciones (en la forma @ & ltnombre de la función& gt) no deben usarse como valores de cadena literales.


Métodos de autenticación multifactor

Los usuarios deben autenticarse con todos o algunos de estos factores:

Autenticación de factor único: La autenticación basada en contraseña es el tipo de autenticación de factor único más utilizado, que se considera inseguro ya que la mayoría de las personas utilizan contraseñas que se comprometen fácilmente o son débiles.

Autenticación de dos factores (2FA): Se requiere que un usuario proporcione un código a través de una aplicación de autenticación en su teléfono inteligente, un token de software (& ldquosoft token & rdquo), un token de hardware (& ldquohard token & rdquo) o un mensaje SMS. Este código puede ser una contraseña única que caduca después de 30 segundos, o una contraseña HMAC única, o HOTP, que no caduca hasta que se utiliza.

Autenticación de tres factores (3FA): Los factores biométricos se tienen en cuenta con la autenticación de tres factores y ndash, algo que es único para cada usuario. Estos incluyen reconocimiento facial, geometría de la mano, geometría de los dedos, escaneos de huellas de la palma y escaneos de huellas dactilares.

Autenticación de cuatro factores (4FA): En la autenticación de cuatro factores, se tiene en cuenta la ubicación geográfica de un usuario y rsquos, junto con la cantidad de tiempo que tardó el usuario en llegar a su destino. Estos tipos de medidas de autenticación normalmente requieren que los usuarios validen esta información. En cambio, se ejecuta en segundo plano mientras se determina su riesgo de autenticación. La autenticación de cuatro factores confirma la geolocalización y velocidad geográfica del usuario y rsquos, así como su ubicación en el mundo. Por ejemplo, habría un problema de seguridad si un usuario se autenticara con un proveedor de autenticación multifactor en Chicago e intentara autenticarse nuevamente en Moscú diez minutos después.

Autenticación de cinco factores (5FA): En la autenticación de cinco factores, un usuario puede generar patrones táctiles o gestos específicos. Por ejemplo, un usuario puede usar el último sistema operativo Windows (8 o posterior) para crear una contraseña de imagen usando una pantalla táctil de computadora, donde puede tocar una imagen o dibujar líneas rectas o círculos para generar una contraseña de gestos única.

BAutenticación iométrica: La biometría de autenticación se está volviendo más común, ya que las capacidades de escaneo de reconocimiento facial y de huellas dactilares en los productos Android y Apple han ganado popularidad entre los consumidores. Existen otros métodos biométricos que aún no han sido autorizados para el uso comercializado de autenticación multifactor, como el reconocimiento de mecanografía, el análisis de firma y escritura a mano, patrones de venas, patrones de marcha, identidad de olor, escáneres de iris, escáneres de retina, impresiones de voz y forma de orejas. reconocimiento.


Esa animación está hecha a medida en Flash, que es una herramienta perfecta y completa para hacer animaciones 2D (e incluso juegos) como esa, pero tiene un costo bastante alto, incluso para la licencia de estudiante ($ $ 19.99 $ por mes o $ $ 199 $ por año).

GeoGebra es una herramienta gratuita para dibujar geometría, tiene capacidades para interactuar, animaciones y tiene una versión beta para gráficos 3D. Es una de las mejores herramientas para hacer geometría, tiene una versión online bastante limitada y una versión completamente instalada (ambas gratuitas).

Cut-the-knot Geoboard es un subprograma fantástico (interacciones, animaciones), y tiene un costo muy cómodo ($ $ 5 $) para una licencia.

Para realizar gráficos arquitectónicos en 3D de una manera rápida y sencilla, recomiendo la versión gratuita de Sketchup, que además es una herramienta fantástica para hacer diseños y también es posible hacer animaciones en 3D.

Como herramienta de diseño 3D en toda regla, realmente recomiendo Blender, un programa gratuito para hacer 3D en el que tienes una capacidad artística total, pero por eso no tiene herramientas específicas de geometría.

Para geometría 2D:

Hay http://www.robocompass.com/app - gratis y fantástico y anima construcciones geométricas.

Para geometría 3D:

He desarrollado un programa de dibujo 3D gratuito que se ejecuta en navegadores modernos (html5). La versión en inglés se llama Geoservant 3D, la versión original en alemán se llama Geoknecht 3D.

Allí puede especificar formas 3D (geometría sólida) y animarlas cambiando sus valores. Consulte la descripción en la página sobre cómo cambiar los valores: mantenga presionada la tecla MAYÚS CTRL ALT y use las teclas del cursor. También eche un vistazo a la galería en la página para ver qué cosas puede lograr con este programa usando una sintaxis simple. Para cada objeto, el programa calcula valores básicos (área, volumen, vectores).

Para diferentes colores, simplemente agregue el código de color HEX al final, p. Ej. esfera (3 | 3 | 1 2) . Si agrega un valor entre 0 y 1 entre paréntesis, puede definir una transparencia para cada objeto, p. Ej. esfera (3 | 3 | 1 2)

Además, se pueden definir variables. Ejemplo:

Se puede crear un hemisferio como esa esfera (0 | 0 | 4 4 0 6.28 1.55 3.14), ¡vea la galería en el sitio!

Idiomas: Geoservant ahora está disponible en alemán, inglés, español, francés, ruso, sueco y chino, y acepta los objetos especificados en esos idiomas, por ejemplo: 长方体 (4 | 1 | 1 2 | 4 | 1)

Alternativa

Más compleja y poderosa: Encuentro la herramienta de geometría 3D tinkercard.com bastante asombrosa.

Otro programa es Geodrafter 2D (inglés) o Geozeichner 2D (alemán) para formas bidimensionales.


Sintaxis de secuencias de comandos

NetCDFConvert2DVariableToArcGISRaster_GeoEco (inputFile, outputRaster, variableName, xLowerLeftCorner, yLowerLeftCorner, cellSize, nodataValue, transpose, mirror, flip, swapHemispheres, coordenadaSystem, proyectadoCoordinateSystem, resumir, transposición, espejo, voltear, intercambiar

En el momento en que se desarrolló esta herramienta, se estaban desarrollando dos versiones principales de netCDF: NetCDF-3, que estaba en la versión 3.6.2, y NetCDF-4, que estaba en la versión 4.0 beta 1. La implementación actual de esta herramienta puede procesar Archivos NetCDF-3 pero no archivos NetCDF-4. La compatibilidad con NetCDF-4 se agregará en un futuro próximo.

Si proporciona un archivo comprimido en un formato de compresión compatible, se descomprimirá automáticamente. Si es un archivo (por ejemplo, .zip o .tar), debe contener exactamente un archivo, que no debe estar en un subdirectorio.

Nombre de una variable en el archivo netCDF. La variable debe tener dos dimensiones y un tipo de datos entero o de punto flotante.

Si no conoce el nombre de la variable, proporcione su mejor estimación. Si una variable no existe con ese nombre, se generará un ValueError. El mensaje de error enumerará los nombres de todas las variables 2D presentes en el archivo de entrada. Seleccione uno de estos nombres y vuelva a intentarlo.

También puede utilizar la herramienta Extraer encabezado NetCDF para volcar el encabezado netCDF en un archivo de texto. Puede encontrar los nombres de las variables buscando en el archivo la sección que comienza con "variables:". Cada variable se enumera como:

donde tipo de datos es el tipo de datos de la variable (p. ej., "flotante"), nombre es el nombre de la variable, y dimensiones es la lista de dimensiones de la variable. Por ejemplo, el siguiente extracto de un Aviso DT-MADT sobre corrientes geostróficas netCDF enumera dos variables que tienen dos dimensiones, Grid_0001 y Grid_0002:

Esta herramienta extrae la variable netCDF a una cuadrícula ArcInfo ASCII y la convierte en un ráster mediante la herramienta de geoprocesamiento ArcGIS ASCII to Raster. Esa herramienta y el formato ráster de ArcGIS tienen varias limitaciones que restringen los tipos de variables netCDF que se pueden convertir correctamente al formato ráster de ArcGIS.

Las variables NetCDF que utilizan el tipo de datos flotante o doble no deben contener valores "infinito" (INF) o "no un número" (NAN). Se generará un ValueError si se descubren estos valores.

El formato ráster de ArcGIS admite el tipo de datos flotantes de 32 bits, pero no el tipo de datos dobles de 64 bits. Si proporciona una variable netCDF con el tipo de datos doble, se convertirá en un ráster flotante de 32 bits mediante la herramienta de geoprocesamiento ArcGIS ASCII to Raster. El comportamiento de la herramienta en esta situación no está documentado. En ArcGIS 9.1 parece ser:

Los valores en los que el exponente oscila entre -38 y +38 se representan correctamente en el ráster flotante de 32 bits resultante, aunque se pierde algo de precisión debido a la mantisa más pequeña del tipo de datos flotantes de 32 bits.

Los valores donde el exponente es menor que -38 (por ejemplo, -39, -40, etc.) se convierten a 0.

Los valores donde el exponente es mayor que +38 se convierten a -INF o + INF, según el signo del valor (por ejemplo, -5.3083635279597874e-212 aparece como -1. # INF en la GUI de ArcCatalog, mientras que 2.5502286890301497e + 084 aparece como 1. # INF).

La herramienta ArcGIS 9.1 ASCII a ráster también presenta algunas peculiaridades al convertir rásteres enteros:

Para datos enteros de 8 bits, la herramienta creará un ráster de 16 bits si aparece el valor -128, a menos que -128 se designe como valor NODATA. Especificar un valor NODATA diferente, como 0, aún produce un ráster de 16 bits si aparece -128.

De manera similar, para datos enteros de 16 bits, la herramienta creará un ráster de 32 bits si el valor -32768 aparece en el archivo ASCII, a menos que se designe como valor NODATA.

Peor aún, para datos enteros de 32 bits, la herramienta informará un error si aparece el valor -2147483648 a menos que se designe como valor NODATA. Aún más extraño, el valor -2147483647 siempre se traduce a NODATA, pase lo que pase.

Para todos los tipos de datos enteros, la herramienta produce un comportamiento extraño cuando especifica un valor NODATA que no es el valor más pequeño posible para el tipo de datos. Por ejemplo, si el rango de datos va de 0 a 255 y 0 se designa como valor NODATA, la herramienta produce un ráster de salida de 8 bits sin signo. Pero si se designa 1 como valor NODATA, produce un ráster de salida firmado de 16 bits, y ArcCatalog muestra en Propiedades del dataset ráster que el valor NoData es -32768, aunque la herramienta Identificar muestra que las celdas que tenían el valor 1 son en realidad NODATA. Se pueden obtener resultados extraños similares para rásteres enteros de otros tipos de datos, cuando se designa un NODATA que no es el valor más pequeño posible.

Coordenada X de la esquina inferior izquierda del ráster.

La coordenada es para la esquina de la celda inferior izquierda, no el centro de esa celda. Por ejemplo, si el ráster es una proyección geográfica de toda la Tierra, la coordenada de la esquina inferior izquierda sería -180,0, correspondiente a una longitud de 180 grados Oeste.

Coordenada Y de la esquina inferior izquierda del ráster.

La coordenada es para la esquina de la celda inferior izquierda, no el centro de esa celda. Por ejemplo, si el ráster es una proyección geográfica de toda la Tierra, la coordenada de la esquina inferior izquierda sería -90,0, correspondiente a una latitud de 90 grados Sur.

Por ejemplo, si el ráster es una proyección geográfica de toda la Tierra, con 720 columnas y 360 filas, tendría un tamaño de celda de 0.5, correspondiente a 1/2 de un grado geográfico.

El formato de datos subyacente requiere que las celdas sean cuadradas. No es posible especificar un tamaño de celda para cada dimensión.

Valor que indica que una celda no tiene datos.

Si es True, la imagen se transpondrá (se volteará sobre el eje diagonal) antes de la conversión. Utilice esta opción para corregir una imagen que tenga el eje este / oeste hacia arriba y hacia abajo en lugar de hacia la izquierda y hacia la derecha.

Si es Verdadero, la imagen se volteará sobre el eje vertical antes de la conversión. Utilice esta opción para corregir una imagen que sea la "imagen reflejada" de lo que se supone que es.

Si es Verdadero, la imagen se volteará sobre el eje horizontal antes de la conversión. Utilice esta opción para arreglar una imagen que está al revés.

Si es Verdadero, se intercambiarán los hemisferios este y oeste de la imagen. Utilice esta opción para cambiar la orientación de una imagen global de una orientación de 0 a 360 centrada en el océano Pacífico a una orientación de -180 a +180 centrada en el océano Atlántico, o viceversa.

Sistema de coordenadas a definir para el ráster de salida. Si no se proporciona un valor, el sistema de coordenadas del ráster de salida permanecerá sin definir.

Nuevo sistema de coordenadas para proyectar el ráster de salida.

El ráster solo se puede proyectar en un nuevo sistema de coordenadas si se define la proyección original. Se generará un error si especifica un nuevo sistema de coordenadas sin definir el sistema de coordenadas original.

La herramienta ArcGIS Project Raster se utiliza para realizar la proyección. La documentación de esa herramienta recomienda que también especifique un tamaño de celda para el nuevo sistema de coordenadas.

He notado que para ciertos sistemas de coordenadas, la herramienta ArcGIS 9.2 Project Raster parece recortar el ráster proyectado en una extensión arbitraria que es demasiado pequeña. Por ejemplo, al proyectar una imagen global de clorofila MODIS Aqua de 4 km en coordenadas geográficas a Lambert_Azimuthal_Equal_Area con meridiano central de -60 y latitud de origen de -63, la imagen resultante se recorta para mostrar solo una cuarta parte del planeta. Este problema no ocurre cuando Project Raster se invoca de forma interactiva desde la interfaz de usuario de ArcGIS; solo ocurre cuando la herramienta se invoca mediante programación (el método ProjectRaster_management del geoprocesador). Por lo tanto, es posible que no lo vea cuando use Project Raster usted mismo, pero puede suceder cuando use herramientas MGET que invocan Project Raster como parte de sus operaciones de geoprocesamiento.

Si encuentra este problema, puede solucionarlo de esta manera:

Primero, ejecute esta herramienta sin especificar un nuevo sistema de coordenadas, para obtener el ráster de salida en el sistema de coordenadas original.

En ArcCatalog, use la herramienta Proyectar ráster para proyectar el ráster al nuevo sistema de coordenadas. Verifique que esté presente todo el ráster, que no se haya recortado en una extensión demasiado pequeña.

En ArcCatalog, busque la extensión del ráster proyectado haciendo clic derecho sobre él en el árbol del catálogo, seleccionando Propiedades y desplazándose hacia abajo hasta Extensión.

Ahora, antes de ejecutar la herramienta MGET que proyecta el ráster, establezca la configuración del entorno de Extensión en los valores que buscó. Si está invocando la herramienta MGET de forma interactiva desde ArcCatalog o ArcMap, haga clic en el botón Entornos en el cuadro de diálogo de la herramienta, abra Configuración general, cambie el menú desplegable Extensión a "Como se especifica a continuación" y escriba los valores que buscó. Si lo está invocando desde un modelo de geoprocesamiento, haga clic con el botón derecho en la herramienta en el modelo, seleccione Hacer variable, Desde el entorno, Configuración general, Extensión. Esto colocará Extent como una variable en su modelo, adjunta a la herramienta MGET. Abra la variable Extensión, cámbiela a "Como se especifica a continuación" y escriba los valores que buscó. Si está invocando la herramienta MGET mediante programación, debe establecer la propiedad Extensión del geoprocesador en los valores que buscó. Consulte la documentación de ArcGIS para obtener más información sobre esto y la configuración del entorno en general.

Ejecute la herramienta MGET. La extensión del ráster de salida debería tener ahora el tamaño adecuado.

Un método de transformación utilizado para convertir entre el sistema de coordenadas original y el nuevo sistema de coordenadas.

Este parámetro es una nueva opción introducida por ArcGIS 9.2. Debe tener ArcGIS 9.2 para usar este parámetro.

Este parámetro solo es necesario cuando especifica que el ráster debe proyectarse a un nuevo sistema de coordenadas y ese nuevo sistema usa un datum diferente al sistema de coordenadas original, o existe alguna otra diferencia entre los dos sistemas de coordenadas que requiere una transformación. Para determinar si se necesita una transformación, recomiendo el siguiente procedimiento:

Primero, ejecute esta herramienta sin especificar un nuevo sistema de coordenadas, para obtener el ráster de salida en el sistema de coordenadas original.

A continuación, use la herramienta ArcGIS 9.2 Project Raster en el ráster de salida para proyectarlo en el sistema de coordenadas deseado. Si se necesita una transformación geográfica, esa herramienta le pedirá una. Anote el nombre exacto de la transformación que utilizó.

Finalmente, si se necesitaba una transformación, escriba el nombre exacto en esta herramienta, vuelva a ejecutarla y verifique que el ráster de salida se proyectó como deseaba.

El algoritmo de remuestreo que se utilizará para proyectar el ráster original en un nuevo sistema de coordenadas. La herramienta ArcGIS Project Raster se utiliza para realizar la proyección y acepta los siguientes valores:

MÁS CERCANO: interpolación del vecino más cercano

BILINEAR - interpolación bilineal

Debe especificar uno de estos algoritmos para proyectar a un nuevo sistema de coordenadas. Se generará un error si especifica un nuevo sistema de coordenadas sin seleccionar un algoritmo.

El tamaño de celda del sistema de coordenadas proyectadas. Aunque este parámetro es opcional, para recibir los mejores resultados, la documentación de ArcGIS recomienda que siempre lo especifique al proyectar a un nuevo sistema de coordenadas.

Las coordenadas xey (en el espacio de salida) utilizadas para la alineación de píxeles.

Este parámetro es una nueva opción introducida por ArcGIS 9.2. Debe tener ArcGIS 9.2 para usar este parámetro. Se ignora si no especifica que el ráster debe proyectarse en un nuevo sistema de coordenadas.

Rectángulo al que se debe recortar el ráster.

Si se especificó un sistema de coordenadas proyectadas, el recorte se realiza después de la proyección y las coordenadas del rectángulo deben especificarse en el nuevo sistema de coordenadas. Si no se especificó un sistema de coordenadas proyectadas, las coordenadas deben especificarse en el sistema de coordenadas original.

La herramienta ArcGIS Clip se utiliza para realizar el clip. El rectángulo de recorte debe pasarse a esta herramienta como una cadena de cuatro números separados por espacios. La interfaz de usuario de ArcGIS formatea automáticamente la cadena correctamente al invocar esta herramienta desde la interfaz de usuario de ArcGIS, no necesita preocuparse por el formato. Pero al invocarlo mediante programación, tenga cuidado de proporcionar una cadena con el formato adecuado. Los números están ordenados IZQUIERDA, ABAJO, DERECHA, ARRIBA. Por ejemplo, si el ráster está en un sistema de coordenadas geográficas, puede recortarse a 10 W, 15 S, 20 E y 25 N con la cadena:

Pueden proporcionarse números enteros o decimales.

Expresión de álgebra de mapas para ejecutar en el ráster de salida.

ADVERTENCIA: El Generador de modelos de geoprocesamiento de ArcGIS puede eliminar de forma aleatoria y silenciosa el valor de este parámetro. Este es un error en ArcGIS. Antes de ejecutar un modelo que ha guardado, abra esta herramienta y valide que el valor del parámetro todavía existe.

La expresión se ejecuta después de proyectar y recortar el ráster convertido (si se especifican esas opciones). Utilice la cadena inputRaster que distingue entre mayúsculas y minúsculas para representar el ráster sobre el que ahora desea realizar el álgebra de mapas. Por ejemplo, para convertir el ráster en un ráster entero y agregar 1 a todas las celdas, use esta expresión:

La cadena inputRaster distingue entre mayúsculas y minúsculas. Antes de ejecutar la expresión de álgebra de mapas, la cadena se reemplaza con la ruta a un ráster temporal que representa el ráster de salida que se está generando. La expresión final debe tener menos de 4000 caracteres o ArcGIS informará un error.

La herramienta Álgebra de mapas de salida única de ArcGIS se utiliza para ejecutar la expresión de álgebra de mapas. Debe tener una licencia para la extensión ArcGIS Spatial Analyst para poder realizar álgebra de mapas.

La sintaxis del álgebra de mapas puede ser muy delicada. Aquí hay algunos consejos que le ayudarán a tener éxito con esta herramienta:

Before using this tool, construct and test out your map algebra expression using the ArcGIS Single Output Map Algebra tool. Then paste the expression into this tool and edit it to use the inputRaster variable rather than the test value you used with Single Output Map Algebra.

If you do develop your expression directly in this tool, start with a very simple expression. Verify that it works properly, add a little to it, and verify again. Repeat this process until you have built up the complete expression.

Always separate mathematical operators from raster paths using spaces. In the example above, the / operator contains a space on either side. Follow this pattern. In some circumstances, ArcGIS will fail to process raster algebra expressions that do not separate raster paths from operators using spaces. The reported error message usually does not indicate that this is the problem, and tracking it down can be very frustrating.

If True, pyramids will be built for the output raster, which will improve its display speed in the ArcGIS user interface. This is the last step performed in post-conversion processing.


Can we use tf.spectral fourier functions in keras?

Let us start with an input that is a simple time series and try to build an autoencoder that simply fourier transforms then untransforms our data in keras.

Then the third line throws an error when entered:

You see, the output of tf.spectral.irfft is float32 but it looks like Lambda thinks it is complex64?? (Complex64 is the input x from the previous step)

We can fix that error at model entry time with:

This is accepted at input time but then when we try to build the model:

Which I guess is reasonable and was the reason I didn't want to cast it in the first place.

Main question: how do I successfully wrap the tf.spectral.irfft function which outputs float32 ?

More general question for learning: Let's assume I actually want to do something between the rfft and the irfft, how can I cast those imaginary numbers into absolute values without breaking keras so I can apply various convolutions and the like?


Ver el vídeo: Proyectar y Reproyectar un shapefile ArcGIS - (Septiembre 2021).