Tag Info

New answers tagged

2

Riffing on @mkennedy's comment above, here's a small script to parse and create a spatial reference if the datum is NAD83. import arcpy inputString = r'PROJCS["GCS North American 1983 UTM Zone 10S (Calculated)",GEOGCS[ETCETCETC]' editString = inputString.replace(" ", "_").upper() if "NORTH_AMERICAN_1983" in editString and "UTM_ZONE" in editString: ...


1

I've come up with a solution that's maybe not the prettiest but it works. Below is code that will parse out the needed information from a feature class field after the Calculate UTM Zone tool is run, and then use it to create the spatial reference object. >>> cursor = arcpy.da.SearchCursor ("test", "utm_zone") >>> for utmStr, in cursor: ...


0

What are you actually doing? If you get your coordinates in to 4326 using Proj4, then you can add those points straight to the map, or form geojson points and add them. I'm suspicious you're overthinking it though! Check out the top example here: http://leafletjs.com/examples/geojson.html - it should be obvious from there how to go from lat/lon (WGS, ...


2

I'm not sure what NOAA thinks are the right coordinates, but I have no problem loading the file into QGIS, or reprojecting it to WGS84 with gdalwarp: QGIS uses this custom projection string: +proj=lcc +lat_1=25 +lat_2=25 +lat_0=25 +lon_0=265 +x_0=0 +y_0=0 +a=6371229 +b=6371229 +units=m +no_defs You can use the same string with pyproj.Proj(). where you ...



Top 50 recent answers are included