Generate a grid of latitude values

ToddMcDonnell

I need to convert a 30m elevation grid to a grid that gives values for the cell latitude rather than elevation.

I tried converting the GRID to point, adding x,y data to the points, then converting back to GRID using the y-coordinate as the value, however, the large region I am working with is preventing the processes from completing successfully.

Any advice will certainly be appreciated.

Todd

LuizVianna

This worked for me, for southern hemisphere (negative lat and lon). It is necessary to set the environmental options before running the script.

import numpy as np
import arcpy
from arcpy.sa import *  
from arcpy import env as E

#Create variables from environment
latmin = E.extent.YMin
latmax = E.extent.YMax
lonmin = E.extent.XMin
lonmax = E.extent.XMax
cellsize = float(E.cellSize)
lowerleft = arcpy.Point(lonmin, latmin)

#Generating inverse meshgrid for latitude grid
lat1 = np.arange(float(-latmax), float(-latmin), float(cellsize))
lon1 = np.arange(float(-lonmax), float(-lonmin), float(cellsize))
lon1, lat1 = np.meshgrid(lon1, lat1)

lat = arcpy.NumPyArrayToRaster(-lat1, lowerleft, float(cellsize), float(cellsize))

#Generating meshgrid for longitude grid
lat2 = np.arange(float(latmin), float(latmax), float(cellsize))
lon2 = np.arange(float(lonmin), float(lonmax), float(cellsize))
lon2, lat2 = np.meshgrid(lon2, lat2)

lon = arcpy.NumPyArrayToRaster(lon2, lowerleft, float(cellsize), float(cellsize))‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍